高精度之階乘和

Maosghoul發表於2019-05-11

之前只會求高精度的階乘，不知道他們怎麼求和，經過別人的知道（特地感謝某部落格園的博主的指導），然後知道了這個新的求階乘和的方法
首先
先看一下之前我是怎麼求階乘的
就是類似於我之前求斐波拉契數列的時候求的方法
然而這種方法，不好來求階乘和

#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("
");
    }
    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("
");
 
    }
    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("
");
    }
    return 0;
}

於是這個問題就解決了啊，程式碼應該不是很複雜biu~~

相關文章