UVALive - 8144 Sacred Scarecrows DP + FMT 未解决

简略题意：$R*C$的庄稼地，有些地方已经种了庄稼了。现在需要放置一些稻草人，使得满足以下两个条件：
$1$. 所有行都包含稻草人。
$2$. 相邻的两列至少包含两个稻草人。

这题目前还没有$AC$，但是值得记录。

注意到行很少，因此可以状压进行$DP$。
令$dp[i][0/1][j]$代表，前i列的稻草人存放了那些行，当前列是否有稻草人。

用$**$代表集合并卷积，那么转移有如下：
$1. dp[i][0] = dp[i-1][1] ** {g[0]}$
$2. dp[i][1] = dp[i-1][0] ** {g/g[0]} + dp[i-1][0] ** {g/g[0]} = (dp[i-1][0] + dp[i-1][1]) ** {g/g[0]}$
障碍物的处理直接把对应卷上的$g$给置零即可，所以可以用$FMT$来优化到$R * (2^R) * C$。
经过用$assert$得出测试数据至少有$20$组…上述算法无法通过。

qls给出了如下解法：
对行容斥，可以知道每列有多少空地，然后 $O(c)$的$dp$.
容斥用 $dfs$ 枚举集合可以 $2^r$，可以做到 $O(2^R * C)$.
暂时还没太理解，先想一想。

#define others
#ifdef poj
#include <iostream>
#include <cstring>
#include <cmath>
#include <cstdio>
#include <algorithm>
#include <vector>
#include <string>
#include <map>
#include <set>
#endif // poj
#ifdef others
#include <bits/stdc++.h>
#endif // others
//#define file
#define all(x) x.begin(), x.end()
using namespace std;
#define eps 1e-8
const double pi = acos(-1.0);

typedef long long LL;
typedef unsigned long long ULL;
void umax(int &a, int b) {
    a = max(a, b);
}
void umin(int &a, int b) {
    a = min(a, b);
}
int dcmp(double x) {
    return fabs(x) <= eps?0:(x > 0?1:-1);
}
void file() {
    freopen("data_in.txt", "r", stdin);
    freopen("data_out.txt", "w", stdout);
}

const LL mod = 1e9+7;

int r, c, lim;
LL dp[2][2][20000], f[20000], g[20000], h[20000];
char G[16][1100];

LL add(LL a, LL b) {
    a += b;
    if(a >= mod) a -= mod;
    return a;
}

LL mul(LL a, LL b) {
    LL sum = 1ll * a * b;
    if(sum >= mod) a %= mod;
    return sum;
}

LL sub(LL a, LL b) {
    a -= b;
    if(a < 0) a += mod;
    return a;
}

void FMT(int col, int on) {
    for(int i = 0; i < lim; i++) g[i] = 0;
    if(on) {
        for(int i = 1; i < lim; i++) {
            bool isok = 1;
            for(int j = 0; j < r; j++) {
                if((i & (1 << j)) && G[j][col] == 'v') isok = 0;
            }
            if(isok) g[i] = 1;
        }
    } else {
        g[0] = 1;
    }
    for(int i = 0; i < r; i++)
        for(int j = 0; j < lim; j++)
            if((j >> i) & 1)
                f[j] = add(f[j], f[j^(1<<i)]);
    for(int i = 0; i < r; i++)
        for(int j = 0; j < lim; j++)
            if((j >> i) & 1)
                g[j] = add(g[j], g[j^(1<<i)]);
    for(int i = 0; i < lim; i++)
        h[i] = mul(f[i], g[i]);
    for(int i = 0; i < r; i++)
        for(int j = 0; j < lim; j++)
            if((j>>i)&1)
                h[j] = sub(h[j], h[j^(1<<i)]);
}

void cpy(LL x[], LL y[]) {
    for(int i = 0; i < lim; i++)
        x[i] = y[i];
}

int main() {
//    file();
int ff = 0;
    while(~scanf("%d%d", &r, &c)) {
        ff++;
//    assert(ff <= 30);
        lim = 1 << r;
        for(int i = 0; i < r; i++)
            for(int j = 1; j <= c; j++)
                scanf(" %c", &G[i][j]);
        for(int i = 0; i < 2; i++)
            for(int j = 0; j < 2; j++)
                for(int k = 0; k < lim; k++) dp[i][j][k] = 0;
        LL ans = 0;
        dp[0][1][0] = 1;
        for(int i = 1; i <= c; i++) {
            int now = i & 1, pre = now ^ 1;
            for(int j = 0; j < lim; j++) dp[now][0][j] = dp[now][1][j] = 0;
            for(int j = 0; j < lim; j++)
                dp[now][0][j] = dp[pre][1][j], dp[pre][1][j] = add(dp[pre][1][j], dp[pre][0][j]);

            cpy(f, dp[pre][1]);
            FMT(i, 1);
            cpy(dp[now][1], h);
        }
        cout<<(dp[c&1][0][(1<<r)-1] + dp[c&1][1][(1<<r)-1])%mod<<'\n';
    }
    return 0;
}