Parity of Tuples (Easy)(贡献+状压dp)

Parity of Tuples (Easy)

题意：

输入$n(1e4),m(10),k(30)$
接下来$n$行，第$i$行输入$a_{i1},a_{i2},\dots ,a_{im}(0\le a_{ij}<2^k)$
求${\sum }_{x=0}^{2^k-1}count(x)3^x\mathrm{m}\mathrm{o}\mathrm{d}10^9+7$
其中$count(x)$表示$i$的个数，对于任意的$1\le j\le m$使$a_{ij}\&x$中位为1的个数为奇数。

题解：

显然可以分析出，不同的$i$,(即不同的行)，它的贡献是相互独立的。于是便对每一行算贡献。
而对于每一行，肯定要从位的角度考虑问题，于是就引入状压$dp$
设$dp[i][j]$表示$0,1,\dots ,i-1$位，状态为$j$(位为1表示奇数)，的贡献。
转移：
$dp[i][j]=dp[i][j]+dp[i-1][j]$不算该位贡献；
$dp[i][j]=dp[i][j]+dp[i-1][j⨁st]\times 3^{2^{i-1}}$算该位的贡献。
最后返回$dp[k][2^m-1]$

代码：

#include<bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
int n,m,k,a[19];
int dp[34][1<<11];
int base[39];
int solve(){
    for(int i=1;i<=m;i++)cin>>a[i];
    for(int i=0;i<=k;i++)for(int j=0;j<(1<<m);j++)dp[i][j]=0;
    dp[0][0]=1;
    for(int i=1;i<=k;i++){
        int zt=0;
        for(int j=1;j<=m;j++)if(a[j]&(1<<(i-1)))zt|=(1<<j-1);
        for(int j=0;j<(1<<m);j++){
            dp[i][j]=(dp[i][j]+dp[i-1][j])%mod;
            dp[i][j]=(dp[i][j]+1ll*dp[i-1][j^zt]*base[i-1])%mod;
        }
    }
    return dp[k][(1<<m)-1];
}
int main(){
    //freopen("tt.in","r",stdin),freopen("tt.out","w",stdout);
    base[0]=3;
    for(int i=1;i<=31;i++)base[i]=1ll*base[i-1]*base[i-1]%mod;
    while(cin>>n>>m>>k){
        int ans=0;
        while(n--)ans=(ans+solve())%mod;
        cout<<ans<<endl;
    }
    return 0;
}