VPF::[Pascal] Разложение числа на простые делители

Савин

Дата 6.12.2006, 21:47 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 3
Регистрация: 10.9.2006

Репутация: нет
Всего: нет

Добрый вечер.
Хотелось бы проконсультироваться насчет того, как в Паскале написать программу, которая для любого, введенного с клавиатуры числа, находила все его простые делители. Тоесть делала каноническое разложение: K=P(1)*P(2)*..*P(n), где Р - простые числа.
Спасибо.

volvo877

Дата 6.12.2006, 22:40 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Комодератор
Сообщений: 2073
Регистрация: 15.11.2004

Репутация: 34
Всего: 116

Код

procedure Factorization(x: word);
var i: word;

 procedure DivX;
 begin
      while (x>1) and (x mod i = 0) do
      begin
           write(i:4);
           x:= x div i;
      end;
 end;

begin
     i:=2;
     DivX;
     i:=3;
     while (i < x div 2) do
     begin
           DivX;
           inc(i,2);
     end;
     if x>1 then writeln(x:4);
end;

begin
  factorization(1200); { <--- Со вводом числа справишься? }
end.

(С) Romtek

Савин

Дата 6.12.2006, 23:49 (ссылка)

(нет голосов)

Загрузка ...

Новичок

Профиль
Группа: Участник
Сообщений: 3
Регистрация: 10.9.2006

Репутация: нет
Всего: нет

Спасибо, мне вывело 2 2 2 2 3 5 5, а можно как-то эти числа теперь в формулу загнать:
F= k(1-1/p(1))*(1-1/p(2))*..(1-1/p(n)), где к - это наше вводимое число, а р(1)..p(n) - простые числа, которые мы получили в этой программе. И помогите, пожалуйста, с вводом числа с клавиатуры, я не очень понимаю как это сделать.
Большое спасибо.

volvo877

Дата 7.12.2006, 10:20 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Комодератор
Сообщений: 2073
Регистрация: 15.11.2004

Репутация: 34
Всего: 116

Так ?

Код

function Factorization(x: word): real;
var i: word;
    p: real;

 function Mult(prev, next: real): real;
 begin
      Mult := prev * (1 - 1 / next);
 end;

 procedure DivX;
 begin
      while (x>1) and (x mod i = 0) do
      begin
           p := Mult(p, i);
           write(i:4);
           x:= x div i;
      end;
 end;
begin
     p := 1;

     i:=2;
     DivX;
     i:=3;
     while (i < x div 2) do
     begin
           DivX;
           inc(i,2);
     end;
     if x>1 then
     begin
           p := Mult(p, x);
           writeln(x:4);
     end;

     Factorization := p;
end;


var n: integer;
    coeff: double;
begin
  write('n = '); readln(n);
  coeff := factorization(n);
  writeln;
  writeln('F = ', (n * coeff):0:0);
end.

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Центр помощи \| Следующая тема »