通过这篇短文想说明两个道理：

   

   

1.    
2. 看似高大上、神秘兮兮的算法，都是paper tiger；
3.    
4. 计算机的计算方式（动辄几Ghz的主频）简直就是为蒙特卡洛度身定做；
5.    

 

一、 准备

1.1 生成任意区间任意大小的伪随机数

 

C语言中的rand()更深远的意义在于其对应于数学（概率论）中的均匀分布（uniformed distributed）。

C语言生成伪随机数的函数：

int rand(void);
  

• 1

该函数随机生成0~RAND_MAX之间内的整数：

#define RAND_MAX 0x7fff  // 0x7fff == 32767
  

• 1

产生随机数需要设置种子：

void srand(unsigned int _Seed);
  

• 1

这两个函数所在的头文件是stdlib.h或者cstdlib，后者又被包含在iostream头文件中。
有了rand()这个可以生成0-RAND_MAX随机数（整数）的函数，经过一定的四则运算和取模运算，便可很容易地得到任意区间的随机数。
以生成(2, 5)之间的随机数（可整可小）为例：

double x = 3*(double(rand())/RAND_MAX)+2;
  

• 1

先通过double(rand())/RAND_MAX使随机数区间转换为(0, 1)，再通过一定的伸缩平移实现对任意区间的仿真，这里的double类型转换不可省略，否则整数之间的除法运算得到的结果仍是整数。

• double vs float

两者的区别在于对浮点数表达的精度不同，double是双精度，float为单精度。

sizeof(double) == 8;sizeof(float) == 4;
  

• 1
• 2

c语言中的浮点数默认是double类型的，除非显式声明为l（或者L）

float x = 1l;
  

• 1

• long vs int

 

16位系统：long是4字节，int是2字节
  32位系统：long是4字节，int是4字节
  64位系统：long是8字节，int是4字节
  更详细的讨论见long vs int。

更详细的内容参见之前的一篇博客C++伪（pseudo）随机数生成及简单应用。

1.2 什么是蒙特卡洛算法？

 

这部分内容会比较枯燥，如果读不下去，可先看后边的实验，再读这部分内容会很容易理解和接受。

蒙特卡罗方法（Monte Carlo method），也称统计模拟方法，是一种数值计算方法。是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明，而被提出的一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数（或更常见的伪随机数）来解决很多计算问题的方法。

通常蒙特卡罗方法：因所求解的问题本身具有内在的随机性，借助计算机的运算能力可以直接模拟这种随机的过程。

更详尽的解释，参见Monte Carlo method, 蒙特卡洛算法在机器学习中的应用可参见增强学习-蒙特卡洛方法。

二、蒙特卡洛算法的实现

2.1 π

#include <stdio.h>#include <stdlib.h>#include <time.h>#include <math.h>int main(int, char**){    long c = 0, d = 0, N = 100000;    double x = 0, y = 0, pi = 0;    srand(unsigned(time(NULL)));    for (unsigned i = 0; i < N; ++i)    {        d += 1;        x = double(rand()) / RAND_MAX;        y = double(rand()) / RAND_MAX;        if (sqrt(x*x + y*y) > 1 && sqrt((x - 1)*(x - 1) + y*y) > 1)            c += 1;    }    printf("s = %f\n", double(c) / d);    return 0;}
   

• 1
• 2
• 3
• 4
• 5
• 6
• 7
• 8
• 9
• 10
• 11
• 12
• 13
• 14
• 15
• 16
• 17
• 18
• 19
• 20
• 21