資料結構
網路
RDBMS
作業系統
Java
MS Excel
iOS
HTML
CSS
Android
Python
C 程式設計
C++
C#
MongoDB
MySQL
Javascript
PHP
物理
化學
生物
數學
英語
經濟學
心理學
社會研究
時裝研究
法學使用 Park-Miller 隨機數生成演算法實現 C++ 程式
Park-Miller 隨機數生成演算法是生成隨機數的另一種方法。
這種型別的隨機數生成器 (RNG) 的通用公式為：X_{k+1} = g X(k) mod n
其中模數 n 是素數或素數的冪，乘法器 g 是模 n 的高乘法階元素，種子 X0 與 n 互素。
演算法
Begin
Declare variables n, a, b, c and seed
Read variables n, a, b, c and seed
Uniform()
Declare variable hi, lo, t
hi=seed divided by b
lo = seed - b * hi
t = a * lo - c * hi
if (t > 0)
seed = t;
else
seed = t + n;
return seed;
Done
For i =0 to n
Call the function random
Done
End示例程式碼
#include <iostream>
using namespace std;
const long n = 2145678965L;
const long a = 763214L;
const long b = 88844L;
const long c = 7766L; i
static long seed = 12345678L;
double uniform() {
long hi = seed / b;
long lo = seed - b * hi;
long t = a * lo - c * hi;
if (t > 0)
seed = t;
else
seed = t + n;
return seed;
}
int main(int argc, char **argv) {
double A[10];
for (int i = 0; i < 10; i++)
A[i] = uniform();
cout << "Random numbers are:\n";
for (int i = 0; i < 10; i++)
cout << A[i] << endl;
}輸出
Random numbers are: 6.50293e+10 4.27187e+10 2.1539e+10 4.62058e+10 1.70792e+10 8.24569e+09 5.93381e+10 3.63839e+10 4.81931e+10 8.91007e+09
更新於： 30-Jul-2019
401 次瀏覽
廣告