欧拉函数模板

来源:互联网 发布:uujuly淘宝女装店 编辑:程序博客网 时间:2024/05/22 04:46
 5 LL eulr(LL n)  6 {  7     LL ans=n;  8     for(LL i=2;i*i<=n;i++)  9     {    10         if(n%i==0) 11         {    12             ans=ans/i*(i-1); 13             while(n%i==0) 14                 n/=i; 15         } 16     } 17     if(n>1) 18         ans=ans/n*(n-1); 19     return ans; 20 }

void Euler(){    phi[1] = 1;    for(int i = 2; i < N; i ++){        if(!phi[i]){            for(int j = i; j < N; j += i){                if(!phi[j]) phi[j] = j;                phi[j] = phi[j] / i * (i-1);            }        }    }}

更快

#include<cstdio>using namespace std;const int N = 1e6+10 ;int phi[N], prime[N];int tot;//tot计数,表示prime[N]中有多少质数 void Euler(){    phi[1] = 1;    for(int i = 2; i < N; i ++){        if(!phi[i]){            phi[i] = i-1;            prime[tot ++] = i;        }        for(int j = 0; j < tot && 1ll*i*prime[j] < N; j ++){            if(i % prime[j]) phi[i * prime[j]] = phi[i] * (prime[j]-1);            else{                phi[i * prime[j] ] = phi[i] * prime[j];                break;            }        }    }}
原创粉丝点击