hdu 4901 The Romantic Hero

http://acm.hdu.edu.cn/showproblem.php?pid=4901

dp1[i][j]是i参与，异或值为j的个数，x1[i][j]是以i位置向前到1的位置的异或值为j的个数，dp2[i][j]是i参与，&值为j的个数，x2[i][j]是以i位置向前到n的位置的&值为j的个数。

#include <cstdio>
#include <cstring>
#include <algorithm>
#define ll __int64
using namespace std;
const int mod=1000000007;

ll dp1[1100][1100];
ll dp2[1100][1100];
ll x1[1100][1100];
ll x2[1100][1100];
int t,n;
int c[1100];

int main()
{
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i=1; i<=n; i++)
        {
            scanf("%d",&c[i]);
        }
        memset(dp1,0,sizeof(dp1));
        memset(dp2,0,sizeof(dp2));
        memset(x1,0,sizeof(x1));
        memset(x2,0,sizeof(x2));
        for(int i=1; i<=n; i++)
        {
            for(int j=0; j<=1024; j++)
            {
                int m=j^c[i];
                dp1[i][m]=(dp1[i][m]+x1[i-1][j])%mod;
            }
            dp1[i][c[i]]=(dp1[i][c[i]]+1)%mod;
            for(int j=0; j<=1024; j++)
            {
                x1[i][j]=(x1[i-1][j]+dp1[i][j])%mod;
            }
        }
        for(int i=n; i>=1; i--)
        {
            for(int j=0; j<=1024; j++)
            {
                int m=j&c[i];
                dp2[i][m]=(dp2[i][m]+x2[i+1][j])%mod;
            }
            dp2[i][c[i]]=(dp2[i][c[i]]+1)%mod;
            for(int j=0; j<=1024; j++)
            {
                x2[i][j]=(x2[i+1][j]+dp2[i][j])%mod;
            }
        }
        ll ans=0,sum;
        for(int i=1; i<=n-1; i++)
        {
            for(int j=0; j<=1024; j++)
            {
                sum=(dp1[i][j]*x2[i+1][j])%mod;
                ans=(ans+sum)%mod;
            }
        }
        printf("%I64d
",ans);
    }
    return 0;
}