标签：函数 筛法 int res vis maxn include 欧拉

公式法求解单个欧拉函数

//https://blog.csdn.net/elbadaernu/article/details/76945418
ll euler(ll n)
{ //返回euler(n)   
       ll res=n,a=n;  
       for(ll i=2;i*i<=a;i++){  
           if(a%i==0){  
               res=res/i*(i-1);//先进行除法是为了防止中间数据的溢出   
               while(a%i==0) a/=i;  
           }  
       }  
      if(a>1) res=res/a*(a-1);  
      return res;  
}

公式法打表

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

欧拉筛法求欧拉函数

//#include<bits/stdc++.h>
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int maxn=4e6+10;
bool vis[maxn];
long long a[maxn];
long long p[maxn];
int main()
{
    int n,m;
    n=3e6;
//  n=12;
    memset(vis,0,sizeof(vis));
    int tot=0;
    a[1]=0;
    for (int i=2;i<=n;i++)
    {
        if (!vis[i])
        {
            p[++tot]=i;
            a[i]=i-1;
        }
        for (int j=1;j<=tot;j++)
        {
            if (p[j]*i>n) break;
            vis[p[j]*i]=1;
            if (i%p[j]==0)
            {
                a[i*p[j]]=a[i]*p[j];
                break;
            }
            else a[i*p[j]]=a[i]*a[p[j]];
        }
    }
}

欧拉筛法求素数

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int maxn=16010;
bool vis[maxn];
int a[maxn];
int main()
{
    int n=16000;
    int tot=0;
    memset(vis,0,sizeof(vis));
    for (int i=2;i<=n;i++)
    {
        if (!vis[i]) a[++tot]=i;
        for (int j=1;j<=tot;j++) if (a[j]*i<=n)
        {           
            vis[a[j]*i]=1;
            if (i%a[j]==0) break;
        }
    }
}

标签：函数,筛法,int,res,vis,maxn,include,欧拉
来源： https://www.cnblogs.com/reshuffle/p/12251576.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。