hdu 3698 Let the light guide us【线段树+dp】 2010 Asia Fuzhou Regional Contest

原题连接：http://acm.hdu.edu.cn/showproblem.php?pid=3698

如果不考虑数据规模的话，这个题用dp是可以解决的，应该不难理解，但是这个题的数据很大

然后的工作就是怎样优化的问题了

线段树优化，我第一次遇到这种题，但是理解起来不是太难

我在做的时候，线段树的延迟标记没有弄好，超时了好几次

只说一点：题目中的要求 |j-k|≤f(i,j)+f(i+1,k)可以理解成 j点：以j为中心的半径f(i,j)的区间和以k为中心的半径f(i+1,k)的区间，如果有交叉，那么就符合条件，这就是code中的查找和查询的依据。

code:

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define fmid (l + r) >> 1;
#define N 5002
#define inf 999999999
int sum[4*N],cnt[4*N];

void PushDown(int rt)//先将这一段负值为一个更小的值
{
    if(cnt[rt] == -1) return;
    if(cnt[rt<<1] == -1) cnt[rt<<1] = cnt[rt];
    else cnt[rt<<1] = min(cnt[rt<<1],cnt[rt]);

    if(cnt[rt<<1|1] == -1) cnt[rt<<1|1] = cnt[rt];
    else cnt[rt<<1|1] = min(cnt[rt<<1|1],cnt[rt]);

    sum[rt<<1] = min(sum[rt<<1],cnt[rt<<1]);
    sum[rt<<1|1] = min(sum[rt<<1|1],cnt[rt<<1|1]);

    cnt[rt] = -1;
}
void build(int l,int r,int rt)
{
    sum[rt] = inf;
    cnt[rt] = -1;
    if(l == r)   return;
    int mid = fmid;
    build(lson);
    build(rson);
}

void update(int L,int R,int data,int l,int r,int rt)
{
    if(L <= l && R >= r)
    {
        if(cnt[rt] != -1) cnt[rt] = min(cnt[rt],data);
        else cnt[rt] = data;
        sum[rt] = min(cnt[rt],sum[rt]);
        return;
    }
    PushDown(rt);
    int mid = fmid;
    if(L <= mid) update(L,R,data,lson);
    if(R > mid) update(L,R,data,rson);

    sum[rt] = min(sum[rt<<1],sum[rt<<1|1]);
}
int query(int L,int R,int l,int r,int rt)
{
    if(L <= l && R >= r)
    {
        return sum[rt];
    }
    PushDown(rt);
    int mid = fmid;
    int re = inf;
    if(L <= mid) re = min(re,query(L,R,lson));
    if(R > mid) re = min(re,query(L,R,rson));
    return re;
}


int n,m;
int t[110][N],f[110][N];

int dp[110][N];
inline int getInt()
{
    char c;
    while(1)
    {
        c = getchar();
        if(c >= '0' && c <= '9') break;
    }
    int re = 0;
    while(c >= '0' && c <= '9')
    {
        re = re*10 + c-'0';
        c = getchar();
    }
    return re;
}
int main()
{
    while(scanf("%d%d",&n,&m) && (n != 0 || m != 0))
    {
        for(int i = 0;i < n;i++)
            for(int j = 0;j < m;j++)
                //t[i][j] = getInt();
                scanf("%d",&t[i][j]);
        for(int i = 0;i < n;i++)
            for(int j= 0;j < m;j++)
                //f[i][j] = getInt();
                scanf("%d",&f[i][j]);


        for(int i = 0;i < m;i++)    dp[0][i] = t[0][i];
        for(int i = 0;i < n-1;i++)
        {
            build(0,m-1,1);
            for(int j = 0;j < m;j++)
            {
                update(j-f[i][j],j+f[i][j],dp[i][j],0,m-1,1);
            }
            for(int j = 0;j < m;j++)
            {
                dp[i+1][j] = query(j-f[i+1][j], j+f[i+1][j],0,m-1,1) + t[i+1][j];
            }
        }
        int re = inf;
        for(int i = 0;i < m;i++)    re = min(re,dp[n-1][i]);
        printf("%d\n",re);
    }
    return 0;
}