Программа печати всех делителей натурального числа A

[8Observer8]

Привет всем!

Правильно ли написана вот эта программа: http://acmp.ru/article.asp?id_text=141

Я для удобства её приведу:

{ Составить программу печати всех делителей натурального числа A }
var a,n,c,d:word;
begin { основная программа }
    readln( a );
    n:=1;
    while ( n <= sqrt(a) ) do begin
       c:=a mod n;
       d:=a div n;
       if c = 0 then begin
          writeln( n );
          if n <> d then writeln( d );
       end;
       inc( n );
    end;
end.

Меня смущает строка:

while ( n <= sqrt(a) ) do begin

А именно, меньше или равно. Функция sqrt возвращает вещественный тип и может возникнуть такая ситуация: sqrt(49) = 7.0000000000000001 Я прав? Ведь n будет приведён к вещественному. Сравнивать на равенство два вещественных типа - некорректно.

[Serge_Bliznykov]

во-первых, напишите так:

while ( sqr(n) <= a )  do begin

во-вторых, можно использовать приведение к целому типу через trunc: n <= Trunc(sqrt(N)+0.1)

в-третьих, имхо, и так, как у Вас написано, будет корректно работать (хотя, Ваша подозрительность, на мой взгляд, вполне оправдана!)

[Poma][a]

Уж обсуждали
тыц
тыц

[8Observer8]

Спасибо большое, парни!) Как посоветовали так и сделал. Только я на С++:

#include <iostream>
#include <cmath>
#include <cstdlib>
using namespace std;

int main(int argc, char** argv) {
    int a;
    cin >> a;

    int b = trunc(sqrt(a)+0.1);

    int n = 1;
    div_t divresult;

    while (n <= b) {
        divresult = div(a, n);
        if (divresult.rem == 0) {
            cout << n << endl;
            if (n != divresult.quot) {
                cout << divresult.quot << endl;
            }
        }
        ++n;
    }
    
    return 0;
}