Question

C语言新手提问：怎么计算1+1/3+1*2/(3*5)+1*2*3/(3*5*7)+...+1*2*...*n/(3*5*7*...*(2n+1))？

这是我写出来的，但是运行不正确，求大神指出错误，无限感激
void main()
{ double pi2;
int n;
int a=1;//这里的变量类型应该用？
int b=1; //同上
int i;
scanf("%d",&n);
pi2=1;
for(i=1;i<=n;i++)
{a=a*i;
b=b*(2*i+1);
pi2=pi2+(double)a/(double)b;
printf("%d %d %i\n",a,b,i); //输出a,b,i,在循环时的值
}
printf("%lf\n",pi2);
getch();
}
最后运行结果是这样的，疑似溢出。：

如果有更好的算法的话也麻烦附上吧

Answer 1

#include <stdio.h>

int main()
{  double pi2; 
	int n;   
	double a=1;
    double b=1;
    int i;
    scanf("%d",&n);
    pi2=1;
    for(i=1;i<=n;i++)
    {a*=i;
    b*=(2*i+1);
    pi2+=a/b;
    //printf("%lg  %lg  %i\n",a,b,i);
    }
    printf("%lg\n",pi2*2);
}

追问

谢谢！如果用强制类型转换的话定义a,b 为long int型再转为double可以吗？

追答

定义为 long 话，如果 n 很大，还是会溢出的。

Answer 2

整型位数不够，要用Long类型