BitComet 旗下网站

转到日志
相关贴吧:
acm of zju

zju1007

楼主 发表于:2008-06-09 11:19:23 [回复]

/*2008.6.5 07:55:08 Accepted 1007 C++ 00:00.62 392K 天将降大任于我
Numerical Summation of a Series http://acm.zju.edu.cn/show_problem.php?pid=1007
  在n->inf的时候,1/(k * (k + 1) * (k + x))和1/(k * k * (k + 1))相等,当n=10000时,两者差已经
小于10^-15,故大于10000的部分我们可以不算了,将hint中给的那个式子积分,
当k->inf时,原式=0.5×(x^(-2)) 于是就得到了程序。

 */
#include<stdio.h>
int main()
{
 double a,s,k;
 for(a=0.000;a<=2.000;a+=0.001)
 {
  s=0;
  for(k=1;k<10000;k++)
  s+=(1-a)/(k*(k+1)*(k+a));   //经推导可得 f(x)-f(1)=(1-x)/(k*(k+1)*(k+x))
  s+=(1-a)/(10000*10000*2)+1.0;   //由于f(1)=1.0,故还应加上1.0,不过公式是怎么来的,我还是没弄明白
  printf("%5.3f %16.12f\n",a,s);
 }
 return 0;
}

心难泰,世风坏,旧时正气今何在?正义寡,人情薄,闻道虽多,茅塞不开。怪!怪!怪! 空等待,几多载,冲出重围人心快!暴雨打,狂风袭,任他折磨,此志难改。耐!耐!耐!

 

您现在还没有登录,请在登录后发贴