HDU5275 数论灵性题

这道题非常常见的数论分块

但是会超时

然后会发现

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;
    }
}