根据我的更广泛的答案here,下面是一个C++实现钛PRNG的:
#include
#include
using namespace std;
long mod1 = 2147483563;
long mod2 = 2147483399;
long mult1 = 40014;
long mult2 = 40692;
long seed1,seed2;
void Seed(int n){
if(n<0) //Perform an abs
n = -n;
if(n==0){
seed1 = 12345;
seed2 = 67890;
} else {
seed1 = (mult1*n)%mod1;
seed2 = n%mod2;
}
}
double Uniform(){
double result;
seed1 = (seed1*mult1)%mod1;
seed2 = (seed2*mult2)%mod2;
result = (double)(seed1-seed2)/(double)mod1;
if(result<0)
result = result+1;
return result;
}
double RandInt(int min, int max){
if(min
return min+int((max-min+1)*Uniform());
else
return max+int((min-max+1)*Uniform());
}
int main(){
Seed(0);
for(int i=0;i<10;i++)
cout<
}
我用了一个TI-83模拟器生成值RandInt(5,97):
这与我上面的程序的输出相匹配:
$ ./a.out
92
89
18
52
42
73
9
36
97
23