codeforces 1114C Trailing Loves (or L'oeufs?) （数论）

https://codeforces.com/contest/1114/problem/C

题意

输入一个n和b。$1\le n\le 10^{18}$，$2\le b\le 10^{12}$。
输出n! 在b进制下，末尾连续0的个数。
例如5 2
120 = 1111000
末尾连续0的个数为3。

题解

$n!=x_k\ast b^{y_k}+x_{k-1}\ast b^{y_{k-1}}+\cdots +x_1\ast b^{y_1}$。
转换成的b进制数为$x_k{x}_{k-1}\dots x_1$。
当$n!mod{b}^k==0$时，$x_k$就为0。
所以只要知道何时不整除$b^k$就可以了。
但是n的范围特别大，不可能算出n的阶乘。
所以考虑质因子分解。
$n!=p_1^{m1}{p}_2^{m2}\dots p_r^{mr}$
$b=p_1^{k1}{p}_2^{k2}\dots p_r^{kr}$
算出b中的每个质因子p在$n!$中的贡献，例如2在6!中出现的次数tmp是6/2+6/4=4。即因子2作为2出现的次数加上作为4出现的次数。以此类推。
然后tmp/k1即使$p^{k1}$出现的次数。
最后取所有的出现次数的最小值便是$b^r$中的r。

代码

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int main() {
	ll n,b;
	cin >> n >> b;
	ll ans = 1ll<<60;
	for(ll i = 2; i <= b; ++i) {
		if(1ll*i*i > b) i = b;
		// cout << i << endl;
		int cnt = 0;
		if(b%i == 0) {
			while(b%i == 0) {
				b /= i;
				cnt++; // p^cnt
			}
			ll mul = 1;
			ll tmp = 0;
			while(n/i >= mul) {
				mul*=i;
				tmp += n/mul;
				
			}
			ans = min(ans, tmp/cnt);
		}
	}
	cout << ans << endl;
	return 0;
}