hdu3641 Treasure Hunting（灵活题，阶乘的素因子个数求法）

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