SDOI2016 储能表

传送门

要求出答案，首先要知道异或值大于$k$的异或和，还要知道方案数，因为每个要减去$k$。
那么用$pair$记录答案，$first$为方案数，$second$为异或和。
数位$DP$枚举$i$和$j$的每一位，需要满足$i⩽n,j⩽m$且异或出的数大等于$k$。
$f[pos][mx\_n][mx\_m][mn\_k]$表示从$pos$开始到最后一位的答案，后面三维记录三个边界。转移就是当前位置从下一个位置得到转移。

#include<bits/stdc++.h>
#define cs const
#define fi first
#define se second
#define re register
#define ll long long
#define mp std::make_pair
#define pii std::pair<int,int>
cs int L=63;
namespace IO{
	cs int Rlen=1<<22|1;
	char buf[Rlen],*p1,*p2;
	inline char gc(){return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;}
	template<typename T>
	inline T get(){
		char ch;T x;
		while(!isdigit(ch=gc()));x=ch^48;
		while(isdigit(ch=gc())) x=((x+(x<<2))<<1)+(ch^48);
		return x;
	}
	inline int gi(){return get<int>();}
	inline ll gl(){return get<ll>();}
}	
using namespace IO;
int T,mod,N[L],M[L],K[L],power[L];ll n,m,k;
inline int add(int x,int y){return (x+y>=mod)?x+y-mod:x+y;}
inline int dec(int x,int y){return (x-y<0)?x-y+mod:x-y;}
inline int mul(int x,int y){return 1ll*x*y%mod;}
inline void Add(int &x,int y){x=(x+y>=mod)?x+y-mod:x+y;}
inline void Dec(int &x,int y){x=(x-y<0)?x-y+mod:x-y;}
inline void Mul(int &x,int y){x=1ll*x*y%mod;}
int vis[L][2][2][2];pii f[L][2][2][2],ans;
inline pii dp(int pos,int mx_n,int mx_m,int mn_k){
	if(pos==L) return mp(1,0);
	if(vis[pos][mx_n][mx_m][mn_k]) return f[pos][mx_n][mx_m][mn_k];
	vis[pos][mx_n][mx_m][mn_k]=1;
	int up_i=mx_n?N[pos]:1,up_j=mx_m?M[pos]:1,dn_k=mn_k?K[pos]:0;
	for(int re i=0;i<=up_i;++i){
		for(int re j=0;j<=up_j;++j){
			if(mn_k&&((i^j)<dn_k)) continue;
			pii F=dp(pos+1,mx_n&&(i==up_i),mx_m&&(j==up_j),mn_k&&((i^j)==dn_k));
			Add(f[pos][mx_n][mx_m][mn_k].fi,F.fi);
			Add(f[pos][mx_n][mx_m][mn_k].se,add(F.se,(i^j)?mul(power[L-1-pos],F.fi):0));
		}
	}return f[pos][mx_n][mx_m][mn_k];
}
inline void init(){
	power[0]=1;
	for(int re i=1;i<L;++i) power[i]=mul(power[i-1],2);
	for(int re i=0;i<L;++i)
		N[i]=(n&(1ll<<(L-i-1)))?1:0,
		M[i]=(m&(1ll<<(L-i-1)))?1:0,
		K[i]=(k&(1ll<<(L-i-1)))?1:0;
	memset(vis,0,sizeof(vis));
	memset(f,0,sizeof f);
}
int main(){
//	freopen("2536.in","r",stdin);
	T=gi();while(T--){
		n=gl()-1,m=gl()-1,k=gl(),mod=gi(),init(),ans=dp(0,1,1,1);
		printf("%d\n",dec(ans.se,mul(k%mod,ans.fi)));
	}
}