高精度计算PI值

实验1.2：高精度计算PI值

所用公式：

#include <stdio.h>
#include <stdlib.h>

typedef struct list{
    int data;
    struct list *next;
    struct list *pre;
}list;

list* Initlist()
{
    list* head;
    head=(list*)malloc(sizeof(list));
    head->next=head->pre=head;
    return head;
}
list* Creatlist(list*head)
{
    int i;
    list *p;
    p=head;
    for(i=0;i<1000;i++)
    {
        list *q=(list*)malloc(sizeof(list));
        q->data=0;
        q->pre=p;
        p->next=q;
        q->next=head;//第一次循环时此时的head和p是一个东西，目的为把链表画成一个圈
        head->pre=q;
        p=p->next;
    }
    return head;
}
int main()
{
    int n,i,a,b;
    scanf("%d",&n);
    list *number,*sum;
    list *p,*q,*x;
    number=Initlist();
    sum=Initlist();
    number=Creatlist(number);
    sum=Creatlist(sum);
    number->next->data=2;//第一位储存2，即2*R(1)=2
    sum->next->data=2;//与上同理
    a=0,b=0;//分别是用来暂时存储进位和余数
    for(i=1;i<2000;i++)//循环两千次，确保精确度
    {


        p=number->pre;//做大数乘法时从链表的后方开始
        while(p!=number)//大于10则把10位的数字给b，个位数字放入data域中。
        {
            a=p->data*i+b;
            p->data=a%10;
            b=a/10;
            p=p->pre;
        }
        b=0;//清空b,为除法做准备
        p=p->next;//大数除法从链表的前方开始
        while(p!=number)//若计算出的数字为自然数，则直接放入data域；若等于0,或为小数，则要计算余数并给b。
        {
            a=p->data+b*10;
            p->data=a/(2*i+1);
            b=a%(2*i+1);
            p=p->next;
        }
        b=0;//清零
        p=number->pre;//大数加法均从末尾开始
        q=sum->pre;
        while(p!=number)//大于10进位，并储存个位数，进位数字赋给b。
        {
            a=p->data+q->data+b;
            q->data=a%10;
            b=a/10;
            p=p->pre;
            q=q->pre;
        }
    }
    printf("3.");
    x=sum->next->next;//从小数开始输出。
    for(i=0;i<n;i++)
    {
        printf("%d",x->data);
        x=x->next;
    }
    return 0;
}

当n=5时，PI=3.14159
当n=500时 ,
PI=3.14159265358979323846264338327950288419716939937510582097494459230781640628620899862803482534211706798214808651328230664709384460955058223172535940812848111745028410270193852110555964462294895493038196442881097566593344612847564823378678316527120190914564856692346034861045432664821339360726024914127372458700660631558817488152092096282925409171536436789259036001133053054882046652138414695194151160943305727036575959195309218611738193261179310511854807446237996274956735188575272489122793818301194912