洛谷--P2158 仪杖队--欧拉函数

这题可以用莫比乌斯反演来写，相关式子：
$ans={\sum }_{i=1}^{n-1}{\sum }_{j=1}^{n-1}[\gcd (i,j)==1]+2(n>=2)$
$当n=1,ans=0$
但我们不用反演去写，我们思考欧拉函数

欧拉函数

$\varphi (n)=n\prod (1-\frac{1}{p})$,p为n的质因数
用和式的写法：$\varphi (n)={\sum }_{i=1}^n[\gcd (i,n)==1]$
分析:
$ans={\sum }_{i=1}^{n-1}{\sum }_{j=1}^{n-1}[\gcd (i,j)==1]+2$
$ans={\sum }_{i=1}^{n-1}{\sum }_{j=1}^i[\gcd (i,j)==1]+{\sum }_{i=1}^{n-1}{\sum }_{j=i}^{n-1}[\gcd (i,j)==1]+1$
考虑欧拉函数：
$ans=2{\sum }_{i=1}^{n-1}\varphi (i)+1(n>=2)$
做完前缀和，查询是O(1)感觉比莫比乌斯还快

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
template <typename T>
void out(T x) { cout << x << endl; }
const int N = 4e4 + 4;
ll fast_pow(ll a, ll b, ll p) {ll c = 1; while(b) { if(b & 1) c = c * a % p; a = a * a % p; b >>= 1;} return c;}
ll phi[N];
void get_phi()
{
    memset(phi, 0, sizeof(phi));
    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);
            }
    }
    for(int i = 1; i < N; i ++)
        phi[i] += phi[i - 1];
}
int main()
{
    ios::sync_with_stdio(false);
    get_phi();
    int n;
    cin >> n;
    cout << (n != 1 ? 2 * phi[n - 1] + 1 : 0) << endl;
}