HDU 4901 多校4 经典计数DP

RT

最近不想写博客，累积了一周多的题目，今天趁着周日放假，全部补上吧

dp[i][j]表示前i个数，操作后的值为j的总个数

注意取或不取，有种完全背包的意味。因为数字小于1024，所以异或的结果也绝对不会超过1024，在循环第二维的时候到1024就行了，不要循环多了，反而会错，循环多了 异或值会超，结果爆了数组，产生了奇怪的结果，一开始我就这么莫名其妙的错，还不知道是为什么，就是爆了数组 他也不会提醒你的

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define LL __int64
using namespace std;
const LL M=1000000000+7;
const int N=1050;
const int maxn=1010;
LL dp1[maxn][N];
LL dp2[maxn][N];
LL s1[maxn][N],s2[maxn][N];
int A[maxn];
int main()
{
    int t,n;
    scanf("%d",&t);
    while (t--)
    {
        scanf("%d",&n);
        for (int i=1;i<=n;i++){
            scanf("%d",&A[i]);
        }
        for (int i=0;i<=n+1;i++){
            for (int j=0;j<=N;j++){
                dp1[i][j]=dp2[i][j]=0;
                s1[i][j]=s2[i][j]=0;
            }
        }
        dp1[0][0]=1;
        for (int i=1;i<n;i++){
            //dp1[i][A[i]]=1;
            for (int j=0;j<=1024;j++){
                dp1[i][j]+=dp1[i-1][j];
                dp1[i][j^A[i]]+=dp1[i-1][j];
                if (dp1[i][j]>=M) dp1[i][j]%=M;
                if (dp1[i][j^A[i]]>=M) dp1[i][j^A[i]]%=M;
            }
            for (int j=0;j<=1024;j++){
               s1[i][j^A[i]]=dp1[i-1][j];
               s1[i][j^A[i]]%=M;
            }
        }
        LL ans=0;
        for (int i=n;i>1;i--){
            dp2[i][A[i]]=1;
            for (int j=1024;j>=0;j--){
                dp2[i][j&A[i]]+=dp2[i+1][j];
                dp2[i][j]+=dp2[i+1][j];
                if (dp2[i][j&A[i]]>=M) dp2[i][j&A[i]]%=M;
                if (dp2[i][j]>=M) dp2[i][j]%=M;
            }
        }
        for (int i=1;i<n;i++){
            for (int j=0;j<=1024;j++){
                ans+=s1[i][j]*dp2[i+1][j];
                ans%=M;
            }
        }
        printf("%I64d
",ans%M);
    }
    return 0;
}