过河 NOI 蓝桥杯 动态规划题 题解

描述

在河上有一座独木桥，一只青蛙想沿着独木桥从河的一侧跳到另一侧。在桥上有一些石子，青蛙很讨厌踩在这些石子上。由于桥的长度和青蛙一次跳过的距离都是正整数，我们可以把独木桥上青蛙可能到达的点看成数轴上的一串整点：0，1，……，L（其中L是桥的长度）。坐标为0的点表示桥的起点，坐标为L的点表示桥的终点。青蛙从桥的起点开始，不停的向终点方向跳跃。一次跳跃的距离是S到T之间的任意正整数（包括S,T）。当青蛙跳到或跳过坐标为L的点时，就算青蛙已经跳出了独木桥。

题目给出独木桥的长度L，青蛙跳跃的距离范围S,T，桥上石子的位置。你的任务是确定青蛙要想过河，最少需要踩到的石子数。

对于30%的数据，L <= 10000；
对于全部的数据，L <= 10^9。

格式

输入格式

输入的第一行有一个正整数L（1 <= L <= 10^9），表示独木桥的长度。第二行有三个正整数S，T，M，分别表示青蛙一次跳跃的最小距离，最大距离，及桥上石子的个数，其中1 <= S <= T <= 10，1 <= M <= 100。第三行有M个不同的正整数分别表示这M个石子在数轴上的位置（数据保证桥的起点和终点处没有石子）。所有相邻的整数之间用一个空格隔开。

输出格式

输出只包括一个整数，表示青蛙过河最少需要踩到的石子数。

样例1

样例输入1

10
2 3 5
2 3 5 6 7

样例输出1

2

限制

1s

来源

NOIp2005 第二题

题解

根据题意列出状态转移方程，青蛙跳到数轴上的i点时，至少要踩到的石头的数量为：
dp[i]=min(dp[i],dp[i-j]+stone[i])
其中i-j的值为青蛙上一跳的位置；
j的取值范围为S<=j<=T，表示青蛙的跳跃距离范围；
从i=0遍历到最后一个石头所在的位置，遍历时，进行此状态转移方程，最终得到的就是青蛙跳到尽头至少需要踩到的石头数。

但是本题的桥长L可能会很大(L <= 10^9),这么长的数组可能存不完，遍历起来也很慢。
所以我们需要用到状态压缩的知识，把桥的长度压缩为较短但等效的长度。

因为石头的数量相对于桥长而言很少(1 <= M <= 100)，两个石头之间的距离很大，当两个石头之间的距离大到一定程度时，青蛙从前一个石头过去，在不同的跳跃距离组合累加之后，可以跳到下一个石头附近的任意一点（该石头的附近具体指距离该石头小于青蛙最大跳跃距离T的点位的集合）。

所以，当两个石头间的距离大于一定程度X时，我们都可以把它们视为距离X,因为大于X的情况与等于X的情况等效。

那么这个X至少为多少才合适呢？

我们先考虑最坏的情况S==T-1，此时青蛙只有两种跳法，而且两种跳法跳跃的距离都比较远，第一次跳就要跳出去很远，由于两种跳法相距不大，只有1，所以需要X很大的情况下，才能保证能跳到X以后的任意一点。

此时X应该为$(T-1)\ast T$或者说$T\ast S$，因为任何大于等于$T\ast S$的数，你都可以先找到最邻近它且大于的数$n\ast T$(n为正整数），或者说写成n个T的和。此时，如果我们把n个T中的m个T替换为T-1,那么累加后得到的和就是$(n-m)\ast T+m\ast (T-1)==n\ast T-m$，因为m可以取任意正整数，所以$n\ast T-m$可以取到$n\ast T$的附近任何一个距离其小于等于T的数。所以青蛙可以跳到$T\ast S==T\ast (T-1)$(最坏的情况为S==T-1时)之后的任意点位。

所以我们的X要至少大于等于$T\ast S==T\ast (T-1)$(最坏的情况为S==T-1时)。

当然还有一种特殊情况，那就是T==S,此时，任何是T的倍数的点位上的石头都必须被踩到，如果该石头的坐标值能被T整除，那么就一定会被踩到。

#include<iostream>
#include<algorithm>
#include<cstring>
#include<vector>
using namespace std;

#define MAX_SIZE 1000000
#define MAX 0x3f3f3f3f

int L,S,T,M,ans=0;
int stone[MAX_SIZE];
int dp[MAX_SIZE];
vector<int> stone_pos;
int main()
{
	cin>>L>>S>>T>>M;
	if(S==T)
	{
		int current_stone;
		for(int i=0;i<M;i++)
		{
			cin>>current_stone;
			if(current_stone%S==0)ans++;
		}
		cout<<ans;
		return 0;
	}
	else
	{
		ans=MAX;
		int last_stone=0,current_stone,sum_dis=0;
		memset(stone,0,sizeof(stone));
		for(int i=0;i<M;++i)
		{
			int j;
			cin>>j;
			stone_pos.push_back(j);
		}
		sort(stone_pos.begin(),stone_pos.end());
		int last_pos=0;
		for(int i=0;i<stone_pos.size();i++)
		{
			int delta_dis=stone_pos[i]-last_pos;
			sum_dis+=min(90,delta_dis);
			last_pos=stone_pos[i];
			stone[sum_dis]=1;
		}
		memset(dp,MAX,sizeof(dp));
		dp[0]=0;
		for(int i=S;i<=sum_dis+T;i++)
		{
			for(int j=S;j<=T;j++)
			{
				if(i>=j)
				{
					dp[i]=min(dp[i],dp[i-j]+stone[i]);
				}
			}
		}
		for(int i=sum_dis;i<=sum_dis+T;i++)
		{
			ans=min(dp[i],ans);
		}
		cout<<ans;
	}
}