素数（线性）筛 + 合数分解

//素数筛 + 合数分解
// O(n) 
#include <bits/stdc++.h>
using namespace std;
 
const int MAXN=10000;
 int prime[MAXN+1];
 int is_prime[MAXN+1];
 void getPrime() {
	 memset(is_prime,0,sizeof(is_prime));
	 for(int i=2; i<=MAXN; i++) {
		 if(!is_prime[i])prime[++prime[0]]=i;//第prime[0]个素数是i
		 for(int j=1; j<=prime[0]&&prime[j]<=MAXN/i; j++) {
			 is_prime[prime[j]*i]=1;
			 if(i%prime[j]==0) break;//保证每个合数只会被它的最小质数筛去
//prime[]数组中的素数是递增的,当i能整除prime[j]，那么i*prime[j+1]这
//个合数肯定被prime[j]乘以某个数筛掉。   因为i中含有prime[j],
//prime[j]比prime[j+1]小，即i=k*prime[j]，
//那么i*prime[j+1]=(k*prime[j])*prime[j+1]=k’*prime[j]，
//接下去的素数同理。所以不用筛下去了。
//因此，在满足i%prime[j]==0这个条件之前以及第一次   
//满足改条件时,prime[j]必定是prime[j]*i的最小因子。
		}
		
	}
	
}

 long long factor[100][2];
 int fatCnt;
 int getFactors(long long x) {
	 fatCnt=0;
	 long long tmp=x;
	 for(int i=1; prime[i]<=tmp/prime[i]; i++) {
		 factor[fatCnt][1]=0;
		 if(tmp%prime[i]==0) {
			 factor[fatCnt][0]=prime[i];
			 while(tmp%prime[i]==0) {
				 factor[fatCnt][1]++;
				 tmp/=prime[i];				
			}
			 fatCnt++;	
		}	
	}
	 if(tmp!=1) {
		 factor[fatCnt][0]=tmp;
		 factor[fatCnt++][1]=1;	
	}
	
	 return fatCnt;
	
}
int main()
{
	getPrime();
	cout<<getFactors(1000);
	
}