2112: 聪明的美食家
Time Limit: 1 Sec Memory Limit: 128 MB
Submit: 279 Solved: 52
[Submit][Status][Web Board]
Description
如果有人认为吃东西只需要嘴巴,那就错了。
都知道舌头有这么一个特性,“由简入奢易,由奢如简难”(据好事者考究,此规律也适合许多其他情况)。具体而言,如果是甜食,当你吃的食物不如前面刚吃过的东西甜,就很不爽了。
大宝是一个聪明的美食家,当然深谙此道。一次他来到某小吃一条街,准备从街的一头吃到另一头。为了吃得爽,他大费周章,得到了各种食物的“美味度”。他拒绝不爽的经历,不走回头路而且还要爽歪歪(爽的次数尽量多)。
Input
两行数据。
第一行为一个整数n,表示小吃街上小吃的数量
第二行为n个整数,分别表示n种食物的“美味度”
Output
一个整数,表示吃得爽的次数
Sample Input
10 3 18 7 14 10 12 23 41 16 24
Sample Output
6
HINT
美味度为0到10000的整数
n<200000
题意:
这个实际上就是大宝从一端走到另一端,然后所经过的美味度必须是单调递增的,其实相邻的美味度也是可以相等的。
思路:
这个实际上就是动态规划中的最长上升子序列的问题,解决这种问题一般有两种方法模板,
一种就是dp的朴素算法,还有一种就是在原来dp基础上加上用二分查找可以降低运算复杂度。
下面来讲一下这两个模板:
先给出一个例子:
【题目描述】
给定N个数,求这N个数的最长上升子序列的长度。
【样例输入】
7
2 5 3 4 1 7 6
【样例输出】
4
什么是最长上升子序列? 就是给你一个序列,请你在其中求出一段不断严格上升的部分,它不一定要连续。
就像这样:2,3,4,7和2,3,4,6就是序列2 5 3 4 1 7 6的两种选取方案。最长的长度是4.
那么,怎么求出它的最大上升子序列长度为4呢?这里介绍两种方法,都是以动态规划为基础的。
首先,我们先介绍较慢(O())的方法。我们记num为到这个数为止,最长上升子序列的长度。
这种方法就是每一次寻找“可以接下去的”,换句话说,设原序列为a,则
当且时,。
对于每一个数,他都是在“可以接下去”的中,从前面的最优值+1转移而来。
因此,这个算法是可以求出正确答案的。复杂度很明显,外层i枚举每个数,内层j枚举目前i的最优值,即O()。
那么,有没有更快的方法呢?当然有。这回要用到二分。
我们回想一下,在上面O()的程序中,哪些地方看起来比较费时?
没错,就是内层用于更新i的循环。因为每一次他都要查找一遍,效率并不高。
回到题目,我们发现,他只要我们求长度,所以?
我们可以模拟一个栈。
所以每遇到一个比栈顶元素大的数,就放进栈里,遇到比栈顶元素小的就二分查找前边的元素,找到一个“最应该被换掉的元素”,用新数去更新前边的元素。
这个算法不难证明也是正确的。因为前面每一次的枚举都换成了二分,内层的复杂度从降到了,外层不变。所以总的复杂度是O()。
接下来,我先给出朴素算法的代码:
#include<cstdio>
const int MAX=1001;
int a[MAX];
int lis(int x)
{
int num[MAX];
for(int i=0;i<x;i++)
{
num[i]=1;
for(int j=0;j<i;j++)
{
if(a[j]<a[i]&&num[j]+1>num[i])
num[i]=num[j]+1;
}
}
int maxx=0;
for(int i=0;i<x;i++)
if(maxx<num[i])
maxx=num[i];
return maxx;
}
int main()
{
int n;
scanf("%d",&n);
for(int i=0;i<n;i++)
scanf("%d",&a[i]);
return !printf("%d\n",lis(n));
}
这个则是二分算法的代码:
#include<cstdio>
#include<algorithm>
const int MAXN=200001;
int a[MAXN];
int d[MAXN];
int main()
{
int n;
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
d[1]=a[1];
int len=1;
for(int i=2;i<=n;i++)
{
if(a[i]>d[len])
d[++len]=a[i];
else
{
int j=std::lower_bound(d+1,d+len+1,a[i])-d;
d[j]=a[i];
}
}
printf("%d\n",len);
return 0;
}
下面给出我们这一道题目的AC代码:
#include<cstdio>
#include<algorithm>
int a[200005];//a用来存放所有的美味度
int d[200005];//b用来存放最长的单调递增的序列
int main()
{
int n;
while(scanf("%d",&n)!=EOF)//输入所有美味度
{
for(int i=1;i<=n;i++)
{
scanf("%d",&a[i]);
}
d[1]=a[1];//先把a【1】放到d【1】中
int len=1;
for(int i=2;i<=n;i++)//从2开始遍历
{
if(a[i]>=d[len])
d[++len]=a[i];
else
{
int j=std::upper_bound(d+1,d+len+1,a[i])-d;//用了二分查找
d[j]=a[i];
}
}
printf("%d\n",len);
}
return 0;
}