ZOJ-1007

来源:互联网 发布:滴胶印花t恤 知乎 编辑:程序博客网 时间:2024/06/06 03:35

ZOJ-1007-Numerical Summation of a Series

题解中的代码,解释一下数学部分。

顺便感谢做题解的大神。

#include<stdio.h>int main(){    double sum,a,k;    for(a=0.000;a<=2.000;a=a+0.001){        sum=0;for(k=1;k< 10000;k++){sum=sum+(1-a)/(k*(k+1)*(k+a));}sum=sum+(1-a)/(2*10000*10000)+1.0; printf("%5.3f %16.12f\n", a, sum );}return 0;}
1. sum的表达式

实际上表达的是f(x)-f(1),精度见下。


2. for中k<10000

题目要求error less than 0.5e-12,因为分母中有k的立方,所以取了10000。(公式求的不一样,会出现不一样的k)


3. 对sum的精度补充

对导出的公式放缩(k,k+1,k+2均换为k),得到的求和式再用hint里的公式,求积分,解决。精度补充完成,最后记得加上f(1)。

0 0
原创粉丝点击