hdu--4681--dp

这个dp应该算lcs加上一些模拟处理吧

lcs我们大家都会求 很简单。

这边的话 就是先算出 a b两个字符串从前往后 与 从后往前的lcs

因为 我们在算d的长度的时候 是向两边进行拓展的。

然后 拿c分别于a和b去暴力匹配 分别找出以a b中的字符串中的每个位置起始点 可以包含c这个字符串的子序列的最左边位置

这边 一定是最左边的 这样就会尽可能地留出多余位置 去和另一个字符串进行 lcs 这样就是使拓展长度可能地增加。

这题 我在work函数 那边花了超长时间 很容易出错= = 。。   

#include <iostream>
#include <cstring>
#include <vector>
#include <utility>
#include <algorithm>
using namespace std;

const int size = 1010;
char a[size] , b[size] , c[size];
int preDp[size][size] , nextDp[size][size];
int lenA , lenB , lenC;
vector< pair<int,int> > va , vb;

void preLcs( )
{
    for( int i = 1 ; i<=lenA ; i++ )
    {
        for( int j = 1 ; j<=lenB ; j++ )
        {
            if( a[i]==b[j] )
                preDp[i][j] = preDp[i-1][j-1] + 1;
            else
                preDp[i][j] = max( preDp[i][j-1] , preDp[i-1][j] );
        }
    }
}

void nextLcs( )
{
    for( int i = lenA ; i>=1 ; i-- )
    {
        for( int j = lenB ; j>=1 ; j-- )
        {
            if( a[i]==b[j] )
                nextDp[i][j] = nextDp[i+1][j+1] + 1;
            else
                nextDp[i][j] = max( nextDp[i][j+1] , nextDp[i+1][j] );
        }
    }
}

void solve( char* str , int num , vector< pair<int,int> >& ve )
{
    int cnt , i , k;
    if( lenC==1 )
    {
        for( i = 0 ; i<num ; i++ )
        {
            if( str[i]==c[1] )
            {
                ve.push_back( make_pair(i+1,i+1) );
            }
        }
    }
    else
    {
        for( i = 0 ; i<num ; i++ )
        {
            cnt = 1;    
            if( str[i]==c[cnt] )
            {
                ++ cnt;
                for( k = i+1 ; k<num ; k++ )
                {
                    if( str[k]==c[cnt] )
                    {
                        ++ cnt;
                    }    
                    if( cnt==lenC+1 )
                    {
                        ve.push_back( make_pair(i+1,k+1) );
                        break;
                    }
                }
            }    
        }
    }
}

int work( )
{
    int vAsize = va.size();
    int vBsize = vb.size();
    int ans = 0;
    for( int i = 0 ; i<vAsize ; i++ )
    {
        
        for( int j = 0 ; j<vBsize ; j++ )
        {
            ans = max( ans , preDp[ va[i].first-1 ][ vb[j].first-1 ]+nextDp[ va[i].second+1 ][ vb[j].second+1 ]+lenC );
        }
    }
    return ans;    
}

int main()
{
    cin.sync_with_stdio(false);
    int t;
    cin >> t;
    for( int k = 1 ; k<=t ; k++ )
    {
        cin >> (a+1) >> (b+1) >> (c+1);
        lenA = strlen( a+1 );
        lenB = strlen( b+1 );
        lenC = strlen( c+1 );
        memset( preDp , 0 , sizeof(preDp) );
        memset( nextDp , 0 , sizeof(nextDp) );
        va.clear();
        vb.clear();
        preLcs( );
        nextLcs( );
        solve( (a+1) , lenA , va );
        solve( (b+1) , lenB , vb );
        cout << "Case #" << k << ": " << work() << endl;
    }
    return 0;
}