【NOIP2016提高A组模拟8.19】公约数

escription

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

Input

输入共一行，一个正整数n。

Output

输出共一行，一个正整数表示答案。

Sample Input

3

Sample Output

1
解释：只有(2,3)满足要求

Data Constraint

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

题目の大意

并不需要啊QvQ

比赛时の想法

好难啊，看到$10^7$就以为是一个O(n)的递推，然后就开始找规律，然后一直都找不到啊QAQ，最厚就打了一发表水了60％QAQ

正解

现在需要知道很多关于a，b，gcd(a,b)=a xor b是否成立，那么显然，当a=b时，这条东西是不成立的（因为a xor b=0）现在我们考虑$a<b$的情况
首先要证明两个东西，①：gcd(a,b)<=a-b,②:a xor b>=a-b
对①的证明，由辗转相除法得：gcd(a,b)=gcd(b,a-b) 如果b=a-b，那么gcd(a,b)=b=a-b,否则要继续辗转下去，那么gcd(a,b)肯定是小于a-b的
对②的证明：我们可以分类讨论（接下来的一段都是用二进制的，贡献是按照二进制的第k位来说的~）
1：现在a的最高位为1，b的最高位为1，此时这一位对a xor b的贡献为0，对a-b贡献为0
2：现在a的最高位为0，b的最高位为1，因为a>b，所以如果有这种情况，肯定有一个更高的位出现3的情况，而如果出现情况3后面的比较已经没有意义了，所以这种情况不用讨论QvQ~
3：现在a的最高位为1，b的最高位为0，此时这一位对a xor b的贡献为1，对a-b贡献为0
4：现在a的最高位为0，b的最高位为0，此时这一位对a xor b的贡献为0，对a-b贡献为0
由上面4种情况我们可以看出：a xor b获得的贡献不小于a-b获得的贡献，即a xor b>=a-b
所以如果要对最终答案有贡献，那么gcd(a,b)=a xor b=a-b
现在我们枚举一个c=a-b，设a=i*c，那么$gcd(a,a-c)=gcd(i*c,(i-1)*c)$，因为i和i-1互质，所以gcd(a,a-c)=c。那么当gcd(a,c)=a-c（也就是gcd(a,b)=b,因为在证明①时我们发现只有当b=a-b时gcd(a,b)=a-b）时就可以把答案加一啦~(≧▽≦)/~

贴代码

好久没有打过这么纯粹的数学problem了QwQ

var
    n,a,b,c,i,j:longint;
    ans:int64;
begin
    readln(n);
    for c:=1 to n do
    begin
        i:=2;
        while i*c<=n do
        begin
            if (i*c) xor c=(i-1)*c then inc(ans);
            inc(i);
        end;
    end;
    writeln(ans);
end.