ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

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

2019-03-21 10:38:44  阅读:277  来源: 互联网

标签:const Arrays 题解 Codeforces long int ans miu define


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 + 7;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = 1000000007;
const double eps = 1e-6;
const double PI = acos(-1);

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

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

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

/*
*/

 

标签:const,Arrays,题解,Codeforces,long,int,ans,miu,define
来源: https://www.cnblogs.com/CJLHY/p/10570014.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有