WikiSort.ru - Не сортированное

ПОИСК ПО САЙТУ | о проекте

Линейный конгруэнтный метод — один из методов генерации псевдослучайных чисел. Применяется в простых случаях и не обладает криптографической стойкостью. Входит в стандартные библиотеки различных компиляторов.

Описание

Линейный конгруэнтный метод был предложен Д. Г. Лемером в 1949 году.[1] Суть метода заключается в вычислении последовательности случайных чисел , полагая

где  — модуль (натуральное число, относительно которого вычисляет остаток от деления; ),  — множитель ( ),  — приращение ( ),  — начальное значение ( ).

Эта последовательность называется линейной конгруэнтной последовательностью. Например, для получим последовательность [2]

Свойства

Линейная конгруэнтная последовательность, определенная числами , , и периодична с периодом, не превышающим . При этом длина периода равна тогда и только тогда, когда[3]:

  1. Числа и взаимно простые;
  2. кратно для каждого простого , являющегося делителем ;
  3. кратно , если кратно .

Наличие этого свойства для случая , где  — число битов в машинном слове, было доказано М. Гринбергом (англ. M. Greenberg).[4] Наличие этого свойства для общего случая и достаточность условий были доказаны Т. Е. Халлом (англ. T. E. Hull) и А. Р. Добеллом (англ. A. R. Dobell).[5]

Метод генерации линейной конгруэнтной последовательности при называют мультипликативным конгруэнтным методом, а при  — смешанным конгруэнтным методом. При генерируемые числа будут иметь меньший период, чем при , но при определенных условиях можно получить период длиной , если  — простое число. Тот факт, что условие может приводить к появлению более длинных периодов, был установлен В. Е. Томсоном (англ. W. T. Thomson) и независимо от него А. Ротенбергом (англ. A. Rotenberg).[2] Чтобы гарантировать максимальность цикла повторения последовательности при , необходимо в качестве значения параметра выбирать простое число. Самым известным генератором подобного рода является так называемый минимальный стандартный генератор случайных чисел, предложенный Стивеном Парком (англ. Stephen Park) и Кейтом Миллером (англ. Keith Miller) в 1988 году. Для него , а .[6][7]

Наиболее часто практикуемым методом генерации последовательностей псевдослучайных чисел является смешанный конгруэнтный метод.[источник не указан 582 дня]

Часто используемые параметры

При выборе числа необходимо учитывать следующие условия:

1) число должно быть довольно большим, так как период не может иметь больше элементов;

2) значение числа должно быть таким, чтобы вычислялось быстро.

На практике при реализации метода исходя из указанных условий чаще всего выбирают , где  — число битов в машинном слове. При этом стоит учитывать, что младшие двоичные разряды сгенерированных таким образом случайных чисел демонстрируют поведение, далёкое от случайного, поэтому рекомендуется использовать только старшие разряды. Подобная ситуация не возникает, когда , где  — длина машинного слова. В таком случае младшие разряды ведут себя так же случайно, как и старшие.[2] Выбор множителя и приращения в основном обусловлен необходимостью выполнения условия достижения периода максимальной длины.

Печально известен «неудачный» (с точки зрения качества выходной последовательности) алгоритм RANDU, на протяжении многих десятилетий использовавшийся в самых разных компиляторах.

Для улучшения статистических свойств числовой последовательности во многих генераторах псевдослучайных чисел используется только часть битов результата. Например, в стандарте ISO/IEC 9899 на язык Си приведен (но не указан в качестве обязательного) пример функции rand(), принудительно отбрасывающей младшие 16 и один старший разряд.

#define RAND_MAX 32767 

static unsigned long int next = 1;

int rand(void)
{
  next = next * 1103515245 + 12345;
  return (unsigned int)(next/65536) % (RAND_MAX + 1);
}

void srand(unsigned int seed)
{
  next = seed;
}

Именно в таком виде функция rand() используется в компиляторах Watcom C/C++. Известны числовые параметры иных алгоритмов, применяемых в различных компиляторах и библиотеках.

Sourcemмножитель aслагаемое cиспользуемые биты
Numerical Recipes[9]23216645251013904223
Borland C/C++232226954771bits 30..16 in rand(), 30..0 in lrand()
glibc (used by GCC)[10] 231110351524512345bits 30..0
ANSI C: Watcom, Digital Mars, CodeWarrior, IBM VisualAge C/C++ [11]231110351524512345bits 30..16
C99, C11: Suggestion in the ISO/IEC 9899 [12]232110351524512345bits 30..16
Borland Delphi, Virtual Pascal2321347758131bits 63..32 of (seed * L)
Microsoft Visual/Quick C/C++232214013 (343FD16)2531011 (269EC316)bits 30..16
Microsoft Visual Basic (6 and earlier)[13]2241140671485 (43FD43FD16)12820163 (C39EC316)
RtlUniform from Native API[14]231 − 1 2147483629 (7FFFFFED16)2147483587 (7FFFFFC316)
Apple CarbonLib, C++11's minstd_rand0[15]231 − 1168070see MINSTD
C++11's minstd_rand[15]231 − 1482710see MINSTD
MMIX by Donald Knuth26463641362238467930051442695040888963407
Newlib26463641362238467930051bits 63...32
VAX's MTH$RANDOM,[16] old versions of glibc232690691
Java2482521490391711bits 47...16
Ранее во многих компиляторах:
RANDU231  655390

Возможность использования в криптографии

Хотя линейный конгруэнтный метод порождает статистически хорошую псевдослучайную последовательность чисел, он не является криптографически стойким. Генераторы на основе линейного конгруэнтного метода являются предсказуемыми, поэтому их нельзя использовать в криптографии. Впервые генераторы на основе линейного конгруэнтного метода были взломаны Джимом Ридсом (Jim Reeds), а затем Джоан Бояр (Joan Boyar). Ей удалось также вскрыть квадратические и кубические генераторы. Другие исследователи расширили идеи Бояр, разработав способы вскрытия любого полиномиального генератора. Таким образом, была доказана бесполезность генераторов на основе конгруэнтных методов для криптографии. Однако генераторы на основе линейного конгруэнтного метода сохраняют свою полезность для некриптографических приложений, например, для моделирования. Они эффективны и в большинстве используемых эмпирических тестах демонстрируют хорошие статистические характеристики[8].

См. также

Примечания

  1. D. H. Lehmer, Mathematical methods in large-scale computing units, Proceedings of a Second Symposium on Large-Scale Digital Calculating Machinery, 1949, Harvard University Press, Cambridge, Mass., 1951, pp. 141—146. MR 0044899 (13,495f)
  2. 1 2 3 Дональд Кнут. Том 2. Получисленные методы // Искусство программирования. Указ. соч. — С. 21—37.
  3. Кнут Д. Э., Искусство программирования. Том 2. Получисленные методы — Вильямс. 2001. с.21-37
  4. M. Greenberger, Method in randomness, Comm. ACM 8 (1965), 177—179.
  5. T.E. Hull and A.R. Dobell «Random Number Generators»,SIAM Review 4-3(1962),230-254
  6. "Бакнелл Д. М. Фундаментальные алгоритмы и структуры данных в Delphi. Библиотека программиста. 2002 год. журнал Delphi Informant Magazine. Глава 6.
  7. Stephen K. Park and Keith W. Miller (1988). Random Number Generators: Good Ones Are Hard To Find. Communications of the ACM 31 (10): 1192—1201
  8. 1 2 Брюс Шнайер. Глава 16. // Прикладная криптография.Триумф.2002. Указ. соч. — С. 275. Архивная копия от 28 февраля 2014 на Wayback Machine
  9. Numerical recipies in C. The art of scientific computing. 2-nd edition. - Cambridge University Press, 1992. - 925 pp.
  10. The GNU C library's rand() in stdlib.h uses a simple (single state) linear congruential generator only in case that the state is declared as 8 bytes. If the state is larger (an array), the generator becomes an additive feedback generator and the period increases. See the simplified code that reproduces the random sequence from this library.
  11. A collection of selected pseudorandom number generators with linear structures, K. Entacher, 1997. Проверено 16 июня 2012.
  12. Last public Committee Draft from April 12, 2011, page 346f. Проверено 21 декабря 2014.
  13. How Visual Basic Generates Pseudo-Random Numbers for the RND Function. Microsoft Support. Microsoft. Проверено 17 июня 2011.
  14. In spite of documentation on MSDN, RtlUniform uses LCG, and not Lehmer's algorithm, implementations before Windows Vista are flawed, because the result of multiplication is cut to 32 bits, before modulo is applied
  15. 1 2 ISO/IEC 14882:2011. ISO (2 September 2011). Проверено 3 сентября 2011.
  16. GNU Scientific Library: Other random number generators

Литература

Ссылки

Данная страница на сайте WikiSort.ru содержит текст со страницы сайта "Википедия".

Если Вы хотите её отредактировать, то можете сделать это на странице редактирования в Википедии.

Если сделанные Вами правки не будут кем-нибудь удалены, то через несколько дней они появятся на сайте WikiSort.ru .




Текст в блоке "Читать" взят с сайта "Википедия" и доступен по лицензии Creative Commons Attribution-ShareAlike; в отдельных случаях могут действовать дополнительные условия.

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

2019-2024
WikiSort.ru - проект по пересортировке и дополнению контента Википедии