圆周率（一）

2022-08-07 10:30:14 阅读：134 来源： 互联网

标签：__ int 圆周率 quadmath 1e8 类型 float128

使用float128类型以提高精度。
float128不是标准库类型而是实现定义的，包含在quadmath.h库中。

(1)该类型的使用必须包含quadmath.h头文件，并且要在编译选项中指定：-std=gnu++11 -lquadmath；
(2)该类型的字节大小为16位；（在本机使用sizeof测试的结果）
(3)该类型不能直接使用cout输出，而需要使用quadmath_snprintf函数（用法同sprintf，对于float128的格式符用%Qf），来转换为字符串后输出；

使用尼拉坎塔级数，计算得到的结果

#include<stdio.h>
#include<quadmath.h>

const int N=1e8;

__float128 f(__float128 i){
  __float128 res=(__float128)1/(__float128)(i*(i+1)*(2*i+1));
  return res;
}
int main(){
  __float128 ans=3;
  for(int i=1;i<=N;i++){
    if(i%2==1)ans+=f(i);
    else ans-=f(i);
  }
  char s[100];
  quadmath_snprintf(s,sizeof(s),"%.40Qf",ans);
  puts(s);
}

当N=1e7时，精确到21位；
当N=1e8时，精确到24位；

3.1415926535897932384626431332797386533518

当N=1e9甚至更加以上的时候，耗时更长，此处不列举结果。

注：
如果只使用double类型，循环次数即使在1e8，位数也只能精确到13-15位左右，这是因为double本身的限制。

标签：__,int,圆周率,quadmath,1e8,类型,float128
来源： https://www.cnblogs.com/jisuanjizhishizatan/p/16558566.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

圆周率（一）