题目: 有一个长为n的数列,求出这个序列中最长上升子序列
n^2算法
首先我们定义dp[i]表示以a[i]为末尾的最长上升子序列的长度,这样我们就可以得到递推关系:对于在i以前的元素,只要满足a[j]<a[i],那么我们就可以尝试把a[i]放在j元素的后面得到一个新的子序列长度的值,即dp[i] = max(dp[i] , dp[j] + 1) (j<i && a[j] < a[i])
#include <iostream>
#include <algorithm>
using namespace std;
int a[1005] = {0};
int dp[1005] = {0};
int main()
{
int n,res = 0; //由于答案的最长上升子序列可能以任意一个元素结束,所以我们需要res来保存最大值
cin >> n;
for( int i = 0 ; i < n ; i++ )
{
cin >> a[i];
}
for( int i = 0 ; i < n ; i++ )
{
dp[i] = 1; //每一个元素至少自己都能算一个递增子序列,长度为1
for( int j = 0 ; j < i ; j++ )
{
if( a[j] < a[i] ) dp[i] = max(dp[i],dp[j]+1);
}
res = max(res,dp[i]);
}
cout << res << endl;
return 0;
}
nlogn算法
我们改变一个dp数组的含义,就可以利用二分优化到nlogn了
dp[i]表示长度为i+1的最长上升子序列中末尾元素的最小值,对于相同长度的最长递增子序列,末尾元素越小更加有优势,所以对于任意的a[j],只要在dp数组中找到第一个大于等于a[j]的元素dp[k],a[j]取而代之即可,这样我们就成功更新了长度为k+1的最长递增子序列末尾元素的最小值了。我们只要维护这个dp数组的大小就可以算出答案了。
#include <iostream>
#include <vector>
#include <algorithm>
using namespace std;
int a[1005] = {0};
vector<int> dp;
int main()
{
int n;
cin >> n;
for( int i = 0 ; i < n ; i++ )
{
cin >> a[i];
}
for( int i = 0 ; i < n ; i++ )
{
if( dp.size() == 0 ) dp.push_back(a[i]);
else if( a[i] > dp[dp.size()-1] ) dp.push_back(a[i]);
else *lower_bound(dp.begin(),dp.end(),a[i]) = a[i];
}
cout << dp.size() << endl;
return 0;
}