Again Prime? No Time.（数论）

传送门

题意：问 n!%m^k=0, 求最大的k
分析：先对m进行素因子分解，然后将n!分成单个数k(1<=k<=n),看k与m有多少相同的素因子
     如果n!的每个与m共有的素因子的个数都比对应m的多，就说明 n!%m^k=0

#include <set>
#include <map>
#include <list>
#include <cmath>
#include <stack>
#include <queue>
#include <string>
#include <bitset>
#include <vector>
#include<cstring>
#include <stdio.h>
#include <iostream>
#include <algorithm>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
inline int read(){int s=0,w=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9') s=s*10+ch-'0',ch=getchar();return s*w;}
const int maxn = 1e5+5;
int tot;
struct node{
    int num,cnt;
}pi[maxn];
void Pi(int n)
{
    for(int i=2;i*i<=n;i++){
        bool flag=false;
        int cnt=0;
        while (n%i==0){
            n/=i;
            if(!flag) pi[++tot].num=i,flag=true;
            cnt++;
        }
        if(flag) pi[tot].cnt=cnt;
    }
    if(n>1) pi[++tot]={n,1};
}
int main()
{
    int t=read(),cas=0;
    while (t--)
    {
        memset(pi,0,sizeof(pi));
        int m=read(),n=read();
        int ans=0;
        tot=0;
        Pi(m);
        int a[1005]={0},l=0;
        for(int i=n;i>=1;i--)
        {
            int k=i;
            for(int j=1;j<=tot;j++)
                while (k%pi[j].num==0) k/=pi[j].num,a[j]++;

        }
        bool flag=false;
        int Min=INF;
        for(int i=1;i<=tot;i++)
            if(a[i]<pi[i].cnt!=0) flag=true;

        if(!flag) {
            for(int i=1;i<=tot;i++)
                Min=min(Min,a[i]/pi[i].cnt);
        }
        ans=Min;
        printf("Case %d:\n",++cas);
        if(ans==INF) puts("Impossible to divide");
        else printf("%d\n",ans);
    }
    return 0;
}
/*
100
5000 10000
Case 1:
624

6 10
Case 2:
4

10 20
Case 3:
4

12 30
Case 4:
13

 */