演算法-尤拉函式
尤拉函式簡介——摘自百度百科
注意以下三個特殊性質
程式設計實現
利用尤拉函式和它本身不同質因數的關係,用篩法計算出某個範圍內所有數的尤拉函式值。
//直接求解尤拉函式
#include<cstdio>
int euler(int n){ //返回euler(n)
int res=n,a=n;
for(int i=2;i*i<=a;i++){//從小到大嘗試n的質因數
if(a%i==0){//如果i是n的質因數
res=res/i*(i-1);//提了一個1/i出來,先進行除法是為了防止中間資料的溢位
while(a%i==0) a/=i;//尤拉函式只記算一種質因數
}
}
if(a>1) res=res/a*(a-1);//如果最後還剩因子
return res;
}
int main(){
int x;
scanf("%d",&x);
printf("%d",euler(x));
return 0;
}
//篩選法打尤拉函式表
#include<cstdio>
#define Max 1000001
int euler[Max];
void Init(){
euler[1 ]=1;
for(int i=2;i<Max;i++)
euler[i]=i;
for(int i=2;i<Max;i++)
if(euler[i]==i)//如果i是質數
for(int j=i;j<Max;j+=i)
euler[j]=euler[j]/i*(i-1);//提一個1/i,先進行除法是為了防止中間資料的溢位
return ;
}
int main(){
Init();
for(int i=1;i<=100;i++)
printf ("%d\n",euler[i]);
return 0;
}
相關推薦
演算法-尤拉函式
尤拉函式簡介——摘自百度百科 注意以下三個特殊性質 程式設計實現 利用尤拉函式和它本身不同質因數的關係,用篩法計算出某個範圍內所有數的尤拉函式值。 //直接求解尤拉函式 #include<cstdio&g
Pollard的rho啟發式因子分解演算法 & [CodeVS 4939] 尤拉函式:Miller-Rabin + Pollard-rho 質因數分解
Pollard的rho啟發式因子分解演算法用於給出整數的一個因子。在一定的合理假設下,如果n有一個因子p,可在O(p√)的期望時間內可找出n的一個因子p。 關於其複雜度,Wikipedia是這樣敘述的: If the pseudo random num
求數論求約數和 與 互質和演算法 (分解質因數與尤拉函式)
Description One day, Qz met an easy problem. But after a 5-hout-long contest in CSU, he became very tired and he wanted to call his girl
【演算法】尤拉函式——小於n的數中與n互質數的數目
尤拉函式 簡介 在數論,對正整數n,尤拉函式是小於n的數中與n互質的數的數目。此函式以其首名研究者尤拉命名(Ruler’so totient function),它又稱為Euler’s totient function、φ函式、尤拉商數等。 例如φ(8)=4
數學 尤拉函式相關
尤拉函式相關 1,\(phi(i)\)表示在1到i的數中與i互質的數的個數。 2,\(O(\sqrt{n})\)求\(phi\) 算數基本定理: \[ phi(i)=i*(p_1-1)/p_1*(p_2-1)/p_2*……*(p_k-1)/p_k \] 列舉質因數套公式即可: code:
一類尤拉函式相關的求和式推導
\(\\\) 寫在前面 因為最近做了不少和尤拉函式相關的求和問題,而這一類求和的推導有沒有涉及到反演和卷積,所以單獨寫一寫。 給出的題目順序與難度大致無關,是按照個人做題的順序安排的。 再次宣告尤拉函式的定義:\(\varphi(x)\) 表示 \([1,x]\) 裡的所有整數中,與 \(x\)
POJ3090 Visible Lattice Points (數論:尤拉函式模板)
題目連結:傳送門 思路: 所有gcd(x, y) = 1的數對都滿足題意,然後還有(1, 0) 和 (0, 1)。 #include <iostream> #include <cstring> using namespace std; const in
[BZOJ4026]dC Loves Number Theory 尤拉函式+線段樹
連結 題意:給定長度為 \(n\) 的序列 A,每次求區間 \([l,r]\) 的乘積的尤拉函式 題解 考慮離線怎麼搞,將詢問按右端點排序,然後按順序掃這個序列 對於每個 \(A_i\) ,列舉它的質因數,由於不同的質因數只算一次,所以我們只關心每個質數它最後一次出現的位置,開一棵線段樹維護
hdu5528(積性函式+尤拉函式)
題意:設(題目已把f(6)的表給出),,給定n(n<=1e9),求g(n) 最主要的是求f(m),考慮到模數大於0的情況比較多,所以考慮考慮求ij mod n==0的情況。。 然後其實從題目給的表中看出對一個a來說,他0的個數和gcd(a,n)有關,其實也很容易證明,對一個數a來說
尤拉函式 線段樹 狀壓 奇數國
讓我們一起來%forever_shi神犇 題意:求區間積的 ϕ \phi ϕ值。 題解:
UVA11426 GCD - Extreme (II) (尤拉函式/莫比烏斯反演)
UVA11426 GCD - Extreme (II) 題目描述 PDF 輸入輸出格式 輸入格式: 輸出格式: 輸入輸出樣例 輸入樣例#1: 10 100 200000 0 輸出樣例#1: 67 13015 143295493160 Solution 這道題我用莫比烏斯反演和尤拉函式都寫了一遍,發現
【hdu 5728 PowMod】【數論】【尤拉函式】【尤拉降冪遞迴取模】【尤拉積性函式】
【連結】 http://acm.hdu.edu.cn/showproblem.php?pid=5728 【題意】 n是無平方因子的數 定義k=∑mi=1φ(i∗n) mod 1000000007,求K^k^k^k......%p 【思路】 先尤拉性質求出k
尤拉篩尤拉函式
尤拉函式 phi[i]表示 1~i 內與 i 互質的個數 通式:phi[i]=x∏(1-pi) pi表示 i 的質因數 是積性函式 phi[i]*phi[j]=phi[i*j] 做法:一般用尤拉篩 先貼一份程式碼: 1 #include<cstdio>
2015 ICPC瀋陽現場賽 F. Frogs (尤拉函式)
題目連結 m個石頭圍成一圈,一群青蛙從0開始跳,第i個青蛙每步跳ai距離,求所有能被跳到的石頭的編號之和。 容易推出石頭x被第i只青蛙跳到的充要條件是,顯然這與下面的命題是等價的: 石頭x被跳到的充要條件是存在一個i,使得。 gcd(x, m)顯然是m的因數,那麼我們就可以列舉m的因
[SDOI2008]沙拉公主的困惑 線性篩_尤拉函式_逆元_快速冪
Code: #include<cstdio> using namespace std; typedef long long ll; const int maxn=10000000+1; long long mod; ll fac[maxn]; ll inv[maxn];
[SDOI2008]儀仗隊 尤拉函式
顯然,橫縱座標要互質才能被看到。 處理出與橫座標 i i i 互質的縱座標的個數,求一遍字首和即可。 Co
GCD - Extreme (II) UVA - 11426 尤拉函式_數學推導
Code: #include<cstdio> using namespace std; const int maxn=4000005; const int R=4000002; const int N=4000002; long long sumv
BZOJ4869 六省聯考2017相逢是問候(線段樹+尤拉函式)
由擴充套件尤拉定理,a^(a^(a^(……^x)))%p中x作為指數的模數應該是φ(φ(φ(φ(……p)))),而p取log次φ就會變為1,也即每個位置一旦被修改一定次數後就會變為定值。線段樹維護區間剩餘修改次數的最大值,暴力修改即可。 可以預處理出每個位置進行k次操作後的值。直接計算是log^3的
HDU 5514.Frogs-尤拉函式 or 容斥原理
Frogs Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/65536 K (Java/Others)Total Submission(s): 4904 &n
尤拉函式(模板)
#include <iostream> #include <cmath> using namespace std; int oula(int n) { int ans=n; int i; for(i=2;i<=sqrt(n);i++) {