华华给月月出题

题目链接

(a*b)^c = (a^c)*(b^c)，我们将合数素数幂分解，然后我们只要处理每个素数的幂，就可以处理出所以合数的幂。

代码：

#include<bits/stdc++.h>
using namespace std;
using ll = long long;
#define fi first
#define se second
const ll mod=1e9+7;
const int N=2e5+10;
//#define int ll
const double eps=1e-6;


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;
}
const int Maxn=2e7+10;
 
int n;
//int phi[Maxn];//记录数的约数个数(欧拉函数) 
bool vis[Maxn];//记录数字是否访问 
int prime[Maxn];//保存素数
int f[Maxn];
void init(){
	vis[1]=1;//1不是素数 
	//phi[1]=1;
	f[1]=1;
	for(int i=2;i<=Maxn;i++)
	{
		if(!vis[i])//没有被访问，也就是没有被筛掉,说明是素数 
		{
			vis[i]=!vis[i];
			prime[++prime[0]]=i;
            //phi[i]=i-1;
            f[i]=qpow(i,n);
		}
		for(int j=1;j<=prime[0]&&i*prime[j]<=Maxn;j++)
		{
			vis[i*prime[j]]=true;
			f[i*prime[j]]=f[i]*f[prime[j]]%mod;
			if(i%prime[j]==0)//a%b==0，那么phi[a*b]=b*phi[a] 
			{
				//phi[i*prime[j]]=phi[i]*prime[j];
				break;
			}
			else ;//phi[i*prime[j]]=phi[i]*(prime[j]-1);//两者互素 
		}
	}
} 

void solve(){
	cin>>n;
	init();
	ll ans=0;
	for(int i=1;i<=n;i++){
		//cout<<f[i]<<endl;
		ans^=f[i];
	}	
	cout<<ans;
}

signed main(){
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    int t=1;
    //cin>>t;
	while (t--){
		solve();
	}
	return 0;
}