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的个数。
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
1 #include <cstdio> 2 #include <cstring> 3 #include <algorithm> 4 #define ll __int64 5 using namespace std; 6 const int mod=1000000007; 7 8 ll dp1[1100][1100]; 9 ll dp2[1100][1100]; 10 ll x1[1100][1100]; 11 ll x2[1100][1100]; 12 int t,n; 13 int c[1100]; 14 15 int main() 16 { 17 scanf("%d",&t); 18 while(t--) 19 { 20 scanf("%d",&n); 21 for(int i=1; i<=n; i++) 22 { 23 scanf("%d",&c[i]); 24 } 25 memset(dp1,0,sizeof(dp1)); 26 memset(dp2,0,sizeof(dp2)); 27 memset(x1,0,sizeof(x1)); 28 memset(x2,0,sizeof(x2)); 29 for(int i=1; i<=n; i++) 30 { 31 for(int j=0; j<=1024; j++) 32 { 33 int m=j^c[i]; 34 dp1[i][m]=(dp1[i][m]+x1[i-1][j])%mod; 35 } 36 dp1[i][c[i]]=(dp1[i][c[i]]+1)%mod; 37 for(int j=0; j<=1024; j++) 38 { 39 x1[i][j]=(x1[i-1][j]+dp1[i][j])%mod; 40 } 41 } 42 for(int i=n; i>=1; i--) 43 { 44 for(int j=0; j<=1024; j++) 45 { 46 int m=j&c[i]; 47 dp2[i][m]=(dp2[i][m]+x2[i+1][j])%mod; 48 } 49 dp2[i][c[i]]=(dp2[i][c[i]]+1)%mod; 50 for(int j=0; j<=1024; j++) 51 { 52 x2[i][j]=(x2[i+1][j]+dp2[i][j])%mod; 53 } 54 } 55 ll ans=0,sum; 56 for(int i=1; i<=n-1; i++) 57 { 58 for(int j=0; j<=1024; j++) 59 { 60 sum=(dp1[i][j]*x2[i+1][j])%mod; 61 ans=(ans+sum)%mod; 62 } 63 } 64 printf("%I64d\n",ans); 65 } 66 return 0; 67 }