1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
#include <bits/stdc++.h>
using namespace std;

int n, m, k;
string s;
int p[5][2];    //0 - parzyste, 1 - nieparzyste
long long pt[10][5];
long long sil[21];
long long DN[37][9];
bool even, now, evenend;
int a;
long double up = 1, down = 0;
long double res;

long long dn(long long x1, long long x2){
    long long w = 1;
    for(int i = x2+1; i <= x1; i++){
        w *= i;
    }
    for(int i = 1; i <= x1 - x2; i++){
        w /= i;
    }
    return w;
}

int main(){
    ios_base::sync_with_stdio(0);
    cin.tie(0);
    cout.tie(0);

    pt[0][1] = 1;
    pt[0][2] = 1;
    pt[0][3] = 1;
    pt[0][4] = 1;
    for(int i = 1; i < 10; i++){
        pt[i][1] = 1;
        pt[i][2] = pt[i-1][2]*2;
        pt[i][3] = pt[i-1][3]*3;
        pt[i][4] = pt[i-1][4]*4;
    }
    sil[0] = 1;
    for(int i = 1; i < 21; i++)
        sil[i] = sil[i-1]*i;
    for(int i = 0; i < 37; i++){
        for(int j = 0; j < min(9, i+1); j++){
            DN[i][j] = dn(i, j);
        }
    }

    cin >> n >> m;
    for(int i = 1; i <= n; i++){
        cin >> s;
        for(int j = 0; j < m; j++){
            if(s[j] == 'O'){
                k++;
                if((j+i+1)%2 == 1 && even == 1)
                    even = 0;
                else if((j+i+1)%2 == 1 && even == 0)
                    even = 1;
            }
        }
    }
    for(int i = 1; i <= n; i++){
        cin >> s;
        for(int j = 0; j < m; j++){
            now = (j+i+1)%2;
            if(i==1 || i==n){
                if(j==0 || j==m-1)
                    a = 2;
                else
                    a = 3;
            } else{
                if(j==0 || j==m-1)
                    a = 3;
                else
                    a = 4;
            }
            if(n==1 || m==1)
                a--;
            p[a][now]++;

            if(s[j] == 'O'){
                up *= a;
                if((j+i+1)%2 == 1 && evenend == 1)
                    evenend = 0;
                else if((j+i+1)%2 == 1 && evenend == 0)
                    evenend = 1;
            }
        }
    }

    if(even != evenend){
        cout << 0;
        return 0;
    }

    int x = 2;
    if(n==1 || m==1)
        x--;
    // cout << "x=" << x << '\n';
    //
    // cout << p[x][1] << ' ' << p[x+1][1] << ' ' << p[x+2][1] << '\n';
    // cout << p[x][0] << ' ' << p[x+1][0] << ' ' << p[x+2][0] << '\n';
    //
    // cout << "DN(2, 0)" << dn(2,0) << '\n';

    for(int i = even; i <= k; i+=2){
        for(int i1 = 0; i1 <= min(p[x][1], i); i1++){
            for(int i2 = 0; i2 <= min(p[x+1][1], i-i1); i2++){
                for(int i3 = i-i1-i2; i3 <= min(p[x+2][1], i-i1-i2); i3++){
                    for(int i4 = 0; i4 <= min(p[x][0], k-i); i4++){
                        for(int i5 = 0; i5 <= min(p[x+1][0], k-i-i4); i5++){
                            for(int i6 = k-i-i4-i5; i6 <= min(p[x+2][0], k-i-i4-i5); i6++){
                                // cout << i1 << ' ' << i2 << ' ' << i3 << '\n';
                                // cout << i4 << ' ' << i5 << ' ' << i6 << '\n';
                                // cout << down << ' ';
                                down += pt[i1+i4][x] * pt[i2+i5][x+1] * pt[i3+i6][x+2]
                                * DN[p[x][1]][i1] * DN[p[x+1][1]][i2] * DN[p[x+2][1]][i3]
                                * DN[p[x][0]][i4] * DN[p[x+1][0]][i5] * DN[p[x+2][0]][i6];
                                // cout << down << '\n';
                                // cout << "------------------------\n";
                            }
                        }
                    }
                }
            }
        }
    }

    // cout << fixed << setprecision(15) << up << ' ' << down << '\n';
    res = up/down;
    cout << fixed << setprecision(15) << res;

return 0;
}