HDU 6357 Hills And Valleys (DP)

题目链接 给出一个由0-9构成的串，允许将其中的一个子区间翻转，求最长不下降子序列的长度及翻转的区间。

枚举翻转区间的端点显然是不行的，可以换个思路，枚举翻转区间的值域。

假设翻转的区间是[x,y]，那么参与构成最长不下降子序列的元素在翻转之前必然是这样的：

0`,1`,.....(x-1)`,x`,y`,(y-1)`,...(x+1)`,x`,y`,(y+1)`.....8`,9`。其中k`表示有非负个k出现。

枚举[x,y]并构造出上面的这个数列b，用数列a表示原来的串，考虑用dp[i][j]表示a串的第i个位置匹配到b串的第j个位置时的最长不下降子序列的长度。显然有dp[i][j] = max(dp[i][j - 1] , dp[i - 1][j] + (a[i] == b[j]))。可以另开数组来存储i和j匹配时翻转区间的左右端点。

#include <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
#include <map>
#include <vector>
#include <iostream>
#include <algorithm>
#define INF 0x3f3f3f3f
#define eps 1e-8
using namespace std;
typedef unsigned int ul;
const int maxn = 100005;
const int maxm = 22;
typedef long long ll;
const ll mod = (1LL << 30);

int kase, n, a[maxn], dp[maxn][maxm], b[maxn];
int ans, ansl, ansr, l, r, no;
int al[maxn][maxm], ar[maxn][maxm];
char s[maxn];

int solve()
{
    for(int i = 0;i <= no;i++) dp[0][i] = 0;
    for(int i = 1;i <= n;i++)
    {
        for(int j = 1;j <= no;j++)
        {
            dp[i][j] = dp[i - 1][j];
            al[i][j] = al[i - 1][j];
            ar[i][j] = ar[i - 1][j];
            if(a[i] == b[j])
            {
                dp[i][j] = dp[i - 1][j] + 1;
                if(l == j && !al[i][j]) al[i][j] = i;
                if(r == j) ar[i][j] = i;
            }
            if(dp[i][j - 1] > dp[i][j])
            {
                dp[i][j] = dp[i][j - 1];
                al[i][j] = al[i][j - 1];
                ar[i][j] = ar[i][j - 1];
            }
        }
    }
    return dp[n][no];
}

int main()
{
    scanf("%d", &kase);
    while(kase--)
    {
        scanf("%d%s", &n, s + 1);
        for(int i = 1;i <= n;i++)
            a[i] = s[i] - '0';
        no = 1;
        for(int i = 0;i <= 9;i++) b[no++] = i;
        no--;
        ansl = ansr = l = r = 1;
        ans = solve();
        for(int i = 0;i <= 9;i++)
        {
            for(int j = i + 1;j <= 9;j++)
            {
                no = 1;
                for(int k = 0;k <= i;k++) b[no++] = k;
                l = no;
                for(int k = j;k >= i;k--) b[no++] = k;
                r = no - 1;
                for(int k = j;k <= 9;k++) b[no++] = k;
                no--;
                //memset(al, 0, sizeof(al));
                //memset(ar, 0, sizeof(ar));
                int tmp = solve();
                if(ans < tmp && al[n][no] && ar[n][no])
                {
                    ans = tmp;
                    ansl = al[n][no], ansr = ar[n][no];
                }
            }
        }
        printf("%d %d %d\n", ans, ansl, ansr);
    }
    return 0;
}