n的阶乘在m进制下结尾0的个数

推导

首先举几个例子：

• $8$在二进制下可以表示为$1\times 2^3$

• $400$在十进制下可以表示为$4\times 10^2$

• $12$在$6$进制下表示为$2\times 6^1$

因此不难发现如果要求一个数$n$在$m$下结尾$0$的个数，实际上就是求出满足$n=t\ast m^k$最大的$k$。

但是如果是$n$的阶乘呢？需要对$n!$质因数分解吗？先看下面⬇️

n!的所有质因数

阶乘数很大，显然没有办法求出阶乘然后质因数分解。考虑$n!$的所有质因数都不会超过$n$，显然可用的质因数就是$n$以内的所有因数。考虑如何求出$n!$所有的质因数的个数。先分析一个简单的问题，求出$n$以内质数$p$的倍数的个数，即$\lfloor \frac{n}{p}\rfloor$，但是$n$的阶乘是前$n$个数相乘，而且对于某个质数，我们只知道它的倍数是没有用的，需要知道每个数上该质数的幂次是多少。而$n$以内$p$的倍数设分解为$p^k\ast q$，显然$k=1,2,3,...$。那么$\lfloor \frac{n}{p}\rfloor$求出的只是$k=1$的数的个数，还有$k=2,3,...$的$p$的倍数，这时我们如果再将上述结果除以$p$，即$\lfloor \frac{\lfloor \frac{n}{p}\rfloor }{p}\rfloor$，也相当于是$\lfloor \frac{n}{p^2}\rfloor$，这求出的是$k=2$时的倍数个数，以此类推。将所有结果相加，得到的恰好是$n!$中质数$p$的个数。

例：求在$8!$中$2$的个数:

• 首先我们先计算出$2$的倍数的个数：$8/2=4$

• 其次我们计算出$4$的倍数的个数：$8/4=2$（上面一个式子求出了第一层，现在求第二层)

• 最后我们解出$8$的倍数的个数：$8/8=1$

$4+2+1=7$，一共$7$个$2$出现了

即：$cnt(x)=[n/(x^1)]+[n/(x^2)]+[n/(x^3)]+...(直到x的次方大于n)$

typedef long long ll
ll cal(int n, int x) {
    if(n==0 || n==1) return 0;
    ll ans = 0;
    while(n) {
        ans += n / x;
        n /= x;
    }
    return ans;
}

总结

因为任何一个数都可以分解为若干质因数的乘积，再从质因数分解的角度考虑，在保证有解的情况下可以按以下方式推导：

设$n!=t\ast m^k=t\ast (p_1^{q_1},p_1^{q_2},...,,p_n^{q_n}{)}^k$，可以求出$n!$中$p_1,p_2,...,p_n$的出现次数，假设分别为$x_1,x_2,...,x_n$，那么有$x_1=k_1{q}_1，x_2=k_2{q}_2，...，x_n=k_n{q}_n$。显然最终的答案 k = m i n { k 1 , k 2 , . . . , k n } k = min\{ k_1,k_2,...,k_n \} k=min{k1​,k2​,...,kn​}。

简单模板

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
typedef long long ll;
const int maxn = 1e4 + 10;
int prime[maxn];     //存所有质数
int a[maxn], b[maxn];  //a存所求数的所有质因子，b存每个质因子的个数
bool isprime[maxn];
int num, cnt;  //num是maxn范围内素数个数

void init() {
    memset(isprime, true, sizeof(isprime));
    for(int i = 2; i <= maxn; i++) {
        if(isprime[i]) {
            prime[num++] = i;
            for(int j = i * i; j <= maxn; j += i)
                isprime[j] = false;
        }
    }
}

void decom(int n) {
    cnt = 0;       //cnt是所有质因数的种数
    memset(b, 0, sizeof b);
    for(int i = 0; i < num; i++) {
        if(prime[i] > n) break;
        int d = n;
        while(d % prime[i] == 0) {
            a[cnt] = prime[i];
            b[cnt]++;
            d /= prime[i];
        }
        if(d < n) cnt++;
    }
}

ll cal(ll n, ll x) {
    if(n == 0 || n == 1) return 0;
    ll ans = 0;
    while(n) {
        ans += n / x;
        n /= x;
    }
    return ans;
}

ll solve(ll n, int m) {
    decom(m);
    ll ans = INF;
    for (int i = 0; i < cnt; i++) {
        ll t = cal(n, a[i]);
        ans = min(ans, t / b[i]);
    }
    return ans;
}

int main(){
    ll n;
    int m;
    cin >> n >> m;
    init();
    cout << solve(n, m) << endl;
    return 0;
}