n!的素因子分解

#include<iostream>
using namespace std;
#define N 1000
bool pr[N+1]={true,true};
int p[N>>1],plen=0;
int main()
{
    int i,j,n,tmp;
    for(i=2;i*i<=N;i++)if(!pr[i])for(j=i;j*i<=N;j++)pr[i*j]=true;
    for(i=2;i<=N;i++)if(!pr[i])p[plen++]=i;
    while(cin>>n)
    {
        for(i=0;i<plen&&p[i]<=n;i++)
        {
            tmp=n;
            while(tmp>=p[i])
            {
                for(j=1;j<=tmp/p[i];j++)cout<<p[i]<<' ';
                tmp/=p[i];
            }
        }
        cout<<endl;
    }
    return 0;
}