翻转序列（动态规划）

题目

描述

你有一个长度为 $n$ 的序列 $a_1,a_2,...,a_n$ 。你现在可以选择一个区间 $(L,R)$，然后翻转这个区间，求操作后的序列非递减子序列的长度最大。

一个非递减子序列是一个索引为$p_1,p2,...,p_k$ 的序列，使得 $p_1<p_2<...<p_k$和$a_{p_1}\le a_{p_2}\le ...\le a_{p_k}$ 。子序列的长度为 $k$。

输入

第一行包含一个整数 $n（1\le n\le 2000）$ ，表示原始序列的长度。

第二行包含 $n$ 个以空格分隔的整数，描述原始序列$a_1,a_2,...,a_n(1⩽a_i⩽2)$ 。

输出

打印一个整数，这意味着新序列的最长非递减子序列的最大可能长度。

样例

输入

4
1 2 1 2

输出

4

提示

【样例解释】

你可以选择翻转 (2，3)，新序列变成 [1，1，2，2] ，答案是 44 。

【数据说明】

分析

由于翻转某块后使得存在最长的非递减子序列，假设翻转的左右区间为翻转（L，R），那么我们可以将这个子序列分为三段，为（1，L-1）、（L，R），（R+1，n），第一段为全1，第二段在翻转前为递减子序列，第二段为全2。
这时有人会问为什么第一段为全1呢，因为假如第一段为111…1122…222的形式后面的2222一定可以并入第二段的递减子序列。
所以第一段和第三段可以用前缀数组和后缀数组轻易得出。
第二段可以用动态规划解决

dp[i][j][0]为从i到j以2结尾的最长递减序列，由于2是这个序列中最大的数
字，那么即为i到j中2出现的次数。
dp[i][j][1]代表从i到j以1结尾的最长递减序列.
dp[i][j][0]=dp[i][j-1][0]+(a[i]==2);
dp[i][j][1]=max(dp[i][j-1][0],dp[i][j-1][1])+(a[i]==1);

代码

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const ll mod=1e9+7;

ll a[2009];
ll dp[2009][2009][2];
ll pre[2009],suf[2009];
void solve()
{
	ll n,i,j;
	cin>>n; 
	for(i=1;i<=n;i++) cin>>a[i];
	for(i=1;i<=n;i++) pre[i]=pre[i-1]+(a[i]==1);
	for(i=n;i>=1;i--) suf[i]=suf[i+1]+(a[i]==2);
	ll ans=0;
	for(i=1;i<=n;i++)
	{
		for(j=i;j<=n;j++)
		{
			dp[i][j][0]=dp[i][j-1][0]+(a[j]==2);
			dp[i][j][1]=max(dp[i][j-1][0],dp[i][j-1][1])+(a[j]==1);
			ans=max(ans,dp[i][j][0]+pre[i-1]+suf[j+1]);
			ans=max(ans,dp[i][j][1]+pre[i-1]+suf[j+1]);
		}
	}
	cout<<ans<<endl;
}
int main(){
    ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);

    int t=1;
    //cin>>t;
    while(t--)
    {
        solve();
    }
    return 0;
}