使用Park& Miller RNG制作大样本?

时间:2018-06-13 12:25:22

标签: c pointers random malloc

我正在尝试使用Park& Miller RNG ran0 1 2 3(3暂时不起作用)生成伪随机数来自 C中的数字食谱“。 这些生成器正常工作,因为它使得样本在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);}

2 个答案:

答案 0 :(得分:0)

scanf("%lo",&j);

这会将您输入的数字解释为八进制整数,即只有数字0到7的数字。

用u。

替换模式中的o
scanf("%lu",&j);

答案 1 :(得分:0)

您的问题是您的格式说明符。说明符&#34;%lo&#34;告诉scanf读取无符号八进制整数的字符串表示形式,并将其存储在相应参数指向的long对象中。

这意味着当您输入&#34; 1000&#34;时,其解释为八进制,而1000 8 与512 10 相同。然后,您为该数字添加2,这解释了为什么您会看到输入的#514; 1000&#34;。

使用&#34; d&#34;而不是&#34; o&#34;在您的格式说明符中表示十进制整数:

scanf("%ld", &j);