唯一质因数分解定理:
任意一个合数a仅能以一种方式,写成如下的乘积形式:
\(a\) =$ p1^{e1}\times p2^{e2}\times ...\times pr^{er}$
\(a\)的因子数= \((e1+1)\times (e2+1)\times ....\times (er+1)\)
const int N = (int)2e5 + 7;
int noprime[N], pcnt, p[N / 2];
int nump[N / 2], yinzi[N / 2];
int n, m, top;
void getprime(){
pcnt = 0;
memset(noprime, 0, sizeof(noprime));
noprime[0] = noprime[1] = 1;
for(int i = 2; i < N; ++i) {
if(!noprime[i])p[pcnt++] = i;
for(int j = 0; j < pcnt && i * p[j] < N; ++j) {
noprime[i * p[j]] = 1;
if(i % p[j] == 0)break;
}
}
}
void cal(int t){
memset(nump, 0, sizeof(nump));
top = -1;
int tmp = (int)sqrt(t*1.0);
for(int i = 0; i < pcnt && p[i] <= tmp; ++i) {
if(t % p[i] == 0) {
yinzi[++top] = p[i];
while (t % p[i] == 0) {
nump[top] ++;
t /= p[i];
}
}
if(t == 1)break;
}
if(t > 1) {
yinzi[++top] = t;
nump[top] ++;
}
}
欧拉筛解释:
当\(i\)是\(prime[j]\)的整数倍时(\(i \% prime[j]\;==\;0\)),跳出循环。
因为\(i\)可以看做\(prime[j]*某个数\), \(i*prime[j+1]\)就可以看做\(prime[j]*某个数*prime[j+1]\),\(prime[j]\) 又小于 \(prime[j+1]\),
所以 \(i*prime[j+1]\) 将来必定会被 \(prime[j]*另一个i\) 给筛掉,这里就不用再做了。
同时我们可以发现在满足程序里的两个条件的时候,\(prime[j]\)必定是\(prime[j]*i\)的最小质因子。