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

题目

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

样例输入：
输入共一行，一个正整数n。
3

样例输出：
输出共一行，一个正整数表示答案。
1

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

---

剖解题目

被虐，不说了QAQ~~

---

思路

懵逼+_+|||。

---

解法

30%：暴力枚举两个点，判断。时间：$O(n^2 log_n)$;
60%：对于两个数a，b，如果有gcd(a,b)=a xor b，我们设c=gcd(a,b)=a xor b。因为a xor b=c等价于a xor c=b，所以有gcd(a,a xor c)=c，而c是(a,b)的gcd，所以c是a的约数。
我们只需要枚举一个c和i，让a=i*c。在a<=n的情况里，然后判断gcd(a,c)?=a xor c就行了。时间：$O(n log^2_n)$.
100%:继续上面思路。因为a=b时，是对答案没有贡献的，我们设a>b。其实我们可以证得gcd(a,b)<=a-b,a xor b>=a-b。如果gcd(a,b)=a xor b=c，那么c=a-b。转化一下就是有b=a-c,又b=a xor c,所以a-c=a xor c，所以我们还是枚举c和i，然后判断a-c?=a xor c,记录答案即可。时间：$O(n log_n)$。

---

代码

好短啊，很少见了，300b都不到。

#include<cstdio>
#include<algorithm>
#include<cstdlib>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)

using namespace std;

int n,ans;

int main()
{
    scanf("%d",&n);
    fo(c,1,n-1){
        int i=2;
        while (i*c<=n){
            int a=i*c;
            if ((a^c)==(a-c)) ++ans;
            ++i;
        }
    }
    printf("%d",ans);
}