UVA 11404 Palindromic Subsequence

逆序后求LCS，然后找出i+j == len情况的最长最小的LCS(子序列在1~i中出现，且其逆序列在i+1到len中出现，对应于回文串偶数情况)， (i+j == len-1情况的最长最小的LCS+str[i+1])(子序列在1~i 中出现，且其逆序在i+2到len中出现，注意其中包括子序列为空的情况，即回文串中只有一个字符，对应于回文串奇数情况)，从这俩种情况中找出最优解

#include <cstring>
#include <cstdio>
#include <iostream>
#include <string>
#include <algorithm>
using namespace std;

const int maxn = 1010;
struct pp
{
    int len;
    string s;
    bool operator <(const pp a)const
    {
        if(len==a.len) return s>a.s;
        return len<a.len;
    }
}dp[maxn][maxn],ans;
char s[maxn],rs[maxn];
int main()
{

    while(~scanf("%s",s+1))
    {
        int len=strlen(s+1);
        s[len+1]='*';
        //rs="";
        //cout<<s+1<<endl;
        for(int i=1;i<=len;i++)
            rs[i]=s[len-i+1];
        rs[len+1]='*';
        for(int i=0;i<=len;i++)
        {
            dp[i][0].len=dp[0][i].len=0;
            dp[i][0].s="";
            dp[0][i].s="";

        }
        for(int i=1;i<=len;i++)
        {
            for(int j=1;j<=len;j++)
            {
                //cout<<"?"<<endl;
                if(s[i]==rs[j])
                {
                   // cout<<"1"<<endl;
                    dp[i][j].len=dp[i-1][j-1].len+1;
                    dp[i][j].s=dp[i-1][j-1].s+s[i];
                }
                else
                {
                    if(dp[i][j-1]<dp[i-1][j])
                    {
                       // cout<<"2"<<endl;
                        dp[i][j].len=dp[i-1][j].len;
                        dp[i][j].s=dp[i-1][j].s;
                    }
                    else
                    {
                       // cout<<"3"<<endl;
                        dp[i][j].len=dp[i][j-1].len;
                        dp[i][j].s=dp[i][j-1].s;
                    }
                }
               // cout<<dp[i][j].s<<endl;
            }

        }
        ans.len=0;
        ans.s="";
        for(int i=1;i<=len;i++)
        {
            if(ans<dp[i][len-i]) ans=dp[i][len-i];
        }
        ans.len*=2;
        for(int i=0;i<len;i++)
        {
            pp tmp=dp[i][len-i-1];
            tmp.len=tmp.len*2+1;
            tmp.s+=s[i+1];
            if(ans<tmp) ans=tmp;
        }
        //cout<<ans.len<<endl;
        if(ans.len&1)
        {
            cout<<ans.s;
            for(int i=ans.len/2-1;i>=0;i--)
                cout<<ans.s[i];
        }
        else
        {
            cout<<ans.s;
            for(int i=ans.len/2-1;i>=0;i--)
                cout<<ans.s[i];
        }
        cout<<endl;
    }
}