c++ - 蒙特卡洛法计算圆周率

Question

我正在尝试使用蒙特卡洛方法计算圆周率。但我总是得到零，我不知道为什么。这是我的代码

   #include <tchar.h>
   #include <Windows.h>
   #include <omp.h>
   #include <iostream>
   #include<math.h>
   using namespace std;
   int main(int argc, char *argv[]){

       int N = 1000, n = 0;
       double x = 0, y = 0;
       double answer; 

       for (int i = 0; i < N; i++){
         x = (double)rand() / (double)RAND_MAX;
         y = (double)rand() / (double)RAND_MAX;

         if (((x*x) + (y*y)) < 1)
                ++n;
       }
       //cout << "n = " <<n << endl; 

       answer = n / N; 
       cout << answer*4.0 << endl;

       //system("pause");
  }

Answer 1

answer 计算中的整数除法:

answer = n / N;

'nuff 说。

编辑 1:
今天是星期五，所以我将添加一些解释。

变量n 和N 声明为整数。

除法优先于任何转换或赋值。除法作为两个整数执行，然后截断小数部分。剩余的值转换为 double，然后分配给变量 answer。

请不要区分标识符。 n 和N 应该是不同的字母。这有助于作者和审稿人避免拼写错误。