#include using namespace std; typedef long long ll; typedef vector vi; typedef vector vl; typedef complex P; typedef pair pii; #define REP(i,n) for(ll i=0;i &vr,int from,int to,Real rate){ Real weight = vr[from]*rate; vr[to] += weight; vr[from] -= weight; if(weight > EPS && to == 4)return true; return false; } int main(){ int h,w; cin>>h>>w; int p[h][w],s[h][w]; REP(i,h)REP(j,w)cin>>p[i][j]; REP(i,h)REP(j,w)cin>>s[i][j]; Real result = 0.0; Real before[w]; vector< vector > now(w,vector(5)), tmp(w,vector(5)); vector top(w); REP(i,w)before[i]=0.0; REP(i,h){ // initialize REP(j,w){ // init REPR(k,5)now[j][k]=0.0; now[j][0]=1.0; // from S Real upp = (Real)p[i][j]/100.0; int uph = 4-s[i][j]; cut(now[j],0,uph,upp); // from before upp = before[j]; REPR(_k,5){ int k = 4-_k; cut(now[j],k,k+1,upp); } // for interacting phase before[j] = now[j][4]; top[j] = now[j][4]; now[j][4] = 0.0; } // interacting Real over = 0.0; // toright REP(j,w){ REPR(_k,5){ int k = 4-_k; cut(now[j],k,k+1,over); } over = top[j]; before[j] += now[j][4]; now[j][4] = 0.0; } over = 0.0; // toleft REP(_j,w){ int j = w-1-_j; REPR(_k,5){ int k = 4-_k; cut(now[j],k,k+1,over); } over = top[j]; before[j] += now[j][4]; now[j][4] = 0.0; } // accumulate REP(j,w){ result += before[j]; } } // bool flag = true; // while(flag){ // flag = false; // // toright // REPR(j,w){ // REPR(_k,5){ // int k = 4-_k; // cut(now[j],k,k+1,top[j-1]); // } // } // // toleft // REPR(_j,w){ // int j = w-1-_j; // REPR(_k,5){ // int k = 4-_k; // cut(now[j],k,k+1,top[j+1]); // } // } // // copying // REP(j,w){ // before[j] += top[j]; // top[j] = now[j][4]; // now[j][4] = 0.0; // if(top[j]>EPS)flag = true; // } // } // REP(j,w){ // result += before[j]; // } // } cout.precision(15); cout << fixed << result << endl; return 0; }