UVA12716 GCD等于XOR GCD XOR

UVA12716 GCD等于XOR GCD XOR

TITLE

思路

∵ gcd ⁡ ( i , j ) ≤ i − j ≤ i ⊕ j \because \gcd(i,j) \leq i-j \leq i \oplus j gcd(i,j)ijij
∵ gcd ⁡ ( i , j ) = i ⊕ j \because \gcd(i,j)=i \oplus j gcd(i,j)=ij
∴ gcd ⁡ ( i , j ) = i − j = i ⊕ j \therefore \gcd(i,j)=i-j=i\oplus j gcd(i,j)=ij=ij
gcd ⁡ ( i , j ) = gcd ⁡ ( i , i − j ) = i − j \gcd(i,j)=\gcd(i,i-j)=i-j gcd(i,j)=gcd(i,ij)=ij
∴ i − j ∣ i \therefore i-j|i iji
枚举 j , i ≥ j ∗ 2 , j ∣ i j,i\ge j*2,j|i jij2,ji
i ⊕ j = i − j , a n s [ j ] + 1 i \oplus j=i-j,ans[j]+1 ij=ij,ans[j]+1
a n s i 放 的 是 i 的 答 案 ans_i放的是i的答案 ansii
a n s i = a n s i + a n s i − 1 , i ≥ 2 ans_i=ans_i+ans_{i-1},i\ge2 ansi=ansi+ansi1,i2
前 缀 和 , a n s i 是 1 至 i 的 答 案 前缀和,ans_i是1至i 的答案 ansi1i

CODE

#include<iostream> 
#include<cstdio>
using namespace std;
const int maxn=30000000,half=15000000;
int ans[maxn+10];
int main()
{
    
    
	int n,x;
	for(int i=1;i<=half;++i)
		for(int j=i<<1;j<=maxn;j+=i)
			if((i^j)==j-i)++ans[j];
	for(int i=2;i<=maxn;++i)ans[i]+=ans[i-1];
	scanf("%d",&n);
	for(int i=1;i<=n;++i)
		scanf("%d",&x),printf("Case %d: %d\n",i,ans[x]);
	return 0;
}

猜你喜欢

转载自blog.csdn.net/weixin_46975572/article/details/116939604
xor