求最长上升子序列

DP O(N^2)

解法

$dp[i]$指以 $S[i]$ 结尾的最长上升子序列长度。

$dp[i]=max(dp[i],dp[k]+1)S[k]<S[i]，1\le k\le i-1$

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn = 2e5 + 10;
const int INF = 1e9 + 5;

int s[maxn], dp[maxn];
int n;

int main() {
	cin >> n;
	int res = 0;
	for (int i = 1; i <= n; i++)cin >> s[i];
	s[0] = -INF;
	for (int i = 1; i <= n; i++) {
		for (int j = i - 1; j >= 0; j--) {
		  //若求非严格应加上=，若求下降可对s[i]取负值
			if (s[j] < s[i])dp[i] = max(dp[i], dp[j] + 1);
			res = max(res, dp[i]);
		}
	}
	cout << res << "\n";
	return 0;
}

贪心+二分 O(N*logN)

解法

$dp[i]$ 指长度为 $i$ 的最长上升子序列的最小的结尾元素。

令 $cnt$ 为最长上升子序列长度。

遍历序列

• 若 $s[i]$大于当前的 $dp[cnt]$，说明$s[i]$可以插入长度为 $cnt$ 的序列，因此$cnt++$，并记录长度为 $cnt+1$的结尾元素。

• 否则，则可以查找并更新长度为$1-cnt$中大于 $S[i]$ 的第一个上升序列的结尾元素

  贪心性质，每个上升序列的结尾元素都应该尽量小。

  比如 $(1,2,1,2,4,5,3)$，当遍历到 $3$ 时

  $dp[1]=1,dp[2]=2,dp[3]=4,dp[4]=5$

  $value$已经比 $3$ 小的不能更新，比 $4$ 大的由于前面存在 $4$ ,因此他的 $value$ 必须比 $4$ 大，也不能更新。因此只能更新 $dp[3]=3$。

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn = 2e5 + 10;
const int INF = 1e9 + 5;

int s[maxn], dp[maxn];
int n;

int main() {
	cin >> n;
	int cnt = 0;
	for (int i = 1; i <= n; i++)cin >> s[i];
	dp[0] = -INF;
	for (int i = 1; i <= n; i++) {
		//若求非严格应加上=，若求下降可对s[i]取负值
		if (s[i] > dp[cnt])dp[++cnt] = s[i];
		else {
			//若求非严格应改成upper_bound
			int k = lower_bound(dp + 1, dp + 1 + cnt, s[i]) - dp;
			dp[k] = s[i];
		}
	}
	cout << cnt << "\n";
	return 0;
}

例题

AtCoder Regular Contest 133 B - Dividing Subsequence

题意

对于两个序列 $A，B$，找到最长的相同长度的子序列$a，b$，使得 $a[i]$ 是 $b[i]$ 的因子。

输出子序列长度。

题解

首先，我们先找到序列$A$中每个元素的倍数在序列$B$中的下标，同一个元素的下标属于同一个集合。

那么，问题就可以转化成求不属于同一个集合的下标的最长上升子序列。

我们在取一个下标集合时，将下标从大到小排序，即可保证每个集合自能取得一个下标。

最后，我们根据序列$A$的顺序，将排好序的下标集合存入，求其最长上升子序列即可。

代码

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int maxn = 2e5 + 10;
const int INF = 1e9 + 5;

int n, a[maxn], b[maxn], v[maxn];
vector<int>vec;
int dp[maxn];

int main() {
	cin >> n;
	for (int i = 1; i <= n; i++)cin >> a[i];
	for (int i = 1; i <= n; i++)cin >> b[i], v[b[i]] = i;
	for (int i = 1; i <= n; i++) {
		vector<int>pre;
		for (int j = 1; j * a[i] <= n; j++) pre.push_back(v[j * a[i]]);
		sort(pre.begin(), pre.end());
		int len = pre.size();
		for (int j = len - 1; j >= 0; j--)vec.push_back(pre[j]);
	}
	int len = vec.size();
	int cnt = 0; dp[0] = vec[0];
	for (int i = 1; i < len; i++) {
		if (vec[i] > dp[cnt])dp[++cnt] = vec[i];
		else {
			int k = lower_bound(dp, dp + cnt, vec[i]) - dp;
			dp[k] = vec[i];
		}
	}
	cout << cnt + 1;

	return 0;
}