http://acm.hdu.edu.cn/showproblem.php?pid=2824
以下是2种求欧拉函数的算法
利用素数筛选:
明显第一种的编程复杂度要低很多
所以,一般情况下(N不是很大),采用第一种即可;
贴在这里供以后复习
定义: 对于正整数n,φ(n)是小于或等于n的正整数中,与n互质的数的数目;
例如: φ( 8 ) = 4 , 因为1, 3 , 5 ,7均和8互质。
性质: 1 . 若p是质数,φ(p) = p - 1 .
2 . 若n是质数p的k次幂,φ(n) = (p - 1 )p ^ (k - 1 )
因为除了p的倍数都与n互质
3 . 欧拉函数是积性函数,若m,n互质,φ(mn) = φ(m)φ(n)
根据这3条性质我们就可以退出一个整数的欧拉函数的公式,因为一个数总可以一些质数的乘积的形式。
E(k) = (p1 - 1 )(p2 - 1 )…(pi - 1 ) * (p1 ^ (a1 - 1 ))(p2 ^ (a2 - 1 ))…(pi ^ (ai - 1 ))
= k * (p1 - 1 )(p2 - 1 )…(pi - 1 ) / (p1 * p2 * …pi)
= k * ( 1 - 1 / p1) * ( 1 - 1 / p2)…( 1 - 1 / pk)
在程序中利用欧拉函数如下性质,可以快速求出欧拉函数的值(a为N的质因素)
若(N % a == 0 && (N / a) % a == 0 ) 则有:E(N) = E(N / a) * a;
若(N % a == 0 && (N / a) % a != 0 ) 则有:E(N) = E(N / a) * (a - 1 );
例如: φ( 8 ) = 4 , 因为1, 3 , 5 ,7均和8互质。
性质: 1 . 若p是质数,φ(p) = p - 1 .
2 . 若n是质数p的k次幂,φ(n) = (p - 1 )p ^ (k - 1 )
因为除了p的倍数都与n互质
3 . 欧拉函数是积性函数,若m,n互质,φ(mn) = φ(m)φ(n)
根据这3条性质我们就可以退出一个整数的欧拉函数的公式,因为一个数总可以一些质数的乘积的形式。
E(k) = (p1 - 1 )(p2 - 1 )…(pi - 1 ) * (p1 ^ (a1 - 1 ))(p2 ^ (a2 - 1 ))…(pi ^ (ai - 1 ))
= k * (p1 - 1 )(p2 - 1 )…(pi - 1 ) / (p1 * p2 * …pi)
= k * ( 1 - 1 / p1) * ( 1 - 1 / p2)…( 1 - 1 / pk)
在程序中利用欧拉函数如下性质,可以快速求出欧拉函数的值(a为N的质因素)
若(N % a == 0 && (N / a) % a == 0 ) 则有:E(N) = E(N / a) * a;
若(N % a == 0 && (N / a) % a != 0 ) 则有:E(N) = E(N / a) * (a - 1 );
以下是2种求欧拉函数的算法
1
void
init()
2 {
3 __int64 i,j;
4 e[ 1 ] = 1 ;
5 for (i = 2 ;i <= N;i ++ )
6 if ( ! e[i])
7 {
8 for (j = i; j <= N; j += i)
9 {
10 if ( ! e[j])
11 e[j] = j;
12 e[j] = e[j] / i * (i - 1 );
13 }
14 }
15 }
2 {
3 __int64 i,j;
4 e[ 1 ] = 1 ;
5 for (i = 2 ;i <= N;i ++ )
6 if ( ! e[i])
7 {
8 for (j = i; j <= N; j += i)
9 {
10 if ( ! e[j])
11 e[j] = j;
12 e[j] = e[j] / i * (i - 1 );
13 }
14 }
15 }
利用素数筛选:
void
init()
{
__int64 i, j;
p[ 0 ] = 1 ; // 记录素数个数
p[ 1 ] = 2 ;
for (i = 3 ; i < N; i += 2 )
{
if (hash[i])
continue ;
p[ ++ p[ 0 ]] = i;
for (j = i * i; j < N; j += i)
hash[j] = true ;
} // 筛素数
e[ 1 ] = 1 ;
for (i = 1 ; i <= p[ 0 ]; i ++ )
e[p[i]] = p[i] - 1 ; // 初始化素数的phi
for (i = 2 ; i < N; i ++ )
{
if ( ! e[i])
{
for (j = 1 ; j <= p[ 0 ]; j ++ )
if (i % p[j] == 0 )
{
if (i / p[j] % p[j])
e[i] = e[i / p[j]] * e[p[j]];
else
e[i] = e[i / p[j] ] * p[j];
break ;
} // 利用上述性质求解
}
}
return ;
}
{
__int64 i, j;
p[ 0 ] = 1 ; // 记录素数个数
p[ 1 ] = 2 ;
for (i = 3 ; i < N; i += 2 )
{
if (hash[i])
continue ;
p[ ++ p[ 0 ]] = i;
for (j = i * i; j < N; j += i)
hash[j] = true ;
} // 筛素数
e[ 1 ] = 1 ;
for (i = 1 ; i <= p[ 0 ]; i ++ )
e[p[i]] = p[i] - 1 ; // 初始化素数的phi
for (i = 2 ; i < N; i ++ )
{
if ( ! e[i])
{
for (j = 1 ; j <= p[ 0 ]; j ++ )
if (i % p[j] == 0 )
{
if (i / p[j] % p[j])
e[i] = e[i / p[j]] * e[p[j]];
else
e[i] = e[i / p[j] ] * p[j];
break ;
} // 利用上述性质求解
}
}
return ;
}
明显第一种的编程复杂度要低很多
所以,一般情况下(N不是很大),采用第一种即可;
贴在这里供以后复习
![](https://i-blog.csdnimg.cn/blog_migrate/ea4429455bb133d049ee0bc6aa94bb1e.gif)