С помощью процедуры привести дробь к несократимому виду Опции

[marwell]

с помощью процедуры Socr привести дробь 1+1/2+1/3+1/4+1/5+...+1/20 к несократимому виду

var a_1,a_2,b_1,b_2,p,q:integer;
function NOD(a:integer; b:integer):integer;
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;
function NOK (A: integer;  B: integer): integer;
begin
     NOK := (a * b) div NOD (a, b);
end;
procedure Socr(a,b:integer; var p,q:integer);
begin
if a=b then begin
p:=1;
q:=1;
end
       else if NOD (a, b)<>1 then
       begin
       p:=a div NOD(a, b);
       q:=b div NOD(a, b);
       end
       else Writeln ('дробь несократима');
       end;
begin
a_1:=1;
a_2:=1;
b_1:=1;
for b_2:=2 to 21 do
begin
a_1:=(NOK(b_1, b_2) div b_1)*a_1;
a_1:=a_1+((NOK(b_1, b_2) div b_2)*a_2);
b_1:=NOK(b_1, b_2);
end;
writeln (a_1,'/',b_1);
Socr(a_1,b_1,p,q);
Writeln (p,'/',q);
readln;
end.

посмотрите, please, правильно?

[volvo]

А ты не копируй КУСОК. Я дал ссылку на работоспособную программу. Почему надо ее сначала исковеркать, а потом удивляться и возмущаться, что она не работает? Ratio - это мой тип, я его там описывал...

[marwell]

А ты не копируй КУСОК. Я дал ссылку на работоспособную программу. Почему надо ее сначала исковеркать, а потом удивляться и возмущаться, что она не работает? Ratio - это мой тип, я его там описывал...

function NOD(A, B: longint): integer;
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;

function NOK(A, B: longint): integer;
begin
  NOK := a * b div NOD (a, b)
end;

type
  ratio = record
    num, denom: longint;
  end;

{ Процедура сокращения дробей (заодно корректируем знак) }
procedure Socr(var res: ratio);
var the_nod: longint;
begin
  repeat
    the_nod := NOD(abs(res.num), abs(res.denom));
    if the_nod > 1 then begin
      res.num := res.num div the_nod;
      res.denom := res.denom div the_nod;
    end;
  until the_nod = 1;

  if res.denom < 0 then begin
    res.num := res.num * (-1);
    res.denom := abs(res.denom);
  end;
end;

{ Сложение дробей }
procedure add_ratio(var res: ratio; curr: ratio);
var new_denom: longint;
begin
  Socr(curr);

  new_denom := NOK(res.denom, curr.denom);
  res.num :=
    res.num * (new_denom div res.denom) +
    curr.num * (new_denom div curr.denom);
  res.denom := new_denom;

  Socr(res);
end;

var
  curr: ratio;
  i: integer;
const
  res: ratio = (num:0; denom:1);

begin
  for i := 1 to 20 do
  begin
    curr.num := 1; curr.denom := i;
    add_ratio(res, curr);
  end;
  writeln(res.num, '/', res.denom);
readln;
end.

а так правильно? результат 6229/23780=0,2619428. Но ведь 1+1/2+1/3+1/4+...+1/20=3.59773966(примерно)

Сообщение отредактировано: marwell - 4.01.2010 20:24

[Lapp]

а так правильно? результат 6229/23780=0,2619428. Но ведь 1+1/2+1/3+1/4+...+1/20=3.59773966(примерно)

Исправь тип функций с integer на LongInt.
И никогда не отключай range check до окончания отладки.