随机投点法-----------7-1 椭圆的面积 (100 分)

一、小明在完成一道求椭圆面积的题目时，忘记了椭圆面积的计算公式，你能编一个程序，利用随机投点法帮他求出椭圆面积的近似值吗？

椭圆的公式为：
a ^2 / x ^2​ + b ^2 / y ^2 ​ =1
提示：求[0,a]之间的随机数为：double x = double(rand())/RAND_MAX * a;

输入格式:
两个整数，椭圆的长半轴a和短半轴b，以空格分隔，0<b<a<10

输出格式:
椭圆的面积，只保留整数部分

输入样例:

6 3

输出样例:

56

二、代码

#include<iostream>
#include<cstdlib>  // 随机数头文件 
#include<ctime>    // 时间头文件 
using namespace std;
int a,b;

// 算法思路: 椭圆面积 / 长方形的面积 = k / n  
double Darts(int n){
	
	srand(time(0)); // 更新随机数 
	int k=0; // 统计点数 
	for(int i = 1;i <= n;i++){
		
		double x = (double)rand() / RAND_MAX * a; // 生成 [ 0, a ] 的随机数 
		double y = (double)rand() / RAND_MAX * b; // 生成 [ 0, b ] 的随机数 
		
		if((x * x) / (a * a) + (y * y) / (b * b) <= 1) // 统计落在椭圆内的点的个数 
			k++;
	}
	return 4 * ( a * b * k) / n; // 这里注意长方形的面积是 2a * 2b = 4ab
}

int main(){
	
	cin>>a>>b;
	cout<< Darts(100000) <<endl; // 提供的点数n尽量多点，这里取100000个 
	
	return 0;
}