公约数

Description

给定一个正整数n，在[1,n]的范围内，求出有多少个无序数对(a,b)满足gcd(a,b)=a xor b。

Data Constraint

对于30%的数据满足n<=1000
对于60%的数据满足n<=10^5
对于100%的数据满足n<=10^7

Solution

这道题是结论题。。。
首先我们要知道一些结论：

1. a xor b=c 等价于 a xor c=b 这点很显然啊
2. gcd(a,b)<=a-b (a>b) 这点你可以想一下辗转相除，或者更相减损术，就懂了
3. a xor b>=a-b (a>b) 1 xor 1=0,1-1=0，相等；1 xor 0=1,1-0=1，相等；0 xor 1=1，0-1退位，异或打；0 xor 0=0,0-0=0，相等。所以异或要大于等于减法。

然后根据第一点，我们就可以拿60分。O(n)枚举c（a和b的gcd），O(log n)枚举i，a=c*i，b=a xor c，然后O(log n)判断a和b最大公约数是不是c就好了。这样是O(n $log^2$n)的
再看看第二点和第三点结论，很显然可以得到a-b=c (c=gcd(a,b))，所以跟60分一样枚举c，再枚举i，a=c * i，b=a-c，很明显gcd(a,b)=c（自己想一想就好了，a=i * c，b=(i-1)*c），所以只用判断a xor b是否等于c就好了，这样就少掉了一个log。

Code

Var
    n:Longint;
    Ans:Int64;

Function Gcd(x,y:Longint):Longint;
Begin
    If y=0 Then Exit(x);
    Gcd:=Gcd(y,x Mod y);
End;

Procedure Work();
Var
    a,b,g,i:Longint;
Begin
    For g:=1 To n Do
        For i:=3 To n Div g Do
        Begin
            a:=g*i; b:=a-g;
            If a Xor b=g Then Inc(Ans);
        End;
End;

Begin
    Readln(n);
    Work();
    Writeln(Ans);
End.