题意:
有一列n个数(n<=1000,每个数的范围是[0,1024]),现在要从这列数中取出两个集合S与T,其中S中每个元素要在T的左边,现在对S的所有数做异或操作,对T中所有数做与操作,若算出了相同的数,那么当前的S,T就是所需要的,问有多少组S,T符合要求,答案对10^9+7取摸。
题解:
关键是T中的与操作,所知两个数做与操作得到的数不会比这两个数大,所以可以限制了数据的范围,首先用dp[x][y]代表从x开始往后这些数中取一部分做与操作,能得到y的个数.这个过程的复杂度是O(n*1024).再从头开始做异或操作,注意的是两个数异或得到的数是可能比这两个数大的,若计算中算出超过1024可以直接舍弃.那么取值的时候为了防止重复统计,每次都是保证了前半部分会取第i位后能得到某个异或值x的个数乘上从i+1位(不一定取i+1位)往后能通过与操作得到值x的个数.
Program:
#include<iostream>
#include<stdio.h>
#include<math.h>
#include<string.h>
#include<algorithm>
#define MAXN 1050
#define ll long long
#define oo 1000000007
using namespace std;
ll a[MAXN],s[MAXN],p[MAXN],dp[MAXN][MAXN];
int main()
{
ll cases,n,i,x,ans;
freopen("input.txt","r",stdin);
freopen("output.txt","w",stdout);
scanf("%I64d",&cases);
while (cases--)
{
scanf("%I64d",&n);
for (i=1;i<=n;i++) scanf("%I64d",&a[i]);
memset(dp,0,sizeof(dp));
for (i=n;i>=1;i--)
{
memcpy(dp[i],dp[i+1],sizeof(dp[i]));
for (x=0;x<=1024;x++)
dp[i][x&a[i]]=(dp[i][x&a[i]]+dp[i+1][x])%oo;
dp[i][a[i]]=(dp[i][a[i]]+1)%oo;
}
memset(p,0,sizeof(p));
ans=0;
for (i=1;i<=n;i++)
{
memset(s,0,sizeof(s));
for (x=0;x<=1024;x++)
if ((x^a[i])<=1024)
s[x^a[i]]=(s[x^a[i]]+p[x])%oo;
s[a[i]]=(s[a[i]]+1)%oo;
for (x=0;x<=1024;x++)
ans=(ans+s[x]*dp[i+1][x])%oo;
for (x=0;x<=1024;x++) p[x]=(p[x]+s[x])%oo;
}
printf("%I64d\n",ans);
}
return 0;
}