LIS定义
LIS(Longest Increasing Subsequence)最长上升子序列
一个数的序列bi,当b1 < b2 < … < bS的时候,我们称这个序列是上升的。对于给定的一个序列(a1, a2, …, aN),我们可以得到一些上升的子序列(ai1, ai2, …, aiK),这里1 <= i1 < i2 < … < iK <= N。
比如,对于序列(1, 7, 3, 5, 9, 4, 8),有它的一些上升子序列,如(1, 7), (3, 4, 8)等等。这些子序列中最长的长度是4,比如子序列(1, 3, 5, 8).
你的任务,就是对于给定的序列,求出最长上升子序列的长度。
两种做法
O(N^2)做法:dp动态规划
状态设计:dp[i]代表以a[i]结尾的LIS的长度
状态转移:dp[i]=max(dp[i], dp[j]+1) (0<=j< i, a[j]< a[i])
边界处理:dp[i]=1 (0<=j< n)
时间复杂度:O(N^2)
举例: 对于序列(1, 7, 3, 5, 9, 4, 8),dp的变化过程如下
dp[i] 初始值 j=0 j=1 j=2 j=3 j=4 j=5
dp[0] 1
dp[1] 1 2
dp[2] 1 2 2
dp[3] 1 2 2 3
dp[4] 1 2 3 3 4
dp[5] 1 2 2 3 3 3
dp[6] 1 2 3 3 4 4 4
求完dp数组后,取其中的最大值就是LIS的长度。【注意答案不是dp[n-1],这个样例只是巧合】
#include <bits/stdc++.h>
using namespace std;
const int MAXX=100000+5;
const int INF=INT_MAX;
int a[MAXX],dp[MAXX]; // a数组为数据,dp[i]表示以a[i]结尾的最长递增子序列长度
int main()
{
int n;
while(cin>>n)
{
for(int i=0; i<n; i++)
{
cin>>a[i];
dp[i]=1; // 初始化为1,长度最短为自身
}
int ans=1;
for(int i=1; i<n; i++)
{
for(int j=0; j<i; j++)
{
if(a[i]>a[j])
{
dp[i]=max(dp[i],dp[j]+1); // 状态转移
}
}
ans=max(ans,dp[i]); // 比较每一个dp[i],最大值为答案
}
cout<<ans<<endl;
}
return 0;
}
O(NlogN)做法:贪心+二分
a[i]表示第i个数据。
dp[i]表示表示长度为i+1的LIS结尾元素的最小值。
利用贪心的思想,对于一个上升子序列,显然当前最后一个元素越小,越有利于添加新的元素,这样LIS长度自然更长。
因此,我们只需要维护dp数组,其表示的就是长度为i+1的LIS结尾元素的最小值,保证每一位都是最小值,这样子dp数组的长度就是LIS的长度。
dp数组具体维护过程同样举例讲解更为清晰。
同样对于序列 a(1, 7, 3, 5, 9, 4, 8),dp的变化过程如下:
- dp[0] = a[0] = 1,长度为1的LIS结尾元素的最小值自然没得挑,就是第一个数。 (dp = {1})
- 对于a[1]=7,a[1]>dp[0],因此直接添加到dp尾,dp[1]=a[1]。(dp = {1, 7})
- 对于a[2]=3,dp[0]< a[2]< dp[1],因此a[2]替换dp[1],令dp[1]=a[2],因为长度为2的LIS,结尾元素自然是3好过于7,因为越小这样有利于后续添加新元素。 (dp = {1, 3})
- 对于a[3]=5,a[3]>dp[1],因此直接添加到dp尾,dp[2]=a[3]。 (dp = {1, 3, 5})
- 对于a[4]=9,a[4]>dp[2],因此同样直接添加到dp尾,dp[3]=a[9]。 (dp = {1, 3, 5, 9})
- 对于a[5]=4,dp[1]< a[5]< dp[2],因此a[5]替换值为5的dp[2],其实这步操作对于求LIS长度没什么意义,有意义的是第三步将dp[1]=7换成dp[1]=3,因为无论是1 3 5 9还是1 2 4 9都不影响最终的LIS长度,当然具体题目具体分析,如hdu1257导弹拦截系统中,每个系统的最高拦截高度进行更新更好理解(虽然题目求的是最少需要多少套拦截系统,万一类似的题目要求输出每个拦截系统的最后高度情况呢==)(dp = {1, 3, 4, 9})
- 对于a[6]=8,dp[2]< a[6]< dp[3],同理a[6]替换值为9的dp[3],道理你懂。 (dp = {1, 3, 5, 8})
ok,这样子dp数组就维护完毕,所求LIS长度就是dp数组长度4。
通过上述求解,可以发现dp数组是单调递增的,因此对于每一个a[i],先判断是否可以直接插入到dp数组尾部,即比较其与dp数组的最大值即最后一位;如果不可以,则找出dp中第一个大于等于a[i]的位置,用a[i]替换之。
这个过程可以利用二分查找,因此查找时间复杂度为O(logN),所以总的时间复杂度为O(NlogN)
#include <bits/stdc++.h>
using namespace std;
const int MAXX=100000+5;
const int INF=INT_MAX;
int a[MAXX],dp[MAXX]; // a数组为数据,dp[i]表示长度为i+1的LIS结尾元素的最小值
int main()
{
int n;
while(cin>>n)
{
for(int i=0; i<n; i++)
{
cin>>a[i];
dp[i]=INF; // 初始化为无限大
}
int pos=0; // 记录dp当前最后一位的下标
dp[0]=a[0]; // dp[0]值显然为a[0]
for(int i=1; i<n; i++)
{
if(a[i]>dp[pos]) // 若a[i]大于dp数组最大值,则直接添加
dp[++pos] = a[i];
else // 否则找到dp中第一个大于等于a[i]的位置,用a[i]替换之。
dp[lower_bound(dp,dp+pos+1,a[i])-dp]=a[i]; // 二分查找
//如果题目中要求的是严格递增的话 dp[upper_bound(dp,dp+pos+1,a[i])-dp]=a[i];
//找到dp中第一个大于a[i]的位置
}
cout<<pos+1<<endl;
}
return 0;
}
一道例题
洛谷【p1439】排列LCS问题
题目描述:
给出1-n的两个排列P1和P2,求它们的最长公共子序列。
输入格式:
第一行是一个数n,
接下来两行,每行为n个数,为自然数1-n的一个排列。
输出格式:
一个数,即最长公共子序列的长度
输入样例:
5
3 2 1 4 5
1 2 3 4 5
输出样例:
3
【数据规模】
对于50%的数据,n≤1000
对于100%的数据,n≤100000
求解:
显然这是一道LCS最长公共子序列的问题,但是N特别大,对于LCS的O(N^2)做法,时间肯定得爆吧,由于题目中的两个序列都是1~n的一个排列,噢很巧哦,由于这个条件LCS的问题可以用LIS来解决。
若其中一个排列是1,2,3…n,那么他们的LCS(最长公共子序列)就是就是另一个序列的LIS(最长上升子序列)。如果两个序列的排列都不是1,2,3…n,那么我们可以认为其中一个序列是1,2,3..n,然后把第一个序列的a[1]映射到1,a[2]映射到2,a[n]映射到n,对b序列也按照a序列的映射规则处理,这样再求b序列的LIS即可。
#include <bits/stdc++.h>
using namespace std;
const int MAXX=100000+5;
const int INF=INT_MAX;
int dp[MAXX];
int a[MAXX],b[MAXX];
int main()
{
int n;
while(cin>>n)
{
int v;
for(int i=0; i<n; i++)
{
cin>>v;
a[v]=i;//把a[i]映射到i
}
for(int i=0; i<n; i++)
{
cin>>v;
b[i]=a[v];//把b数组按照a数组的映射规则进行映射
}
for(int i=0; i<n; i++)
dp[i]=INF; //初始化
int pos=0; // 记录dp当前最后一位的下标
dp[0]=b[0];
for(int i=1; i<n; i++)
{
if(b[i]>=dp[pos])
dp[++pos]=b[i];
else
dp[lower_bound(dp,dp+pos+1,b[i])-dp]=b[i];
}
cout<<pos+1<<endl;
}
return 0;
}
HDU1257
最少拦截系统
Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 65050 Accepted Submission(s): 25333
Problem Description
某国为了防御敌国的导弹袭击,发展出一种导弹拦截系统.但是这种导弹拦截系统有一个缺陷:虽然它的第一发炮弹能够到达任意的高度,但是以后每一发炮弹都不能超过前一发的高度.某天,雷达捕捉到敌国的导弹来袭.由于该系统还在试用阶段,所以只有一套系统,因此有可能不能拦截所有的导弹.
怎么办呢?多搞几套系统呗!你说说倒蛮容易,成本呢?成本是个大问题啊.所以俺就到这里来求救了,请帮助计算一下最少需要多少套拦截系统.
Input
输入若干组数据.每组数据包括:导弹总个数(正整数),导弹依此飞来的高度(雷达给出的高度数据是不大于30000的正整数,用空格分隔)
Output
对应每组数据输出拦截所有导弹最少要配备多少套这种导弹拦截系统.
Sample Input
8 389 207 155 300 299 170 158 65
Sample Output
2
Source
AC代码:
#include<bits/stdc++.h>
using namespace std;
int main()
{
int n,i,j,k,dp[10000],a[10000];
while(cin>>n)
{
memset(dp,0,sizeof(dp));
memset(a,0,sizeof(a));
for(i=0;i<n;i++)
cin>>a[i];
dp[1]=a[0];
k=1;
for(i=1;i<n;i++)
{
if(a[i]>dp[k]){
dp[++k]=a[i];
}else{
int pos=lower_bound(dp,dp+k+1,a[i])-dp;
dp[pos]=a[i];
}
}
cout<<k<<endl;
}
}
---------------------
作者:wbin233
来源:CSDN
原文:https://blog.csdn.net/wbin233/article/details/77570070
版权声明:本文为博主原创文章,转载请附上博文链接!