欧拉函数
2015-09-18 21:25
169 查看
欧拉函数是指:对于一个正整数n,小于n且和n互质的正整数(包括1)的个数,记作φ(n) 。
通式:φ(x)=x*(1-1/p1)*(1-1/p2)*(1-1/p3)*(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数,x是不为0的整数。φ(1)=1(唯一和1互质的数就是1本身)。
对于质数p,φ(p) = p - 1。注意φ(1)=1.
欧拉定理:对于互质的正整数a和n,有aφ(n) ≡ 1 mod n。
欧拉函数是积性函数——若m,n互质,φ(mn)=φ(m)φ(n)。
若n是质数p的k次幂,φ(n)=p^k-p^(k-1)=(p-1)p^(k-1),因为除了p的倍数外,其他数都跟n互质。
特殊性质:当n为奇数时,φ(2n)=φ(n)
欧拉函数还有这样的性质:
设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)。
筛法实现
#include <iostream>
using namespace std;
int main() {
int eu[110];
for (int i = 1; i<110; i++) eu[i] = i;
for (int i = 2; i<110; i++) {
if (eu[i] == i) {
for (int j = i; j<110; j+=i) eu[j] = eu[j]/i*(i-1);
}
}
for (int i = 1; i<101; i++) cout << i<< " " << eu[i] << endl;
return 0;
}
通式:φ(x)=x*(1-1/p1)*(1-1/p2)*(1-1/p3)*(1-1/p4)…..(1-1/pn),其中p1, p2……pn为x的所有质因数,x是不为0的整数。φ(1)=1(唯一和1互质的数就是1本身)。
对于质数p,φ(p) = p - 1。注意φ(1)=1.
欧拉定理:对于互质的正整数a和n,有aφ(n) ≡ 1 mod n。
欧拉函数是积性函数——若m,n互质,φ(mn)=φ(m)φ(n)。
若n是质数p的k次幂,φ(n)=p^k-p^(k-1)=(p-1)p^(k-1),因为除了p的倍数外,其他数都跟n互质。
特殊性质:当n为奇数时,φ(2n)=φ(n)
欧拉函数还有这样的性质:
设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)。
int eular(int n) { //欧拉函数实现 int ans = n; for (int i = 2; i*i<=n; i++) { if (n%i == 0) { ans -= ans/i; while (n%i == 0) n = n/i; } } if (n > 1) ans -= ans/n; return ans; }
筛法实现
#include <iostream>
using namespace std;
int main() {
int eu[110];
for (int i = 1; i<110; i++) eu[i] = i;
for (int i = 2; i<110; i++) {
if (eu[i] == i) {
for (int j = i; j<110; j+=i) eu[j] = eu[j]/i*(i-1);
}
}
for (int i = 1; i<101; i++) cout << i<< " " << eu[i] << endl;
return 0;
}
相关文章推荐
- php连接sqlserver
- LEETCODE-Climbing Stairs
- 京东2015在线笔试 (1)
- Java HashMap的工作原理
- poj 2479 最大连续子段和 dp算法
- 欧拉函数
- Effective Java 英文 第二版 读书笔记 Item 7:Avoid finalizers
- 京东2015在线笔试 (1)
- mysql和oracle的一个汉字占几个字符
- 安卓sqlite数据库的使用
- VMware.Workstation.11.1.0安装rhel-server-5.5-x86_64-dvd.ISO安装
- 题目:1<n<9, n个猴子,第一个猴子扔一个,然后平均分为n份,拿一份; 第二个猴子扔一个,然后平均分n份,拿一份。。。依次到最后一个猴子也是扔一个,均分n份,拿一份。
- Android绘图:Canvas绘制联系人姓氏侧滑列表
- Aactivity和Service之间的通信
- Java记录 -36- 冒泡排序
- SVC和PEDSVC
- 在 Ubuntu 14.04 上裝好 CUDA 7.0
- 分糖果
- SVN - 基础知识
- Hadoop安装准备(2)