#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] = 0.0; top[j] = now[j][4]; now[j][4] = 0.0; } // interacting 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; }