Codeforces Round #577 D Treasure Hunting (dp+贪心）

大致题意

给一个NxM的平面，有k个宝藏的位置(x,y)，你的初始位置在(1,1), 因为一些特殊原因，你的行动方式只有3种，向上，向左，向右，并且只能在给定的q列位置向上移动。数据范围都是2e5。

思路

每一层的状态之和前一层的到达方式有关，可以考虑dp。如果无脑dp的话可以写成dp[i][j][k]，第i层，从第j列上来，第k列出去，肯定是暴毙的。然后我想到可以贪心考虑每一行最边缘的4个位置，状态可以只用
dp[i][0/1]表示，从左和从右走即可。因为总归要取完整个一行的宝藏。但是我又死在怎么转移上面了。
然后我看了这篇博客：https://www.cnblogs.com/birchtree/p/11301444.html
然后终于明白了，方法是 不要考虑每一行的回头之类的东西，每一行只要考虑从左到右，和从右到左的路程。将折返的路程全都算在转移里面。状态dp[i][0]表示第i行从右走到最左端的答案，dp[i][1]表示从最左端走到最右端的答案。
然后转移就只有简单的4种，即从pre行的左走到now行的左，还有就是左到右，右到左，右到右。对于每一种转移，计算其距离dis的时候，只需要枚举列作为通道的答案，贪心的考虑每个端点的左右两个最近的可行列，这样只要枚举4个列即可，然后加上行的距离即可。

代码

#include<bits/stdc++.h>
using namespace std;
#define maxn 200005
#define maxm 1000006
#define ll long long int
#define INF 0x3f3f3f3f
#define inc(i,l,r) for(int i=l;i<=r;i++)
#define dec(i,r,l) for(int i=r;i>=l;i--)
#define mem(a) memset(a,0,sizeof(a))
#define sqr(x) (x*x)
#define inf (ll)2e18+1
int read(){
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}
     while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
     return f*x;
}
int n,m,k,q,b[maxn],x,y,l[maxn],r[maxn];
ll dp[maxn][2];
vector<int>a[maxn];
ll calc(int y1,int y2,int line){return (ll)abs(y1-line)+abs(y2-line);}
ll dis(int x1,int y1,int x2,int y2)
{
    ll d=(ll)abs(x1-x2);
    return d+(ll)(min(min(calc(y1,y2,l[y1]),calc(y1,y2,r[y1])),min(calc(y1,y2,l[y2]),calc(y1,y2,r[y2]))));
}
int main()
{
    n=read();m=read();k=read();q=read();
    inc(i,1,k){x=read();y=read();a[x].push_back(y);}
    inc(i,1,q)b[i]=read();
    b[0]=-INF;b[q+1]=INF;
    sort(b+1,b+q+1);
    inc(i,1,m){
        l[i]=b[lower_bound(b+1,b+q+1,i)-b-1];
        r[i]=b[lower_bound(b+1,b+q+1,i)-b];
    }
    sort(a[1].begin(),a[1].end());
    if(!a[1].empty()){
        dp[1][0]=dp[1][1]=a[1].back()-1;
        a[1][0]=a[1].back();
    }
    else{
        dp[1][0]=dp[1][1]=0;
        a[1].push_back(1);
    }
    int pre=1;
    inc(i,2,n){
        if(!a[i].empty()){
            sort(a[i].begin(),a[i].end());
            ll l_l=dp[pre][0]+dis(pre,a[pre].front(),i,a[i].front());
            ll l_r=dp[pre][0]+dis(pre,a[pre].front(),i,a[i].back());
            ll r_l=dp[pre][1]+dis(pre,a[pre].back(),i,a[i].front());
            ll r_r=dp[pre][1]+dis(pre,a[pre].back(),i,a[i].back());
            ll tmp=a[i].back()-a[i].front();
            dp[i][0]=tmp+min(l_r,r_r);
            dp[i][1]=tmp+min(l_l,r_l);
            pre=i;
        }
    }
    printf("%lld\n",min(dp[pre][0],dp[pre][1]));
    return 0;
}