poj 2649 素因子分解 ------ n！的素因子

n!的素因子分解中的素数p的个数为

n/p+n/p^2+n/p^3+n/p^4+......

原理：n/p 计算出了1～n中含有素因子p的个数

            n/p^2计算出了1～n中含有因子p^2的个数

            .........

           求和就是n！中素数p的个数

注意：n或m为0时，都不能整除

                                                               
//poj2649
#include<iostream>
#include<vector>
#include<cstring>
using namespace std;
const int N=1000000;
bool prime[N+1];
long long p[N+1];
int k=0;
void getprime()
{
    memset(prime,0,sizeof(prime));
    prime[0]=1;
    prime[1]=1;
    for(int i=2;i<=N;i++)
    {
        if(!prime[i]) p[k++]=i;
        for(int j=0;j<k;j++)
        {
            if(p[j]*i>N) break;
            prime[p[j]*i]=1;
            if(i%p[j]==0) break;
        }
    }
}
struct factor//a：素因子，k素因子的个数
{
    int a,k;
};
vector<factor> v;
void getfen(int t)//对t分解质因子，结果存在v中
{
    v.clear();
    for(int i=0;p[i]*p[i]<=t;i++)
    {
        if(t%p[i]==0)
        {
            factor a;
            a.a=p[i];
            a.k=0;
            while(t%p[i]==0)
            {
                a.k++;
                t=t/p[i];
            }
            v.push_back(a);
        }
    }
    if(t&&t!=1)
    {
        factor a;
        a.a=t;
        a.k=1;
        v.push_back(a);
    }
}
int main()
{
    getprime();
    int n,m;
    while(cin>>m>>n)
    {
        getfen(n);
        int flag=1;
        for(int i=0;i<v.size();i++)
        {
            int kk=0,l=v[i].a;//kk代表m中素因子v[i].a的个数
            while(m/l)
            {
                kk+=m/l;
                l=l*v[i].a;
            }
            if(kk<v[i].k)
            {
                flag=0;
                break;
            }
        }
        if(m==0||n==0) flag=0;
        if(flag)
            cout<<n<<" divides "<<m<<"!"<<endl;
        else cout<<n<<" does not divide "<<m<<"!"<<endl;
    }
}