USACO 3.2 Factorials(高精度乘法）

水过啊。。直接模拟O（n^2)的算法，有O（n*logn)的算法，研究。

浙大关于求N！最后非0位的阶乘模版。2012.11.22

#include <cstdio>
#include <cstring>
#include <string>
#include <cmath>
#include <queue>
using namespace std;
int lastdigit(char* buf)
{
    const int mod[20] = {1,1,2,6,4,2,2,4,2,8,4,4,8,4,6,8,8,6,8,2};
    int len = strlen(buf),a[10000],i,c,ret = 1;
    if(len == 1)
    return mod[buf[0]-'0'];
    for(i = 0;i < len;i ++)
    {
        a[i] = buf[len-1-i]-'0';
    }
    for(;len;len-=!a[len-1])
    {
        ret = ret*mod[a[1]%2*10+a[0]]%5;//注意模版敲对了。。
        for(c = 0,i = len-1;i >= 0;i --)
        {
            c = c*10+a[i],a[i]=c/5,c%=5;
        }
    }
    return ret+ret%2*5;
}
int main()
{
    char str[501];
    while(scanf("%s",str)!=EOF)
    {
        printf("%d\n",lastdigit(str));
    }
    return 0;
}

 

/*
   ID: cuizhe
   LANG: C++
   TASK: fact4
*/
#include<stdio.h>
#include<string.h>
int p[300000];
int main()
{
    int a,b,i,j;
    freopen("fact4.in","r",stdin);
    freopen("fact4.out","w",stdout);
    scanf("%d",&a);
    p[0]=1;
    j=0;
    for(b=1; b<=a; b++)
    {
        for(i=0; i<=j; i++)
        {
            p[i]=p[i]*b;
        }
        for(i=0; i<=j; i++)
        {
            if(p[i]>9)
            {
                p[i+1]+=p[i]/10;
                p[i]=p[i]%10;
                if(i+1>j)
                    j=i+1;
            }
        }
    }
    for(i = 0; i <= j; i++)
    {
        if(p[i] != 0)
        {
            printf("%d\n",p[i]);
            break;
        }
    }
    return 0;
}