51Nod 1301 集合异或和 —— 异或DP

题目：http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1301

参考博客：https://blog.csdn.net/qq_33229466/article/details/78418414

原来是DP。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
int const maxn=,mod=1e9+;
int n,m,w,bin[];
ll ans,f[][][];
int get(int x){int ret=; while(x)ret++,x>>=; return ret;}
void init()
{
    bin[]=;
    for(int i=;i<;i++)bin[i]=(bin[i-]<<);
}
int main()
{
    init();
    scanf("%d%d",&n,&m); w=get(max(n,m));
    for(int i=,p,q;i<=w;i++)
    {
        memset(f,,sizeof f);
        p=; f[][][]=;
        int lim=;
        for(int j=;j<=max(n,m);j++)//枚举数
        {
            p^=; q=!p;
            if(j>=bin[lim])lim++;
//            int tmp=(j&bin[i]); //j&bin[i]不仅是0或1！！！
            bool tmp=(j&bin[i]);
            for(int s=;s<bin[lim];s++)//<
            {
                f[p][s][]=f[q][s][]; f[p][s][]=f[q][s][];
                for(int k=;k<=;k++)
                {
                    if(j<=m)(f[p][s][k]+=f[q][s^j][k^tmp])%=mod;//给B
                    if(j<=n)(f[p][s][k]+=f[q][s^j][k])%=mod;//给A
                }
            }
        }
        for(int s=bin[i];s<bin[i+];s++)ans=(ans+f[p][s][])%mod;
    }
    printf("%lld\n",ans);
    return ;
}