SSLOJ 1461.异或

题目

题目描述
SarvaTathagata是个神仙,一天他在研究数论时,书上有这么一个问题:求不超过n两两的数的gcd。
SarvaTathagata这么神仙的人当然觉得这个是sb题啦。学习之余,他还发现gcd的某一个特别好的性质:如果有两个数i,j满足gcd(i,j)=i ^ j(这里的^为c++中的异或)的话,那么这两个数组成的数对(i,j)就是一个nb的数对(这里认为(i,j)和(j,i)为相同的,并不需要计算2次)。
当然,SarvaTathagata并不会只满足于判断一个数对是否nb,他还想知道满足两个数都是不超过n并且nb的数对有多少个。
由于SarvaTathagata实在是太神仙了,他认为这种题实在是太简单了。于是他找到了你,看看你是否能解决这个问题。

输入 共一行一个整数n,含义如题所述。 输出 一行一个整数,表示nb的数对的个数。 输入样例

样例输入1:
12
样例输入2:
123456

样例输出1:
8
样例输出2:
214394

说明 提示
样例1中共有八对,分别是:

分析

注意,在这里a>b, “ ^ ”表示异或(xor)

容易证明: 当gcd(a,b) = a ^ b时,gcd(a,b) = a-b = a ^ b:
1.gcd(a,b) <= a - b(应该很容易理解):
设gcd(a,b) = k , a' = a / k , b' = b / k
a - b = a' k + b' k = k(a' - b')
又∵ a > b
∴ a' > b'
又∵ a , b ∈ Z
∴a' - b' >= 1
∴k <= a-b
即gcd(a,b) <= a-b
2.a ^ b >= a - b
首先,异或是按位的,没有进退位
对比一下:

0 xor 1 = 1 0 - 1 = 1 (需要从前面借1 , 小于异或)
1 xor 0 = 1 1 - 0 = 0(相等于异或)
1 xor 1 = 0 1-1 = 0(相等于异或)
0 xor 0 = 0 0 - 0 = 0(相等于异或)
所以,a ^ b >= a - b
综上,gcd(a,b) <= a - b <= a ^ b
又∵gcd(a,b) = a ^ b
∴gcd(a,b) = a-b = a ^ b

=============== 分割线 ==========================
设c = a-b 则b = a-c
gcd(a,b) = gcd(a,a-b) = gcd(a,c)(详见更相减损术

重点来了:
当a-c = a^c时
b = a^c
b ^ b = a ^ b ^ c
0 = a ^ b ^ c
即a ^ b = c
又∵c = a - b = gcd(a,b)
∴gcd(a,b) = a - b = a ^ b

因此,只要满足a-c = a ^ c,就有gcd(a,b) = a ^ b,我们枚举a,c即可
又∵gcd(a,b) = gcd(a,c) = a - b = c
∴a是c的整数倍

代码

#include <iostream>
using namespace std;
int main(){
	int n , ans = 0;
	cin >> n;
	for(int c = 1 ; c <= n ; c++)//枚举c
		for(int a = (c << 1) ; a <=n ; a += c)//枚举c的倍数,即a
			if(a - c == (a ^ c)){
				ans += 1;
			}
	cout << ans;
	return 0;
}
posted @ 2020-11-11 10:38  追梦人1024  阅读(62)  评论(0编辑  收藏  举报