因子【Wannafly挑战赛25 A】

题目链接

  思路：

遇到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;
}