【洛谷P3607】【USACO17JAN】Subsequence Reversal （DP）

       链接：洛谷P3607

       这题很明显就是$DP$，只是可以翻转一个子序列有些头疼，其实我们可以转化一下，翻转子序列其实是从内向外两两交换，这样就可以在转移的同时进行了。

       设$dp_{i,j,k,l}$表示下标在$[i,j]$内，最小值不小于$k$，最大值不大于$l$的最多进行一次翻转得到的最长的不下降子序列的长度。那么转移很显然，无非就是每次判断是否需要交换$a_i$和$a_j$，分别更新，然后再用$max(max(dp_{i+1,j,k,l},dp_{i,j-1,k,l}),max(dp_{i,j,k+1,l},dp_{i,j,k,l-1}))$更新$dp_{i,j,k,l}$即可。

       最终答案即为$dp_{1,n,1,50}$。

       时间复杂度$\Theta(50^2n^2)$。

       蒟蒻的代码：

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
using namespace std;
int n;
int dp[51][51][51][51];
int a[51];
inline void getmax(int &x,int y){x=max(x,y);}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)scanf("%d",a+i);
    for(int i=1;i<=n;i++)dp[i][i][a[i]][a[i]]=1;
    for(int l=2;l<=n;l++)
    {
        for(int i=1;i+l-1<=n;i++)
        {
            int j=i+l-1;
            for(int ln=1;ln<=50;ln++)
            {
                for(int L=1;L+ln-1<=50;L++)
                {
                    int R=L+ln-1;
                    getmax(dp[i][j][L][R],max(dp[i+1][j][L][R],dp[i][j-1][L][R]));
                    getmax(dp[i][j][L][R],max(dp[i][j][L+1][R],dp[i][j][L][R-1]));
                    getmax(dp[i][j][min(L,a[i])][R],dp[i+1][j][L][R]+(a[i]<=L));
                    getmax(dp[i][j][L][max(R,a[j])],dp[i][j-1][L][R]+(a[j]>=R));
                    getmax(dp[i][j][min(L,a[j])][R],dp[i+1][j-1][L][R]+(a[j]<=L));
                    getmax(dp[i][j][L][max(R,a[i])],dp[i+1][j-1][L][R]+(a[i]>=R));
                    getmax(dp[i][j][min(L,a[j])][max(R,a[i])],dp[i+1][j-1][L][R]+(a[j]<=L)+(a[i]>=R));
                }
            }
        }
    }
    printf("%d",dp[1][n][1][50]);
    return 0;
}