我是靠谱客的博主 冷酷心情,最近开发中收集的这篇文章主要介绍华华给月月出题题意:思路:,觉得挺不错的,现在分享给大家,希望可以做个参考。

概述

传送门

题意:

求题中所给式子的答案。

思路:

题目唯一的难点在于想到 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且xleq n 并且 yleq 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;
}

最后

以上就是冷酷心情为你收集整理的华华给月月出题题意:思路:的全部内容,希望文章能够帮你解决华华给月月出题题意:思路:所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供,作为学习参考使用,或来自网络收集整理,版权属于原作者所有。
点赞(42)

评论列表共有 0 条评论

立即
投稿
返回
顶部