欧拉函数以及线性筛法求欧拉函数

 

比如：求6的欧拉函数6=2*3，所以其欧拉函数为6*（2-1）/2*（3-1）/3；
原理：1~N中与N互素的个数等于N - 与N不互质的数，与N不互质的数由N的因子的倍数筛掉，但有的会被多个因子筛掉会多筛，所以还要加上几个因子乘积的倍数（简述一下原理，在此不做证明）

#include<iostream>
using namespace std;
typedef long long ll;
ll a, n;

int main()
{
    cin >> n;
    while (n--)
    {
        cin >> a;
        ll res = a;
        for (int i = 2; i <= a / i; i++)
        {
            if (a % i == 0)res =res* (i-1)/i;
            while (a % i == 0)a /= i;
        }
        if (a > 1)res =res* (a-1)/a;
        cout << res << endl;
    }
}

前置知识：线性筛法求质数

原理：合数都会被其最小质因数筛掉，剩余的为素数

//prime线性筛法求素数，保证每一个合数被它的最小质因数筛掉，prime数组记录素数，v数组记录每个数的最小质因数（同时也起标记作用，标记一个数是否被筛掉）
//如 2 3 4 5 6 7 8 9 10 11 12这一段数（结合下方代码食用）
2没被用过，2会加入prime数组，2遍历一遍prime数组将4筛掉，v[4]==2代表其最小质因数为2。 
3没被用过，3会把6 9筛掉
4被标记为2，不加入prime数组，遍历prime数组[2,3]，将8筛掉，但12不会被筛掉，因为12的最小质因数为2，12 应该在i==6时被6*2筛掉，不是3所以就跳出循环。

#include<iostream>
#include<vector>
using namespace std;
const int N = 1000010;
int n, cnt = 0;//cnt记录素数的个数
int prime[N], v[N];
int main()
{
    cin >> n;
    for (int i = 2; i <= n; i++)
    {
        if (v[i] == 0) {
            v[i] = i;
            prime[++cnt] = i;
        }
        for (int j = 1; j <= cnt; j++)   
        {//如果一个数的最小质因子小于prime中的质数就break，保证每个数被其最小质因子删掉      
                                                         
            if (prime[j] > v[i] || prime[j] * i > n)break;
            v[i * prime[j]] = prime[j];
        }
    }
    for (int i = 1; i <= cnt; i++)
        cout << prime[i] << ' ';

}

线性筛法求欧拉函数（欧拉函数是积性函数）

下方代码：phi[i * prime[j]] = phi[i] * prime[j]//prime[j]是i的因子

原理:{
因为prime[j]是i和i*prime[j]的共同因子，有欧拉函数ϕ(N) = N×(p1−1)/p1×(p2−1)/p2×…×(pm−1)/pm可得求phi[i*prime[j]]的欧拉函数只需要将phi[i]的基数N扩大prime[j]倍
}
phi[i * prime[j]] = phi[i] * (prime[j] - 1);//prime[j]不是i的因子
因为prime[j]不是i的因子但是是prime[j]*i的因子，首先将的phi[i]的基数扩大prime[j]倍，由于新加了一个prime[j]的因子
所以还要乘以（prime[j]-1)/prime[j]最后化简可以得到phi[i * prime[j]] = phi[i] * (prime[j] - 1)

#include<iostream>
using namespace std;
typedef long long ll;
const int N = 1000010;
ll a, n, cnt = 0;
ll prime[N], v[N], phi[N];

ll get_eulors()
{
    phi[1]=1;
    for (int i = 2; i <= n; i++)
    {
        if (v[i] == 0) {
            prime[++cnt] = i;
            phi[i] = i - 1;//i为质数，与其互素的个数为i-1
        }
        for (int j = 1; j <= cnt,prime[j] * i <= n; j++)
        {
            v[prime[j] * i] = prime[j];
            if (i % prime[j] == 0) {
                phi[i * prime[j]] = phi[i] * prime[j];
                break;
            }
            phi[i * prime[j]] = phi[i] * (prime[j] - 1);
        }
    }
    ll res = 0;

    for (int i = 1; i <= n; i++)res += phi[i];
    return res;
}

int main()
{
    cin >> n;
    cout<<get_eulors();
}