碰撞圆周率

最新推荐文章于 2022-01-23 21:04:02 发布

爱蹄子的羊头

最新推荐文章于 2022-01-23 21:04:02 发布

阅读量2k

点赞数 1

分类专栏： C++

本文链接：https://blog.csdn.net/warrially/article/details/89352334

版权

这篇博客分享了一段有趣的代码，用于计算圆周率，展现了编程中的数学之美。

摘要由CSDN通过智能技术生成


#include "stdafx.h"
#include <iostream>

// 碰撞后速度
void txpz(double m1, double m2, double & v1,  double & v2)
{
	double v1_ = ((m1-m2)* v1 + 2 * m2 * v2) / (m1 + m2);
	double v2_ = ((m2-m1) * v2 + 2 * m1 * v1) / (m1 +m2);

	v1 = v1_;
	v2 = v2_;
}

// 完全弹性碰撞( 返回碰撞次数 )
int cishu(double m1, double m2) 
{
	double v1 = -1.0;
	double v2 = 0;
	int nTime = 0;
	while (true) {
		if ((v1 >= 0) && (v2 < v1)) 
		{
			break;
		}

		txpz(m1, m2, v1, v2);
		nTime ++;

		if (v2 < 0)
		{
			v2 = -v2;
			nTime ++;
		}
	}
	using namespace std;
	cout<<m1<< "\t"<<m2<< "\t"<< nTime<<endl;
	return nTime;
}

int main()
{
	for (int i = 1; i < 8; i ++)
	{
		cishu(pow(100,i) * 1.0, 1.0);
	}
	
    return 0;
}

看一下一个有趣的代码

// 下面再放一份LUA版本