编程求pi
A. 编程计算pi值
1、对于精度设计,一般与时间成反比,算法和时间也是成正比的
2、时间和CPU和内存又是成比例的。
以下非本人算法:
#include
<cstdio>
using
namespace
std;
//某年Obfuscated
C
Contest佳作
long
a=10000,b,c=2800,d,e,f[2801],g;
main(){for(;b-c;)f[b++]=a/5;
for(;d=0,g=c*2;c-=14,printf(
"%.4d
",e+d/a),e=d%a)
for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);}
===========================
B. C语言编程用公式法求π,应该怎么做
给你提个思路:
在这个公式中有几个关键点:
每一项中有一个规律变化数字1,3,5,7,后面应该是9,11........;可以用i=1; i+=2;来表示。
每一项中都有规律变化分数相乘1/2 × 3/4 ×......这个乘法的截止方法为分母<i。
可以根据以上方法列出计算循环代码,最后一项小于10-6为跳出循环条件;
C. 怎么用C语言编程计算出π的值
一个足以让你吐血的计算400位pi的C程序:
int a=10000, b, c=2800, d, e, f[2801], g;
main() {
for(; b-c; f[b++]=a/5);
for(; d=0,g=c*2; c -=14,printf("%.4d",e+d/a),e=d%a)
for(b=c; d+=f[b]*a,f[b]=d%--g,d/=g--,--b; d*=b);
}
理论上只要用大数算法扩了就能一直算下去。我在Palm上面实现的那个就是基于这个算法,没有扩,算到4680位(受到Palm平台32k数组的长度限制)。
这个程序的特点是:你抱着想看看算pi原理的希望来读这个程序,结果发现就是看了也还是看不懂~~
转的,但是网络知道不允许贴地址哦
第二种:
用C语言编程π*π/6=1/1*1+1/2*2+1/3*3+。。。+1/n*n
#include<stdio.h>
#include<math.h>
main()
{
int n,i;
double s=0;
printf("输入n:");
scanf("%d",&n);
for(i=1;i<=n;i++)
{s=s+1.0/(i*i);
}
printf("pai=%f",sqrt(s*6));
}
第三种:
求pi=(1+1/(1*3))*(1+1/(3*5))*(1+1/(5*7))…………直到最后一项的分数小于10的负6次方为止。
#include<stdio.h>
#include<math.h>
#include<stdlib.h>
void main()
{
double pi,i,j,sigh,term;
i=1,j=3,sigh=1,term=1;
while(sigh>1e-6)
{
sigh=1/(i*j);
term=(1+sigh)*term;
i=i+2;
j=j+2;
}
pi=term;
printf("pi=%f\n",pi);
}
D. 编程求pi的值
利用“正多边形逼近”的方法求出π值在很早以前就存在,我们的先人祖冲之就是用这种方法在世界上第一个得到精确度达小数点后第6位的π值的。
利用圆内接正六边形边长等于半径的特点将边数翻番,作出正十二边形,求出边长,重复这一过程,就可获得所需精度的π的近似值。
假设单位圆内接多边形的边长为2b,边数为i,则边数加倍后新的正多边形的边长为:
周长为:
y=2 * i * x i:为加倍前的正多边形的边数
*程序与程序注释
#include<stdio.h>
#include<math.h>
void main()
{
double e=0.1,b=0.5,c,d;
long int i; //i: 正多边形边数
for(i=6;;i*=2) //正多边形边数加倍
{
d=1.0-sqrt(1.0-b*b); //计算圆内接正多边形的边长
b=0.5*sqrt(b*b+d*d);
if(2*i*b-i*e<1e-15) break; //精度达1e-15则停止计算
e=b; //保存本次正多边形的边长作为下一次精度控制的依据
}
printf("pai=%.15lf\n",2*i*b); //输出π值和正多边形的边数
printf("The number of edges of required polygon:%ld\n",i);
}
*运行结果
pi=3.141592653589794
The number of edges of required polygon:100663296
来源:考试大-计算机二级考试
E. 用C++编程计算π
#include<iostream>
#include<conio.h>
using namespace std;
void main()
{
double dPI;//π
int iFenMu;//分母
int iFenZi;//分子
int iCount;//计数
int iSign;//记录正负号
//初始化各变量
iFenZi = 4;
iFenMu = 1;
iSign = 1;
dPI = 0;
for(iCount=1;iCount<=1000;iCount++)//用for语句实现
{
//公式π=4 – 4/3 + 4/5 – 4/7 + 4/9 -4/11 + … …
dPI += ((iSign)*(((double)iFenZi)/((double)iFenMu)));//π的累加,相当于dPI=dPI+(((iSign)*(((double)iFenZi)/((double)iFenMu))));
cout<<"PI["<<iCount<<"]="<<dPI<<" ";//输出第iCount项的值
if(iCount%4==0)
{
cout<<endl;
}
iFenMu+=2;//分母增加2
iSign=-iSign;//正负号交替变化
}
cout<<"Press any key to continue"<<endl;
getch();
iFenZi = 4;
iFenMu = 1;
iSign = 1;
dPI = 0;
iCount =1;
while(iCount<=1000)//用while语句实现
{
dPI += ((iSign)*(((double)iFenZi)/((double)iFenMu)));
cout<<"PI["<<iCount<<"]="<<dPI<<" ";
if(iCount%4==0)
{
cout<<endl;
}
iFenMu+=2;
iSign=-iSign;
iCount++;
}
cout<<"Press any key to continue"<<endl;
getch();
iFenZi = 4;
iFenMu = 1;
iSign = 1;
dPI = 0;
iCount =1;
do//用do while语句实现
{
dPI += ((iSign)*(((double)iFenZi)/((double)iFenMu)));
cout<<"PI["<<iCount<<"]="<<dPI<<" ";
if(iCount%4==0)
{
cout<<endl;
}
iFenMu+=2;
iSign=-iSign;
iCount++;
}while(iCount!=1001);
}
F. C语言编程计算出π的值 最好有三种方法!! 急求 好人一生平安啊
我只说说π是怎么近似计算出来的,大一应该学过高数,可以知道:tan π/ 4 =1,从而π=4*arctan1,应用泰勒公式将arctan(x)展开,就可以得到
把这展开式展开,直到最后一项的绝对值小于1e-6为止。
下面是C语言递归实现:
#include<stdio.h>
#include<math.h>
void main()
{
float fuc(int n);
float s,pi;
for(int i=1;i<10000;i++)
{s=1.0/(2*i-1);
if(s<1e-6)
break;}
pi=fuc(i);
pi=pi*4;
printf ("pi=%10.6f ",pi);
}
float fuc(int n)
{
float pi;
if(n==1) pi=1;
else if(n%2==0)
pi=fuc(n-1)+(-1.0)/(2*n-1);
else pi=f(n-1)+1.0/(2*n-1);
return(pi);
}
程序运行结果pi=3.141498.
希望能明白明白π的实现原理,增进你的理解。
G. C语言编程用公式法求π,应该怎么做
给你提个思路:
在这个公式中有几个关键点:
每一项中有一个规律变化数字1,3,5,7,后面应该是9,11........;可以用i=1;
i+=2;来表示。
每一项中都有规律变化分数相乘1/2
×
3/4
×......这个乘法的截止方法为分母
追问:
能写个核心的代码吗,可以提高悬赏
追答:
午休的时候给你写吧!现在有事。
评论
0
0
0
加载更多
H. 用C语言编程计算π的值
#include<stdio.h>
#include<math.h>//开根号我用了数学库函数
int main()
{
double Pi=2,x=0,t;//我把第一项2作为Pi的初始值
while(1)
{
x=sqrt(x+2);//x是分母
t=2/x;//t是现在要乘的数
Pi=Pi*t;
if(t-1<1e-14) break;//因为t一定大于1,就这么写了;精确到小数点后14位
}
printf("%.14lf\n",Pi);//小数点后14位
return 0;
}