数据结构
核心思路
通过最小素因子
来消除。
代码
int prime[maxn];
int is_prime[maxn];
int Euler(int n)
{
int cnt=0;
memset(prime,0,sizeof(prime));
for (int i = 0; i <= n; ++i)
is_prime[i] = 1;
for (int i = 2;i <= n; i++) //i不仅代表着现在搜索的数,也代表着“倍数”
{
if (is_prime[i])//找到新的素数并储存
prime[++cnt] = i;
for (int j = 1; j <=cnt && i*prime[j] <= n; j++)
{
is_prime[i*prime[j]] = 0;
if (i % prime[j] == 0)
break;
//对于 i%prime[j] == 0 就break的解释 :
//当 i是prime[j]的倍数时,i = kprime[j],如果继续运算 j+1,
//i * prime[j+1] = prime[j] * k prime[j+1],
//这里prime[j]是最小的素因子,当i = k * prime[j+1]时会重复,所以才跳出循环。
}
}
return cnt;
}