Let the light guide us HDU

//inline 快读 交C++ 加速
//懒标记的不同记录
//DP+线段树 题目中的要求 |j-k|≤f(i,j)+f(i+1,k)可以理解成
// j点：以j为中心的半径f(i,j)的区间和以k为中心的半径f(i+1,k)的区间
//一个是插入区间 ，一个是查询最小的区间
//代码参考http://blog.csdn.net/wukonwukon/article/details/7020590
#include<stdio.h>
#include<string.h>
#define LL(x) (x<<1)
#define RR(x) ((x<<1)|1)

int n,m;
const int M1=105;
const int M2=5005;
const int INF=10000005;
int T[M1][M2];
int F[M1][M2];
int dp[M1][M2];
int flag;
struct Seg{
    int l,r;
    int min;
    int flag;
}tree[M2*4];

inline void build(int l,int r,int k){
    tree[k].flag=-1;
    tree[k].min=INF;
    tree[k].l=l;
    tree[k].r=r;
    if(l==r)return ;
    int mid=(l+r)>>1;
    build(l,mid,LL(k));
    build(mid+1,r,RR(k));
}

inline int Min(int a,int b){
    return (a<b?a:b);
}

inline void push_down(int k){
    if(tree[k].flag!=-1){
        int l=LL(k),r=RR(k);
		if(tree[l].flag==-1)tree[l].flag=tree[k].flag;
		else tree[l].flag=Min(tree[l].flag,tree[k].flag);

		if(tree[r].flag==-1)tree[r].flag=tree[k].flag;
		else tree[r].flag=Min(tree[r].flag,tree[k].flag);
			
		if(tree[l].min>tree[k].flag){
			tree[l].min=tree[k].flag;
		}
		if(tree[r].min>tree[k].flag){
			tree[r].min=tree[k].flag;
		} 
        tree[k].flag=-1;
    }
}

inline void update(int l,int r,int val,int k){
    if(l<=tree[k].l && r>=tree[k].r){
        if(tree[k].min>val){    
            tree[k].min=val;
        }
		if(tree[k].flag==-1)tree[k].flag=val;
		else tree[k].flag=Min(tree[k].flag,val);
        return ;
    }
    push_down(k);
	int L=LL(k);
	int R=RR(k);
    int mid=(tree[k].l+tree[k].r)>>1;
    if(l<=mid) update(l,r,val,L);
    if(r>mid) update(l,r,val,R);
    tree[k].min=Min(tree[L].min,tree[R].min);

}

inline int query(int l,int r,int k){
    if(l<=tree[k].l && tree[k].r<=r){
        return tree[k].min;
    }
    push_down(k);
    int mid=(tree[k].l+tree[k].r)>>1;
	int ret=INF;
    if(l<=mid) ret=Min(ret, query(l,r,LL(k)));
    if(r>mid) ret=Min(ret, query(l,r,RR(k)));
	return ret;
}
int get_val() {
    int res;
    char c;
    while ((c=getchar()) < '0' || c > '9');
    res = c - '0';
    while ((c=getchar()) >= '0' && c <= '9')
        res = res * 10 + c - '0';
    return res;
}

int main(){
    int i,j,k;
	int l,r;
    while(scanf("%d%d",&n,&m)!=EOF && n+m){
        for(i=0;i<n;i++)
            for(j=0;j<m;j++)
               // T[i][j]=get_val() ;
				scanf("%d",&T[i][j]);
        for(i=0;i<n;i++)
            for(j=0;j<m;j++)
				scanf("%d",&F[i][j]);
                //F[i][j]=get_val() ;
        for(i=0;i<m;i++)
			dp[0][i]=0;
        for(i=0;i<n-1;i++){
            build(0,m,1);
            for(j=0;j<m;j++){
                 l=j-F[i][j];
                 r=j+F[i][j];
                if(l<0)l=0;
                if(r>=m)r=m-1;
                update(l,r,dp[i][j]+T[i][j],1);
            }
            for(j=0;j<m;j++){
                 l=j-F[i+1][j];
                 r=j+F[i+1][j];
                if(l<0)l=0;
                if(r>=m)r=m-1;
                dp[i+1][j]=query(l,r,1);
            }
        }
        int ans=INF;
        for(i=0;i<m;i++){
            ans=Min(ans,dp[n-1][i]+T[n-1][i]);
		}
        printf("%d\n",ans);
    }
    return 0;
}