欧拉函数定义
φ(n)=n(1-1/p1)(1-1/p2)(1-1/p3)*(1-1/p4)……(1-1/pn)*
其中p1, p2……pn为n的所有质因数,n是不为0的整数。φ(1)=1(唯一和1互质的数就是1本身)。
通俗一点,eular函数提供了小于等于n的质因数数量
欧拉函数的性质
- 欧拉函数是积性函数——若m,n互质,φ(mn)=φ(m)φ(n)。
- 欧拉函数是积性函数,但不是完全积性函数。若m,n互质,则φ(m∗n)=φ(m)∗φ(n)φ(m*n)=φ(m)φ(n)φ(m∗n)=φ(m)∗φ(n)。特殊的,当m=2,n为奇数时,φ(2n)=φ(n)。
对于其性质这一部分,由于本人才疏学浅,并没有做到与欧拉函数性质相关的题目,所以先在这里给我挖一个坑
欧拉函数的写法
单点写法
ll eular(ll n){
ans = n;
for(int i=2;i<=sqrt(n);i++){
//寻找质因数,注意只需要从最小的开始找就可以了,找到的自动就是质因数
//因为是在找质因数,所以你就不用找过这个
if(n%i == 0){
//ans = ans*(1-1/i);
//这里千万不要用这个公式,因为1/i不是你想象中的答案
ans = ans/i*(i-1);
while(n%i == 0) n/=i;
}
}
//注意:如果n只由一个质因数的一次方组成,那么这个数就会被留到最后,因而还需要加一步特判来进行
if(n > 1) ans = ans/n*(n-1);
return ans;
}
尤其是最后一步判断,n是否为1,很多博客都没有写这是为什么,其实这是因为,我们每一步都在改变n的值,如果n的因数里面有一个为质因数,而且只含有1次,那么就会失效,比如说:33,42
欧拉函数打表
const ll m = 3000000;
bool vis[m];
int eu[m];
void eular(ll n){
eu[1] = 1;
for(int i=2;i<m;i++){
if(!eu[i]){
for(int j = i;j<m;j+=i){
if(!eu[j]){
eu[j] = j;
}
//相当于把每一个包含因子i的函数φ(n)都赋了初值 “n ”
eu[j] = eu[j]*(i-1)/i;
//还是那个道题,从小到大来将已经访问过的包括其倍数都赋值
//那么自然就是找到了每个n对应的质因数
//那么对于没有访问过的质因数对应的倍数 j 如果包括某个值的话就
//会执行 eu[j] = eu [j] *(i-1) /i;
}
}
}
}
其中的注释可能有点难懂,但相比其他博主,这个多看一会就能懂
递推求欧拉函数
for(i=1; i<=maxn; i++)
p[i]=i;
for(i=2; i<=maxn; i+=2)
p[i]/=2;
for(i=3; i<=maxn; i+=2)
if(p[i]==i)
{
for(j=i; j<=maxn; j+=i)
p[j]=p[j]/i*(i-1);
}
递推由于并不是很常用,所以有需要的话可以看一下其他博主的文章