杭电1042

题意：计算N!（0<N<=10000)

Analyse:

大数乘小数问题。还是用比较原始的模拟。这道题调试了好一段时间，真受不鸟，高人指点才发现是个低级错误，d的值没有先用temp储存而错了。

#include<stdio.h>
#define UNIT 100000
typedef struct
{
    __int64 key[10000];//每个key存五位
    int digit;//每个key算作一个digit
}bigint;
main()
{
    bigint carry;//用来保存相乘后的进位
    bigint num;//每次大数相乘的结果
    int n;
    __int64 d,temp;//d为大数相加的进位
    int i,j;//j为大数key的索引，i为因数

    while(scanf("%d",&n)!=EOF)
    {
        num.key[0]=1;
        num.digit=1;
        for(i=1;i<=n;i++)
        {
            //设定carry的位数为1
            carry.digit=1;
            //进到"0"位的有0
            carry.key[0]=0;
            //相乘，进位存到carry[j+1]
            for(j=0;j<num.digit;j++)
            {
                carry.key[j+1] = (i*num.key[j])/UNIT;
                //carry的位数增一
                carry.digit++;
                num.key[j] = (i*num.key[j])%UNIT;
            }
            //初始化进位为0
            d=0;
            //carry各个位加上去num的各个位
            for(j=0;j<num.digit && j<carry.digit;j++)
            {
                temp=d;
                d=(carry.key[j]+num.key[j]+d)/UNIT;
                num.key[j]=(carry.key[j]+num.key[j]+temp)%UNIT;
            }
            //假若num比carry位数多，而且有进位就继续加
            for(;j<num.digit && d;j++)
            {
                temp=d;
                d=(d+num.key[j])/UNIT;
                num.key[j]=(temp+num.key[j])%UNIT;
            }
            //假若carry比num位数多
            for(;j<carry.digit;j++)
            {
                temp=d;
                d=(carry.key[j]+d)/UNIT;
                num.key[j]=(carry.key[j]+temp)%UNIT;
                if(num.key[j]!=0)
                    num.digit++;
            }
            if(d>0)
            {
                num.key[j]=d;
                num.digit++;
            }
        }
        for(j=num.digit-1;j>=0;j--)
            printf(j==num.digit-1?"%I64d":"%05I64d",num.key[j]);
        printf("\n");
    }
}