华华给月月出题

108 阅读1分钟

本文已参与「新人创作礼」活动,一起开启掘金创作之路。 传送门

题意:

求题中所给式子的答案。

思路:

题目唯一的难点在于想到f(x)=xnf(x) = x^n是一个完全积性函数,下面简单证明一下: 设有f(x)=xn,f(y)=ynxn并且yn,则:f(xy)=(xy)n=xnyn=f(x)f(y)设有f(x) = x^n,f(y) = y^n且x\leq n 并且 y\leq n,则:f(x*y) = (x*y)^n = x^n*y^n = f(x)*f(y)。 证毕! 然后用欧拉筛求积性函数就好啦。

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define endl '\n'
const int mod = 1e9+7;
ll qpow(ll a,ll b)
{
	ll res = 1;
	while(b)
	{
		if(b&1)res = res*a%mod;
		a = a*a%mod;
		b>>=1;
	}
	return res;
}
ll vis[13000010],pri[13000010];
int cnt;
void oula(int n)
{
	for(int i = 2; i <= n; i++)
	{
		if(!vis[i])pri[++cnt] = i,vis[i] = qpow(i,n);
		for(int j = 1; j <= cnt && pri[j]*i <= n; j++)
		{
			vis[pri[j]*i] = vis[pri[j]]*vis[i]%mod;
			if(pri[j]%i == 0)
			{
				break;
			}
		}
	}
}
int main()
{
	ll n;
	cin>>n;
	ll ans = 0;
	vis[1] = 1;
	oula(n);
	for(int i = 1; i <= n; i++)ans^=vis[i];
	cout<<ans;
}