欧拉函数和最大公约数

分析:如果两个数的最大公约数是一个质数p,那么这两个数都除以p,得到的两个数的最大公约数一定是1.

反证法:如果得到的两个数的最大公约数不是1,那么把此时的最大公约数乘以上边的最大公约数,得到的一定比上述的最大公约数大,那么上述的最大公约数就不是最大那两个数的最大公约数,所以结论错误。即得到的两个数的最大公约数一定是1.

由于发现两个数都除以p之后,得到的数的最大公约数是1,那么我们可以想到欧拉函数,此时就可以先处理欧拉函数和欧拉函数的前缀和,然后枚举1~n的所有质数,每次求1~n/p(下取整)中与n/p(下取整)互质的个数,由于(1,2),(2,1)属于两个那么还需要乘以2,(1,1)(1,1)属于1个,最后还得减去1.

cpp 复制代码
#include<bits/stdc++.h>

using namespace std;

const int N = 1e7 + 10;

int hpi[N];
int primes[N],cnt;
bool st[N];
int n;
long long s[N];

void init()
{
    hpi[1]=1;
    
    for(int i=2;i<=n;i++)
    {
        if(!st[i]) 
        {
            primes[cnt++]=i;
            hpi[i]=i-1;
        }
        
        for(int j=0;primes[j]<=n/i;j++)
        {
            st[primes[j]*i]=true;
            if(i%primes[j]==0)
            {
                hpi[primes[j]*i]=primes[j]*hpi[i];
                break;
            }
            hpi[i*primes[j]]=hpi[i]*(primes[j]-1);
        }
    }
    
    for(int i=1;i<=n;i++) s[i]=s[i-1]+hpi[i];
}
int main()
{
    cin>>n;
    
    init();
    long long res=0;
    for(int i=0;i<cnt;i++)
    {
        int p=primes[i];
        res+=(2*s[n/p]-1);
    }
    cout<<res<<endl;
    return 0;
}
相关推荐
??tobenewyorker14 分钟前
力扣打卡第二十一天 中后遍历+中前遍历 构造二叉树
数据结构·c++·算法·leetcode
贾全35 分钟前
第十章:HIL-SERL 真实机器人训练实战
人工智能·深度学习·算法·机器学习·机器人
GIS小天1 小时前
AI+预测3D新模型百十个定位预测+胆码预测+去和尾2025年7月4日第128弹
人工智能·算法·机器学习·彩票
满分观察网友z1 小时前
开发者的“右”眼:一个树问题如何拯救我的UI设计(199. 二叉树的右视图)
算法
森焱森2 小时前
无人机三轴稳定化控制(1)____飞机的稳定控制逻辑
c语言·单片机·算法·无人机
循环过三天2 小时前
3-1 PID算法改进(积分部分)
笔记·stm32·单片机·学习·算法·pid
闪电麦坤953 小时前
数据结构:二维数组(2D Arrays)
数据结构·算法
凌肖战3 小时前
力扣网C语言编程题:快慢指针来解决 “寻找重复数”
c语言·算法·leetcode
埃菲尔铁塔_CV算法3 小时前
基于 TOF 图像高频信息恢复 RGB 图像的原理、应用与实现
人工智能·深度学习·数码相机·算法·目标检测·计算机视觉
NAGNIP4 小时前
一文搞懂FlashAttention怎么提升速度的?
人工智能·算法