第十八届同济大学程序设计竞赛H.三阳开泰（数位dp+二进制）

题目传送门

题意

给你$T(1\le T\le 100)$组询问，每次询问四个数$A,B,C,X(1\le A,B,C,X\le 10^{18}$，问你有多少个三元组$(a,b,c)$满足$a$ xor $b\le X$,$a$ xor $c\le X$,$b$ xor $c\le X$。

思路

这个数据范围很容易想到数位dp，要是加上下限，问你三个小区间内有多少就更像数位dp了。

我们用$f[pos][la][lb][lc][lab][lac][lbc]$表示状态数组，$pos$表示当前搜到了第$pos$位，$la,lb,lc$分别表示当前搜索到的$a,b,c$是否被$A,B,C$限制，$lab,lac,lbc$则表示$a$ xor $b$,$a$ xor $c$,$b$ xor $c$是否被$X$限制。然后套模板就行了。

$Code$

// Author : ACfunhsl
// Time : 2021/5/22 16:42:57
#define int long long
const int N = 1e6+50;
const int inf = 0x3f3f3f3f;
const int mod = 1e9+7;

int aa[N],bb[N],cc[N],xx[N],cnta,cntb,cntc,cntx;
int f[100][2][2][2][2][2][2];

int dfs(int pos,int la,int lb,int lc,int lab,int lac,int lbc)
{
	if(pos==-1) return 1;
	int &ans = f[pos][la][lb][lc][lab][lac][lbc];
	if(ans!=-1) return ans;
	ans = 0;
	int upa = la?aa[pos]:1;
	int upb = lb?bb[pos]:1;
	int upc = lc?cc[pos]:1;
	int upx = xx[pos];
	for(int i=0;i<=upa;i++)
	{
		for(int j=0;j<=upb;j++)
		{
			for(int k=0;k<=upc;k++)
			{
				if(lab && (i^j)>upx) continue;
				if(lac && (i^k)>upx) continue;
				if(lbc && (j^k)>upx) continue;
				ans = (ans + dfs(pos-1,la&&(i==upa),lb&&(j==upb),lc&&(k==upc),lab&&(i^j==upx),lac&&(i^k==upx),lbc&&(j^k==upx))%mod)%mod;
			}
		}
	}
	return ans;
}

int solve(int a,int b,int c,int x)
{
	cnta = cntb = cntc = cntx = 0;
	mem(aa,0);
	mem(bb,0);
	mem(cc,0);
	mem(xx,0);
	while(a)
		aa[cnta++] = a&1,a>>=1;
	while(b)
		bb[cntb++] = b&1,b>>=1;
	while(c)
		cc[cntc++] = c&1,c>>=1;
	while(x)
		xx[cntx++] = x&1,x>>=1;
	return dfs(63,1,1,1,1,1,1);
}

signed main()
{
	int t;
	cin>>t;
	while(t--)
	{
		int a,b,c,x;
		cin>>a>>b>>c>>x;
		memset(f,-1,sizeof f);
		cout<<solve(a,b,c,x)<<endl;
	}



	return 0;
}