c - 使用 Park&Miller RNG 制作大样本?
问题描述
我正在尝试使用“C 中的数字食谱”中的 C++ 中的 Park&Miller RNG ran0 1 2 3(3 暂时不起作用)生成伪随机数。这些生成器正常工作,因为它使样本均匀分布在 0 和 1 之间。我的目标是拥有一个巨大的样本(至少 10^10),所以我使用 malloc() 函数来分配内存。我的问题是样本总是很小(我想生成“j”个随机数,例如:如果我初始化j=1000 我将有 514 个随机数而不是 1000,为了测试我做了不同的“for”循环)。我是 C++(和编程)的初学者,我担心指针问题,但我不明白如何解决它。如果有人有解决方案或只是一个建议,这将对我有很大帮助,因为这个问题严重限制了我的实习(请宽容我的近似英语:'))。这是代码:
int main(){
FILE* simulationr0 = fopen("E:/Simulation_stage/SimulationRan0","w");
FILE* simulationr1 = fopen("E:/Simulation_stage/SimulationRan1","w");
FILE* simulationr2 = fopen("E:/Simulation_stage/SimulationRan2","w");
FILE* simulationr3 = fopen("E:/Simulation_stage/SimulationRan3","w");
float d;
long j;
long seed;
printf("tapez 0 pour utiliser ran0\ntapez 1 pour utiliser ran1\ntapez 2 pour utiliser ran2\ntapez 3 pour utiliser ran3\ntapez 4 pour les utiliser tous");
scanf("%f",&d);
if( d!=0 & d!=1 & d!= 2 & d!= 3 & d!=4){
printf("Erreur: valeur incorrecte//incorrect value\n");
exit(0);}
printf("Combien de nombres pseudo-aleatoires ?// How many random numbers ?\n");
scanf("%lo",&j);
printf("Quelle graine ? Which seed ?\n");
scanf("%lo", &seed);
long *pseed=&seed;
if(d==4){
float * ranp0;
float * ranp1;
float * ranp2;
float * ranp3;
ranp0 = (float*) malloc (j*sizeof(float));
ranp1 = (float*) malloc (j*sizeof(float));
ranp2 = (float*) malloc (j*sizeof(float));
ranp3 = (float*) malloc (j*sizeof(float));
for (int i=0;i<j+2;i+=1){
ranp0[i]=ran0(pseed);
ranp1[i]=ran1(pseed);
ranp2[i]=ran2(pseed);
ranp3[i]=ran3(pseed);
fprintf(simulationr0,"%f\n", ranp0[i]);
fprintf(simulationr1,"%f\n", ranp1[i]);
fprintf(simulationr2,"%f\n", ranp2[i]);
fprintf(simulationr3,"%f\n", ranp3[i]);
}
fclose(simulationr0);
fclose(simulationr1);
fclose(simulationr2);
fclose(simulationr3);
system("PAUSE");
free(ranp0);
free(ranp1);
free(ranp2);
free(ranp3);}
if(d==0){
float * ranp0;
ranp0 = (float*) malloc (j*sizeof(float));
for (int i=0;i<j+2;i++){
ranp0[i]=ran0(pseed);
fprintf(simulationr0,"%f\n", ranp0[i]);}
fclose(simulationr0);
system("PAUSE");
free(ranp0);}
if(d==1){
float * ranp1;
ranp1 = (float*) malloc (j*sizeof(float));
for (int i=0;i<j+2;i++){
ranp1[i]=ran1(pseed);
fprintf(simulationr1,"%f\n", ranp1[i]);}
fclose(simulationr1);
system("PAUSE");
free(ranp1);}
if(d==2){
float * ranp2;
ranp2 = (float*) malloc(j*sizeof(long));
for (int i=0;i!=j;i=i+1){
ranp2[i]=ran2(pseed);
fprintf(simulationr2,"%f\n", ranp2[i]);}
fclose(simulationr2);
system("PAUSE");
free(ranp2);}
if(d==3){
float * ranp3;
ranp3 = (float*) malloc(j*sizeof(float));
for (int i=0;i<j+2;i++){
ranp3[i]=ran3(pseed);
fprintf(simulationr3,"%f\n", ranp3[i]);}
fclose(simulationr3);
system("PAUSE");
free(ranp3);}
system("PAUSE");
exit(1);}
解决方案
scanf("%lo",&j);
这会将您输入的数字解释为八进制整数,即只有数字 0 到 7 的数字。
将模式中的 o 替换为 u。
scanf("%lu",&j);
推荐阅读
- sql - ORACLE SQL 根据“-”后整数的长度用零替换/填充“-”
- django - Django:在设计模型之间的关系时需要帮助
- python - 硒网络驱动程序:
- sql-server - 输出每个订单的最大联系日期,其中联系日期大于订单日期
- laravel-4 - 通过 Sendgrid 发送的电子邮件延迟交付给客户
- asp.net-mvc - 无法在 web.config 文件中使用两个连接字符串。一个是EntityConnection string Database first 一个是normal DBConnection
- apache-spark - 在 spark/pyspark (2.2*, 2.3.0) 中处理高维数据帧
- ansible - 运行安全的 linux ad-hoc 命令时,ansible 结果显示已更改
- java - 如何在 Flux.subscribe 之后收集数据以将其作为 json 数组发送?
- javascript - Moment.js - 从输入表单创建“时间”唯一变量