ZOJ - 3621——Factorial Problem in Base K

题意：给出一个k进制的数s，求s!里一共有多少个0。
思路：把k分解质因数，s的阶乘里每个0都代表着因式中存在一组k的质因数。计算s!包含k的质因数各多少个，其中最小的一组就是答案。对于k的每个质因数m，可以这样想，1*2……*s中，每隔m个数就有一个可以被m整除，所以将s/m就是可以被m整除一次的数的个数。s/m/m就是可以被m整除两次的数的个数……以此类推。
代码如下：

#include<cstdio>
#include<cstdlib>
#include<math.h>
#include<iostream>
#include<cstring>
#include<map>
#include<set>
#include<list>
#include<stack>
#include<algorithm>
#include<queue>
#include<vector>
#include<time.h>
#include<iomanip>
using namespace std;

typedef unsigned long long ull;

ull a[70];
ull b[70];
int factor(ull n){
    ull temp,i,now;
    int tot = 0;
    temp = (ull)((double)sqrt(n)+1);
    now = n;
    for(i = 2; i<= temp; ++i)if(now%i == 0){
        a[++tot] = i;
        b[tot] = 0;
        while(now%i==0){
            ++b[tot];
            now /= i;
        }
    }
    if(now != 1){
        a[++tot] = now;
        b[tot] = 1;
    }
    return tot;
}

char s[100];
int k;

ull fpow(ull x,int n)
{
    ull ret = 1;
    while(n){
        if(n&1)ret *= x;
        x*=x;
        n>>=1;
    }
    return ret;
}

ull cnt()
{
    ull ret = 0;
    int len = strlen(s);
    for(int i = 0;i < len;++i){
        ull tmp;
        if(s[i] >='0' && s[i] <= '9'){
            tmp = s[i] - '0';
        } else if(s[i]>='A'&&s[i]<='Z'){
            tmp = s[i] - 'A' + 10;
        } else if(s[i]>='a' && s[i] <= 'z'){
            tmp = s[i] - 'a' + 36;
        }
        ret += tmp * fpow(k,len-i-1);
    }
    return ret;
}

ull sav[70];

int main()
{
//    freopen("data.txt","r",stdin);
    while(scanf("%s%d",s,&k)!=EOF){
        ull num = cnt();
        int tot = factor(k);
        memset(sav,0,sizeof(sav));
        ull ans = 0;
        for(int i =1;i<=tot;++i){
            ull tmp = num;
            while(tmp){
                sav[i] += tmp/a[i];
                tmp /=a[i];
            }
            ans = max(ans,sav[i]);
        }

        for(int i = 1;i<=tot;++i){
            ans = min(ans,sav[i]/b[i]);
        }
        cout<<ans<<endl;
    }
    return 0;
}