51nod 1301 集合异或和——异或dp

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

好题！看了TJ才会。

因为是不可重集合，所以当然有前 i 个表示A和B都考虑的前 i 个，新加一个讨论放A、放B、不放。

A<B在异或上看就是有一位，它前面的A和B都一样，该位A是0、B是1。该位可以枚举。然后就能dp了。

注意边界细节……和标程对拍真愉快……

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=,M=,mod=1e9+;
int n,m,mi,mj,lm,dp[][M][],ans,bin[];
int calc(int a)
{
    int ret=; while(a) a>>=,ret++; return ret;
}
void ad(int &x,int y)
{
    x=(x+y>=mod?x+y-mod:x+y);
}
int main()
{
    scanf("%d%d",&n,&m);
    mi=max(n,m);  lm=calc(m);  int tmp=calc(mi);
    bin[]=;  for(int i=;i<=tmp;i++) bin[i]=(bin[i-]<<);//tmp not lm!!!
    bin[tmp+]=(bin[tmp]<<);
    mj=bin[tmp+]-;
    for(int t=;t<=lm;t++)
    {
        memset(dp[],,sizeof dp[]);///not dp[0] if mj isn't gu ding
        dp[][][]=;
//        mj=1;
        for(int i=,u,v;i<=mi;i++)
        {
            u=(i&);v=!u;
//            if(i>=bin[mj])mj++;
//            for(int j=0;j<=bin[mj];j++)
            for(int j=;j<=mj;j++)
            {
                dp[u][j][]=dp[v][j][];
                dp[u][j][]=dp[v][j][];
                if(i<=n)//A
                {
                    ad(dp[u][j][],dp[v][j^i][]);
                    ad(dp[u][j][],dp[v][j^i][]);
                }
                if(i<=m)//B
                {
                    bool d=(i&bin[t]);
                    ad(dp[u][j][],dp[v][j^i][^d]);
                    ad(dp[u][j][],dp[v][j^i][^d]);
                }
            }
        }
        int d=(mi&);// mi not n!!!
        for(int i=bin[t];i<bin[t+];i++) ad(ans,dp[d][i][]);
    }
    printf("%d\n",ans);
    return ;
}