luoguP5495：Dirichlet 前缀和

题意：给定数组a[]的生成方式，然后b[i]=∑a[j]  ，(i%j==0)，求所有b[i]的异或和。所有运算%2^32;

思路：高维前缀和的思想，先筛出所有素数，然后把每个素数当成一维，那么分开考虑即可。复杂度O(NloglogN)；

如果有这一维就加进去就可以了~神奇。

#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
const int maxn=;
#define uint unsigned int
uint seed;
inline uint getnext(){
    seed^=seed<<;
    seed^=seed>>;
    seed^=seed<<;
    return seed;
}
uint a[maxn],ans;
bool vis[maxn];int p[maxn/],cnt,N;
void solve()
{
    rep(i,,N){
        if(!vis[i]) p[++cnt]=i;
        for(int j=;j<=cnt&&i*p[j]<=N;j++){
            vis[i*p[j]]=;
            if(i%p[j]==) break;
        }
    }
    rep(i,,cnt)
     for(int j=;p[i]*j<=N;j++) a[p[i]*j]+=a[j];
    ans=a[];
    rep(i,,N) ans^=(a[i]+a[]);
}
int main()
{
    cin>>N>>seed;
    rep(i,,N) a[i]=getnext();
    solve();
    cout<<ans<<endl;
    return ;
}