64,266
社区成员
发帖
与我相关
我的任务
分享
//如果要很高很高的精度,可以通过计算,以字符串来保存。
//bit为精确到小数点后的位数,datas为需要输出到的字符串的指针(datas至少要有bit+3个元素)
//最后几位可能是错的(如需准确的前n位数,可以带入n+15位去运算,最后的字符串截取至n位)。
void mult(int *a,int b,int *s,int bit);
void divi(int *a,int b,int *s,int bit);
void incr(int *a,int *b,int *s,int bit);
bool eqs(int *a,int *b,int bit);
void get_pi(int bit,char* datas)
{
int *pi=new int[bit+1];
int *ls=new int[bit+1];
int *sl=new int[bit+1];
int *p=new int[bit+1];
for(int i=0;i<=bit;i++)
*(pi+i)=*(ls+i)=*(sl+i)=*(p+i)=0;
memset(pi,0,sizeof(pi));
memset(ls,0,sizeof(ls));
memset(sl,0,sizeof(sl));
memset(p,0,sizeof(p));
*pi=*ls=*sl=1;
for(int is=1;;is++)
{
mult(ls,is,sl,bit);
divi(sl,2*is+1,ls,bit);
incr(pi,ls,p,bit);
if(eqs(pi,p,bit))
break;
int *t;
t=p;
p=pi;
pi=t;
}
mult(p,2,pi,bit);
datas[0]=*pi+'0';
datas[1]='.';
for(int in=1;in<=bit;in++)
datas[in+1]=*(pi+in)+'0';
datas[bit+2]='\0';
delete[] pi;
delete[] ls;
delete[] sl;
delete[] p;
}
//以下4个函数是get_pi函数内调用的函数
void mult(int *a,int b,int *s,int bit)
{
for(int i=bit,c=0;i>=0;i--)
{
int y=(*(a+i))*b+c;
c=y/10;
*(s+i)=y%10;
}
}
void divi(int *a,int b,int *s,int bit)
{
for(int i=0,c=0;i<=bit;i++)
{
int y=(*(a+i))+c*10;
c=y%b;
*(s+i)=y/b;
}
}
void incr(int *a,int *b,int *s,int bit)
{
for(int i=bit,c=0;i>=0;i--)
{
int y=(*(a+i))+(*(b+i))+c;
c=y/10;
*(s+i)=y%10;
}
}
bool eqs(int *a,int *b,int bit)
{
int i=0;
while(((*(a+i))==(*(b+i)))&&(i<=bit))
i++;
return i>bit;
}
const double E=2.71828182845904523536029;
const double Pi=3.14159265358979323846264338328;
#include <stdio.h>
long a=10000;
long b;
long c=2800;
long d;
long e;
long f[2801];
long g;
int main() {
for(;b-c;) f[b++]=a/5;
// while (1) {
// if (0==b-c) break;
// f[b]=a/5;
// b++;
// }
//f[0 - 2800] = 10000/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);
// while (1) {
// d=0;
// g=c*2;
// if (0==g) break;
// b=c;
// while (1) {
// d+=f[b]*a;
// f[b]=d%--g;
// d/=g--;
// --b;
// if (0==b) break;
// d*=b;
// }
// c-=14;
// printf("%.4d",e+d/a);
// e=d%a;
// }
return 0;
}