牛客小白月赛12 华华给月月出题

最新推荐文章于 2022-03-02 16:43:37 发布

kmlver

最新推荐文章于 2022-03-02 16:43:37 发布

阅读量352

点赞数 4

本文链接：https://blog.csdn.net/qq_42671946/article/details/88427563

版权

题意：给一个N，求1到N i^n的异或和。

N<=1.3e7，显然只能O（n)，因为i^n是积性函数，所以可以线性筛O（n)求出1到N的n次方，然而比赛的时候我写了个假的线性筛。。。。

#include <bits/stdc++.h>
using namespace std;
const int mod = 1e9+7;
int n;
bool isprim[13000004];
long long num[13000004];
vector<long long>prim;
long long qk(long long a,long long b)
{
    long long ans=1;
    while(b)
    {
        if(b&1)ans=ans*a%mod;
        b>>=1;
        a=a*a%mod;
    }
    return ans;
}
int main (){
    cin>>n;
    memset(isprim,true, sizeof(isprim));
    for (int i = 2; i <= n; ++i) {
        if(isprim[i]){
            num[i]=qk(i,n);
            prim.push_back(i);
        }
        for (int j = 0; j < prim.size()&&i*prim[j]<=n; ++j) {
            num[prim[j]*i]=num[prim[j]]*num[i]%mod;
            isprim[prim[j]*i]=0;
            if(i%prim[j]==0)break;
        }
    }
    long long ans=1;
    for (int i = 2; i <= n; ++i) {
        ans^=num[i];
    }
    printf("%lld\n",ans);
    return 0;
}