Question

求问一道C++题目

用蒙特卡洛方法计算阴影部分面积C++，产生随机数，答案是19.5039

Answer 1

#include <iostream>
#include <random>
#include <ctime>
#include <cmath>
using namespace std;

int main()
{
    srand(time(0));
    double pointx;
    double pointy;
    double countnum=0;                //countnum为随机生成的10000000个点中在阴影部分的点
    for(int i=0;i<10000000;i++){     //随机生成10000000个点
        pointx=rand()%21;
        pointy=rand()%11;
        if(pointx/pointy>2.0){        //当点位于对角线下方时
                if(pointx<10&&pointx>=5){    //当点的x大于5小于10时
                    if(sqrt(pow(pointx-5,2)+pow((5-pointy),2))>=5)  //当点到左边第一个圆的距离大于半径时：
                        countnum++;
                }
                else  {
                        if(sqrt(pow(abs(pointx-15),2)+pow(abs(5-pointy),2))>=5)//当点到右边第一个圆的距离大于半径时：
                         countnum++;

                }

        }
    }

    cout<< countnum/10000000*100;  //求出出现在隐形部分的概率 乘于100为面积
    return 0;
}

追问

不可以运行哦。出错了
补了
#include 也不对……

追答

我用code::block 运行正常 应该是缺少某个头文件 贴一下你的报错信息

追问

我知道啦。是缺少

#include
可以运行了，但是和答案有偏差，这是为什么呢？

追答

不好意思 刚才的代码有点小问题 正确代码如下 误差已经解决

刚才点的取值没有取值到小数，导致精度缺失。 最好得到的概率应该乘于长方体的面积 得到阴影的面积：

#include <iostream>
#include <ctime>
#include <cmath>
#include <random>
using namespace std;

int main()
{
    srand(time(0));
    double pointx;
    double pointy;
    double countnum=0;

    for(int i=0;i<100000000;i++){
        pointx=rand()%20+double(rand()%100)/100;
        pointy=rand()%10+double(rand()%100)/100;
       // cout<<pointx<<" "<<pointy<<endl;
        if(pointy==0.0||pointx/pointy>2){
                if(pointx<10&&pointx>5){
                    if(sqrt(pow(pointx-5,2)+pow((5-pointy),2))>5)
                        countnum++;
                }
                else if(pointx<20&&pointx>=10) {
                        if(sqrt(pow(abs(pointx-15),2)+pow(abs(5-pointy),2))>5)
                         countnum++;

                }

        }
    }

    cout<< countnum/100000000*200;
        return 0;
}

追问

其实还是有点问题，但是启发了我思路，采纳啦，辛苦！

Answer 2

百度作业帮比较好吧！