Codeforces 915G Coprime Arrays 莫比乌斯反演 (看题解)

Coprime Arrays

啊，我感觉我更本不会莫比乌斯啊啊啊， 感觉每次都学不会， 我好菜啊。

#include<bits/stdc++.h>
#define LL long long
#define fi first
#define se second
#define mk make_pair
#define PLL pair<LL, LL>
#define PLI pair<LL, int>
#define PII pair<int, int>
#define SZ(x) ((int)x.size())
#define ull unsigned long long

using namespace std;

const int N = 2e6 + ;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = ;
const double eps = 1e-;
const double PI = acos(-);

int n, k, miu[N], Pown[N], sum[N];

int Power(int a, int b) {
    int ans = ;
    while(b) {
        if(b & ) ans = 1LL * ans * a % mod;
        a = 1LL * a * a % mod; b >>= ;
    }
    return ans;
}

int main() {
    miu[] = ;
    for(int i = ; i < N; i++)
        for(int j = i + i; j < N; j += i)
            miu[j] -= miu[i];
    scanf("%d%d", &n, &k);
    for(int i = ; i <= k; i++) Pown[i] = Power(i, n);
    int ans = , tmp = ;
    for(int i = ; i <= k; i++) {
        for(int j = i; j <= k; j += i) {
            sum[j] = ((sum[j] + 1LL * miu[i] * (Pown[j / i] - Pown[j / i - ])) % mod + mod) % mod;
        }
        tmp = (tmp + sum[i]) % mod;
        ans = (ans + (tmp ^ i)) % mod;
    }
    printf("%d\n", ans);
    return ;
}

/*
*/