题目链接
思路:
遇到N!这样的大数很显然是没办法直接去处理的,题目中告诉我们的已知是:(N!)%P^k==0与(N!)%P^(k+1)!=0。
怎么处理N!是一个很复杂的事情,那我们从P开始考虑,尝试着将P拆成几个质因子的乘积形式,例如12可以拆成2*2*3的形式,那么如何与N!关联起来?不妨可以去看到N!中最多有几个这样的质因子,只要大于这个临界值就肯定不会被整除,于是满足条件,因为被除的数可能大小不一,被拆成的个数也不同,所以要结合所有质因子来判断。
求N!中m的个数:
ll solve(int m, ll n)
{
ll sum=0;
while(n)
{
if(n/m) sum+=n/m;
n=n/m;
if(N==0) break;
}
return sum;
}
主要求解方程:
ll ans=N;
for(int i=2; i<=P; i++)
{
if(cnt[i]) ans=min(ans, solve(i, N)/cnt[i]);
}
完整代码:
#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
using namespace std;
typedef long long ll;
const int maxN=10005;
ll N, P;
bool sushu[maxN];
int cnt[maxN];
void susai()
{
memset(sushu, true, sizeof(sushu));
for(int i=2; i<=maxN-1; i++)
{
if(sushu[i])
{
for(int j=2*i; j<maxN; j+=i) sushu[j]=false;
}
}
}
ll solve(int m, ll n)
{
ll sum=0;
while(n)
{
if(n/m) sum+=n/m;
n=n/m;
if(N==0) break;
}
return sum;
}
int main()
{
susai();
while(scanf("%lld%lld", &N, &P)!=EOF)
{
ll pp=P;
memset(cnt, 0, sizeof(cnt));
for(int i=2; i<=pp; i++)
{
if(sushu[i])
{
if(pp%i==0)
{
while(pp%i==0)
{
cnt[i]++;
pp/=i;
}
}
}
}
ll ans=N;
for(int i=2; i<=P; i++)
{
if(cnt[i]) ans=min(ans, solve(i, N)/cnt[i]);
}
printf("%lld\n", ans);
}
return 0;
}