NOIP2017模拟赛 龙珠（dp+单调队列优化）

问题描述

你得到了一个龙珠雷达，它会告诉你龙珠出现的时间和地点。 　　
龙珠雷达的画面是一条水平的数轴，每一个窗口时间，数轴的某些点上会出现同一种龙珠，每当你获得其中一颗龙珠，其它龙珠就会消失。下一个窗口时间，数轴上又会出现另一种龙珠。总共有n个窗口时间，也就是总共有n种龙珠。
　　假设你会瞬间移动，你从数轴的x点移动到y点，耗时0秒，但是需要耗费|x-y|的体力。同时，挖出一颗龙珠也需要耗费一定的体力。请问，最少耗费多少体力，就可以收集齐所有种类的龙珠。

输入格式

第一行，三个整数n,m,x,表示共有n个窗口时间，每个窗口时间会出现m个龙珠，x是一开始你所处的位置。 接下来有两个n*m的矩阵。
对于第一个矩阵，坐标为(i,j)的数字表示第i个窗口时间，第j个龙珠的位置。
对于第二个矩阵，坐标为(i,j)的数字表示第i个窗口时间，挖取第j个龙珠所需的体力。 输出格式 一个整数，表示所需最小体力

样例输入

3 2 5
2 3
4 1
1 3
1 1
1 3
4 2

样例输出

8

题解

f[i][j]=min(f[i-1][k]+abs(wz[i-1][k]-wz[i][j])+cst[i][j]

将此方程式的绝对值展开 可以得到
1.f(i,j)=min（f(i-1,k)-p[i-1][k].wz) || 1<=k<=m ）+p[i][j].c+p[i][j].wz
2.f(i,j)=min（f(i-1,k)+p[i-1][k].wz) || 1<=k<=m ）+p[i][j].c-p[i][j].wz

因为后者的值固定，所以可以考虑单调队列，维护一个wz[i][j]-wz[i-1][k]小于零时时的最小值和wz[i][j]-wz[i-1][k]大于等于零时的最小值

代码

#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdio>
#include<queue>
#include<vector>
using namespace std;
#define inf 1e9
struct node{
    int wz,cst;
}; 
bool cmp(node a,node b){
    return a.wz<b.wz;
}
node a[51][5001];
int n,m,x;
int f[51][5001];
int k;
int main()
{
    int i,j;
    scanf("%d%d%d",&n,&m,&x);
    for(i=1;i<=n;i++){
        for(j=1;j<=m;j++){
            scanf("%d",&a[i][j].wz);
        }
    }
    for(i=1;i<=n;i++){
        for(j=1;j<=m;j++){
            scanf("%d",&a[i][j].cst);
        }
    }
    for(i=1;i<=n;i++) sort(a[i]+1,a[i]+1+m,cmp);
    memset(f,30,sizeof(f));
    for(i=1;i<=m;i++) f[1][i]=abs(a[1][i].wz-x)+a[1][i].cst;

    for(i=2;i<=n;i++){
    int minn=inf;
        for(j=1,k=1;j<=m;j++){

            while(k<=m&&a[i][j].wz>=a[i-1][k].wz){
                minn=min(minn,f[i-1][k]-a[i-1][k].wz);
                k++;
            }
            f[i][j]=minn+a[i][j].cst+a[i][j].wz;
        }
         minn=inf;
        for(j=m,k=m;j>=1;j--){

            while(k>=1&&a[i][j].wz<a[i-1][k].wz){
                minn=min(minn,f[i-1][k]+a[i-1][k].wz);
                k--;
            }
            f[i][j]=min(f[i][j],minn+a[i][j].cst-a[i][j].wz);
        }
    }
    int ans=inf;
    for(i=1;i<=m;i++) ans=min(ans,f[n][i]);
    cout<<ans;
}