高精度计算pi（C语言）

题目

解题入手

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

typedef struct node{
 int date;
 struct node* pre;
 struct node* next;
}node;

void insert(node* list,int date){
//尾插法向链表中添加元素
 node* p = list;
 if(p==NULL){
  return;
 }
 else{
  while(p->next){
   p = p -> next;
  }
   node* pnew=(node*)malloc(sizeof(node));
   pnew -> date = date;
   p -> next = pnew;
   pnew-> pre =p;
   pnew -> next = NULL;
 }
}

void tranverselist(node* list,int n){
//遍历链表并输出指定位数的小数
 node* p = list;
 p=p->next;
 printf("%d.",p->date);
 int i;
 p=p->next;
 for(i=0;i<n;i++){
  printf("%d",p->date);
  if(p->next){
   p=p->next;
  }
 }
 printf("\n");
}

int main()
{
    node *num,*sum;//num为每次相加的 R(n); um最终的值约等于 Π；
    num=(node*)malloc(sizeof(node));
     sum=(node*)malloc(sizeof(node));
    num -> pre = NULL;
   sum -> pre = NULL;
    num -> next = NULL;
    sum -> next = NULL;

   int n;
   scanf("%d",&n); //n为题目要求输出的的Π的小数位数；
    int i;
 for(i=0;i<520;i++){//补够位数
  insert(num,0);
  insert(sum,0);
 }
 node* p1 = num -> next;
 node* p2 = sum -> next;
  //使两个链表除头节点外第一个节点的数据为2
 p1->date=2;
 p2->date=2;
  //指向num,和sum的尾节点
 while(p1->next){
  p1=p1->next;
 }
 node* numtail = p1;
 while(p2->next){
  p2=p2->next;
 }
 node* sumtail = p2;
 //开始运算；
 int temp = 0;//计算过程中产生的临时数据
 int ret = 0; //计算过程中产生的借位or进位
 int t;//计算时有规律出现的除数 (2n+1)
 for(i=1,t=3;i<10000;i++,t=t+2){//根据数学的知识，算的次数越多越精确
   //i同时承担计数和做乘数的任务
   //先计算乘法（即公式中R(n)*n）
   node* p3 = numtail;
   ret=0;
   while(p3) {
    temp = p3 -> date * i + ret;
    p3 -> date = temp%10;
    ret = temp/10;
    p3 = p3 -> pre;
   }
   ret = 0;
   p3 = num -> next;
   //计算除法 (即公式中的R(n)*n/(2*n+1))
   while(p3){
     temp = p3 -> date + ret*10;
     ret = temp % t;
    p3 -> date = temp/t;
    p3=p3->next;
   }
   ret = 0;
   node*  p4 = sumtail;
   p3 = numtail;
   //将得到的R(n)的值加入到结果中
   while(p3&&p4){
    temp = p3 -> date + p4->date +ret ;
    ret = temp/10;
    p4->date=temp%10;
    p3=p3->pre;
    p4=p4->pre;
   }
 }
 tranverselist(sum,n);
 return 0;
}