【Luogu】P4067储能表（数位DP）

　　题目链接

　　好的

　　看到这题之后我一直在想反演，然后想不出来，一度以为自己脑子有问题

　　然后我脑子真的有问题，这题tm根本就不是反演

　　设f[i][j][k][l]表示现在已经DP到从高位往低数的第i位，有没有碰到n的上界，有没有碰到m的上界，有没有碰到k的上界

　　然后记忆化DFS搞一搞，把方案数和异或和都记录下来，最后减掉输出

　　

#include<cstdio>
#include<cstring>
#include<cctype>
#include<cstdlib>
#include<algorithm>
#define maxn 50
using namespace std;
inline long long read(){
    long long num=,f=;
    char ch=getchar();
    while(!isdigit(ch)){
        if(ch=='-')    f=-;
        ch=getchar();
    }
    while(isdigit(ch)){
        num=num*+ch-'';
        ch=getchar();
    }
    return num*f;
}

long long n,m,e,mod;

struct Ans{
    long long xum,sum;
}f[][][][];
bool vis[][][][];

long long Maxi;

void dfs(long long len,long long nown,long long nowm,long long nowk){
    if(len>Maxi){
        f[len][nown][nowm][nowk]=(Ans){,};
        return;
    }
    if(vis[len][nown][nowm][nowk])    return;
    Ans &ret=f[len][nown][nowm][nowk];
    vis[len][nown][nowm][nowk]=;
    long long mn=(n>>(Maxi-len))&,mm=(m>>(Maxi-len))&,mk=(e>>(Maxi-len))&;
    for(long long i=;i<=(nown?mn:);++i)
        for(long long j=;j<=(nowm?mm:);++j){
            if(nowk&&(i^j)<mk)    continue;
            dfs(len+,nown&&(i==mn),nowm&&(j==mm),nowk&&((i^j)==mk));
            Ans now=f[len+][nown&&(i==mn)][nowm&&(j==mm)][nowk&&(i^j)==mk];
            ret.sum=(ret.sum+now.sum)%mod;
            ret.xum=(ret.xum+(1ll<<(Maxi-len))*(i^j)%mod*now.sum+now.xum)%mod;
        }
    return;
}

int main(){
    long long T=read();
    while(T--){
        n=read()-;m=read()-;e=read();mod=read();
        memset(vis,,sizeof(vis));    Maxi=;
        memset(f,,sizeof(f));
        long long nn=n,nm=m,ne=e;long long cnt=;
        while(nn){    cnt++;    nn>>=;}
        Maxi=max(Maxi,cnt);    cnt=;
        while(nm){    cnt++;    nm>>=;}
        Maxi=max(Maxi,cnt);    cnt=;
        while(ne){    cnt++;    ne>>=;}
        Maxi=max(Maxi,cnt);
        dfs(,,,);
        Ans now=f[][][][];
        long long ans=(now.xum%mod-(e%mod)*now.sum%mod+mod)%mod;
        printf("%lld\n",ans);
    }
    return ;
}