先 #include<math.h>
反3角函数有 acos(double),asin(double),atan(double),atan(double,double),返回值 double
cos(double x)
在C/C++中使用pi (π) 值
马青公式:
π=16arctan1/5-4arctan1/239
走进无限美妙的数学世界 http://amuseum.cdstm.cn/AMuseum/math/4/46/4_46_1004.htm
马青公式计算圆周率程序 - CSDN博客 https://blog.csdn.net/MHL_1208980380/article/details/54284631
用马青公式计算圆周率最精简的C++代码(万进制):
#include<stdio.h>
#include<iostream>
using namespace std;
int main(void)
{ //本程序每四位数输出结果,如果请求计算的位数不是4的整数倍,最后输出可能会少1~3位
long a[2]={956,80},b[2]={57121,25},i=0,j,k,p,q,r,s=2,t,u,v,N,M=10000;
printf("%9cMachin%6cpi=16arctan(1/5)-4arctan(1/239)nPlease input a number.n",32,32);
cin>>N,N=N/4+3;
long *pi=new long[N],*e=new long[N];
while(i<N)pi[i++]=0;
while(--s+1)
{
for(*e=a[k=s],i=N;--i;)e[i]=0;
for(q=1;j=i-1,i<N;e[i]?0:++i,q+=2,k=!k)
for(r=v=0;++j<N;pi[j]+=k?u:-u)u=(t=v*M+(e[j]=(p=r*M+e[j])/b[s]))/q,r=p%b[s],v=t%q;
}
while(--i)(pi[i]=(t=pi[i]+s)%M)<0?pi[i]+=M,s=t/M-1:s=t/M;
for(cout<<"3.";++i<N-2;)printf("%04ld",pi[i]);
delete []pi,delete []e,cin.ignore(),cin.ignore();
return 0;
}
原文链接: https://www.cnblogs.com/wxl845235800/p/8894191.html
欢迎关注
微信关注下方公众号,第一时间获取干货硬货;公众号内回复【pdf】免费获取数百本计算机经典书籍
原创文章受到原创版权保护。转载请注明出处:https://www.ccppcoding.com/archives/272648
非原创文章文中已经注明原地址,如有侵权,联系删除
关注公众号【高性能架构探索】,第一时间获取最新文章
转载文章受原作者版权保护。转载请注明原作者出处!