HDU 3698 Let the light guide us dp+线段树维护区间极值

很简单的dp 线段树维护区间最小值 第一题线段树 留个纪念

#include <cstdio>
#include <algorithm>
using namespace std;
const int inf=0x3f3f3f3f;
int n,m,T[105][5005],f[105][5005],dp[5005],t[5005];
struct s {
    int l,r,val,lazy;
    void upd(int v) {
        val=min(val,v);
        lazy=min(lazy,v);
    }
} seg[5005<<2];
void pup(int x) {
    seg[x].val=min(seg[x<<1].val,seg[x<<1|1].val);
}
void init(int x,int l,int r) {
    seg[x].l=l,seg[x].r=r;
    seg[x].val=seg[x].lazy=inf;
    if (l==r)
        return;
    int m=(l+r)>>1;
    init(x<<1,l,m);
    init(x<<1|1,m+1,r);
    pup(x);
}
void pdown(int x) {
    int lv=seg[x].lazy;
    if (lv!=inf) {
        seg[x].lazy=inf;
        seg[x<<1].upd(lv);
        seg[x<<1|1].upd(lv);
    }
}
void upd(int x,int l,int r,int v) {
    int L=seg[x].l,R=seg[x].r;
    if (l<=L&&r>=R) {
        seg[x].upd(v);
        return;
    }
    pdown(x);
    int m=(L+R)>>1;
    if (m>=l)
        upd(x<<1,l,r,v);
    if (m<r)
        upd(x<<1|1,l,r,v);
    pup(x);
}
int qry(int x,int l,int r) {
    int L=seg[x].l,R=seg[x].r;
    if (l<=L&&r>=R)
        return seg[x].val;
    pdown(x);
    int m=(L+R)>>1,ret=inf;
    if (m>=l)
        ret=min(ret,qry(x<<1,l,r));
    if (m<r)
        ret=min(ret,qry(x<<1|1,l,r));
    pup(x);
    return ret;
}
int main()
{
    while (scanf("%d%d",&n,&m)==2&&n+m) {
        for (int i=1;i<=n;++i)
            for (int j=1;j<=m;++j)
                scanf("%d",&T[i][j]);
        for (int i=1;i<=n;++i)
            for (int j=1;j<=m;++j)
                scanf("%d",&f[i][j]);
        init(1,1,m);
        for (int i=1;i<=m;++i)
            upd(1,max(1,i-f[1][i]),min(m,i+f[1][i]),T[1][i]);
        for (int i=2;i<=n;++i) {
            for (int j=1;j<=m;++j)
                t[j]=qry(1,max(1,j-f[i][j]),min(m,j+f[i][j]))+T[i][j];
            init(1,1,m);
            for (int j=1;j<=m;++j)
                upd(1,max(1,j-f[i][j]),min(m,j+f[i][j]),t[j]);
        }
        printf("%d\n",qry(1,1,m));
    }
    return 0;
}