|
ru.algorithms
- RU.ALGORITHMS ----------------------------------------------------------------
From : Serhiy Savychenko 2:5020/400 21 Sep 2001 18:06:39
To : Dmitry Pankov
Subject : Re: Генерирование случайных чисел
--------------------------------------------------------------------------------
Приветствую, "Dmitry Pankov" <Dmitry.Pankov@p27.f58.n5022.z2.fidonet.org>!
Вы сообщили:
> >> DP> Слышал в кратце про метод вычетов, который использует
рекуррентное
> >> DP> соотношение x[i]=b*x[i-1]%M, где % - остаток от деления, а b и M
>
> SS> В новом издании есть и новые методы, типа x[i]=(b*x[i-33]-x[i-97])%M
> SS> или что-то в этом роде. Причем весь смысл именно в этих цифрах 13 и
53.
>
> Что-то я про 13 и 53 не понял. Это b и M должны ими быть?
Очепятка, имелись ввиду конечно 33 и 97. Это метод Фибоначчи.
Важно что этот метод начинает работать, только если уже есть 97 чисел
(думаю понятно почему). Вот тот алгоритм сначала и заполняет массив, а
потом уже начинает генерировать.
> SS> Этот метод хорош для игрушек, или других программ с использованием
ввода
> SS> оператора. Для метода Монте-Карло например он не подходит.
>
> Вот мне как раз ^^^^^^^^^^^^^^^^^^^^^^ и нужен сабж.
Раз тебе нужен сабж для метода Монте-Карло, то очень важны длинный период
и равномерность распределения. Поэтому приведенный алгоритм лучше всего.
> SS> Hиже привожу код на фортране, тут более понятен алогоритм, кроме того
> SS> читай комментарии вначале о качестве метода. Кнут, кстати, хвалит
этот
> SS> алгоритм. Есть реализация на асме, но немного старая. Я проводил года
три
> SS> назад тесты, остался доволен.
>
> [SKIP]
>
> Жуть, как наворочено ;( А нет чего-нибудь попроще, не очень раздутого на
более
> популярном языке (на Си или Паскале например)?
А чем фортран не нравиться? Это же эха алгоритмов, а не языков. В
принципе, как
я уже писал, этот алгоритм в свое время переправил вообще на асм.
Сергей.
--- ifmail v.2.15dev5
* Origin: Digital Generation (2:5020/400)
Вернуться к списку тем, сортированных по: возрастание даты уменьшение даты тема автор
|
