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

题目描述：

有$t\le 10000$组询问，每组询问给出$n,a,b\le 1e9$，求：

$$a\oplus(a+b)\oplus(a+2b)……\oplus(a+(n-1)b)$$

解题思路：

考虑直线下整点个数，那么异或和的最低位就是直线$a+b*i,(0<=i<n)$下的整点个数的奇偶性，一般的，第$k$位为直线$(a+b*i)/2^k,(0<=i<n)$下的整点个数的奇偶性。

直线下整点计数要用到类欧几里得算法：
设$f(a,b,c,n)=\sum\limits_{i=0}^n\lfloor \frac{ai+b}{c}\rfloor$
为了方便，以下的除法皆为整除。

当$a\ge n$或$b\ge n$
$f(a,b,c,n)=f(a\%c,b\%c,c,n)+(a/c)(n(n+1)/2)+(b/c)(n+1)$

当$a,b\le n$，设$m=(an+b)/c$
$f(a,b,c,n)=\sum\limits_{i=0}^n\sum\limits_{j=1}^m[(ai+b)/c\ge j]$
$f(a,b,c,n)=\sum\limits_{i=0}^n\sum\limits_{j=0}^{m-1}[(ai+b)/c\ge j+1]$
$f(a,b,c,n)=\sum\limits_{i=0}^n\sum\limits_{j=0}^{m-1}[ai\ge cj+c-b]$
$f(a,b,c,n)=\sum\limits_{i=0}^n\sum\limits_{j=0}^{m-1}[ai>cj+c-b-1]$
$f(a,b,c,n)=\sum_\limits{i=0}^n\sum_\limits{j=0}^{m-1}[i>(cj+c-b-1)/a]$
$f(a,b,c,n)=\sum_\limits{j=0}^m(n-(cj+c-b-1)/a)$
$f(a,b,c,n)=nm-f(c,c-b-1,a,m-1)$
继续递归处理。
可以发现只看第一和第三维表示状态的话，每次由$(a,c)$变成$(c,a\%c)$，因此复杂度与欧几里得算法一致，所以叫类欧几里得算法。

总时间复杂度为$O(Tlog^21e9)$

还有其他形式的类欧几里得可看这里

#include<bits/stdc++.h>
#define ll long long
using namespace std;
int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();c!='-'&&(c<'0'||c>'9');c=getchar());
    if(c=='-')f=-1,c=getchar();
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}
int f(ll a,ll b,ll c,ll n)
{
    if(a>=c||b>=c)return f(a%c,b%c,c,n)^((((ll)n*(n+1)/2)&1)&(a/c&1))^((n+1&1)&(b/c&1));
    if(n<0||!a)return 0;
    ll m=((ll)a*n+b)/c;
    return ((n&1)&(m&1))^f(c,c-b-1,a,m-1);
}
int main()
{
    //freopen("lx.in","r",stdin);
    for(int T=getint();T;T--)
    {
        ll n=getint()-1,b=getint(),a=getint(),ans=0;
        for(int i=0;i<=60;i++)
            if(f(a,b,1ll<<i,n))ans|=(1ll<<i);
        printf("%lld\n",ans);
    }
}