华华给月月出题

传送门

题意:

求题中所给式子的答案。

思路:

题目唯一的难点在于想到 f ( x ) = x n f(x) = x^n f(x)=xn是一个完全积性函数,下面简单证明一下:
设 有 f ( x ) = x n , f ( y ) = y n 且 x ≤ n 并 且 y ≤ n , 则 : f ( x ∗ y ) = ( x ∗ y ) n = x n ∗ y n = 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)。 f(x)=xn,f(y)=ynxnyn,f(xy)=(xy)n=xnyn=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;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值