传送门
题意:问 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;
}