欧拉函数知识点总结及欧拉函数打表代码(数论)

一、概念:

  在数论,对正整数n,欧拉函数是小于n的正整数中与n互质的数的数目(φ(1)=1)。

  例如φ(8)=4,因为1,3,5,7均和8互质。

  欧拉函数用希腊字母φ表示,φ(N)表示N的欧拉函数.

  对φ(N)的值,我们可以通俗地理解为小于N且与N互质的数的个数(包含1).

 (初学者一定注意:此处的欧拉函数与图论中的欧拉回路、欧拉公式不同,欧拉大神这一辈子著作等身啊)。

对于互质的理解:

互质数为数学中的一种概念,即两个或多个整数的公因数只有1的非零自然数。公因数只有1的两个非零自然数,叫做互质数。 

互质数具有以下定理:
(1)两个数的公因数只有1的两个非零自然数,叫做互质数;举例:2和3,公因数只有1,为互质数;
(2)多个数的若干个最大公因数只有1的正整数,叫做互质数;
(3)两个不同的质数,为互质数;
(4)1和任何自然数互质。两个不同的质数互质。一个质数和一个合数,这两个数不是倍数关系时互质。不含相同质因数的两个合数互质;
(5)任何相邻的两个数互质;(必为一奇一偶)
(6)任取出两个正整数他们互质的概率(最大公约数为一)为6/π^2。
 

二、通式:

 

   

  其中p1, p2……pn为x的所有质因数,x是不为0的整数。

       φ(1) = 1(唯一和1互质的数(小于等于1)就是1本身)。 (注意:每种质因数只一个。比如 12 = 2*2*3 那么      φ(12) = 12 * (1-1/2) * (1-1/3)=4  )

       若 n = p^k  (  p为 质数 ),则 φ(n) = p^k-p^(k-1) = (p-1)p^(k-1),( 除 p 的倍数外,其他数均为 p 的互质数 )。

       若n = p( p 为质数),则  φ(n) = p-p^(1-1) = p-1。

三、性质:

(1)   p^k型欧拉函数:

若N是质数p(即N=p), φ(n)= φ(p)=p-p^(k-1)=p-1。

若N是质数p的k次幂(即N=p^k),φ(n)=p^k-p^(k-1)=(p-1)p^(k-1)。

(2)mn型欧拉函数

设n为正整数,以φ(n)表示不超过n且与n互素的正整数的个数,称为n的欧拉函数值。若m,n互质,φ(mn)=(m-1)(n-1)=φ(m)φ(n)。

(3)特殊性质:

若n为奇数时,φ(2n)=φ(n)。

对于任何两个互质 的正整数a,n(n>2)有:a^φ(n)=1 mod n (恒等于)此公式即 欧拉定理

当n=p 且 a与素数p互质(即:gcd(a,p)=1)则上式有: a^(p-1)=1 mod n (恒等于)此公式即 费马小定理

四、模板

1.直接求小于或等于n,且与n互质的个数:

 1 int  eular(int n)
 2 {
 3     int i,ret=n;
 4     for(i=2; i<=sqrt(n); i++)
 5     {
 6         if(n%i==0)
 7         {
 8             ret=ret/i*(i-1);///这里先使用除法是为了防止溢出,ret=ret*(1-1/p(1))
 9             while(n%i==0)///为了完全消除我们已经除完了刚才得到的那个i因子,确保下一个得到的i是n的素因子
10             {
11                 n/=i;
12             }
13         }
14     }
15     if(n>1)///可能还剩下一个素因子没有除 
16     {
17         ret=ret/n*(n-1);
18     }
19     return ret;
20 }

 

2.筛选模板:求[1,n]之间每个数的质因数的个数

 

如果我们要求的数比较多,如果一个一个求那么很容易就超时,所以我们自然而然就想到——打表。

如果我们依照上述思想,来个最朴素的打表。

 1 void euler()
 2 {
 3     p[1]=1;
 4     for(int i=2; i<=MAXN; i++)
 5     {
 6         int n=i;
 7         p[i]=i;
 8         for(int j=2; j<=sqrt(n); j++)
 9         {
10             if(n%j==0)
11             {
12                 p[i]=p[i]/j*(j-1);
13                 while(n%j==0)
14                 {
15                     n=n/j;
16                 }
17             }
18         }
19         if(n>1)
20         {
21             p[i]=p[i]/n*(n-1);
22         }
23     }
24     for(int i=2; i<MAXN; i++)
25     {
26         p[i]+=p[i-1];
27     }
28 }

但是这种打表方法还是不太理想,效率不是太高,下面给出两种改进后的算法。

 

对于这种打表方法,我开始不是很明白,网络上也没有比较完全的说明,下面给出一些我自己的理解。

之前的打表方式是枚举每个数,再去找每个数的质因子,由于有一些数有相同的质因子,无疑会增大时间复杂度。在这里我们选择去枚举素因子,素因子可是少很多。比如素因子是2,那么我们便能求出那些含有2这个素因子的一些数的欧拉函数,例如2,4,8,10......当然像10,12,20,这样的数不仅仅只含有2这一个素因子,我们先欧拉函数的通式先求出一部分来,之后枚举其他素因子的时候再次计算,也就跳过了if(!euler[j])。同时下次需要枚举的数将会是避开已经求出来或者部分求出来(例如2下的10,12,20)数,作为一个新的素因子再次进入循环。

是不是有点像筛法求素数的思路。。。。

 1 #define size 1000001
 2 int euler[size];
 3 void Init()
 4 {
 5     euler[1]=1;
 6     for(int i=2; i<size; i++)
 7     {
 8         if(!euler[i])
 9         {
10             for(int j=i; j<size; j+=i)
11             {
12                 if(!euler[j])
13                 {
14                     euler[j]=j;
15                 }
16                 euler[j]=euler[j]/i*(i-1);//先进行除法是为了防止中间数据的溢出
17             }
18         }
19     }
20 }

 

 

转载于:https://www.cnblogs.com/wkfvawl/p/9317125.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
### 回答1: 欧拉函数也称为欧拉-费马函数或费马函数,用符号φ(n)表示表示小于等于n且与n互质的正整的个。 以下是一个用MATLAB编写的求欧拉函数代码示例: ```matlab function result = euler_function(n) result = n; % 初始化结果为n p = 2; % 设置初始的素为2 % 计算欧拉函数 while p^2 <= n if mod(n, p) == 0 result = result * (1 - 1 / p); while mod(n,p) == 0 n = n / p; end end p = p + 1; end % 处理剩余的质因 if n > 1 result = result * (1 - 1 / n); end end ``` 使用方法: ```matlab n = 12; % 输入要求解欧拉函数n result = euler_function(n); % 调用欧拉函数的方法 disp(result); % 输出结果 ``` 这个代码通过迭代计算n的所有质因,并根据欧拉函数的公式进行计算,最终得到欧拉函数的值。 ### 回答2: 欧拉函数(Euler's totient function)是一个计算小于等于给定正整n的所有与n互质的正整的个。以下是用Matlab实现欧拉函数代码: ```matlab function euler = euler_function(n) euler = n; % 将结果初始化为n p = 2; % 从最小的素2开始 while p * p <= n % 当p的平方大于n时结束循环 if mod(n, p) == 0 % 若n能整除p euler = euler - euler / p; % 将结果减去n除以p的值(贡献) while mod(n, p) == 0 % 如果n继续能整除p,则继续除以p n = n / p; end end p = p + 1; % 继续查找下一个素 end if n > 1 % 若n不为1,则n是一个大于sqrt(n)的质 euler = euler - euler / n; % 将结果减去n除以n的值(贡献) end end ``` 这个函数接受一个正整n作为输入,并返回n的欧拉函数值。函数首先将结果初始化为n,然后从最小的素2开始,一直遍历到sqrt(n)为止。如果n能被p整除,它将以p的贡献减少结果。然后继续查找下一个素。最后,如果n不等于1,则n是一个大于sqrt(n)的质,将其贡献减少结果。 希望以上的代码可以满足你的需求。 ### 回答3: 欧拉函数数论一个重要函数,用于计算小于n且与n互质的正整的个欧拉函数的公式为:φ(n)=n×(1-1/p1)×(1-1/p2)×⋯×(1-1/pk)其p1,p2,...,pk是n的所有质因。 下面是求欧拉函数的MATLAB代码: ```matlab function result = euler_phi(n) result = n; for i = 2 : round(sqrt(n)) if (n % i == 0) result = result * (1 - 1/i); while (n % i == 0) n = n / i; end end end if (n > 1) result = result * (1 - 1/n); end end ``` 在这段代码,我们首先将结果初始化为n。然后从2到sqrt(n)遍历,如果n能被i整除,则说明i是n的一个质因。这时,我们将结果乘以(1-1/i),然后不断将n除以i直到n不能再被i整除为止。最后,如果n大于1,说明n本身也是一个质因,我们将结果乘以(1-1/n)。 这段代码可以在MATLAB调用,例如: ```matlab n = 12; result = euler_phi(n); disp(result); ``` 这样就可以得到12的欧拉函数值。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值