2.3. Оценка L для алгоритмов вида ...
Для оценки длины отрезка апериодичности, характерного для последовательностей вида (4), можно воспользоваться следующей элементарной вероятностной моделью. Пусть имеется конечное множество, содержащее N различных чисел. Произведем последовательность независимых опытов, в каждом из которых из этого множества извлекается и записывается одно число. Вероятность извлечения каждого числа в каждом из опытов равна (выборка с возвратом). Обозначим через L случайную величину: номер опыта, в котором впервые будет извлечено уже записанное число.
Теорема 3. В указанной модели для любого
Доказательство. Обозначим записанные числа через Легко вычислить вероятность того, что L окажется равным где :
так как
а вероятность того, что совпадет с одним из различных чисел равна
Фиксируем произвольное и пусть . Тогда
Легко видеть, что
Так как то поэтому
и интересующую нас вероятность можно записать в форме
Рассмотрим теперь разбиение интервала точками где (рис. 7), и составим интегральную сумму для функции этому разбиению:
Из неравенства видно, что при
Рис. 7.
. Поэтому последнее слагаемое в этой сумме стремится при к нулю, и
Отсюда и из соотношения (9) следует (8). Теорема доказана.
Следствие. Так как математическое ожидание неотрицательной случайной величины с функцией распределения равно , то из (8) следует, что при
Рассмотренная в настоящем пункте модель предложена в статье [75], там же получена формула (10), а теорема 3 доказана Л. Н. Большевым и Д. И. Голенко. Несмотря на грубость модели, оценки (8) и (10) оказались полезными во многих опытах, связанных со способами «перемешивания», когда теоретико-числовая оценка периода практически невозможна, В качестве N выбирается количество возможных значений функции в конкретной ЭВМ; тогда можно ожидать, что длина отрезка апериодичности окажется порядка
Например, в методе п. 2.2.3 на ЭВМ «Стрела» значение функции Ф( (до нормализации) записывается в форме значит, На практике для ряда подобных алгоритмов были экспериментально обнаружены отрезки апериодичности длиной от до 3,5 105
Замечание. В книге [19] параграф, содержащий теорему 3, называется «Критерий проверки периодичности последовательности псевдослучайных чисел». Однако по этому «критерию» судить о качестве псевдослучайных чисел нельзя. У некоторых вполне удовлетворительных последовательностей, полученных по формуле (б), длина L гораздо больше, чем оценка (10) (фактически ). В то же время некоторые последовательности с не удовлетворяют ни одному из основных тестов § 3.