Problem:
JOI 旗 (JOI Flag)

問題

情報オリンピック日本委員会では，今年の日本情報オリンピック (JOI) を宣伝するために， JOI のロゴをモチーフにした旗を作ることになった．旗は「良い旗」でなければならない．「良い旗」とは，アルファベットの J, O, I のいずれかの文字を，縦 M 行，横 N 列の長方形状に並べたもので， J, O, I が以下の図のように (すなわち，J の右隣に O が，その J の下隣に I が) 並んでいる箇所が旗の少なくとも 1 か所にあるものである．

以下の図に，「良い旗」の例を 2 つ示す．

以下の図に，「良い旗」ではない例を 2 つ示す．

いま M, N の値，および旗の一部の場所について J, O, I のどの文字にするかが決まっており，入力としてその情報が与えられる．考えられる「良い旗」は何通りあるかを計算し，その数を 100000 (=105) で割った余りを出力するプログラムを作成せよ．

入力

入力は 1 + M 行からなる．

1 行目には旗の大きさを表す 2 つの整数 M, N (2 ≤ M ≤ 20, 2 ≤ N ≤ 20) が空白で区切られて書かれている．

1 + i 行目 (1 ≤ i ≤ M) には，N 文字からなる文字列が書かれている．各文字は J, O, I, ? のいずれかであり， j 文字目 (1 ≤ j ≤ N) が J, O, I のいずれかである場合は i 行 j 列の場所の文字がそれぞれ J, O, I に決まっていること， ? である場合はまだ決まっていないことを表す．

出力

考えられる「良い旗」の個数を 100000 (=105) で割った余りを 1 行で出力せよ．

入出力例

入力例 1

2 3
??O
IIJ

出力例 1

4

入力例 1 においては，以下の図の 4 通りの「良い旗」が考えられる．

入力例 2

2 2
??
??

出力例 2

3

入力例 3

3 3
??I
???
O?J

出力例 3

53

入力例 4

5 4
JOI?
????
????
????
?JOI

出力例 4

28218

入力例 4 においては，「良い旗」は 2428218 通り考えられるので，それを 100000 (=105) で割った余りである 28218 を出力する．

問題文と自動審判に使われるデータは、情報オリンピック日本委員会が作成し公開している問題文と採点用テストデータです。

Solution in C++ (judge verdict: Wrong Answer)

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <climits>
#include <iostream>
#include <vector>
#include <string>
#include <sstream>
#include <algorithm>
#include <utility>
#include <set>
#include <map>
#include <stack>
#include <queue>
#include <deque>
#include <functional>
#define FI first
#define SE second
#define PB push_back
#define MT make_tuple
using namespace std;
static const int MBIT = 699000;
static const int mod = 100000;

short N, M;
short fld[20 * 20];
int dp[2][3][MBIT];
vector<int> BIT;
int ans = 0;
int all = 1;
int main() {
    cin.tie(0);
    ios_base::sync_with_stdio(0);
    cin >> N >> M;
    for (int i = 0; i < N; ++i) {
        for (int j = 0; j < M; ++j) {
            char c;
            cin >> c;
            if (c == '?') fld[i*M + j] = 3, all *= 3;
            if (c == 'J')fld[i*M + j] = 0;
            if (c == 'O')fld[i*M + j] = 1;
            if (c == 'I')fld[i*M + j] = 2;
            all %= mod;
        }
    }
    if (fld[0] != 3)dp[0][fld[0]][0] = 1;
    else {
        dp[0][0][0] = 1;
        dp[0][1][0] = 1;
        dp[0][2][0] = 1;
    }
    int maxbit = 0;
    for (int i = 0; i < M; ++i) {
        maxbit += pow(2, i);
    }
    for(int i = 0; i <= maxbit; i++){
        bool isbit=false;
        for(int j = 0 ;j < M-1; ++j){
            int bit1 = 1<<j,bit2 = 1<<(j+1);
            if((i&bit1)&&(i&bit2)){
                isbit=true;
                break;
            }
        }
        if(!isbit)BIT.PB(i);
    }
    int bit = 1 << M;
    int bit1 = 1 << (M - 1);
    int bit2 = 1 << (M - 2);
    for (int i = 0; i < M * N - 1; ++i) {
        int p1 = (i + 1) & 1;
        int p2 = i & 1;
        for (int j = 0; j < 3; ++j) {
            for (int k = 0; k < BIT.size(); ++k) {
                int tk = BIT[k];
                int t1 = tk << 1;
                if (t1&bit)t1 -= bit;
                int t2 = t1 + 1;
                if ((tk&bit1) && j == 2 && (i + 1) % M != 0)continue;
                if (fld[i + 1] == 3) {
                    dp[p1][0][t1] = (dp[p1][0][t1] + dp[p2][j][tk]) % mod;
                    if (!(tk&bit2) || (i + 2) % M == 0)dp[p1][2][t1] = (dp[p1][2][t1] + dp[p2][j][tk]) % mod;
                    if (j == 0 && (i + 1) % M != 0)dp[p1][1][t2] = (dp[p1][1][t2] + dp[p2][j][tk]) % mod;
                    else dp[p1][1][t1] = (dp[p1][1][t1] + dp[p2][j][tk]) % mod;
                }
                else if (fld[i + 1] == 0) {
                    dp[p1][0][t1] = (dp[p1][0][t1] + dp[p2][j][tk]) % mod;
                }
                else if (fld[i + 1] == 1) {
                    if (j == 0 && (i + 1) % M != 0)dp[p1][1][t2] = (dp[p1][1][t2] + dp[p2][j][tk]) % mod;
                    else dp[p1][1][t1] = (dp[p1][1][t1] + dp[p2][j][tk]) % mod;
                }
                else if (fld[i + 1] == 2) {
                    if (!(tk&bit2) || (i + 2) % M == 0)dp[p1][2][t1] = (dp[p1][2][t1] + dp[p2][j][tk]) % mod;
                }
                dp[p2][j][tk] = 0;
            }
        }
    }
    for (int j = 0; j < 3; ++j) {
        for (int k = 0; k < BIT.size(); ++k) {
            ans += dp[(M  * N - 1) & 1][j][BIT[k]];
            ans %= mod;
        }
    }
    cout << (all%mod - ans%mod + mod) % mod << endl;
    return 0;
}