看数据范围,贼像数位DP,但是平时都是写的只有一个数字限制的数位DP,数字变成了三个,那就在dfs里面开一个for循环3重枚举。
将十进制全部转为二进制,进行每位的枚举。这也是方便每一位进行异或处理之后的判断。
从高位往低位思考,如果当前枚举的这位数字两两异或大于x的这一位,并且这两个数的异或有限制,那么就不满足题目的第二个条件。这个 “有限制” 等同于枚举数字时的限制。
记得pos等于每个数二进制最长的,这里直接从70位开始。
代码
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N = 100;
const int mod = 1e9 + 7;
int A,B,C,X;
vector<int> a,b,c,x;
int dp[N][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;
if(~dp[pos][la][lb][lc][lab][lac][lbc]) return dp[pos][la][lb][lc][lab][lac][lbc];
int upa = la ? a[pos] : 1,upb = lb ? b[pos] : 1,upc = lc ? c[pos] : 1;
int res = 0;
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) > x[pos]) continue;
if(lac && (i ^ k) > x[pos]) continue;
if(lbc && (j ^ k) > x[pos]) continue;
int pab = lab & ((i ^ j) == x[pos]);
int pac = lac & ((i ^ k) == x[pos]);
int pbc = lbc & ((j ^ k) == x[pos]);
res = (res + dfs(pos-1,la & (i == upa),lb & (j == upb),lc & (k == upc),pab,pac,pbc)) % mod;
}
return dp[pos][la][lb][lc][lab][lac][lbc] = res;
}
int slove(){
a.clear();b.clear();
c.clear();x.clear();
memset(dp,-1,sizeof(dp));
while(a.size() <= 70) a.push_back(A&1),A>>=1;
while(b.size() <= 70) b.push_back(B&1),B>>=1;
while(c.size() <= 70) c.push_back(C&1),C>>=1;
while(x.size() <= 70) x.push_back(X&1),X>>=1;
return dfs(70,1,1,1,1,1,1);
}
signed main(){
IOS
int tt; cin>>tt;
while(tt --){
cin>>A>>B>>C>>X;
cout<<slove()<<endl;
}
return 0;
}