【高精】阶乘和

咳咳，这是一道高精的题，好久以前十分弱智的题（弱智吗？？？），只是看起来好难，于是一直放弃。。。

可是今天看见了，发现1.6（openjudge）里面竟然就这一个没打钩，于是只好硬着头皮做了。。

先贴上题目：

15:阶乘和

总时间限制: 1000ms 内存限制: 65536kB
描述：

用高精度计算出S=1!+2!+3!+…+n!（n≤50）

其中“!”表示阶乘，例如：5!=5*4*3*2*1。

输入正整数N，输出计算结果S。
输入一个正整数N。输出计算结果S。
样例输入

5

样例输出

153

于是我先潇洒地打完框架，然后先算算50的阶乘有多少位，一算，哇！才六十多位，我就定义了一百个数组（其实是105个==），然后为了节省空间时间，把每一个存进四位，也就是/10000。

Then，我照着上面这道题（求阶乘）的思路，打出了阶乘这部分框架，然后用高精加把他们加起来，就是这俩货：

int a[105]/*阶乘和*/,k[105]/*阶乘*/;

为了进一步节省我就用这个阶乘：

    for(int i=2;i<=n;i++)
        for(int j=100;j>=aa;j--)
        {
            k[j]*=i;
            k[j-1]+=k[j]/10;
            k[j]%=10;
        }
        加起来。。。。。。

是不是很快！
然后我高高兴兴打完代码，才试验几个数就交了，发现是wrong answer！！

#include
int a[105]/*阶乘和*/,k[105]/*阶乘*/,kk=100/*k当前位置*/,aa=100/*a当前位置*/;
void wzk()//计算k位置
{
    int i=1;
    while(!k[i]) i++;
    kk=i;
}
void wza()//计算a位置
{
    int i=1;
    while(!a[i]) i++;
    aa=i;
}
void cc()//把阶乘累加
{
    for(int i=100;i>=kk;i--)//不到a的最高位，而到k的最高位，是因为一个数的阶乘大于等于比它小的阶乘之和
    {
        a[i]+=k[i];
        a[i-1]+=a[i]/10000;
        a[i]%=10000;
    }
    wza();
}
int main()
{
    int n;
    scanf("%d",&n);
    a[100]=1;//从100开始存，ps:50的阶乘只有65位
    k[100]=1;
    for(int i=2;i<=n;i++)
    {
        for(int j=100;j>=kk;j--)
        {
            k[j]*=i;
            k[j-1]+=k[j]/10000;
            k[j]%=10000;
        }
        wzk();
        cc();
    }
    for(int i=100;i>=aa;i--)
        printf("%3d",a[i]);
    return 0;
}

这是半个小时打出来的。。。。
后来发现好多错误：
1.最后的输出，当输出不足三位时会用空格补位。

2.中间的乘法，不能边乘边进位，不然会发现进位的数也跟着乘了的尴尬。。。

这些我改了好几遍才改出来，，，，最后交上去对了，已经用了一个多小时了。。。。。。。

下面是你们最爱的代码时间：：：：

#include
int a[105]/*阶乘和*/,k[105]/*阶乘*/,kk=100/*k当前位置*/,aa=100/*a当前位置*/;
void wzk()//计算k位置
{
    int i=1;
    while(!k[i]) i++;
    kk=i;
}
void wza()//计算a位置
{
    int i=1;
    while(!a[i]) i++;
    aa=i;
}
void cc()//把阶乘累加
{
    for(int i=100;i>=kk;i--)//加法不用加完再进位
    {
        a[i]+=k[i];
        a[i-1]+=a[i]/10;
        a[i]%=10;
    }
    wza();
}
int main()
{
    int n;
    scanf("%d",&n);
    a[100]=1;//从100开始存，ps:50的阶乘只有65位
    k[100]=1;
    for(int i=2;i<=n;i++)
    {
        for(int j=100;j>=aa;j--)
            k[j]*=i;
        for(int j=100;j>=aa;j--)//坑！一定要先乘完再进位。。。
        {
            k[j-1]+=k[j]/10;
            k[j]%=10;
        }
        wzk();
        cc();
    }
    for(int i=aa;i<=100;i++)
        printf("%d",a[i]);
    return 0;
}

不用谢！！！