Модификация многомерных псевдослучайных последовательностей с использованием двойственных lfsr-cns генераторов Текст научной статьи по специальности «Математика»

МОДИФИКАЦИЯ МНОГОМЕРНЫХ ПСЕВДОСЛУЧАЙНЫХ ПОСЛЕДОВАТЕЛЬНОСТЕЙ С ИСПОЛЬЗОВАНИЕМ ДВОЙСТВЕННЫХ LFSR-CNS ГЕНЕРАТОРОВ

А.Н. Калугин

Институт систем обработки изображений РАН Самарский государственный аэрокосмический университет имени академика С.П. Королева

Аннотация

В работе рассматривается новый метод модификации многомерной псевдослучайной последовательности точек, основанный на использовании пары двойственных генераторов ЬРЗЯ-С^. Состояние генератора, восстановленное по элементу многомерной последовательности, интерпретируется как состояние двойственного генератора, что позволяет сгенерировать точку, отличную от точки исходной последовательности. Приводятся сравнительные результаты исследования исходной и модифицированной последовательности с использованием взвешенного спектрального критерия.

Введение

Большинство современных методов решения задач статистической физики определенным образом используют методы Монте-Карло. Методы Монте-Карло при решении многомерных задач (физические задачи являются исходно многомерными), требуют формирования множества многомерных случайных (или псевдослучайных, квазислучайных) точек.

Для формирования многомерной псевдослучайной последовательности за редким исключением [1], [2], в большинстве случае используются последовательности случайных чисел (одна или несколько), которые [3], [4], [5] распределяются по координатам выходной многомерной последовательности. Такие методы «увеличения размерности» основаны как на разбиении одномерной псевдослучайной последовательности, так и использовании нескольких одномерных генераторов различных типов или одного типа с различными параметрами.

Иногда использование методов «увеличения размерности» одномерной последовательности обусловлено необходимостью (в виду специфики задачи или метода ее решения) использования параллельных алгоритмов обработки данных, а, следовательно, именно параллельных версий генераторов случайных точек. Однако для ряда задач [6] требуется именно многомерная последовательность псевдослучайных точек, а не параллельный способ ее формирования.

Известно [7], [8], [9], что метод «увеличения размерности» псевдослучайной последовательности чисел может привести к нежелательным результатам, а именно: к потере случайного характера итоговой последовательности, к несоответствию распределения элементов итоговой последовательности требуемому.

В данной работе предлагается схема модификации многомерной псевдослучайной последовательности, полученной методом «увеличения размерности», в целях улучшения ее многомерных свойств, основанная на использование пары двойственных ЬР8Я-С№ генераторов. С использованием взвешенного спектрального критерия проводится анализ применимости предложенной схемы (сравнение свойств исходной и модифицированной последова-

тельностей) на примере многомерной версии генератора ЯаМи.

Предварительные сведения. Генератор LFSR-CNS

Генератор ЬР8К-СШ [2] основан на использовании линейных рекуррентных

последовательностей в конечных полях канонических систем счисления в решетках. Приведем основные положения этих теорий.

Пусть есть конечное поле из q = р* элементов, р - простое число. Функцию, являющуюся решением линейного рекуррентного соотношения

У(п) = -Ь,-1 у(п -1) -... - ¿оу(п - *), (1)

где

¿о,...,6,-1 е GF(q), Ьо * 0, У = (у(0),...,у(,-1)),

будем называть линейной рекуррентной последовательностью порядка * (или, для краткости, рекуррентной функцией) с начальными условиями (значениями) У = (у(0),..., у(* -1)).

Рекуррентная последовательность (1) с максимально возможным периодом, равным Т = qs -1, называется т -последовательностью (см. [10]) . Последовательность {У (п)} = {У (0), У (1),...},

где

У(1) = (у(0, у0 +1),..., у0 + , - 1))Т (2)

называется «гусеницей последовательности (1)»

Заметим, что последовательность (2) может быть записана в матричном виде:

У (0 = (3)

где все арифметические операции выполняются в поле GF(q).

Матрица G е GF(q)sxs - матрица Фробениуса, сопровождающая матрица [11] характеристического многочлена рекуррентной функции рекуррентного соотношения (1).

Пусть, далее,{Й0,«1,-.,®к-1} - множество линейно независимых векторов пространства . Линей-

ная оболочка с целыми коэффициентами Л векто- где ров {â0,ai,...,âk_i} называется решеткой (lattice) в

Rk с базисом a = {â0,â1,...,âk-1}.

Л = {f 11 = Çoâo +|iâi +... + -iâk-1}, (4)

где ^ e Z , i = 0,1,...,k-1.

Пусть, далее, M - вложение M : Л ^ Л, det(M ) ф 0,

и множество D - конечное подмножество Л , D сЛ , 0 e D.

Тройка объектов (Л, M, D) называется системой счисления [12] (или, иначе, тройка (Л, M, D) обладает свойством единственности представления), если для любого элемента |еЛ существует единственное представление

i (I)

| = £ Mi (âi ), (5)

i=0

где âi e D и

l(1 ) = max(â ф 0), l(|) <œ . (6)

В этом случае, вложение M называется основанием системы счисления, а D - множеством цифр.

Далее, в качестве решетки Л, будем рассматривать решетку Л = Zk.

Система счисления (Zk, M, D) называется канонической, если D образует полную систему вычетов по модулю M и

D = {veo, v = 0,1,...,|detM| -1}, (7)

где eo = 1-ao + 0• a\ +... + 0• a¿_1; a = (5o,ai,...,5¿-i) -базис Zk.

Если (Zk, M, D) - каноническая система счисления, то для любого ZeZk существует единственное представление в виде

Z = ¿Zi (с (eo)), Zi £ {0,1,...,| detM|_1}. (8)

i=0

Вложения M могут быть сконструированы различным способом. Рассмотрим канонические системы счисления, порождаемые многочленами с целыми коэффициентам. Рассмотрим многочлен

f ( x) = ckxk + ck-1 xk-1 + ... + Со ,

(9)

где ск e Z , ck = 1.

В работе [13] показано, что для многочленов

f = xk + c1 x + q , (10)

где -1 < c1 < q - 2 q > 2, q = p ;

f2 = xk + pxk—1 + pxk—2 +... + px + p,

2 < р е N , q = р .

тройка ,Мг,Б), где Б = {0,1,...,q-1}, а М/ -вложение, имеющее в стандартном базисе решетки Zk вид сопровождающей матрицы многочлена (9) /

M f =

Г о 1 о

о

-co

1 —ck—1

Для Z3 многочлены (10), (11), порождающие бинарные канонические системы счисления имеют вид:

g1 = 2 + x3,

g2 = 2 + 2 x + 2 x2 + 2 x3.

(12) (13)

Общая схема генератора ЬР8К-С№ может быть записана в следующем виде:

Этап 1. Формирование гусеницы. Генерация «гусеницы» (2) последовательности (1). Используется т -последовательность, заданная линейным рекуррентным соотношением (1) в конечном поле из q элементов GF(q)

ДО = ^Т(0)^( д),

с ненулевыми начальными условиями У (0). Вектор У (0 в данном случае называется вектором состояния генератора.

Этап 2. Формирование многомерной точки. Векторы У (/') (2) рассматриваются как цифры записи некоторого числа в q -ичной канонической системе счисления.

^ = S Y (i) ;-1(Mj-1e) = ñ [G Y (0)]GF(q),

(14)

j=1

(11)

где Н е Zk, II = (МЧ,М1?,...,М5-1е)

Таким образом, каждому вектору состояния (2) У(/) поставлен в соответствие элемент й, е Zk .

Заметим, что вследствие единственности представления (5), различным состояниям генератора У (/) соответствуют различные элементы й е Zk

Для удобства практического использования далее используются генераторы с рекуррентной функцией (1) порядка

5 = гк, г е N. (15)

Генераторы ЬР8К-С№, использующие одинаковую т -последовательность на первом этапе, но различные системы счисления на втором этапе будем называть двойственными.

Этап 3. Унификация. Заметим, что в общем случае множество U точек й, е Zk , соответствующих всем возможным состояниям генератора У(,), имеет сложную нерегулярную форму [14], [2], однако для многочленов (10) и (11) существует взаимно однозначное инъективное вложение и в единич-

ный куб Ck = ([0,1))

P: U ^ Ck.

(16)

Вложение P называется унификацией фундаментальной области генератора. Пусть

P(u)t = щ , i = 0,1,...,k -1,

(17)

и ui - наименьшии неотрицательный вычет класса Ui (mod qt).

Тогда унификация (16) может быть записана в виде

P(U) = (U).

qq

Обозначив Fqi = t^ , общая схема генератора LFSR-CNS может быть записана в виде

щ= 1.[H[Gii?(0)]GF(q) ](Ft )k . (18)

Модификация существующих схем многомерной генерации

Сопровождающая матрица многочлена f1 при

q = 2 имеет следующий вид:

(о ... -2^ 10 ... о о '•. ;

м л =

V0 ••• 1 0 ...

V /k xk

(19)

Лемма 1. Для матрицы М/ (19) справедливо равенство

М/ +'е = (-2/М/ е , I = 0,1,..., к -1.

Доказательство. Нетрудно видеть, что единственным отличным от нуля главным минором (любого порядка) матрицы М/ является ее определитель, таким образом

Sj =

[0,1 < j < k; [det Mfj, j = k,

(20)

где Sj - сумма главных миноров матрицы М/ ] -го порядка.

Используя разложение по последнему столбцу, получаем:

det Mf = (-1)k • 2. (21)

Характеристическое уравнение матрицы M f

(-X)k + S1(-X)k-1 + S2(-X)k -2 +... +Sk-1(-X) + Sk = 0

с учетом (20) и (21) примет вид: (-X)k + (-1)k • 2 = 0,

или

Xk = -2 .

(22)

Применяя теорему Гамильтона-Кэли [11], из (22) получаем

(Mf )k = (-2) • Ik ,

(23)

где I к - единичная матрица порядка к.

Равенство (23) доказывает утверждение леммы. ■ Таким образом, точка на выходе ЬР8К-С№ генератора (без учета шага 3) может быть представлена в виде

И-1

й, = X У (0 j (М je) = j=0

= £(Mje)]>J(-2)l7(i)k/+j . j=0 /=0

(24)

Равенство (24) свидетельствует, что значения по каждой из координат в базисе (М0е, Ме, М2е,...,Мк-1е) на выходе второго этапа генератора отвечающего многочлену / (q = 2 ) генерируются независимо.

Также из (24) следует, что элементы последовательности Прj й, по каждой из координат получены

в результате интерпретации координат вектора состояния генератора, как цифр в представлении определенного целого числа в негабинарной системе счисления: элементы последовательности Пр j й,

j = 0,1,..., к -1 - суть элементы множества 2

Е = [ u

t-1

= Х (-2)1 ст/, VCTe(0,1}t

(25)

/=0

где 2 - множество целых чисел, имеющих длину представления в «негабинарной» системе счисления, не превышающую t.

Негабинарная система счисления, система счисления с основанием -2 и цифрами {0,1}, была впервые предложена Витторио Грюнвальдом в его работе [15]. Возможность единственного представления любого целого числа в негабинарной системе счисления следует, например, из существования соответствующего многочлена / + 2 = 0 (многочлена

типа /1 ), порождающего бинарную систему счисления. Интересной особенностью «негабинарной» системы счисления является возможность беззнакового представления всех целых чисел. Положительные целые числа имеют нечетную длину представления в негабинарной системе счисления, отрицательные целые числа - четную [16].

Заметим, что на шаге 3 процедуры генерирования ЬР8Я-СМ8 Ек отображается (взаимно однозначно) в куб Ск , более точно

на декартову степень множества

-1 ( п [0, ^)):

Р: Zk ^ Ск , Ск =

-1 ( п [0, qt))

Вид множества (25) обуславливает возможность применения пары генераторов ЬР8Я-СМ8, отвечающих многочленам /1 и /2 для модификации существующих схем многомерной генерации.

Схема модификации многомерной последовательности

1. Исходные данные. Пусть произвольный многомерный генератор (или многомерная версия произвольного одномерного генератора) производит точки последовательности {х(/)} точек многомерного куба Ск, или, после соответствующего масштабирования, точки множества Ск = [^ п [0, qt)]к,

х(') е Ск.

В таком случае, проекция точки х(') на каждую из координат может быть представлена в виде ее разложения в традиционной бинарной системе счисления

г-1

Пр} х(/) = £ 21 ст(/)((1), ст(/)(■/) е {0,1}г. (26)

/=0

2. Конверсия в негабинарную СС. С помощью достаточно простого алгоритма, можно перевести представление ст(') координат точки (26) в «нега-бинарную» систему счисления.

Алгоритм перевода основан на следующих фактах (27)-(30):

-1 • (-2)' = 1 • (-2)г+1 +1 • (-2)'; (27)

2 • (-2)' =-1- (-2)'+1 = 1- (-2)'+2 +1- (-2)'+1; (28)

1- (-2)2' = 1- 22'; (29)

1- 22'+1 =-1- (-2)2'+1. (30)

Блок-схема алгоритма перевода из бинарной в негабинарную систему счисления представлена на рис. 1.

Замечание. «Усечение» и(ст) на последнем шаге алгоритма до г бит корректно в силу единственности представления элемента в негабинарной системе счисления, а также взаимной однозначности унификации (16).

И хотя в общем случае и(ст) Ф ст, тем не менее, справедливо равенство:

х(0 = рр (и(0),

й( г) = (£(-2У и(ст(/)(0))/,

/=0

£(-2Уи(ст(/)(1))/,...,§(-2Уи(ст(/)(0))/)Т , (31)

/=0 /=0

где рР - оператор унификации.

3. Состояние генератора. В результате шагов 1, 2, по заданной точке X (/) выходной последовательности произвольного генератора мы получили состояние

У (')/к +1 =о(ст(/)( 1 ))/ , / = 0,1,...,г -1, 1 = 0,1,...,к -1

(32)

ЬР8Я-СМ8 генератора, соответствующего многочлену /1 .

< 1=0,1 м

нет

Тда

1=0,1,...,1-1

нет

Рис. 1. Блок-схема алгоритма перевода и(ст)

битового представления и е {0,1}г числа из бинарной системы счисления в негабинарную 4. Интерпретация. Интерпретируем состояние (32) как состояние ЬР8Я-СМ8 генератора, соответствующего многочлену /2 .

5. Вычисление модифицированной точки.

Действуя, далее, согласно второму этапу схемы генератора ЬР8К-С№, по состоянию генератора (32) вычислим модифицированную точку х(/) е Ск, или,

Исходный генератор

х(Г)еСк \-

Второй этап модифицированной схемы Таусворта

Рг0х{1) 1

1 1

Рг\х(0 1

о(/)об{0,7}'

после соответствующего масштабирования, точку многомерного единичного куба Ск .

Иллюстрация предложенной схемы модификации представлена на рис. 2. — 1

конверсия

Бинарная СС

конверсия

конверсия

о(о(0о)б{0,7}'

0(0(006 {0,7}'

о(о(/)ы)е{0,7}'

Негабинарная СС

х(¿)еС*

Выход

Состояние генератора ЬРБЯ-СЫБ, отвечающего многочлену /2

интерпретация

Состояние генератора

ЬРБЯ-Ст, отвечающего многочлену /,

те {о,1}'

ы

Сравнение свойств исходной и модифицированной многомерной последовательностей

Покажем, что использование предложенной схемы может привести к значительному улучшению характеристик модифицированного генератора по сравнению с исходным. В качестве «теоретической» оценки качества генератора случайных чисел служит значение определенного критерия (или его оценка), которое дает некоторые основания ожидать определенных характеристик поведения генератора в реальных задачах моделирования. На практике используются различные критерии, применимость которых ограничена типом исследуемого генератора [17], [18]. Относительно недавно, в качестве универсального критерия применимого для генераторов всех типов был предложен взвешенный спектральный критерий (ШарИопу) [19].

Сформулируем схему исследования генератора с использованием взвешенного спектрального критерия. Пусть Н = (Й0,-,Цс-\) еИк и х = (ха,...,Хк-1) е Ск . Пусть, далее тригонометрическая функция, соответствующая мультииндексу Н над Ск , определена равенством

к-1

е^ (х) = П е2*^^ ;

,=0

обозначение для мнимой единицы использовано в целях избежания нежелательных совпадений с индексом произведения.

Если ю = (хп ) -1 - множество точек из Ск, то-

V "/п=0

гда сумма Вейля [19] соответствующая мультиин-дексу Н для ю имеет вид:

Рис. 2. Схема модификации многомерной последовательности

N-1

1 N-1 ^ ^

SN е, ю) = N X е% (Хп), н * 0.

п=0

Функционал взвешенного спектрального критерия FN (ю) для множества точек ю задается соотношением

\1/2

(

FN (ю) =

1

Iх Го^ ^(ен'ю)|

(1 + п2 /3)к-1 г(Н )2

где г(Н) = Птах{1,| н } Н = (Н0,...,Нк-1) е .

¿=0

Последовательность ю = (хп )П=0 распределена

равномерно в Ск если, и только если, limN^да FN (ю) = 0 . Взвешенный спектральный критерий является мерой «нерегулярности» распределения точек последовательности ю .

Существует эффективный алгоритм позволяющий вычислять значения взвешенного спектрального критерия, основанный на следующем факте (лемма 2), доказанном в [20].

Лемма 2. Пусть

2 2

g(х) = 1 + (1 -2х)2, х е [0,1). 6 1

Пусть, далее

к-1

/(Х) = -1+ Пg(х,), Х = (Х0,...,Хк-1) е Ск .

I=0

Тогда для любого множества точек ю = (хп)

N-1 п=0 '

F?(ffl)=

1

r (i) = (r(ki),r(ki + 1),...,(ki + к - 1))r .

(35)

(1 + n2 /3)k -1

1 N-1N-1

• N2 X X f (- ^ mod 1).

N n=0 m=0

(33)

В случае если вычисление величины .2(ю) в

(33) аналитически связано с определенными трудностями, может быть проведено [17] численное исследование (ю) как функции N. Заметим, что вычислительная сложность вычисления .2(ю) по формуле (33) 2). Подобная асимптотика ограничивает величину N, которая может быть использована для исследования.

Как показано в [17], для равномерно распределенной случайной последовательности ю = (Хп ) математическое ожидание случайной величины р2(ю) равно

E ) = -1.

(34)

Свойство (34) позволяет построить следующую схему численного исследования последовательности [17]:

Шаг 1. Для размерности к и каждого из рассматриваемых N (N = N1,N2,■■■Nn ) с использованием тестируемого генератора сформируем К = 20 упорядоченных множеств юп, 0 < п < К (каждое множество юп состоит из N последовательных элементов х' выходной последовательности генератора).

Шаг 2. Для каждого из рассматриваемых N (N = N1,N2,...Nn ), используя (33), вычислим значения величины

N - ^2(юп), 0 < п < К .

Для равномерно распределенной случайной последовательности выборочное среднее N - .2(юп) должно быть близко 1.

Шаг 3. Если значения N - ( юп) близки 1 и не

наблюдается роста значения N • Fy(ro„) с ростом N = N1,N2,...Nn , качество распределения генерируемой последовательности признается приемлемым, иначе - неприемлемым.

Применим взвешенный спектральный критерий для исследования генератора Randu. Randu - линейный конгруэнтный генератор [21] со следующими параметрами: Randu= LCG(231,216 + 3 = 65539,0,1).

Пусть генератор Randu производит последовательность {r(i)}.

Многомерная последовательность {r (i)} для генератора Randu может быть получена по последовательности {r (i)} методом leap-frog [22]

Значения взвешенного спектрального критерия для последовательности (35) и размерности к = 3 приведены в таблице 1.

Таблица 1. Значения N • Fv(®n ) для генератора Randu, к = 3

N Выборочное среднее Максимум Минимум

512 1,02435 1,26149 0,77270

1024 1,04499 1,24142 0,91941

2048 1,03218 1,17598 0,84050

4096 1,07306 1,35123 0,92571

8192 1,07657 1,25191 0,92023

16384 1,13013 1,28906 0,95938

32768 1,26980 1,48367 1,08636

65536 1,57123 1,76214 1,23414

Можно заметить, что значения N - РдКюп) для {г (/)} возрастают с ростом N и имеют заметное отклонение от 1. Качество многомерной последовательности полученной с помощью генератора Randu признается [17], [18] неприемлемой в смысле взвешенного спектрального критерия.

Применив к многомерной последовательности (35), полученной на выходе генератора Randu, рассмотренную схему модификации, мы имеем следующие оценки взвешенного спектрального критерия для модифицированной последовательности (см. таблицу 2).

Таблица 2. Значения N - р2 (юп ) для модифицированной многомерной последовательности, полученной с помощью генератора Ка^й, к = 3

N Выборочное среднее Максимум Минимум

512 0,99924 1,13273 0,87289

1024 1,02247 1,24364 0,85744

2048 1,00165 1,23207 0,89696

4096 0,97197 1,10564 0,80539

8192 0,98350 1,22752 0,72497

16384 0,98558 1,25042 0,77338

32768 0,99881 1,22190 0,82052

65536 0,98671 1,20719 0,79435

Можно заметить, что модифицированная последовательность обладает приемлемыми (в смысле взвешенного спектрального критерия свойствами).

Значение N - .2(юп) не превышает 1 при N > 4096 ,

что свидетельствует об улучшении спектральных свойств.

Полученные экспериментальные данные (см. таблицу 2) свидетельствуют о возможности применения предложенной схемы модификации для улучшения свойств существующих схем многомерной генерации.

Заключение Предложенная схема модификации многомерной последовательности не зависит от специфики методов формирования исходной многомерной последовательности, существование генераторов LFSR-CNS отвечающих многочленам (10), (11) в пространстве любой размерности (начиная с двухмерного) позволяет модификацию практически любой многомерной псевдослучайной последовательности.

Следует отметить, что применение предложенной схемы должно сопровождаться отдельным исследованием свойств видоизмененного генератора для каждой конкретной исходной многомерной последовательности перед применением в реальных задачах.

Благодарности Работа выполнена при поддержке Министерства образования и науки РФ, правительства Самарской области (проект № 295Е2.3Д, дипломная работа), Американского фонда гражданских исследований и развития (CRDF Project SA-014-02) в рамках российско -американской программы «Фундаментальные исследования и высшее образование» (BRHE), а также при поддержке гранта Президента РФ № НШ-1007.2003.01, Российского фонда фундаментальных исследований (проекты №№ 03-0100736, 05-01-96501).

Литература

1. Wolfram S. Random sequence generation by cellular automata // Adv. Appl. Math. 7, 123 (1986).

2. Калугин А.Н. Трехмерное обобщение генератора LFSR случайных точек // Компьютерная оптика, 2005. № 27. С. 131-134

3. Coddington P., Random Number Generators for Parallel Computers, NHSE Review, Second Issue, Northeast Parallel Architectures Center, 1996. http://nhse.cs.rice.edu / NHSEreview / RNG

4. Entacher K. Parallel Streams of Linear Random Numbers in the Spectral Test // ACM Transactions on Modeling and Computer Simulation 9, 1999. № 1. С. 31-44.

5. Entacher K., Uhl A., Wegenkittl S. Parallel Random Number Generation: Long- range Correlations Among Multiple Processors // In P. Zinterhof, M. Vajtersic, and A. Uhl, editors, Parallel Computation, volume 1557 of Lecture Notes in Computer Science, Springer, New York, 1999. P. 107-116.

6. Vattulainen I. Framework for testing random numbers in parallel calculations // Phys. Rev. E, 59, 6, 7200 (1999).

7. Coddington P. Analysis of Random Number Generators Using Monte-Carlo Simulation, Int. J. Mod. Phys. 1994. C. 5. 547 р.

8. Coddington P. Tests of random number generators using Ising model simulations // Int. J. of Mod. Phys., 1996. C. 7(3). P. 295-303.

9. Ferrenberg A.M., Landau D.P. and Wong Y.J. Monte Carlo simulations: Hidden errors from good // Random number generators, Phys. Rev. Lett. 69, 3382 (1992).

10. Golomb S.W. Shift Register sequence // Holden-Day, San Francisco, 1967.

11. Гантмахер Ф.Р. Теория матриц // М.: Наука, 1988. -552 с.

12. Katai I., Kovacs B. Canonical number systems in imaginary quadratic fields // Acta Mathematica Academiae Sci-entarium Hungaricae. 37 (1-3), 1981. Р. 159-164.

13. Kovacs A., Generalized binary number systems // Annales Univ. Sci. Budapest, Sect. Comp. 20, 2001. Р. 195-206.

14. Chernov V.M. Fast uniform distribution of sequences for fractal sets // Proceedings of International Conference on Computer Vision and Graphics, 2004, September 22-24, 2004, Warsaw, Poland, Computational IMAGING AND VISION SERIES, Kluwer Academic Press (accepted for publication)

15. Vittorio Grunwald. Giornale di Matematiche di Battaglini (1885), 203-221, 367

16. Pawlek Z. and Wakulicz A. Bulletin de l'Academie Polonaise des Scienses, Classe III, 5 (1957), 233-236; Serie des sciences techniques 7 (1959), 713-721.

17. Hellekalk P., Niedderreiter H. The Weighted Spectral Test: Diaphony, ACM Trans. on Model. and Comp. Simul., 1998. Vol 8. No. 1 Р. 43-60.

18. Ripley B. Stochasitc Simulation // Wiley, New York, 1987.

19. Айерлэнд К., Роузен М. Классическое введение в современную теорию чисел // М.: Мир, 1987. - 416 с.

20. Zinterhof P. Über einige Abschätzungen bei der Approximation von Funktionen met Gleivhverteilundsmethoden. Sitzungsber. Österr. Akad. Wiss. Math.-Natur. Kl. II 185, 121-132 (1976).

21. Fishman G., Moore L. An exhaustive analysis of multiplicative congruential random number generators with modulus 231-1 // SIAM J. Sci. Statist. Comput. 1986. № 7. Р. 24-45.

22. Celmaster W. and Moriarty K.J.M. A method for vectorized random number generators // J. Comput. Phys. 1986. № 64. 271 р.