合唱队形(LIS) 动态规划

问题

题目：[百练-2711]

思路

动态规划最重要的是第一步和第二步，也就是说一定要定义到合适的状态。并且找到状态转移的函数。这个题目如果只看一边，那就是一个LIS题目，所以问题在于怎么样结合LIS解题。

本题本质是求最长的合唱队形，不妨考虑如下的状态定义。设dp_inc[i]代表以第i个位置的元素结尾的LIS，dp_dec[i]代表以第i个位置的元素开始的LDS。那么对于第i个位置的元素来说，dp_inc[i] + dp_dec[i] - 1。就是以该点为中轴的合唱队形长度。

代码

#include <iostream>
#include <fstream>
#include <cstring>
//#define LOCAL

const int MAXN = 100;

int dp_inc[MAXN + 5];
int dp_dec[MAXN + 5];

int arr[MAXN + 5];

int cal_lis( int* dp, int n );
int cal_lds( int* dp, int n );

int main( void )
{
#ifdef LOCAL
    std::ifstream cin( "input.dat" );
#endif
    int n = 0;
    while( std::cin >> n )
    {
        for( int i = 0; i < n; ++i )
        {
            std::cin >> arr[i];
        }
        cal_lis( dp_inc, n );
        cal_lds( dp_dec, n );

        int max = 0;
        for( int i = 0; i < n; ++i )
        {
            int t = dp_inc[i] + dp_dec[i] - 1;
            max = ( t > max ) ? t : max;
        }
        std::cout << n - max << std::endl;
    }
#ifdef LOCAL
    cin.close();
#endif
    return 0;
}

int cal_lis( int* dp, int n )
{
    std::memset( dp, 0, n*sizeof(int) );

    dp[0] = 1;
    for( int i = 1; i < n; ++i )
    {
        dp[i] = 1;
        for( int j = 0; j < i; ++j )
        {
            if( arr[j] < arr[i] && dp[j] + 1 > dp[i] )
                dp[i] = dp[j] + 1;
        }
    }

    return 0;
}
int cal_lds( int* dp, int n )
{
    std::memset( dp, 0, n*sizeof(int) );

    dp[n-1] = 1;
    for( int i = n-2; i >= 0; --i )
    {
        dp[i] = 1;
        for( int j = n-1; j > i; --j )
        {
            if( arr[j] < arr[i] && dp[j] + 1 > dp[i] )
                dp[i] = dp[j] + 1;
        }
    }

    return 0;
}

总结

其实刚开始做的时候没有想到，明明已经知道LIS。但是还是没有想到。所以，一定要按照科学的方法论来做题。既然DP有着明确的步骤，并且这道题目又知道一定是要结合LIS。那么，按着步骤来。先考虑状态定义，不就很容易想到方法了嘛。所以，一定要冷静下来。按着总结的步骤来，DP就是状态定义和转移函数的事情。