欧拉函数模板

复制代码
void eular()
{
    memset(vis,0,sizeof(vis));
    vis[0]=vis[1]=1;
    for(i=2;i*i<=N;i++)
    {
        if(vis[i]==0)
        {
            for(j=i*i;j<=N;j+=i)
                vis[j]=1;
        }
    } //这段求出了N内的所有素数 
    for(i=1;i<=N;i++)
        phi[i]=i;
    for(i=2;i<=N;i++)
    {
        if(vis[i]==0)
        {
            for(j=i;j<=N;j+=i)//这里从i开始,必定能整除i,其倍数也同理
                phi[j]=phi[j]/i*(i-1); //此处注意先/i再*(i-1),否则范围较大时会溢出
        }
    }
}
复制代码

递归求欧拉函数

for (i = 1; i <= maxn; i++) phi[i] = i;
for (i = 2; i <= maxn; i += 2) phi[i] /= 2;
for (i = 3; i <= maxn; i += 2) if(phi[i] == i) {
for (j = i; j <= maxn; j += i)
phi[j] = phi[j] / i * (i - 1);

单独求欧拉函数

复制代码
unsigned euler(unsigned x)
{
    // 就是公式
    unsigned i, res=x;
    for (i = 2; i < (int)sqrt(x * 1.0) + 1; i++)
        if(x%i==0)
        {
            res = res / i * (i - 1);
            while (x % i == 0) x /= i; // 保证i一定是素数
        }
    if (x > 1) res = res / x * (x - 1);
    return res;
}
复制代码

 

posted @   calmound  阅读(269)  评论(0编辑  收藏  举报
编辑推荐:
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
· 基于Microsoft.Extensions.AI核心库实现RAG应用
· Linux系列:如何用heaptrack跟踪.NET程序的非托管内存泄露
· 开发者必知的日志记录最佳实践
阅读排行:
· TypeScript + Deepseek 打造卜卦网站:技术与玄学的结合
· Manus的开源复刻OpenManus初探
· AI 智能体引爆开源社区「GitHub 热点速览」
· C#/.NET/.NET Core技术前沿周刊 | 第 29 期(2025年3.1-3.9)
· 从HTTP原因短语缺失研究HTTP/2和HTTP/3的设计差异
历史上的今天:
2012-07-14 Havel定理
点击右上角即可分享
微信分享提示