注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

fanhq666的博客

Fan-Fun

 
 
 

日志

 
 

pi  

2008-05-31 11:36:11|  分类: 程序 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
#include <cstdlib>
#include <iostream>
#include <fstream>
#define N 20010
using namespace std;
void mult (int *a,int b,int *s)
{
     for (int i=N,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)
{
     for (int i=0,c=0;i<=N;i++)
     {
         int y=(*(a+i))+c*10;
         c=y%b;
         *(s+i)=y/b;
     }
}
void incr(int *a,int *b,int *s)
{
     for (int i=N,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 i=0;
     while (((*(a+i))==(*(b+i)))&&(i<=N)) i++;
     return i>N;
}
int main(int argc, char *argv[])
{
    int lpi[N+1],lls[N+1],lsl[N+1],lp[N+1];
    
    int *pi=lpi,*ls=lls,*sl=lsl,*p=lp;
    for (int i=0;i<=N;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 i=1;true;i++)
    {
        mult(ls,i,sl);
        divi(sl,2*i+1,ls);
        incr(pi,ls,p);
        if (eqs(pi,p)) break;
        int *t;
        t=p;
        p=pi;
        pi=t;
        if (i>=20010) break;
        if (i%50==0) cout << i << "  ";
    }
    cout << endl;
    mult(p,2,pi);
    ofstream fout("pi.txt");
    fout << *pi << ".";
    for (int i=1;i<=N;i++)
    {
        fout << *(pi+i);
        if (i%10==0) fout << " ";
        if (i%80==0) fout << endl;
    }
    return EXIT_SUCCESS;
}
  评论这张
 
阅读(303)| 评论(0)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017