题目 https://codeforces.com/contest/1114/problem/C
思路
求 n! 转化为 b 进制下末尾有多少个 0.
解题思路:
原题:swjtuOJ 2090
这是一道很好的数论题。
首先转换一下思路:
要求 n! 在 b 进制下有多少个尾 0 就相当于 求 n! % (b^k) == 0 的最大 k。
那么我们现在把 n! 看作一个数 A。问题就是 求 A % (b^k) == 0 的最大 k;
我们知道有素数分解定理: b = p1^a1 * p2^a2 * p3^a3 ...;
那么我们如果可以求得 A 里面 p1^b1 * p2^b2 * p3^b3 ... 的 b1, b2, b3...
那么答案 ans = min(ans, ai/bi ) 了也就是要整除,首先要满足最小的那个能整除。
(1)首先对 b 进行素因子分解,直接暴力(log b), 用一个数组离散化形成该素因子的编号和该素因子的幂的映射 或者 用map存储该素因子的幂,得到所有素因子以及素因子的幂
(2)对于每一个素因子p,计算对应的 A(即 n! ) 中素因子p的幂,两者相除取所有p幂的最小值就是对应的最大整数。
这里求 n! 下 素因子 p 的幂 用累除法,因为存在推论:
n! 下 p 的幂 = [ n/p ] + [ n/(p^2) ] + [ n/(p^3) ] ...
代码
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MAXN = 2e6+1000;
ll N, B;
vector<ll>prime;
map<ll, int>mmp;
void get_p(ll n)
{
ll len = sqrt((double)n);
for(ll i = 2; i <= len; i++){
if(n%i == 0){
prime.push_back(i);
while(n%i == 0){
n/=i;
mmp[i]++;
}
}
}
if(n != 1){
prime.push_back(n);
mmp[n]++;
}
}
ll calc(ll n, ll p)
{
ll res = 0;
while(n){
res += n/p;
n/=p;
}
return res;
}
int main()
{
ll sum = 1ll;
scanf("%lld %lld", &N, &B);
get_p(B);
ll ans = (1ll<<62);
for(ll i = 0; i < prime.size(); i++){
ans = min(ans, calc(N, prime[i])/mmp[prime[i]]);
}
printf("%I64d\n", ans);
return 0;
}