#include "wombats.h"
#include<bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#include <ext/rope>
//#define int long long
#define pb push_back
#define s second
#define f first
#define pf push_front
#define inf 100000000000000000
#define bitebi __builtin_popcountll
#define FOR( i , n ) for( int i = 0 ; i < n ; i ++ )
#define YES cout <<"YES\n"
#define NO cout << "NO\n"
#define debug cout << "Here Fine" << endl ;
#define pr pair < int , int >
#define fbo find_by_order // returns iterator
#define ook order_of_key // returns strictly less numbers than key
using namespace std ;
#pragma GCC optimize("Ofast")
#pragma GCC target("avx,avx2,fma")
using namespace __gnu_pbds;
using namespace __gnu_cxx;
#define ordered_set tree<int, null_type,less<int>, rb_tree_tag,tree_order_statistics_node_update>
const double Pi=acos(-1.0);
const double EPS=1E-8;
const int mod = 1000000007 ;
const int mod1 = 998244353 ;
const int N = 2e5 + 10 ;
mt19937 R(time(0));
map < int , int > ma , ma1 ;
int r , c , h[ 5000 ][ 200 ] , v[ 5000 ][ 200 ] ;
vector < int > dp[ 201 ] ;
vector < int > lst , nw ;
void calculate_dp( int pos ){
nw.clear() ; lst.clear() ;
FOR( i , c ) lst.pb( 0 ) , nw.pb( INT_MAX ) ;
int val = 0 ;
for( int j = pos - 1 ; j >= 0 ; j -- ){
val = val + h[ r - 1 ][ j ] ;
lst[ j ] = val ;
}
val = 0 ;
for( int j = pos ; j < c - 1 ; j ++ ){
val = val + h[ r - 1 ][ j ] ;
lst[ j + 1 ] = val ;
}
dp[ pos ].clear() ;
for( int i = r - 2 ; i >= 0 ; i -- ){
FOR( i , c ) nw[ i ] = INT_MAX ;
int val = INT_MAX ;
for( int j = 0 ; j < c ; j ++ ){
val = min( val , lst[ j ] + v[ i ][ j ] ) ;
nw[ j ] = val ;
if( j != c - 1 )
val = val + h[ i ][ j ] ;
}
val = INT_MAX ;
for( int j = c - 1 ; j >= 0 ; j -- ){
val = min( val , lst[ j ] + v[ i ][ j ] ) ;
nw[ j ] = min( nw[ j ] , val ) ;
if( j != 0 )
val = val + h[ i ][ j - 1 ] ;
}
if( i == 0 ){
for( auto x : nw ) dp[ pos ].pb( x ) ;
}
swap( lst , nw ) ;
}
}
void init(int R, int C, int H[5000][200], int V[5000][200]){
FOR( i , R ){
FOR( j , C - 1 ) h[ i ][ j ] = H[ i ][ j ] ;
}
FOR( i , C ){
FOR( j , R - 1 ) v[ j ][ i ] = V[ j ][ i ] ;
}
r = R , c = C ;
FOR( i , C )
calculate_dp( i ) ;
}
void changeH(int P, int Q, int W){
h[ P ][ Q ] = W ;
FOR( i , c ) calculate_dp( i ) ;
}
void changeV(int P, int Q, int W){
v[ P ][ Q ] = W ;
FOR( i , c ) calculate_dp( i ) ;
}
int escape(int V1, int V2){
return dp[ V2 ][ V1 ] ;
}
//static int H[5000][200];
//static int V[5000][200];
//int main() {
// int R, C, E, P, Q, W, V1, V2, event, i, j;
//
// assert(scanf("%d%d", &R, &C) == 2);
// for (i = 0; i < R; ++i)
// for (j = 0; j < C-1; ++j)
// assert(scanf("%d", &H[i][j]) == 1);
// for (i = 0; i < R-1; ++i)
// for (j = 0; j < C; ++j)
// assert(scanf("%d", &V[i][j]) == 1);
//
// init(R, C, H, V);
//
// std::vector<int> answers;
//
// assert(scanf("%d", &E) == 1);
// for (i = 0; i < E; i++) {
// assert(scanf("%d", &event) == 1);
// if (event == 1) {
// assert(scanf("%d%d%d", &P, &Q, &W) == 3);
// changeH(P, Q, W);
// } else if (event == 2) {
// assert(scanf("%d%d%d", &P, &Q, &W) == 3);
// changeV(P, Q, W);
// } else if (event == 3) {
// assert(scanf("%d%d", &V1, &V2) == 2);
// answers.push_back(escape(V1, V2));
// }
// }
// for (int ans : answers) {
// printf("%d\n", ans);
// }
// return 0;
//}
