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=20000010;
#define uint unsigned int
uint seed;
inline uint getnext(){
    seed^=seed<<13;
    seed^=seed>>17;
    seed^=seed<<5;
    return seed;
}
uint a[maxn],ans;
bool vis[maxn];int p[maxn/10],cnt,N;
void solve()
{
    rep(i,2,N){
        if(!vis[i]) p[++cnt]=i;
        for(int j=1;j<=cnt&&i*p[j]<=N;j++){
            vis[i*p[j]]=1;
            if(i%p[j]==0) break;
        }
    }
    rep(i,1,cnt)
     for(int j=2;p[i]*j<=N;j++) a[p[i]*j]+=a[j];
    ans=a[1];
    rep(i,2,N) ans^=(a[i]+a[1]);
}
int main()
{
    cin>>N>>seed;
    rep(i,1,N) a[i]=getnext();
    solve();
    cout<<ans<<endl;
    return 0;
}