关闭

【N!】【大数】

433人阅读 评论(0) 收藏 举报
分类:

对于求大整数的阶乘,可以采用分段相乘的方法,其理论基础是加法的分配律,乘法的分配律。

for exmaple:123456789*123=123*(123*10^5+56789)=123*123*10^5+123*56789;

如果我们用一个数组result[0]=56789,result[1]=123,那么123456789*123就是123*result[0]^6和result[1]*123+result[0]/10^6构成的一个数,另result[0]=123*result[0]^6,result[1]=result[1]*123+result[0]/10^6,则这个数就是result[1]result[0]。

#include<stdio.h>

int main(){
    int n,j,i,t,p;
    int result[16000];
    while(scanf("%d",&n)!=EOF){
        t=p=0;
        result[0]=1;                //t记录result中所存数的个数
        for(j=1; j<=n; j++){                 //n!=1*2*3*4...*n;
            for(i=0; i<=t; i++){
                result[i]=result[i]*j+p;    //如果sum>99999不成立,则不需进位,进位为0
                p=result[i]/100000;          //进位
                result[i]%=100000;        //数组sum中每位存5位数字
            }
            if(p){         //t在最高位有进位的情况下++,最高位存进位,并让进位为0
                t++;
                result[t]=p;
                p=0;
            }            // 如果进位不为0, 则下一次乘的时候将加上该进位,显然不合理
        }
        printf("%d",result[t]);//先把最高位输出,因为最高位没有前导0
        for(i=t-1; i>=0; i--){
            printf("d",result[i]);//s[i]不足5位,说明需补0
        }
        //也说明result[i]在对100000取余的 时候去掉了前导0,所以这里要加上
        printf("\n");
    }
    return 0;
}

//某位大牛的代码,自己小小改了一下。一时找不到出处了,有空一定给补~

0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:10455次
    • 积分:289
    • 等级:
    • 排名:千里之外
    • 原创:18篇
    • 转载:3篇
    • 译文:0篇
    • 评论:0条
    文章分类
    文章存档