Генерация случайных чисел (Pascal)

[Korlet]

Написать программу, находящую при генерации случайных чисел повторяющиеся значения и подсчитывающую количество шагов между ними. Для генерации чисел использовать линейный конгруэнтный метод.

+ Использовать метод Р. Ковэю.
Заключается он в следующем:
Для случая, когда m представляется степенью двойки, интересный квадратичный метод предложил Р. Ковэю:
X0 mod 4=2, Xn+1=Xn(Xn+1) mod 2e, n0.
Данный метод гарантирует достаточно большой период.

Заранее благодарен,
Студент Фёдоров.

[FPaul]

Сообщение от Korlet

Студент Фёдоров.

А родители знают, что ты ЭТИМ занимаешься?

Исходный код из статьи в Википедии

#define RAND_MAX 32767 
 
static unsigned long int next = 1;
 
int rand(void)
{
next = next * 1103515245 + 12345;
return (unsigned int)(next/65536) % RAND_MAX;
}
 
void srand(unsigned int seed)
{
next = seed;
}

Копируется в Pascal

program primer;
 
{LCG - Linear congruential generator}
const
  LCG_M = 32767;
  LCG_A = 1103515245;
  LCG_C = 12345;
const
  {состояние}
  LCG_Curr: dword = 1;
 
  procedure LCG_Init(a: word);
  begin
    LCG_Curr := a;
  end;
 
  function LCG: word;
  begin
    LCG_Curr := LCG_Curr * LCG_A + LCG_C;
    LCG := hi(LCG_Curr) mod LCG_M;
  end;
 
var
  k: integer;
begin
  LCG_Init(1); {можно инициализировать секундами текущего времени}
  for k := 1 to 5 do
    writeln(LCG);
end.

Это, правда, не Ковэю, но переделка - небольшая.

[FPaul]

  function LCG: word;
  begin
    LCG_Curr := LCG_Curr * (LCG_Curr+1);
    LCG := hi(LCG_Curr) mod LCG_M;
  end;

[Korlet]

Нихрена не работает.

[FPaul]

У чёткого пацана - чёткая диагностика.

[Serge_Bliznykov]

FPaul, согласен.
Умеют же люди нескольким словами чётко и конкретно обрисовать ситуацию.

[FPaul]

Да. И ещё обидно за плохую грамматику - в данном случае "ни" - частица, а не предлог. А ведь интеллигент должен ещё материться уметь... Куда катится мир...

В общем, Korlet, посмотрел ещё раз код. Я действительно некорректно реализовал метод. Так будет корректнее.

  function LCG: word;
  begin
    LCG_Curr := (LCG_Curr * (LCG_Curr + 1)) mod LCG_M;
    LCG := LCG_Curr;
  end;

Там ещё в самой программе остались некорректные вещи. Я надеюсь, что после грамматики, ты найдёшь в интернете описание метода и самостоятельно исправишь две строки. Пойдёт на пользу.