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