计蒜之道 百度AI小课堂-上升子序列

计蒜之道 百度AI小课堂-上升子序列

题目描述

给一个长度为 (n) 的数组 (a) 。试将其划分为两个严格上升子序列，并使其长度差最小。

输入格式

输入包含多组数据。

数据的第一行为一个正整数 (T) ，表示数据组数。

每组数据包括两行:

第一行包括一个正整数 (n)

第二行包括一个长度为 (n) 的数组 (a)。

输出格式

对于每组数据输出一行一个整数，表示两个子序列的最小长度差。若不存在划分方案则输出(-1)

数据范围

(T <= 10)
简单: (n <= 10^5, 0 <= a_i <= 2^{30}), 最多只存在一种合法的划分方案
中等: (n <= 10^3, 0 <= a_i <= 2^{30}), 划分方案不超过(10^{18})
困难: (n <= 10^5, 0 <= a_i <= 2^{30}), 划分方案不超过(10^{18})

思路：直接上题解

如果一个位置(k)满足对于(∀_i∈[1,k], ∀_j∈[k+1,n]， a_i < a_j)恒成立，那么称(k)为一个分界点。找出所有分界点后可将序列分为若干段，每段之间相互独立，同时每段只存在不超过一种合法的方案数。求出每段的划分(O(n))，再进行背包dp((O(nlogK))((K)为可能的划分方案数))

其实就是说对于某个(i),满足(max(1,i) < min(i + 1,n)),则序列([1,i]和[i+1,n])是相互独立的。
同理再来一个(j > i) 满足(max(1,j) < min(j + 1,n)), 则序列([1,j]和[j+1,n])相互独立， 序列([1,i])已经被划分出去，所以左边剩下([i+1,j])。
简单证明这样的每段如果合法,只存在一种划分。
即证明对于合法的一段([i,j](j > i + 1)), 若不存在(k)，满足(i<=k<j, max(i,k) < min(k+1,j)),那么段([i,j])只有一种划分方案。
反证不存在的情况有多种划分成立， 则存在的情况下只有一种划分，存在的情况下可以分成不可再分的几段，每段有多种划分，拼起来也会有多种划分， 与只有一种划分矛盾。
每段只有一种划分， (m)段组合起来的方案数为(2^m = K)， 所以(m = logK)

最后问题变成一个背包，从每段中选择一个序列长度组合起来，最后是否出现长度为(j)的序列((dp[j]))

#include<bits/stdc++.h>
#define LL long long 
using namespace std;
const int INF = 1e9;
const int N = 1e5 + 10;
int a[N];
int mi[N], mx[N];
int dp[2][N];
int T, n, o;
bool gao(int l, int r){
    int first = 0, second = 0, cnt = 0;
    for(int i = l;i <= r;i++){
        if(first < a[i]) first = a[i], cnt++;
        else if(second < a[i]) second = a[i];
        else return false;
    }
    int d[2] = {cnt, r - l + 1 - cnt};
    for(int i = 0;i <= n;i++) {
        dp[o ^ 1][i] = 0;
        for(int k = 0;k < 2;k++){
            if(i >= d[k]) dp[o ^ 1][i] |= dp[o][i - d[k]];
        }
    }
    o ^= 1;
    return true;
}
int main(){
	cin>>T;
	while(T--){
		cin>>n;
		for(int i = 1;i <= n;i++) cin>>a[i];
        if(n == 1){
            cout<<-1<<endl;
            continue;
        }
        mx[0] = -1, mi[n + 1] = (1<<30) + 1;
        for(int i = 1;i <= n;i++) mx[i] = max(mx[i - 1], a[i]);
        for(int i = n;i >= 1;i--) mi[i] = min(mi[i + 1], a[i]);
        memset(dp, 0, sizeof(dp));
        dp[0][0] = 1;
        int last = 1;
        bool flag = true;
        o = 0;
        for(int i = 1;i <= n;i++){
            if(mx[i] < mi[i + 1]){
                if(!gao(last, i)){
                    flag = false;
                    break; 
                }
                last = i + 1;
            }
        }
        int ans = -1;
        if(flag){
            for(int i = n / 2;i >= 0;i--) if(dp[o][i]){
                ans = n - 2 * i;
                break;
            }
        }
        cout<<ans<<endl;
	}
}