zju acm 1007

(转载)

#include<stdio.h> 

int   main() 

 long   k; /*   dummy   summation   variable   */ 
 long   n; /*   number   of   terms   to   sum   */ 
 double   kd; /*   double   version   of   k   */ 
 double   x; /*   the   variable   */ 
 double   sum; /*   used   for   computation   summations   */ 
 double   px,   p1; /*   psi(x)   and   psi(1)   */ 
 double   dpx; /*   value   of   the   first   accelerated   series   at   x   */ 
 double   dp2; /*   value   of   the   first   accelerated   series   at   2   */ 
 double   ddpx; /*   value   of   the   second   accelerates   series   at   x   */ 
 long   i; /*   counter   for   elements   of   table   */ 
 
 p1   =   1.0; 
 dp2   =   0.25; 
 n   =   5000; 

 /*   loop   for   the   2001   table   entries   */ 
 for   (i=0;   i<2001;   i++   ) 
 { 
  x   =   i   *   0.001; 
  
  /*   sum   second   accelerated   series   */ 
  sum   =   0.0; 
  for(   k=1;   k<n;   k++   ) 
  { 
   kd   =   k; 
   sum   +=   1.0/(kd*(kd+1)*(kd+2)*(kd+x)); 
  } 
  ddpx   =   (2.0-x)*sum; 
  
  /*   get   the   value   of   the   first   accelerated   series   */ 
  
  dpx   =   (1.0-x)*(dp2+ddpx); 
  
  /*   get   value   of   original   series   */ 
  
  px   =   p1+dpx; 
  
  printf("%5.3f   %16.12f/n",   x,   px); 
 } 
 return   1; 

=================================================

很遗憾,WA ,不过测试数据都对得上 , 不知道哪里出问题了

阅读更多
文章标签: table 测试
博主设置当前文章不允许评论。

acm zju 额度cn

2009年03月24日 272KB 下载

北大acm 1007

2008年08月21日 744B 下载

DNA排序_北大ACM1007

对象排序

kakingka kakingka

2015-04-07 13:26:47

阅读数:941

浙大zoj月赛解题报告及代码

2011年03月25日 18KB 下载

浙大acm最新模板!

2011年03月25日 643KB 下载

acm新手训练方案新手必备

2011年03月25日 8KB 下载

没有更多推荐了,返回首页

不良信息举报

zju acm 1007

最多只允许输入30个字

加入CSDN,享受更精准的内容推荐,与500万程序员共同成长!
关闭
关闭