#include using namespace std; #define rep(i,n) for(int i=0;i<(n);++i) #define loop for(;;) #define trace(var) cerr<<">>> "<<#var<<" = "<; template inline ostream& operator<<(ostream&os, pair p) { return os << '(' << p.first << ", " << p.second << ')'; } template inline ostream& operator<<(ostream&os, tuple t) { return os << '(' << get<0>(t) << ", " << get<1>(t) << ", " << get<2>(t) << ')'; } template inline ostream& operator<<(ostream&os, set v) { os << "(set"; for (T item: v) os << ' ' << item; os << ")"; return os; } template inline ostream& operator<<(ostream&os, vector v) { if (v.size() == 0) { return os << "(empty)"; } os << v[0]; for (int i=1, len=v.size(); i inline istream& operator>>(istream&is, vector&v) { rep (i, v.size()) is >> v[i]; return is; } using vi = vector; using vvi = vector; using vd = vector; using vvd = vector; using vb = vector; inline vi bit_pattern(int I, int w) { vi r; rep (i, w) r.push_back((I & (1 << i)) ? 1 : 0); return r; } int main() { cerr << fixed << setprecision(5); int h, w; cin >> h >> w; vvi p(h+1, vi(w, 0)); vvi s(h+1, vi(w, 0)); rep (i, h) rep (j, w) cin >> p[i+1][j]; rep (i, h) rep (j, w) cin >> s[i+1][j]; rep (j, w) { p[0][j] = 0; s[0][j] = 4; } vvd Pr(h+1, vd(w, 0)); for (int i = 1; i <= h; ++i) { vvd pt(w, vd(5, 0)); // pt[j][s] = Pr[score(i,j) == s] ただし pt[j][4] = Pr[score>=4] rep (j, w) { pt[j][4 - s[i][j]] = p[i][j]/100.0; pt[j][0] = 1.0 - p[i][j]/100.0; if (s[i][j] == 4) pt[j][0] = 1.0; } // 前の行からの伝播 rep (j, w) { vd memo(6, 0); rep (k, 5) { memo[k+1] += pt[j][k] * Pr[i-1][j]; memo[k] += pt[j][k] * (1.0 - Pr[i-1][j]); } rep (k, 5) pt[j][k] = memo[k]; pt[j][4] += memo[5]; } vector used(w, false); rep (_, w) { int idx = -1; double mp = -1; rep (j, w) { if (not used[j] and pt[j][4] > mp) { mp = pt[j][4]; idx = j; } } used[idx] = true; if (idx > 0) { // to left int j = idx; vd memo(6, 0); rep (k, 5) { memo[k+1] += pt[j-1][k] * mp; memo[k] += pt[j-1][k] * (1.0 - mp); } rep (k, 5) pt[j-1][k] = memo[k]; pt[j-1][4] += memo[5]; } if (idx < w-1) { // to right int j = idx; vd memo(6, 0); rep (k, 5) { memo[k+1] += pt[j+1][k] * mp; memo[k] += pt[j+1][k] * (1.0 - mp); } rep (k, 5) pt[j+1][k] = memo[k]; pt[j+1][4] += memo[5]; } } rep (j, w) Pr[i][j] = pt[j][4]; } //rep (i, h+1) { rep (j, w) cerr << Pr[i][j] << ' '; cerr << endl; } double ans = 0.0; rep (i, h+1) rep (j, w) ans += Pr[i][j]; cout << ans << endl; return 0; }