先上模板:
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
#include<queue>
#include<stack>
#include<vector>
#include<map>
#include<set>
#include<sstream>
#include<cstdlib>
#include<time.h>
#define inf 1000000000
#define mem(a,b) memset((a),b,sizeof(a))
typedef long long ll;
const int N=100010;
using namespace std;
int n,k;
int a[100010];
int dp[100010];
int LIS() //最长上升子序列
{
for(int i=0; i<=n; i++)
dp[i]=inf;
int pos=1;
dp[1]=a[1];
for(int i=2; i<=n; i++)
{
if(a[i]>dp[pos])
dp[++pos] = a[i];
else
dp[lower_bound(dp+1,dp+pos+1,a[i])-dp]=a[i];
}
return pos;
}
//int LIS() //最长不降子序列
//{
// for(int i=0; i<=n; i++)
// dp[i]=inf;
// int pos=1;
// dp[1]=a[1];
// for(int i=2; i<=n; i++)
// {
// if(a[i]>=dp[pos])
// dp[++pos] = a[i];
// else
// dp[upper_bound(dp+1,dp+pos+1,a[i])-dp]=a[i];
// }
// return pos;
//}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
scanf("%d",&n);
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
int ans=LIS();
printf("%d\n",ans);
}
return 0;
}
#include<iostream>
#include<iomanip>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<string>
#include<queue>
#include<stack>
#include<vector>
#include<map>
#include<set>
#include<sstream>
#include<cstdlib>
#define rep(i,a,n) for(int i=a;i<=n;i++)
#define repp(i,n,a) for(int i=n;i>=a;i--)
#define ll long long
#define ull unsigned long long
#define mem(a,x) memset((a),(x),sizeof ((a)))//x只能是0或-1或false或true
#define debug(x) cout<<"X: "<<(x)<<endl
#define de cout<<"************"<<endl
#define lowbit(x) ((x)&(-x))
#define lson rt<<1
#define rson rt<<1|1
#define gcd(a,b) __gcd(a,b)
#define lcm(a,b) a*b/(__gcd(a,b))
#define inf 0x3f3f3f3f//1e9+6e7
#define eps 1e-8
#define mod 1e9+7
#define N 100010
#define maxn 40010
const double pi=acos(-1.0);
using namespace std;
int a[40010];
int dp[maxn];//dp[i]表示长度为i+1的LIS结尾元素的最小值
int n;
int LIS(int a[]) //返回最长上升子序列的长度 O(n*logn)
{
dp[0]=a[0];
int pos=0;
for(int i=1; i<n; i++)
{
if(a[i]>dp[pos])
dp[++pos] = a[i];
else
dp[lower_bound(dp,dp+pos+1,a[i])-dp]=a[i];
}
return pos+1;
}
int main()
{
int T;
cin>>T;
while(T--)
{
scanf("%d",&n);
for(int i=0; i<n; i++)
{
scanf("%d",&a[i]);
dp[i]=inf;
}
int ans=LIS(a);
cout<<ans<<endl;
}
return 0;
}
转自:https://blog.csdn.net/wbin233/article/details/77570070
LIS(Longest Increasing Subsequence)最长上升子序列
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],因此长度为3的LIS,结尾元素为4会比5好,越小越好嘛。(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]; // 二分查找
}
cout<<pos+1<<endl;//输出最长上升子序列的长度
}
return 0;
}