[NOIP2014 提高组] 飞扬的小鸟（线性dp）：类似完全背包问题的通式优化

[NOIP2014 提高组] 飞扬的小鸟

题目背景

NOIP2014 提高组 D1T3

题目描述

Flappy Bird 是一款风靡一时的休闲手机游戏。玩家需要不断控制点击手机屏幕的频率来调节小鸟的飞行高度，让小鸟顺利通过画面右方的管道缝隙。如果小鸟一不小心撞到了水管或者掉在地上的话，便宣告失败。

为了简化问题，我们对游戏规则进行了简化和改编:

游戏界面是一个长为 $n$，高为 $m$ 的二维平面，其中有 $k$ 个管道（忽略管道的宽度）。

小鸟始终在游戏界面内移动。小鸟从游戏界面最左边任意整数高度位置出发，到达游戏界面最右边时，游戏完成。

小鸟每个单位时间沿横坐标方向右移的距离为 $1$，竖直移动的距离由玩家控制。如果点击屏幕，小鸟就会上升一定高度 $x$，每个单位时间可以点击多次，效果叠加；如果不点击屏幕，小鸟就会下降一定高度 $y$。小鸟位于横坐标方向不同位置时，上升的高度 $x$ 和下降的高度 $y$ 可能互不相同。

小鸟高度等于 $0$ 或者小鸟碰到管道时，游戏失败。小鸟高度为 $m$ 时，无法再上升。

现在,请你判断是否可以完成游戏。如果可以，输出最少点击屏幕数；否则，输出小鸟最多可以通过多少个管道缝隙。

输入格式

第 $1$ 行有 $3$ 个整数 $n,m,k$，分别表示游戏界面的长度，高度和水管的数量，每两个整数之间用一个空格隔开；

接下来的 $n$ 行，每行 $2$ 个用一个空格隔开的整数 $x$ 和 $y$，依次表示在横坐标位置 $0\sim n-1$ 上玩家点击屏幕后，小鸟在下一位置上升的高度 $x$，以及在这个位置上玩家不点击屏幕时，小鸟在下一位置下降的高度 $y$。

接下来 $k$ 行，每行 $3$ 个整数 $p,l,h$，每两个整数之间用一个空格隔开。每行表示一个管道，其中 $p$ 表示管道的横坐标，$l$ 表示此管道缝隙的下边沿高度，$h$ 表示管道缝隙上边沿的高度（输入数据保证 $p$ 各不相同，但不保证按照大小顺序给出）。

输出格式

共两行。

第一行，包含一个整数，如果可以成功完成游戏，则输出 $1$，否则输出 $0$。

第二行，包含一个整数，如果第一行为 $1$，则输出成功完成游戏需要最少点击屏幕数，否则，输出小鸟最多可以通过多少个管道缝隙。

样例 #1

样例输入 #1

10 10 6 
3 9  
9 9  
1 2  
1 3  
1 2  
1 1  
2 1  
2 1  
1 6  
2 2  
1 2 7 
5 1 5 
6 3 5 
7 5 8 
8 7 9 
9 1 3

样例输出 #1

1
6

样例 #2

样例输入 #2

10 10 4 
1 2  
3 1  
2 2  
1 8  
1 8  
3 2  
2 1  
2 1  
2 2  
1 2  
1 0 2 
6 7 9 
9 1 4 
3 8 10

样例输出 #2

0
3

提示

【输入输出样例说明】

如下图所示，蓝色直线表示小鸟的飞行轨迹，红色直线表示管道。

【数据范围】

对于 $30\%$ 的数据：$5\le n\le 10,5\le m\le 10,k=0$，保证存在一组最优解使得同一单位时间最多点击屏幕 $3$ 次；

对于 $50\%$ 的数据：$5\le n\le 20,5\le m\le 10$，保证存在一组最优解使得同一单位时间最多点击屏幕 $3$ 次；

对于 $70\%$ 的数据：$5\le n\le 1000,5\le m\le 100$；

对于 $100\%$ 的数据：$5\le n\le 10000$，$5\le m\le 1000$，$0\le k<n$，$0<x,y<m$，$0<p<n$，$0\le l<h\le m$， $l+1<h$。

思路

• 首先，这道题是否上升下降然后怎样上升下降我们可以通过子操作来判定，因此本道题拥有最优子结构且求最优解，因此本道题可以尝试 dp。
• 然后我们的状态表示就是 $f_{i,j}$ 表示所有鸟飞过第 i 个当前管道且当前高度为 j 的所有最短操作次数集合。
• 我们肯定很容易划分集合，一个集合是从下面（这种“/”）转移过来，一个是从上面（这种“\”）转移过来。
• 而且这题还有管道的限制，也就是非法状态。那只要首先先忽略非法状态直接转移，在下一次迭代更新前再把非法状态清空就可以了。
• 状态转移： f i , j = m i n { f i − 1 , j − x i + 1 , f i − 1 , j − 2 x i + 2 , f i − 1 , j − 3 x i + 3..... } f_{i,j}=min\{f_{i-1,j-x_i}+1,f_{i-1,j-2x_i}+2,f_{i-1,j-3x_i}+3.....\} fi,j​=min{fi−1,j−xi​​+1,fi−1,j−2xi​​+2,fi−1,j−3xi​​+3.....}。这是对于从下往上飞，如果从上往下飞，直接就 $f_{i,j}=min(f_{i,j},f_{i-1,j+x_i})$，为什么这样呢？我们规定有升才会有降。
• 之后，我们对第一个状态转移进行化简（类似完全背包问题的那种化简）=>$f_{i,j}=min(f_{i-1,j-x_i}+1,f_{i,j-x_i}+1)$。同时注意一个细节，也就是鸟是飞不出屏幕的，也就是说对于 $j=M$ ，这个情况要特判。

细节

• 管道：管道不是每个地方都有的，他有固定的位置存放。因此先把数组初始化为全部合法，然后再读入。同时题目也要求输出鸟如果失败通过的最多管道数，因此可以用一个前缀和数组统计一下每个点前有多少个管道经过。

• 

• $f_{i,j}$ 我们在状态转移的时候，第一维度为 $i$ ，第二维度为 $j$。

代码

//二刷

//这个状态转移方程跟背包问题极其相似

//f[i][j]=min(f[i-1][j-xi]+1,f[i-1][j-2*xi]+2....)

//优化后变成了
//f[i][j]=min(f[i][j-xi]+1,f[i-1][j-xi]+1)

#include<iostream>
#include<algorithm>
#include<cstring>

using namespace std;

const int N = 1e4+10,M = 1e3+10,INT = 0x3f3f3f3f;

int n,m,k;
int x[N],y[N];
bool st[N];
int l[N],h[N];
int f[N][M];
int cnt[N];
struct E{
    int x,y;
}pipe[N];

int main(){
    cin>>n>>m>>k;
    
    for(int i=0;i<n;i++){
        cin>>x[i]>>y[i];
    }
    
    for(int i=1;i<=n;i++) pipe[i]={0,m+1};
    
    for(int i=1;i<=k;i++){
        int q,x,y;
        cin>>q>>x>>y;
        pipe[q]={x,y};
        cnt[q]++;
    }
    
    for(int i=1;i<=n;i++){
        cnt[i]+=cnt[i-1];
    }
    
    memset(f,0x3f,sizeof f);
    
    for(int i=1;i<=m;i++)f[0][i]=0;
    
    for(int i=1;i<=n;i++){
        int pre_up=x[i-1],pre_down=y[i-1];
        for(int j=1;j<=m;j++){/枚举高度
            if(j>pre_up){
                f[i][j]=min(f[i][j-pre_up]+1,f[i][j]);
                f[i][j]=min(f[i-1][j-pre_up]+1,f[i][j]);
            }
            
            if(j==m){
                //特判
                for(int k=m-pre_up;k<=m;k++){
                    f[i][m]=min(f[i][m],f[i][k]+1);
                    f[i][m]=min(f[i][m],f[i-1][k]+1);
                }
            }
        }
        
        
        for(int j=1;j<=m;j++){
            if(j+pre_down<=m){
                f[i][j]=min(f[i][j],f[i-1][j+pre_down]);
            }
        }
        
        for(int j=0;j<=pipe[i].x;j++)f[i][j]=INT;
        for(int j=pipe[i].y;j<=m;j++)f[i][j]=INT;
    }
    
    
    int res=INT;
    
    for(int i=1;i<=m;i++)res=min(res,f[n][i]);
    
    if(res<INT/2)puts("1"),cout<<res;
    else{
        puts("0");
        for(int i=n;i>=1;i--){
            for(int j=1;j<=m;j++){
                if(f[i][j]<INT/2){
                    cout<<cnt[i];
                    return 0;
                }
            }
        }
    }
    return 0;
}