[luogu] P3202 [HNOI2009]通往城堡之路（贪心）

P3202 [HNOI2009]通往城堡之路

题目描述

听说公主被关押在城堡里，彭大侠下定决心：不管一路上有多少坎坷，不管城堡中的看守有多少厉害，不管救了公主之后公主会不会再被抓走，不管公主是否漂亮、是否会钟情于自己，他将义无反顾地朝着城堡前进。

可是，通往城堡的路上出现了一些情况。抽象地说，假象地图在二维平面的第一象限。在每个横轴的x位置上有一个高为hx的支撑点，如果彭大侠没有跳到支撑点上，那么他就会掉下去，牺牲在路途。

开始时彭大侠在起点（1，h1）处，而城堡的入口在（n，hn）处。彭大侠每次可以从支撑点（x，hx）跳到支撑点（x+1，hx+1）。但是彭大侠每次的跳跃能量只有d，也就是说，每次跳跃必须满足条件|hx+1-hn|<=d。换句话说，如果两个相邻支撑点的纵向落差大于d，那么彭大侠就无法跳跃了！幸运的是，彭大侠还有一个杀手锏。 在起点处，他可以花一个金币，把某个支撑点升高1个单位，或者降低1个单位。但是，起点处和城堡入口处的支撑点高度不能改变，并且一旦离开起点彭大侠就无法使用该杀手锏。

彭大侠被告知100个金币可兑换一单位生命。于是他希望通过少花金币来保存更多单位的生命。

他终于找到了你这位热心的高手，请你帮他规划一下以便耗费尽量少的金币来到达城堡。

输入输出格式

输入格式：

文件第一行包含一个整数m(m<=5)，表示问题求解次数。接下来的2m行依次表示每次求解的输入数据块。每个输入数据块占2行，其中第一行包含两个整数n和d，分别表示从起点到城堡入口处必须经过的支撑点数和每次跳跃允许的最大纵向落差，n和d之间用空格隔开，输入数据保证2<=n<=5000，0<=d<=109；第二行包含用空格隔开的n个非负整数h1、h2、...、hn，其中hi(1<=i<=n)表示第i个支撑点的高度，特别地，h1表示彭大侠出发时所在支撑点的高度，hn表示城堡入口所在支撑点的高度，输入数据保证对所有1<=i<=n有0<=hi<=109。

输出格式：

有m行，第I(1<=I<=m)行表示第I次求解时彭大侠到达城堡必须耗费的最少金币数量。若无论怎样使用杀手锏他都无法到达城堡，则输出impossible。输入数据保证答案在int64范围之内。

输入输出样例

输入样例#1： 复制

3

10 2

4 5 10 6 6 9 4 7 9 8

3 1

6 4 0

4 2

3 0 6 3

输出样例#1： 复制

6

impossible

4

说明

对样例中的第一个输入数据块，d=2，把第三个支撑点降低3个单位，把第六个支撑点降低1个单位，把第七个支撑点升高2个单位，原序列变成：4 5 7 6 6 8 6 7 9 8，这时任意相邻支撑点的纵向落差没有超过2，彭大侠可以到达城堡！

对样例中的第二个输入数据块，d=1，这时不管怎样调节第二个支撑点的高度，都无法使任意相邻支撑点的纵向落差不超过1。

对样例中的第三个输入数据块，d=2，这时，把第二个支撑点升高1个单位，把第三个支撑点降低3个单位就满足条件了。

【数据规模】

 20%     n<=100

 40%     n<=1000

100% n<=5000

题解

这么神仙的题目为什么不是黑题啊qwq

我感觉[HNOI2012]取石头游戏没这个难啊。

这个贪心不看题解，真的有人想得到这么贪心吗？

首先我们造一组伪解，以最坏的情况考虑，定义

\(b[1]=a[1];\)然后\(b[]\)依次减少\(d.\)

因为我们是这样定义的了，所以我们不再需要降低任意一个区间。

为什么？

因为最后是要让\(b[n]=a[n]\)的，

我们考虑依次提高后缀，为什么是提高后缀，因为我们这样\(O(n)\)考虑一遍的话，可以考虑到是否提高每一个点。

又因为设\(s\)为当前后缀内的\(a[i]>b[i]\)，那么当\(s\)为正时，我们提高这个后缀就一定是优的。而提高这个区间的值为这个区间内最小的(\(a[i]-b[i]\))(这时的\(a[i]>b[i]\))，这样可以使这个点提高到正常高度，因此我们最多提高\(n-1\)次，又每次给区间加上提高的值为\(O(n)\)。所以总复杂度为\(O(n^2)\)。

其中因为我们每次只增加最小值，那么我们当前加入的值一定是让答案更优的，如果下降就只会让答案更劣。

Ps:每次加上的提高值为\(min(a[i]-b[i],b[j-1]+d-b[j]);\)(j为后缀的开头，我们要保证提高的合法性，不能跳不到了。)

Code

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=1e6+5;
ll n,d,ans,a[N],b[N];

ll read(){
    ll x=0,w=1;char ch=getchar();
    while(ch>'9'||ch<'0'){if(ch=='-')w=-1;ch=getchar();}
    while(ch>='0'&&ch<='9')x=x*10+ch-'0',ch=getchar();
    return x*w;
}

void solve(){
    n=read();d=read();ans=0;
    for(ll i=1;i<=n;i++)a[i]=read();
    if(abs(a[1]-a[n])>(n-1)*d)
    {cout<<"impossible"<<endl;return ;}
    b[1]=a[1];
    for(ll i=2;i<=n;i++)b[i]=b[i-1]-d;
    while(b[n]!=a[n]){
        ll mx=-1e9,j,mn=1e9,s=0,addx=1e9;
        for(ll i=n;i>=2;i--){
            if(a[i]>b[i])s++,mn=min(mn,a[i]-b[i]);
            else s--;
            if(s>mx&&b[i]<b[i-1]+d){
                mx=s;j=i;addx=mn;
            }
        }
        addx=min(addx,b[j-1]+d-b[j]);
        for(ll i=j;i<=n;i++)b[i]+=addx;
    }
    for(ll i=1;i<=n;i++)ans+=abs(a[i]-b[i]);
    cout<<ans<<endl;
}

int main(){
    ll T=read();
    while(T--)solve();
    return 0;
}