高精度之阶乘和

之前只会求高精度的阶乘，不知道他们怎么求和，经过别人的知道（特地感谢某博客园的博主的指导），然后知道了这个新的求阶乘和的方法
首先
先看一下之前我是怎么求阶乘的
就是类似于我之前求斐波拉契数列的时候求的方法
然而这种方法，不好来求阶乘和

#include <stdio.h>
#include <string.h>
int a[10010][10010];
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        int i,j;
        int d=0,cnt;
        a[1][0]=1;
        for(i=2;i<=n;i++)
        {
            cnt=0;
            for(j=0;j<=d;j++)
            {
                a[i][j]=a[i-1][j]*i+cnt;
                cnt=a[i][j]/10;
                a[i][j]%=10;
            }
            while(cnt!=0)
            {
                a[i][++d]=cnt%10;
                cnt/=10;
            }
        }
        for(i=d;i>=0;i--)
            printf("%d",a[n][i]);
        printf("\n");
    }
    return 0;
}

因此，为了更好的求阶乘和，下面是一种新的方法，求阶乘

#define INF 0x3f3f3f3f
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<string>
using namespace std;
int a[10010];
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        memset(a,0,sizeof(a));
        a[0]=1;
        int i,j,k,d=0//数位;
        for(k=2;k<=n;k++)
        {
            for(j=0;j<=d;j++)
                a[j]*=k;
            for(j=0;j<=d;j++)
            {
                if(a[j]>=10)
                {
                    a[j+1]+=a[j]/10;
                    a[j]=a[j]%10;
                    if(j==d)
                        d++;
                }
            }
 
        }
        for(i=d;i>=0;i--)
            printf("%d",a[i]);
        printf("\n");
 
    }
    return 0;
}

因此，下面就是理所当然的阶乘和的解决方法啦

#define INF 0x3f3f3f3f
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<string>
using namespace std;
int a[10010];
int sum[10010];
int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        memset(a,0,sizeof(a));
        memset(sum,0,sizeof(sum));
        a[0]=1;//阶乘计数
        sum[0]=1//计数这里开始为1计数;
        int i,j,k,d=0//数位;
        for(k=2; k<=n; k++)
        {
            for(j=0; j<=d; j++)
                a[j]*=k;
            for(j=0; j<=d; j++)
            {
                if(a[j]>=10)
                {
                    a[j+1]+=a[j]/10;
                    a[j]=a[j]%10;
                    if(j==d)
                        d++;
                }
            }
            //上面求玩阶乘，下面就是求和的步骤，精度堆叠的计算方法很像
            for(i=0; i<=d; i++)
                sum[i]+=a[i];
            for(i=0; i<=d; i++)
            {
                if(sum[i]>=10)
                {
                    sum[i+1]+=sum[i]/10;
                    sum[i]=sum[i]%10;
                    if(j==d)
                        d++;
                }
            }
        }
        for(i=d; i>=0; i--)
            printf("%d",sum[i]);
        printf("\n");
    }
    return 0;
}

于是这个问题就解决了啊，代码应该不是很复杂biu~~