素因子分解 模板

每个大于1的数都可以写成素数的乘积。

//the p's size = N/lnN;10^5 and 10^6 can use 10
const int N = 1000000;
bool primeTable[N+5];
int p[N/10],tot;
int e[N/10];

void make_primeTable(){
    tot = 0;
    fill(primeTable,primeTable+N,1);
    primeTable[0] = false;
    primeTable[1] = false;

    int maxed = sqrt(N);
    for(int i = 2; i <= maxed; ++i){
        if(primeTable[i] == true){
            p[tot++] = i;
            for(int j = i*i; j <= N; j += i)
                primeTable[j] = false;
        }
    }
    for(int i = maxed+1; i <= N; ++i)
        if(primeTable[i] == true)  p[tot++] = i;
}
void add_factor(int n){
    memset(e,0,sizeof(e));
    if(n < 2)  return;
    for(int i = 0; i < tot; ++i){
        while(n%p[i] == 0){
            e[i]++;
            n /= p[i];
        }
        if(n == 1)  break;
    }
    //if(n > 1){说明此时还有一个超出范围的素数，通常有且仅有一个，例如对10^12素因子分解枚举10^6次内的素数，那么>10^6的数只有一个
}