前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >C++经典算法题-蒙地卡罗法求 PI

C++经典算法题-蒙地卡罗法求 PI

作者头像
cwl_java
发布2022-11-30 08:36:25
2760
发布2022-11-30 08:36:25
举报
文章被收录于专栏:cwl_Java

14.Algorithm Gossip: 蒙地卡罗法求 PI

说明

蒙地卡罗为摩洛哥王国之首都,该国位于法国与义大利国境,以赌博闻名。蒙地卡罗的 基本原理为以乱数配合面积公式来进行解题,这种以机率来解题的方式带有赌博的意味,虽然在精确度上有所疑虑,但其解题的思考方向却是个值得学习的方式。

解法

蒙地卡罗的解法适用于与面积有关的题目,例如求PI值或椭圆面积,这边介绍如何求PI值;假设有一个圆半径为1,所以四分之一圆面积就为PI,而包括此四分之一圆的正方形面积就为1,如下图所示:

在这里插入图片描述
在这里插入图片描述

如果随意的在正方形中投射飞标(点)好了,则这些飞标(点)有些会落于四分之一圆内,假设所投射的飞标(点)有n点,在圆内的飞标(点)有c点,则依比例来算,就会得到上图中最后的公式。

至于如何判断所产生的点落于圆内,很简单,令乱数产生X与Y两个数值,如果X2+Y2等于1 就是落在圆内。

代码语言:javascript
复制
#include <stdio.h>
#include <stdlib.h> 
#include <time.h>
#define N 50000 
    int main(void) {
        int i, sum = 0; double x, y;
        srand(time(NULL)); for(i = 1; i < N; i++) {
            x = (double) rand() / RAND_MAX; y = (double) rand() / RAND_MAX; if((x * x + y * y) < 1)
                sum++;
        }
        printf("PI = %f\n", (double) 4 * sum / N); return 0;
    }
本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2020-01-17,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 14.Algorithm Gossip: 蒙地卡罗法求 PI
    • 说明
      • 解法
      领券
      问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档