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=+;
const int N=;
const int maxn=;
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=;i<=n;i++){
            scanf("%d",&A[i]);
        }
        for (int i=;i<=n+;i++){
            for (int j=;j<=N;j++){
                dp1[i][j]=dp2[i][j]=;
                s1[i][j]=s2[i][j]=;
            }
        }
        dp1[][]=;
        for (int i=;i<n;i++){
            //dp1[i][A[i]]=1;
            for (int j=;j<=;j++){
                dp1[i][j]+=dp1[i-][j];
                dp1[i][j^A[i]]+=dp1[i-][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=;j<=;j++){
               s1[i][j^A[i]]=dp1[i-][j];
               s1[i][j^A[i]]%=M;
            }
        }
        LL ans=;
        for (int i=n;i>;i--){
            dp2[i][A[i]]=;
            for (int j=;j>=;j--){
                dp2[i][j&A[i]]+=dp2[i+][j];
                dp2[i][j]+=dp2[i+][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=;i<n;i++){
            for (int j=;j<=;j++){
                ans+=s1[i][j]*dp2[i+][j];
                ans%=M;
            }
        }
        printf("%I64d\n",ans%M);
    }
    return ;
}