这道题非常常见的数论分块
但是会超时
然后会发现
i小的时候,在等差数列异或和的时候很多计算都是多余的
所以可以从这里入手
前一半暴力
其实这个问题也可以从那个等差数列异或和函数入手
把他的遍历上限改改
这道题就能A
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll M=3e7+7;
ll f(ll a,ll b,ll c,ll n){
if (n==0) return 0;
ll ans=(b/c)*n*(n-1)/2+(a/c)*n;
ans=ans+f((b*n+a)%c,c,b%c,((b%c)*n+(a%c))/c);
return ans;
}
ll F(ll x,ll y,ll z){
ll ans=0;
for (ll u=0;u<40;u++)
ans=ans|((f(x,z,(1ll<<u),(y-x)/z+1)&1)<<u);
return ans;
}
int main(){
int T;
cin>>T;
while(T--){
ll n;
cin>>n;
ll ans=0;
ll c;
for(ll i=1;i<=min(M,n);++i)
ans=ans^(n%i);
for(ll i=M+1;i<=n;i=c+1){
c=n/(n/i);
ans=ans^F(n%c,n%i,n/i);
}
cout<<ans<<endl;
}
}