http://acm.hdu.edu.cn/showproblem.php?pid=3641
题意:有N个ai(ai<=100),bi,M=a1^b1*a2^b2*a3^b3…*an^bn ,求最小的 x 使得 x! % M ==0.
思路:先把M的每个素因子的个数求出来,再二分x,判断M的每个素因子是不是在mid!里个数都比M多,是的话就符合了。
知识点:对于n的阶乘中含有的某个质因子x的个数ans求法:
while(n)
{
n/=x;
ans+=n;
}
觉得自己最近刷题状态不太好……另外吐槽一下hdoj,G++过不了的代码C++居然过了
#include<iostream>
#include<algorithm>
#include<string>
#include<map>//int dx[4]={0,0,-1,1};int dy[4]={-1,1,0,0};
#include<set>//int gcd(int a,int b){return b?gcd(b,a%b):a;}
#include<vector>
#include<cmath>
#include<stack>
#include<string.h>
#include<stdlib.h>
#include<cstdio>
#define mod 1e9+7
#define ll long long
using namespace std;
ll x[115];
int main(){
int t;
cin>>t;
while(t--){
memset(x,0,sizeof(x));
ll n,a,b;
cin>>n;
for(int i=0;i<n;++i){
cin>>a>>b;
for(int j=2;j*j<=a;++j){
int cnt=0;
while(a%j==0){
a/=j;
cnt++;
}
x[j]+=cnt*b;
}
if(a>1)
x[a]+=b;
}
ll l=0,r=1000000000000000000,ans; //这个r也很讲究
while(l<=r){
ll mid=(l+r)/2;
ll w,cnt,u=0;
for(int i=1;i<=100;++i){
if(x[i]>0){
w=mid;cnt=0;
while(w){
w/=i;
cnt+=w; //这个while用到了知识点,不是cnt++而是cnt+=w
}
if(cnt<x[i]){
u=1;
break;
}
}
}
if(u==1)
l=mid+1;
else{
r=mid-1;
ans=mid;
}
}
cout<<ans<<endl;
}
return 0;
}