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

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=2005,M=4100,mod=1e9+7;
int n,m,mi,mj,lm,dp[2][M][2],ans,bin[25];
int calc(int a)
{
    int ret=0; while(a) a>>=1,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[1]=1;  for(int i=2;i<=tmp;i++) bin[i]=(bin[i-1]<<1);//tmp not lm!!!
    bin[tmp+1]=(bin[tmp]<<1);
    mj=bin[tmp+1]-1;
    for(int t=1;t<=lm;t++)
    {
        memset(dp[0],0,sizeof dp[0]);///not dp[0] if mj isn‘t gu ding
        dp[0][0][0]=1;
//        mj=1;
        for(int i=1,u,v;i<=mi;i++)
        {
            u=(i&1);v=!u;
//            if(i>=bin[mj])mj++;
//            for(int j=0;j<=bin[mj];j++)
            for(int j=0;j<=mj;j++)
            {
                dp[u][j][0]=dp[v][j][0];
                dp[u][j][1]=dp[v][j][1];
                if(i<=n)//A
                {
                    ad(dp[u][j][0],dp[v][j^i][0]);
                    ad(dp[u][j][1],dp[v][j^i][1]);
                }
                if(i<=m)//B
                {
                    bool d=(i&bin[t]);
                    ad(dp[u][j][0],dp[v][j^i][0^d]);
                    ad(dp[u][j][1],dp[v][j^i][1^d]);
                }
            }
        }
        int d=(mi&1);// mi not n!!!
        for(int i=bin[t];i<bin[t+1];i++) ad(ans,dp[d][i][1]);
    }
    printf("%d\n",ans);
    return 0;
}