2022 CCPC 威海 G(队友真滴强)

题目
题意: 给定数x，n组查询，每次查询[l,r]中有多少个数满足 gcd(k*x异或x,x) = 1.
思路: 打表。打表发现是有循环的。因为假设x有cnt位，当这cnt位全都移动到左边之后，他们和异或就没有关系了，因为x只有cnt位。只需要k达到2的cnt次方，就会把后cnt位的数都变到左边，之后就重新进行前2的cnt次方的操作，每次k增加1对于kx^x的贡献能在前2的cnt次方操作中找到。简单来说就是循环节是2 ^(logx)，每次kx ^x的增量一样。而后cnt位移到前边恰好是x的倍数，对于gcd来说没有影响。gcd(a,x) = gcd(x,a%x)
代码:

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> PII;
const int N = 2e6+10;
#define int long long
int n,m,k,T;
bool a[N];
int s[N];
ll gcd(ll a,ll b)
{
	if(b==0) return a;
	return gcd(b,a%b);
}
void solve()
{
	int x;
	cin>>x>>n;
	ll base = 0,xx = x;
	while(xx) 
	{
		base++;
		xx /= 2;
	}
	base = (1ll<<base);
	for(int i=1;i<=base;++i) //1e6
	{
		ll t = ((1ll*i*x)^x);
		a[i] = ( gcd(t,x) == 1);
		s[i] = s[i-1]+a[i];
	}
	while(n--)
	{
		ll l,r; cin>>l>>r;
		l--;
		ll ans1 = s[l%base] + 1ll*(l/base)*s[base];
		ll ans2 = s[r%base] + 1ll*(r/base)*s[base];
		ll ans = ans2 - ans1;
		cout<<ans<<"\n";
	}
}
signed main(void)
{
	ios::sync_with_stdio(false),cin.tie(0),cout.tie(0);
	solve();
	return 0; 
}