uva 10139 求n的阶乘能否整除m

n跟m的范围都是0-2^31，只能用质因数分解把m分成若干个质数相乘的形式，然后对每个质因数的数量判断n的阶乘中这个质因数的数量是不是足够，不够的话就不行。

最后m可能是一个质数，没有被2到根号m中的质数分解，于是分解过后要是m还要大于n，这也不能分解，但是当n==0，m==1时是可以的，这种情况要特判。

#include<bitset>
#include<map>
#include<vector>
#include<cstdio>
#include<iostream>
#include<cstring>
#include<string>
#include<algorithm>
#include<cmath>
#include<stack>
#include<queue>
#include<set>
#define inf 0x3f3f3f3f
#define mem(a,x) memset(a,x,sizeof(a))

using namespace std;

typedef long long ll;
typedef pair<int,int> pii;

inline int in()
{
int res=0;
char c;
while((c=getchar())<'0' || c>'9');
while(c>='0' && c<='9')res=res*10+c-'0',c=getchar();
return res;
}
const int N=500000;
bitset <N>vis;
int prime[200000];

int getNum(int i,int m)
{
int res=0;
while(m)
{
res+=m/prime[i];
m/=prime[i];
}
return res;
}

int main()
{
int p=-1;
for(int i=2;i<N;i++)
{
if(!vis[i]) prime[++p]=i;
for(int j=0;j<=p && prime[j]*i<N;j++)
{
vis[prime[j]*i]=1;
if(i%prime[j]==0)break;
}
}

int n,m;
while(~scanf("%d%d",&n,&m))
{
if(n>=m)
{
printf("%d divides %d!\n",m,n);
continue;
}
bool can=1;
int record=m;
for(int i=0;i<=p && can && prime[i]*prime[i]<=m;i++)
{
int tmp=0;
while(m%prime[i]==0)//质因数分解
{
tmp++;            //当前质因数的个数
m/=prime[i];
}
if(tmp){
int tmp2=getNum(i,n); //n！中prime[i]的个数
if(tmp2<tmp) can=0;

}
}
if(n!=0 && m>n)can=0;///
if(can) printf("%d divides %d!\n",record,n);
else printf("%d does not divide %d!\n",record,n);

}
return 0;
}