UVA 11404 Plalidromic Subsquence （回文子序列，LCS）

最长回文子序列可以用求解原串s和反转串rv的LCS来得到，因为要求回文串分奇偶，dp[i][j]保存长度，

要求字典序最小，dp[i][j]应该表示回文子序列的端点，所以边界为单个字符，即i+j=len+1。

这题最麻烦的地方在于字典序，我是写了个比较函数，有点暴力（常数大）。

也可以反着定义，这时结点就要保存那个状态的字符串了（这样定义比较字典序的时候常数小）

#include<bits/stdc++.h>
using namespace std;

#define MP make_pair
#define fi first
#define se second

const int LEN = 1e3+;
char s[LEN],rv[LEN];
int dp[LEN][LEN];
pair<int,int> pre[LEN][LEN];
char val[LEN][LEN];

inline void updata(int i,int j,int v,char c,const pair<int,int> &prv)
{
    dp[i][j] = v;
    pre[i][j] = prv;
    val[i][j] = c;
}

const auto nil = MP(,);

#define dim(x) [x.fi][x.se]
bool cmpLex(pair<int,int> a,pair<int,int> b)
{
    //
    while(a != nil && val[a.fi][a.se] == val[b.fi][b.se]){
        a = pre dim(a); b = pre dim(b);
    }
    return val dim(a) < val dim(b);
}

//#define LOCAL
int main()
{
#ifdef LOCAL
    freopen("in.txt","r",stdin);
#endif

    while(gets(s)){
        int len = strlen(s);
        for(int i = ; i < len; i++){
            rv[len--i] = s[i];
        }
        int hd1,hd2,vl = ;
        for(int i = ; i <= len; i++){
            int j = len+-i;
            dp[i][j] = ; val[i][j] = s[i-]; pre[i][j] = nil;
            if(dp[i][j] > vl || ( dp[i][j] == vl && cmpLex( MP(i,j), MP(hd1,hd2) ) ) ){//
                vl = dp[i][j];
                hd1 = i; hd2 = j;
            }
            for(int k = ; k < j; k++) dp[i][k] = ;
            for(j++; j <= len; j++){
                if(s[i-] == rv[j-]){
                    updata(i,j,dp[i-][j-]+,s[i-],make_pair(i-,j-));
                    if(dp[i][j] > vl || (dp[i][j] == vl && cmpLex(MP(i,j),MP(hd1,hd2)) ) ){//
                        vl = dp[i][j];
                        hd1 = i; hd2 = j;
                    }
                }else {
                    if(dp[i-][j] > dp[i][j-] || (dp[i-][j] == dp[i][j-] && cmpLex(MP(i-,j),MP(i,j-)) ) ){//
                        updata(i,j,dp[i-][j],val[i-][j],pre[i-][j]);
                    }else {
                        updata(i,j,dp[i][j-],val[i][j-],pre[i][j-]);
                    }
                }

            }
        }
        int pv = (vl+)>>,ln = vl;

        auto u = MP(hd1,hd2);
        for(int i = ; i < pv; i++){
            s[i] = val[u.fi][u.se];
            u = pre[u.fi][u.se];
        }

        s[ln] = '\0';
        for(int i = pv; i < ln; i++){
            s[i] = s[ln--i];
        }
        puts(s);
    }
    return ;
}