Мультипликативный метод

@nobo · Регистрация: 31.10.2013

Мультипликативный метод. Широкое распространение получил мультипликативный метод формирования случайной последовательности:
$https://www.cyberforum.ru/cgi-bin/latex.cgi? X(i+1)=aX(i)(mod M), i=0,1$ …..
где X(0)-произвольное*нечетное число, неотрицательное;
a- коэффициент, $https://www.cyberforum.ru/cgi-bin/latex.cgi? a = 8t \pm 3,$
t - любое целое положительное число;
M-значение модуля.

Необходимо сгенерировать до 1000 значений, проблема заключается в с циклом.

# Мультипликативный метод
# Дано

import numpy as np
# Мультипликативный метод
# Дано

x0=754211
c=456789
M=268435455
t=1
a=5
i=x0

"Решение"
Mult_i = np.array([(a * x0) % M for i in range(1000)])

Mult_2 = (a * x0) % M
Mult_3 = (a * Mult_2) % M
Mult_4 = (a * Mult_3) % M
Mult_5 = (a * Mult_4) % M
Mult_6 = (a * Mult_5) % M
Mult_7 = (a * Mult_6) % M
Mult_8 = (a * Mult_7) % M
Mult_9 = (a * Mult_8) % M
Mult_10 = (a * Mult_9) % M
.......
print(x0)
print(Mult_2)
print(Mult_3)
print(Mult_4)
print(Mult_5)
print(Mult_6)
print(Mult_7)
print(Mult_8)
print(Mult_9)
print(Mult_10)
.......

@ioprst · 14.08.2018, 08:48

nobo, у вас в формуле a*х_i, а в цикле a*х₀ => каждый раз одно значение, т.к. а и х₀ - константы. Может быть так?

Python

1
2
3

tmp = [(a * x0) % M] # init first elem
for i in range(1, 1000):
    tmp.append((a * tmp[i-1]) % M)

@nobo · 14.08.2018, 13:32 **[ТС]**

Спасибо большое! долго думал, как написать. Вы помогли!

@nobo 2 / 2 / 0 Регистрация: 31.10.2013 Сообщений: 210
		1
	Мультипликативный метод 13.08.2018, 20:37. Просмотров 2312. Ответов 2 Метки цикл (Все метки) Мультипликативный метод. Широкое распространение получил мультипликативный метод формирования случайной последовательности: $https://www.cyberforum.ru/cgi-bin/latex.cgi? X(i+1)=aX(i)(mod M), i=0,1$ ….. где X(0)-произвольноенечетное число, неотрицательное; a- коэффициент, $https://www.cyberforum.ru/cgi-bin/latex.cgi? a = 8t \pm 3,$ t - любое целое положительное число; M-значение модуля. Необходимо сгенерировать до 1000 значений, проблема заключается в с циклом. # Мультипликативный метод # Дано import numpy as np # Мультипликативный метод # Дано x0=754211 c=456789 M=268435455 t=1 a=5 i=x0 "Решение" Mult_i = np.array([(a x0) % M for i in range(1000)]) Mult_2 = (a * x0) % M Mult_3 = (a * Mult_2) % M Mult_4 = (a * Mult_3) % M Mult_5 = (a * Mult_4) % M Mult_6 = (a * Mult_5) % M Mult_7 = (a * Mult_6) % M Mult_8 = (a * Mult_7) % M Mult_9 = (a * Mult_8) % M Mult_10 = (a * Mult_9) % M ....... print(x0) print(Mult_2) print(Mult_3) print(Mult_4) print(Mult_5) print(Mult_6) print(Mult_7) print(Mult_8) print(Mult_9) print(Mult_10) ....... 0

@nobo 2 / 2 / 0 Регистрация: 31.10.2013 Сообщений: 210
	14.08.2018, 13:32 [ТС]	3
	Спасибо большое! долго думал, как написать. Вы помогли! 0

Опции темы