高精度计算PI值

高精度计算PI值

所用公式：

#include

#include

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

{

printf("%d",x->data);

x=x->next;

}

return 0;

}

当n=5时，PI=3.14159

当n=500时 ,

PI=3.14159265358979323846264338327950288419716939937510582097494459230781640628620899862803482534211706798214808651328230664709384460955058223172535940812848111745028410270193852110555964462294895493038196442881097566593344612847564823378678316527120190914564856692346034861045432664821339360726024914127372458700660631558817488152092096282925409171536436789259036001133053054882046652138414695194151160943305727036575959195309218611738193261179310511854807446237996274956735188575272489122793818301194912

高精度计算PI值相关教程

如何在SQL Server计算XX年第XX周是哪几天

如何在SQL Server计算XX年第XX周是哪几天 SQL Server有一个很强大的函数：DATEPART() 它可以计算2017-10-17是属于今年的 第几天：select datepart(dy,'2017-10-17') 第几周：select datepart(wk,'2017-10-17') 第几月：select datepart(mm,'2017-10-17') 第

地基中的附加应力计算 —— 平面问题

地基中的附加应力计算 —— 平面问题 平面问题成立的条件 荷载面为长条形，宽度有限，长度无限(一般当长度是宽度的5倍时，可以视为平面问题) 长度方向为均布荷载 平面问题之 线荷载—费拉曼解 其结果仅与R1，p 和 β 有关，表明长度垂直的任何平面上的应力

matlab中常用的矩阵计算

matlab中常用的矩阵计算 matlab常用函数链接 Matlab基本语法与指令链接 参考博客网址 x = A\b; 本博客参考来源超链接 附上一些有用的matlab学习的链接 MATLAB多图合一代码模板链接 Matlab基本语法与指令链接 matlab常用函数链接 MATLAB曲线拟合添加链接描述 M

js加减乘除计算器

js加减乘除计算器 1.首先我们需要规划好计算器的结构 计算器的结构可以用好几种方法写：表格、button、列表…这篇博文我是用表格写的 body里先构建4*6的表格并调整好宽、高、边框等属性，需要使用table、tr、td标签以及td标签的colspan属性进行行的合并操作

车流量计数统计算法

车流量计数统计算法 车流量计数统计算法是目前安防领域重要的应用方向，根据实时或历史视频流，实时统计不同类型的车流量 比如： 小客车，客货两用车，出租车，公交车，中客车，大客车，小货车 等类别的个数实时统计 算法目前的实现为：多目标检测+多目标跟

计算机科学——web

计算机科学——web ## HTML文件 一.HTML简介 1.定义：HTML是一种文本类，由浏览器解释执行的标签语言(是sgml下的一个应用)。 2.特性：具有平台无关性。 二.HTML文件 1.定义：使用HTML语言编写的文件称为HTML文件，也叫web页面(是一种纯文本文件)，扩展名

服务计算-程序包开发读简单配置文件 v1

服务计算-程序包开发，读简单配置文件 v1 熟悉程序包的编写习惯(idioms)和风格(convetions) 熟悉 io 库操作 使用测试驱动的方法 简单 Go 程使用 事件通知 读取该配置文件，并保持监听 本次任务要求实现读取ini文件，观察实验要求后发现，代码主要需要实

在树莓派4B上部署intel NCS2神经网络计算棒与yolov3全过程(附代

在树莓派4B上部署intel NCS2神经网络计算棒与yolov3全过程(附代码)(一) 该项目全部经实测具有可行性，且经过比赛检验 不同时间软件版本可能不同，请按文中使用版本为准 软硬件要求 1.树莓派4B与16G及以上的SD卡和读卡器 2.intel NCS 2 神经网络计算棒 3.Win