ZOJ-1007

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)。