Question

C语言；已知E=1+1/1!+1/2!+1/3!+.......,求E的近似值，要求最后一项的值小于10的负6次方。

Answer 1

#include <iostream>
#include <assert.h>
template<int N>
struct fact
{
enum
{
value = N * fact<N - 1>::value
};
};

template<>
struct fact<1>
{
enum
{
value = 1
};
};

template<>
struct fact<0>
{
enum
{
value = 1
};
};

template<int l, int r>
struct greater_
{
enum
{
value = l > r
};
};

template<int i, int miss_divisor, int end>
struct e_loop
{
static const double value = 1.0f / fact<i>::value + e_loop< i + 1, miss_divisor, greater_< fact<i>::value, miss_divisor>::value >::value;
};

template<int i, int miss_divisor>
struct e_loop<i, miss_divisor, 1>
{
static const double value = 1.0f / fact<i>::value;
};

template<int miss_divisor>
struct e
{
static const double value = e_loop<0, miss_divisor, 0>::value;
};

int main(int argc, char* argv[])
{
std::cout << (e<1000000>::value) << std::endl;
//(e<1000000>::value)由编译器运算 exe里运行时间为0哦 由于模板参数不能用float 1e-6 用 1 / 1e +6表示
// 模板参数填除数就可以了所以 这里填1000000 （1e + 6 == 1000000, 1e - 6 = 1 / 1000000）

//以下是验证
double e2 = 1.0 + 1.0
+ 1.0 / fact<2>::value
+ 1.0 / fact<3>::value
+ 1.0 / fact<4>::value
+ 1.0 / fact<5>::value
+ 1.0 / fact<6>::value
+ 1.0 / fact<7>::value
+ 1.0 / fact<8>::value
+ 1.0 / fact<9>::value
+ 1.0 / fact<10>::value;

std::cout << e2 << std::endl;

assert(!(1.0 / fact<9>::value < 1e-6));
assert(1.0 / fact<10>::value < 1e-6);
return 0;
}

真的只要过程？

Answer 2

楼主是想知道思路还是？

追问

只要过程，思路我会自己琢磨\(^o^)/~

追答

等下 我写个大概给你，不过要加点分才行

int jiec(int m) //这个是求阶乘
{
for(int i=0;i1e-6;x++)
{
a=jiec(x);
b=(float)1/a;
E+=b;
}
printf("%f\n",E);
}

Answer 3

#include "stdio.h"
void main(void){
double E,x;
int i,j;
for(E=0,i=1;;i++){
for(j=2,x=1;j<i;j++)
x*=j;
E+=(x=1/x);
if(x<1E-6) break;
}
printf("E=%f.\n",E);
}