由于输出的数太大，定义int范围只在-32768到32767，跟本不够用，float就可以了，它范围是-3.4乘以10的-38次方到3.4乘以10的38次方

刚好写过一个这样的.改N就可以了.基本可以秒出.

#include <stdio.h>   
#include <stdlib.h>   
 
#define N 20    //要计算的N   
#define P 10000  
long s[N]={1,1},a=1,b=0;  
 
void factorial(long s[],long n)  
{  
for(;a<=*s;)   
    {   
     b+=s[a]*n;   
     s[a]=b%P;   
     b/=P;   
     if((*s==a)&&b)   
       (*s)++;   
     ++a;   
    }  
b=0,a=1;  
}   
 
void add(long sum[],long s[])  
{  
 int i;
for(i=1;i<=s[0];++i)  
    if((sum[i]+=s[i])>9999) {sum[i]%=10000,sum[i+1]++;}  
}  
 
int main()  
{   
int i=2;  
long sum[N]={1,1};  
for(;i<=N;factorial(s,i),add(sum,s),++i);  
printf("%d",sum[s[0]]);  
for(i=s[0]-1;i>0;printf("%04d",sum[i]),--i);   
putchar('\n');  
system("pause");   
return 0;   
}

呵呵，你们再试试做100！试试，用链表做

100!

/*求1!+2!+3!+4!+......+20!的值*/
#include<stdio.h>
main()
{
int i,j;
float k,sum;
for(i=1;i<=20;i++)
{
k=1;
for(j=1;j<=i;j++)
{
k=k*j;
}
sum+=k;
}
printf("1!+2!+3!+4!+......+20!=%f\n",sum);
getch();
}

结果：458284170950269045000000000000.000000

用递归函数也可以的，如下：
#include "stdio.h"
double fun(double n)
{
    double f;
    if(n<0)printf("Data Error!");
    else if(n==0||n==1)f=1;
    else
    f=fun(n-1)*n;
    return(f);
     
    }
void main()
{
   
  double k=0,i=1;
  do{
      k+=fun(i);
      i++;
  }while(i<=20);
  printf("%f",k);  
}

貌似sum变量要设为long型吧~~

唉，要有针对性呀！！！

#include<stdio.h>
double f(int n)
{
    if(n==1)return 1;
    else
       return (n*f(n-1));
}
void main()
{
    double  sum=0;
    int i;
     for(i=1;i<21;i++)
        sum+=f(i);
       printf("%f",sum);
       getch();
}这个数太大，可以用以上程序，但结果是实型的，你自己去看看吧！

实数的精度只有16位吧，也不精确的


ms 32楼可行。。。？