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;
}
-
sum的表达式
实际上表达的是f(x)-f(1),精度见下。 -
for中k<10000
题目要求error less than 0.5e-12,因为分母中有k的立方,所以取了10000。(公式求的不一样,会出现不一样的k) -
对sum的精度补充
对导出的公式放缩(k,k+1,k+2均换为k),得到的求和式再用hint里的公式,求积分,解决。精度补充完成,最后记得加上f(1)。
网友评论