题目链接
(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;
}