Как найти общий делитель двух чисел программа

Алгоритм Евклида

Для начала разберемся, что это и как это работает. Алгоритм Евклида позволяет найти нам наибольший общий делитель чисел. Как это работает:
Пусть a = 18, b = 30.
Цикл: a!=0 and b!=0
Если a > b, то a = a % b, если меньше, то b = b % a, таким образом мы сначала находим остаток деления, а потом повторяем действия. У нас a < b, значит, ищем остаток деления b % a (30 % 18) = 12, присваиваем b = 12, повторяем цикл но теперь у нас уже a > b(b = 12)
значит выполняем a % b (18 % 12) = 6? снова переходим к циклу, теперь снова b > a, значит выполняем b % a (30 % 6), остаток от деления 0, на этом мы прекращаем наш цикл и узнаем, что наибольший общий делитель 18 и 30 = 6. и выводим a + b (b = 0, a = 6).

Python

#!/usr/bin/env python
a = 18
b = 30
 
while a!=0 and b!=0:
    if a > b:
        a = a % b
    else:
        b = b % a
 
print (a+b)

Perl:

 sub nod
 {
 return  $_[0] != 0  ?  nod ( ( $_[1] % $_[0] ), $_[0] )  :  $_[1];
 }

C:

 int gcd(int a, int b) {
   int c;
   while (b) {
      c = a % b;
      a = b;
      b = c;        
   }
   return abs(a);
 }

Pascal:

  function nod( a, b: longint): longint; 
  begin
   while (a <> 0) and (b <> 0) do
     if a >= b then 
       a:= a mod b 
     else 
       b:= b mod a;
   nod:= a + b;
  end;

Java:

public class GCD {
    public static int gcd(int a,int b) {
        while (b !=0) {
            int tmp = a%b;
            a = b;
            b = tmp;
        }
        return a;
    }
}

C#:

int gcd(int a, int b)
{
    while (b != 0)
        b = a % (a = b);
    return a;
}

Given two integer numbers, the task is to find count of all common divisors of given numbers?

Examples : 

Input : a = 12, b = 24
Output: 6
// all common divisors are 1, 2, 3, 
// 4, 6 and 12

Input : a = 3, b = 17
Output: 1
// all common divisors are 1

Input : a = 20, b = 36
Output: 3
// all common divisors are 1, 2, 4

It is recommended to refer all divisors of a given number as a prerequisite of this article. 

Naive Solution 
A simple solution is to first find all divisors of first number and store them in an array or hash. Then find common divisors of second number and store them. Finally print common elements of two stored arrays or hash. The key is that the magnitude of powers of prime factors of a divisor should be equal to the minimum power of two prime factors of a and b.

  • Find the prime factors of a using prime factorization.
  • Find the count of each prime factor of a and store it in a Hashmap.
  • Prime factorize b using distinct prime factors of a.
  • Then the total number of divisors would be equal to the product of (count + 1) 
    of each factor.
  • count is the minimum of counts of each prime factors of a and b.
  • This gives the count of all divisors of a and b.

C++

#include <bits/stdc++.h>

using namespace std;

map<int, int> ma;

void primeFactorize(int a)

{

    for(int i = 2; i * i <= a; i += 2)

    {

        int cnt = 0;

        while (a % i == 0)

        {

            cnt++;

            a /= i;

        }

        ma[i] = cnt;

    }

    if (a > 1)

    {

        ma[a] = 1;

    }

}

int commDiv(int a, int b)

{

    primeFactorize(a);

    int res = 1;

    for(auto m = ma.begin();

             m != ma.end(); m++)

    {

        int cnt = 0;

        int key = m->first;

        int value = m->second;

        while (b % key == 0)

        {

            b /= key;

            cnt++;

        }

        res *= (min(cnt, value) + 1);

    }

    return res;

}

int main()

{

    int a = 12, b = 24;

    cout << commDiv(a, b) << endl;

    return 0;

}

Java

import java.util.*;

import java.io.*;

class GFG {

    static HashMap<Integer, Integer> ma = new HashMap<>();

    static void primeFactorize(int a)

    {

        for (int i = 2; i * i <= a; i += 2) {

            int cnt = 0;

            while (a % i == 0) {

                cnt++;

                a /= i;

            }

            ma.put(i, cnt);

        }

        if (a > 1)

            ma.put(a, 1);

    }

    static int commDiv(int a, int b)

    {

        primeFactorize(a);

        int res = 1;

        for (Map.Entry<Integer, Integer> m : ma.entrySet()) {

            int cnt = 0;

            int key = m.getKey();

            int value = m.getValue();

            while (b % key == 0) {

                b /= key;

                cnt++;

            }

            res *= (Math.min(cnt, value) + 1);

        }

        return res;

    }

    public static void main(String args[])

    {

        int a = 12, b = 24;

        System.out.println(commDiv(a, b));

    }

}

Python3

import math

ma = {}

def primeFactorize(a):

    sqt = int(math.sqrt(a))

    for i in range(2, sqt, 2):

        cnt = 0

        while (a % i == 0):

            cnt += 1

            a /= i

        ma[i] = cnt

    if (a > 1):

        ma[a] = 1

def commDiv(a, b):

    primeFactorize(a)

    res = 1

    for key, value in ma.items():

        cnt = 0

        while (b % key == 0):

            b /= key

            cnt += 1

        res *= (min(cnt, value) + 1)

    return res

a = 12

b = 24

print(commDiv(a, b))

C#

using System;

using System.Collections.Generic;

class GFG{

static Dictionary<int,

                  int> ma = new Dictionary<int,

                                           int>();

static void primeFactorize(int a)

{

    for(int i = 2; i * i <= a; i += 2)

    {

        int cnt = 0;

        while (a % i == 0)

        {

            cnt++;

            a /= i;

        }

        ma.Add(i, cnt);

    }

    if (a > 1)

        ma.Add(a, 1);

}

static int commDiv(int a, int b)

{

    primeFactorize(a);

    int res = 1;

    foreach(KeyValuePair<int, int> m in ma)

    {

        int cnt = 0;

        int key = m.Key;

        int value = m.Value;

        while (b % key == 0)

        {

            b /= key;

            cnt++;

        }

        res *= (Math.Min(cnt, value) + 1);

    }

    return res;

}

static void Main()

{

    int a = 12, b = 24;

    Console.WriteLine(commDiv(a, b));

}

}

Javascript

<script>   

      let ma = new Map();

    function primeFactorize(a)

    {

        for(let i = 2; i * i <= a; i += 2)

        {

            let cnt = 0;

            while (a % i == 0)

            {

                cnt++;

                a = parseInt(a / i, 10);

            }

            ma.set(i, cnt);

        }

        if (a > 1)

        {

            ma.set(a, 1);

        }

    }

    function commDiv(a,b)

    {

        primeFactorize(a);

        let res = 1;

        ma.forEach((values,keys)=>{

            let cnt = 0;

            let key = keys;

            let value = values;

            while (b % key == 0)

            {

                b = parseInt(b / key, 10);

                cnt++;

            }

            res *= (Math.min(cnt, value) + 1);

        })

        return res;

    }

    let a = 12, b = 24;

    document.write(commDiv(a, b));

</script>

Output: 

6

Time Complexity: O(√n log n) 
Auxiliary Space: O(n)

Efficient Solution – 
A better solution is to calculate the greatest common divisor (gcd) of given two numbers, and then count divisors of that gcd. 

C++

#include <bits/stdc++.h>

using namespace std;

int gcd(int a, int b)

{

    if (a == 0)

        return b;

    return gcd(b % a, a);

}

int commDiv(int a, int b)

{

    int n = gcd(a, b);

    int result = 0;

    for (int i = 1; i <= sqrt(n); i++) {

        if (n % i == 0) {

            if (n / i == i)

                result += 1;

            else

                result += 2;

        }

    }

    return result;

}

int main()

{

    int a = 12, b = 24;

    cout << commDiv(a, b);

    return 0;

}

Java

class Test {

    static int gcd(int a, int b)

    {

        if (a == 0)

            return b;

        return gcd(b % a, a);

    }

    static int commDiv(int a, int b)

    {

        int n = gcd(a, b);

        int result = 0;

        for (int i = 1; i <= Math.sqrt(n); i++) {

            if (n % i == 0) {

                if (n / i == i)

                    result += 1;

                else

                    result += 2;

            }

        }

        return result;

    }

    public static void main(String args[])

    {

        int a = 12, b = 24;

        System.out.println(commDiv(a, b));

    }

}

Python3

from math import sqrt

def gcd(a, b):

    if a == 0:

        return b

    return gcd(b % a, a)

def commDiv(a, b):

    n = gcd(a, b)

    result = 0

    for i in range(1,int(sqrt(n))+1):

        if n % i == 0:

            if n/i == i:

                result += 1

            else:

                result += 2

    return result

if __name__ == "__main__":

    a = 12

    b = 24;

    print(commDiv(a, b))

C#

using System;

class GFG {

    static int gcd(int a, int b)

    {

        if (a == 0)

            return b;

        return gcd(b % a, a);

    }

    static int commDiv(int a, int b)

    {

        int n = gcd(a, b);

        int result = 0;

        for (int i = 1; i <= Math.Sqrt(n); i++) {

            if (n % i == 0) {

                if (n / i == i)

                    result += 1;

                else

                    result += 2;

            }

        }

        return result;

    }

    public static void Main(String[] args)

    {

        int a = 12, b = 24;

        Console.Write(commDiv(a, b));

    }

}

PHP

<?php

function gcd($a, $b)

{

    if ($a == 0)

        return $b;

    return gcd($b % $a, $a);

}

function commDiv($a, $b)

{

    $n = gcd($a, $b);

    $result = 0;

    for ($i = 1; $i <= sqrt($n);

                 $i++)

    {

        if ($n % $i == 0)

        {

            if ($n / $i == $i)

                $result += 1;

            else

                $result += 2;

        }

    }

    return $result;

}

$a = 12; $b = 24;

echo(commDiv($a, $b));

?>

Javascript

<script>

    function gcd(a, b)

    {

        if (a == 0)

            return b;

        return gcd(b % a, a);

    }

    function commDiv(a, b)

    {

        let n = gcd(a, b);

        let result = 0;

        for (let i = 1; i <= Math.sqrt(n); i++) {

            if (n % i == 0) {

                if (n / i == i)

                    result += 1;

                else

                    result += 2;

            }

        }

        return result;

    }

    let a = 12, b = 24;

    document.write(commDiv(a, b));

</script>

Output :  

6

Time complexity: O(n1/2) where n is the gcd of two numbers.
Auxiliary Space: O(1)

This article is contributed by Aarti_Rathi and Shashank Mishra ( Gullu ). If you like GeeksforGeeks and would like to contribute, you can also write an article using write.geeksforgeeks.org or mail your article to review-team@geeksforgeeks.org. See your article appearing on the GeeksforGeeks main page and help other Geeks.

Please write comments if you find anything incorrect, or you want to share more information about the topic discussed above.

Another Approach:

1. Define a function “gcd” that takes two integers “a” and “b” and returns their greatest common divisor (GCD) using the Euclidean algorithm.
2. Define a function “count_common_divisors” that takes two integers “a” and “b” and counts the number of common divisors of “a” and “b” using their GCD.
3. Calculate the GCD of “a” and “b” using the “gcd” function.
4. Initialize a counter “count” to 0.
5. Loop through all possible divisors of the GCD of “a” and “b” from 1 to the square root of the GCD.
6. If the current divisor divides the GCD evenly, increment the counter by 2 (because both “a” and “b” are divisible by the divisor).
7. If the square of the current divisor equals the GCD, decrement the counter by 1 (because we’ve already counted this divisor once).
8. Return the final count of common divisors.
9. In the main function, define two integers “a” and “b” and call the “count_common_divisors” function with these integers.
10. Print the number of common divisors of “a” and “b” using the printf function.

C

#include <stdio.h>

int gcd(int a, int b) {

    if(b == 0) {

        return a;

    }

    return gcd(b, a % b);

}

int count_common_divisors(int a, int b) {

    int gcd_ab = gcd(a, b);

    int count = 0;

    for(int i = 1; i * i <= gcd_ab; i++) {

        if(gcd_ab % i == 0) {

            count += 2;

            if(i * i == gcd_ab) {

                count--;

            }

        }

    }

    return count;

}

int main() {

    int a = 12;

    int b = 18;

    int common_divisors = count_common_divisors(a, b);

    printf("The number of common divisors of %d and %d is %d.n", a, b, common_divisors);

    return 0;

}

C++

#include <bits/stdc++.h>

using namespace std;

int gcd(int a, int b) {

    if(b == 0) {

        return a;

    }

    return gcd(b, a % b);

}

int count_common_divisors(int a, int b) {

    int gcd_ab = gcd(a, b);

    int count = 0;

    for(int i = 1; i * i <= gcd_ab; i++) {

        if(gcd_ab % i == 0) {

            count += 2;

            if(i * i == gcd_ab) {

                count--;

            }

        }

    }

    return count;

}

int main() {

    int a = 12;

    int b = 18;

    int common_divisors = count_common_divisors(a, b);

    cout<<"The number of common divisors of "<<a<<" and "<<b<<" is "<<common_divisors<<"."<<endl;

    return 0;

}

Java

import java.util.*;

public class Main {

  public static int gcd(int a, int b) {

    if(b == 0) {

      return a;

    }

    return gcd(b, a % b);

  }

  public static int countCommonDivisors(int a, int b) {

    int gcd_ab = gcd(a, b);

    int count = 0;

    for(int i = 1; i * i <= gcd_ab; i++) {

      if(gcd_ab % i == 0) {

        count += 2;

        if(i * i == gcd_ab) {

          count--;

        }

      }

    }

    return count;

  }

  public static void main(String[] args) {

    int a = 12;

    int b = 18;

    int commonDivisors = countCommonDivisors(a, b);

    System.out.println("The number of common divisors of " + a + " and " + b + " is " + commonDivisors + ".");

  }

}

Python3

import math

def gcd(a, b):

    if b == 0:

        return a

    return gcd(b, a % b)

def count_common_divisors(a, b):

    gcd_ab = gcd(a, b)

    count = 0

    for i in range(1, int(math.sqrt(gcd_ab)) + 1):

        if gcd_ab % i == 0:

            count += 2

            if i * i == gcd_ab:

                count -= 1

    return count

a = 12

b = 18

common_divisors = count_common_divisors(a, b)

print("The number of common divisors of", a, "and", b, "is", common_divisors, ".")

C#

using System;

public class MainClass

{

    public static int GCD(int a, int b)

    {

        if (b == 0)

        {

            return a;

        }

        return GCD(b, a % b);

    }

    public static int CountCommonDivisors(int a, int b)

    {

        int gcd_ab = GCD(a, b);

        int count = 0;

        for (int i = 1; i * i <= gcd_ab; i++)

        {

            if (gcd_ab % i == 0)

            {

                count += 2;

                if (i * i == gcd_ab)

                {

                    count--;

                }

            }

        }

        return count;

    }

    public static void Main()

    {

        int a = 12;

        int b = 18;

        int commonDivisors = CountCommonDivisors(a, b);

        Console.WriteLine("The number of common divisors of {0} and {1} is {2}.", a, b, commonDivisors);

    }

}

Javascript

function gcd(a, b) {

    if(b === 0) {

        return a;

    }

    return gcd(b, a % b);

}

function count_common_divisors(a, b) {

    let gcd_ab = gcd(a, b);

    let count = 0;

    for(let i = 1; i * i <= gcd_ab; i++) {

        if(gcd_ab % i === 0) {

            count += 2;

            if(i * i === gcd_ab) {

                count--;

            }

        }

    }

    return count;

}

let a = 12;

let b = 18;

let common_divisors = count_common_divisors(a, b);

console.log(`The number of common divisors of ${a} and ${b} is ${common_divisors}.`);

Output

The number of common divisors of 12 and 18 is 4.

The time complexity of the gcd() function is O(log(min(a, b))), as it uses Euclid’s algorithm which takes logarithmic time with respect to the smaller of the two numbers.

The time complexity of the count_common_divisors() function is O(sqrt(gcd(a, b))), as it iterates up to the square root of the gcd of the two numbers.

The space complexity of both functions is O(1), as they only use a constant amount of memory regardless of the input size.

Last Updated :
13 Apr, 2023

Like Article

Save Article

Вариант 1

var a,b:longint;
 
function NOD(x,y:longint):longint; { функция поиска наиб. общ. делителя }
begin
   if x<>0 then NOD:=NOD(y mod x,x) else NOD:=y;
end;
 
function NOK(x,y:longint):longint; { функция поиска наим. общ. кратного }
begin
   NOK:=( x div NOD(x,y) ) * y;
end;
 
begin { основная программа }
    readln(a,b);
    writeln( 'НОД этих чисел = ', NOD(a,b) );
    writeln( 'НОК этих чисел = ', NOK(a,b) );
end.

Вариант 2 Переборный алгоритм

var a, b, d: integer;
begin 
    write('Введите два числа: ');
    readln(a, b);
    if a < b then d := a + 1 else d := b + 1;
    {так как мы используем цикл с постусловием, необходимо минимальное значение увеличить на один, 
    иначе цикл repeat, в силу своих конструктивных 
    особенностей, не учтет это минимальное число и 
    не сделает его кандидатом в НОД. Например, 5 и 25.}
    repeat d := d - 1 
    until (a mod d = 0) and (b mod d = 0); 
write('NOD = ', d) 
end.

Вариант 3

var
m,n,r:integer;
label lb;
begin
write('Введите первое число:');readln(m);
write('Введите второе число:');readln(n);
lb:r:=m mod n;
if r=0 then writeln('НОД = ',n)
else
 begin
  m:=n;
  n:=r;
  goto lb;
 end;
end.

Вариант 4 Алгоритм Евклида с вычитанием

Пусть a и b — целые числа, тогда верны следующие утверждения:

Все общие делители пары a и b являются также общими делителями пары a — b, b;

И наоборот, все общие делители пары a — b и b являются также общими делителями пары a и b; НОД(A,  B) = НОД(A — B, B), если A > B; НОД(A, 0) = A.

Доказательство:

Если t — произвольный общий делитель a и b, то он делит и разность a — b. Действительно, из a = t * u и b = t * v следует, что a — b = t * u — t * v = t * (u — v). То есть t — также общий делитель а — b и b. Обратно, если t — произвольный делитель общий делитель a — b и b, то он делит и их сумму a — b + b = a. Это можно доказать аналогично предыдущему. Поэтому t — также общий делитель a и b. Делаем вывод, что множество общих делителей a и b совпадает с множеством делителей a — b и b. В частности, совпадают и наибольшие общие делители этих пар. Наибольшее целое, на которое делится число a, есть само число а. Число 0 делится на любое число. Отсюда наибольший общий делитель а и 0 равен а. Доказанная формула(3) позволяет свести вычисление наибольшего делителя одной пары к вычислению наибольшего общего делителя другой пары, в которой числа уже меньше. Очевидная же формула (4) дает нам понять, когда надо остановиться.

var a, b: integer;
begin 
    write('a = ');
    readln(a);
    write('b = ');
    readln(b);
    while a <> b 
        do if a > b then a := a - b else b := b - a;
    writeln('NOD = ', a);
end.

Вариант 5 Алгоритм Евклида с делением

Пусть a и b — целые числа, а r — остаток от деления a на b. Тогда НОД(a, b) = НОД(b, r). Эта формула также позволяет свести вычисление наибольшего общего делителя одной пары чисел к вычислению наибольшего обшего делителя другой пары чисел.

var a, b: integer;
begin
    write('a = ');
    readln(a);
    write('b = ');
    readln(b);
    while (a <> 0) and (b <> 0) 
        do if a >= b then a := a mod b else b := b mod a;
    write(a + b)
end.

Вариант № 6

Program test2(input,output);
Const N = 5;
Var
       С: array[1..5] of integer;
       A,B:integer;
function HOК (A, В:integer):integer;
begin
        HOK:=A*B/ HOD(A,B);
end;
function НОD(А, В:integer):integer;
var
      X,Y:integer;
begin
X:= A; Y: = В;
1:IF X = Y THEN HOD:=X;
IF X > Y THEN begin
                           X:= X – Y;goto 1;
                           end;
IF Y > X THEN begin
                           Y:= Y – X;goto 1;
                           end;
end;
Begin
FOR i= 1 ТО N READ (C[i]);
A:= С ([l])
FOR i = 1 TO N–1 begin B:=С[i + 1];
                                          A:= HOK(A,B);
                               end;
writeln ("HOK="; A);
end.

Вариант 7

Program N_O_D (Input, Output); 
Var 
A, B: LongInt; 
NOD : LongInt; 
 
Begin
 
WriteLn ('PASCAL: Нахождение Н.О.Д. двух заданных чисел.');
Writeln ('Введите числа, для которых ищется НОД:'); 
Write('Первое число: ');ReadLn (A);
Write('Второе число: ');ReadLn (B); 
 
If (A < B)ThenNOD := A Else NOD := B;
 
While Not( (A mod NOD = 0) and (B mod NOD = 0) ) do 
NOD := NOD - 1;
 
WriteLn ('НОД = ',NOD);
 
ReadLn;
End. 

Program N_O_D (Input, Output); 
Var 
A, B: LongInt; 
NOK, NOD : LongInt; 
 
Begin
 
WriteLn ('PASCAL: Нахождение Н.О.К. двух заданных чисел.');
WriteLn ('Введите числа, для которых ищется НОК:'); 
Write ('Первое число: ');ReadLn (A); 
Write ('Второе число: ');ReadLn (B);
 
If (A < B)ThenNOD := A Else NOD := B;
 
While Not ( (A Mod NOD = 0) And (B Mod NOD = 0) ) Do 
NOD := NOD - 1;
 
A := A Div NOD;
B := B Div NOD; 
NOK := A * B * NOD; 
WriteLn ('НОК = ', NOK); 
 
ReadLn;
End. 

Алгоритм Евклида

Приветствуем читателей и посетителей нашего сайта! Сегодня на learnpascal.ru открывается новая рубрика — Алгоритмы. В этой рубрике мы с вами будем разбирать различные алгоритмы, а также их реализацию на Паскале.

Для освоения материала сегодняшнего урока вам понадобится знание циклов и ветвлений.

Сегодня мы рассмотрим три алгоритма(из пяти) на нахождение наибольшего общего делителя двух целых чисел, два из которых непосредственно связывают с именем Евклида. Еще два мы рассмотрим в следующей части.
Наибольший общий делитель (НОД) двух чисел a и b — наибольшее целое число, которое делит их оба.
Пример: НОД(25, 5) = 5; НОД(12, 18) = 6.

Переборный алгоритм

Начинаем перебор с d — наименьшего из двух чисел. Это первый, очевидный кандидат на роль их наибольшего общего делителя. А затем, пока d не делит оба числа, уменьшаем его на единицу. Как только такое деление будет обеспечено, останавливаем уменьшение d.

var
  a, b, d: integer;

begin
  write('Введите два числа: ');
  readln(a, b);
  if a < b then d := a + 1 else d := b + 1; 
  {так как мы используем цикл с постусловием, 
  необходимо минимальное значение увеличить на один,
  иначе цикл repeat, в силу своих конструктивных особенностей,
  не учтет это минимальное число  
  и не сделает его кандидатом в НОД. Например, 5 и 25.}
  repeat d := d - 1
  until (a mod d = 0) and (b mod d = 0);
  write('NOD = ', d)
end.

Обратимся к этой программе, например, с числами 30 и 18. Тогда на пути к ответу(числу 6) ей придется перебрать числа: 18, 17, 16, 15, 14, 13, 12, 11, 10, 9, 8, 7, 6.

Алгоритм Евклида «с вычитанием»

Пусть a и b — целые числа, тогда верны следующие утверждения:

  1. Все общие делители пары a и b являются также общими делителями пары a — b, b;
  2. И наоборот, все общие делители пары a — b и b являются также общими делителями пары a и b;
  3. НОД(A,  B) = НОД(A — B, B), если A > B;
  4. НОД(A, 0) = A.

Доказательство:

  1. Если t — произвольный общий делитель a и b, то он делит и разность a — b. Действительно, из a = t * u и b = t * v следует, что a — b = t * u — t * v = t * (u — v). То есть t — также общий делитель а — b и b.
  2. Обратно, если t — произвольный делитель общий делитель a — b и b, то он делит и их сумму a — b + b = a. Это можно доказать аналгично предыдущему. Поэтому t — также общий делитель a и b.
  3. Делаем вывод, что множество общих делителей a и b совпадает с множеством делителей a — b и b. В частности, совпадают и наибольшие общие делители этих пар.
  4. Наибольшее целое, на которое делится число a, есть само число а. Число 0 делится на любое число. Отсюда наибольший общий делитель а и 0 равен а.

Доказанная формула(3) позволяет свести вычисление наибольшего делителя одной пары к вычислению наибольшего общего делителя другой пары, в которой числа уже меньше. Очевидная же формула (4) дает нам понять, когда надо остановиться.

Вкратце алгоритм Евклида «с вычитанием» будет таким. Вычитаем из большего числа меньшее и заменяем большее на разность до тех пор, пока одно из чисел не обратится в нуль. Тогда оставшееся ненулевое число — наибольший общий делитель.

Пример. Пусть а = 82 и b = 60. НОД(82, 60) = НОД(22, 60) = НОД(22, 38) = НОД(22, 16) = НОД(6, 16) = НОД(6, 10) = НОД(6, 4) = НОД(2, 4) = НОД(2, 2) = НОД(2, 0) = 2.

На предпоследнем шаге алгоритма, перед появлением 0, оба числа равны, иначе не мог возникнуть 0. Поэтому мы будем извлекать НОД именно в этот момент.

Блок — схема алгоритма Евклида «с вычитанием»

Алгоритм Евклида

Программа

var
  a, b: integer;

begin
  write('a = '); 
  readln(a);
  write('b = ');
  readln(b);
  while a <> b do 
    if a > b then 
      a := a - b 
    else 
      b := b - a;
  writeln('NOD = ', a);
end.

Алгоритм Евклида с «делением»

Пусть a и b — целые числа, а r — остаток от деления a на b. Тогда НОД(a, b) = НОД(b, r).

Эта формула также позволяет свести вычисление наибольшего общего делителя одной пары чисел к вычислению наибольшего обшего делителя другой пары чисел.

Пример. НОД(82, 60) = НОД(22, 60) = НОД(22, 16) = НОД(6, 16) = НОД(6, 4) = НОД(2, 4) = НОД(0, 2) = 2.

var
  a, b: integer;

begin
  write('a = ');
  readln(a);
  write('b = ');
  readln(b);
  while (a <> 0) and (b <> 0) do
    if a >= b 
    then 
      a := a mod b 
    else 
      b := b mod a;
  write(a + b)
end.

На сегодня все! Еще несколько модификаций алгоритма Евклида и способов нахождения НОД вы узнаете на следующих уроках.

Программа для нахождения НОД чисел

Исходник программы, задача которой – нахождение наибольшего общего делителя (НОД или Алгоритм Евклида) для двух, введённых с клавиатуры чисел. Используется цикл WHILE, есть пояснительные комментарии ко всем важным строкам программы. Скачать исходник просмотреть исходный код программы можно ниже.

Исходный код программы:

var a,b,c: integer; //Описание переменных
begin //Начало программы
  writeln('Введите a,b: '); //Диалог с пользователем
  read(a,b);  //Чистывание чисел
  while b<>0 do //Вход в цикл while, пока b не равно 0
  begin
    c := a mod b; //Присваивание с остатка деления a/b
    a := b; 
    b := c;
  end;
  writeln('Наибольший Общий Делитель = ',a); //Вывод делителя
end.//конец программы

Скачать:

nod.pas

Дата: 2012-04-07 18:52:48   Просмотров: 46713

Теги: Паскаль исходник Pascal WHILE циклы

Добавить комментарий