ZOJ3621 Factorial Problem in Base K

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3621

题意就是给一个k进制数s，求s!在k进制下末尾有多少个0。

对k进行质因数分解，统计每个质因数及其个数p。然后再计算每个质因数在s!中出现的个数q，则min{q/p}即为最后答案。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <map>
#include <set>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;

long long a[100],b[100],ge;

void init(int kk)
{
    int k=kk;
    ge=0;
    int i,j;

    for (i=2;i<=kk;i++)
    if (k%i==0 && k>0)
    {
        ge++;
        a[ge]=i;

        j=0;
        while (k%i==0 && k>0)
        {
            j++;
            k/=i;
        }

        b[ge]=j;
    }
}

int main()
{
    freopen("a","r",stdin);

    char s[1000];
    long long k,i,a1,len,n;

    while(scanf("%s %lld\n",s,&k)!=EOF)
    {
        len=strlen(s);
        init(k);

        a1=1;
        n=0;
        for(i=len-1;i>=0;i--)
        {
            if(('0'<=s[i])&&(s[i]<='9'))n+=(s[i]-'0')*a1;
            if(('A'<=s[i])&&(s[i]<='Z'))n+=(s[i]-'A'+10)*a1;
            if(('a'<=s[i])&&(s[i]<='z'))n+=(s[i]-'a'+36)*a1;
            a1*=k;
        }

        long long min=-1,m,sum;
        for (i=1;i<=ge;i++)
        {
            sum=0;
            m=n;

            while (m>=a[i])
            {
                sum+=m/a[i];
                m/=a[i];//而不是用个变量不停乘a[i]因为可能溢出
            }

            if (min==-1) min=(sum/b[i]);//初始赋-1然后判断的方法比较好，maxint可能不够
            else if ((sum/b[i])<min) min=(sum/b[i]);
        }

        printf("%lld\n",min);
    }

    return 0;
}