JZOJ5243【GDOI2018模拟8.8】超级绵羊异或 类欧几里得算法

好像没什么人去改这题啊。。。
题意：求$a xor(a+b)xor(a+b*2).....xor (a+b*(n-1))$

考虑计算答案的第x位是否为1
那么对于a+bi，判断$(a+bi)/(1<<k)$（整除）的奇偶性，表示看a+bi的第k位是否为1。最终答案的第k位就看$sigma((a+bi)/(1<<k))i=0..n-1$的奇偶性。
是类欧几里得的模板题。
关于类欧请看http://blog.csdn.net/werkeytom_ftd/article/details/53812718

#include<cstdio>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
typedef long long ll;
const int mo=2;
int n,m;
ll a,b;
typedef long long ll;
inline ll f(ll a,ll b,ll c,ll n)
{
    if (!a)return (n+1)*(b/c)%mo;
    if (a>=c||b>=c)
    {
        return ((a/c)*n*(n+1)/2+(b/c)*(n+1)+f(a%c,b%c,c,n))%mo;
    }
    else
    {
        ll m=1ll*(a*n+b)/c;
        return (1ll*m*n-1ll*f(c,c-b-1,a,m-1)+mo)%mo;
    }
}
int main()
{
    freopen("shxor.in","r",stdin);
    freopen("shxor.out","w",stdout);
    int cas;
    scanf("%d",&cas);
    while(cas--)
    {
        scanf("%d%lld%lld",&n,&a,&b);
        ll ans=0;
        fd(x,60,0)
        {
            ll c=1ll<<x;
            if (f(b,a,c,n-1))ans+=c;
        }
        printf("%lld\n",ans);
    } 
}