[已解决]蒙特卡洛算法

一世轻尘

xdm，为啥我用蒙特卡洛算法计算pi的值在VScode和DEV中的运算结果不同啊？？？

1. #include<stdio.h>
   
2. #include<stdlib.h>
   
3. #include<time.h>
   
4. #include<math.h>
   
5. int main()
   
6. {
   
7.     float x,y,pi;
   
8.     long long i,D=0,C=0;
   
9.     srand((unsigned)time(NULL));
   
10.     for(i=0;i<100000000;i++)
    
11.     {
    
12.         C++;
    
13.         x=(float)rand()/RAND_MAX;
    
14.         y=(float)rand()/RAND_MAX;
    
15.         if(sqrt(x*x+y*y)<=1)
    
16.             D++;
    
17.     }
    
18.     pi=4.0*(float)D/C;
    
19.     printf("pi的值是%f\n",pi);
    
20.     return 0;
    
21. }

复制代码

最佳答案

月排行榜 / 总排行榜

jackz007

2020-12-28 10:33:50

        楼主的代码稍作改动，使用的也是 gcc（和楼主的编译器相同），无论怎么运行也不会是 4. 多啊

1. #include<stdio.h>
   
2. #include<stdlib.h>
   
3. #include<time.h>
   
4. 
   
5. int main()
   
6. {
   
7.         double x , y , pi                     ;
   
8.         unsigned long long D , C              ;
   
9.         srand((unsigned) time(NULL))          ;
   
10.         for(C = D = 0 ; C < 100000000 ; C ++) {
    
11.                 x = rand() / (double)RAND_MAX ;
    
12.                 y = rand() / (double)RAND_MAX ;
    
13.                 if(x * x + y * y <= 1) D ++   ;
    
14.         }
    
15.         pi = 4.0 * D / C                      ;
    
16.         printf("pi 的值是 %lf\n" , pi)        ;
    
17. }

复制代码

        编译、运行实况

1. D:\0002.Exercise\C>g++ -o pi pi.c
   
2. 
   
3. D:\0002.Exercise\C>pi
   
4. pi 的值是 3.141310
   
5. 
   
6. 
   
7. D:\0002.Exercise\C>pi
   
8. pi 的值是 3.141593
   
9. 
   
10. D:\0002.Exercise\C>pi
    
11. pi 的值是 3.141441
    
12. 
    
13. D:\0002.Exercise\C>

复制代码

跳转到最佳答案楼层

jackz007

        楼主的代码稍作改动，使用的也是 gcc（和楼主的编译器相同），无论怎么运行也不会是 4. 多啊

1. #include<stdio.h>
   
2. #include<stdlib.h>
   
3. #include<time.h>
   
4. 
   
5. int main()
   
6. {
   
7.         double x , y , pi                     ;
   
8.         unsigned long long D , C              ;
   
9.         srand((unsigned) time(NULL))          ;
   
10.         for(C = D = 0 ; C < 100000000 ; C ++) {
    
11.                 x = rand() / (double)RAND_MAX ;
    
12.                 y = rand() / (double)RAND_MAX ;
    
13.                 if(x * x + y * y <= 1) D ++   ;
    
14.         }
    
15.         pi = 4.0 * D / C                      ;
    
16.         printf("pi 的值是 %lf\n" , pi)        ;
    
17. }

复制代码

        编译、运行实况

1. D:\0002.Exercise\C>g++ -o pi pi.c
   
2. 
   
3. D:\0002.Exercise\C>pi
   
4. pi 的值是 3.141310
   
5. 
   
6. 
   
7. D:\0002.Exercise\C>pi
   
8. pi 的值是 3.141593
   
9. 
   
10. D:\0002.Exercise\C>pi
    
11. pi 的值是 3.141441
    
12. 
    
13. D:\0002.Exercise\C>

复制代码

一世轻尘

jackz007 发表于 2020-12-28 10:33
楼主的代码稍作改动，使用的也是 gcc（和楼主的编译器相同），无论怎么运行也不会是 4. 多啊
  
...

现在好了，可能是我之前忘记保存了，刚开始用VScode