题目
UVa 11404 Palindromic Subsequence
题解
刚开始想的是倒过来求LCS,但是不知道如何保证字典序,于是看了看题解发现基本上都是用string直接记录路径,因为这个时候string虽然慢点但确实比char[]方便太多,但是有个问题是LCS不一定是答案,比如别人给的这组数据
Input
kfclbckibbibjccbej
(反向之后为jebccjbibbikcblcfk)Output
bcibbibc
然后发现前半部分或后半部分一定是一个可行解的半部分,至于为什么会这样我不是很清楚,所以未采用这种写法,若有人知道还请指教。
用dp[i][j].ans表示i~j的最大回文串长度,dp[i][j].s表示回文串(诶好像ans是不用保存的,反正长度也能测出来,不过可能要慢些),然后就分情况讨论啦,这部分比较简单看代码就懂。
代码
//QWsin
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=1000+10;
char s[maxn];
struct DP{
int ans;string s;DP(){}
DP(int ans,string s):ans(ans),s(s){}
}dp[maxn][maxn];
inline string make_string(char *s,int l,int r)//由于蒟蒻并不会使用string的构造函数,只好写了这么个鬼畜的东西
{
string ret;
for(int i=l;i<=r;++i) ret.push_back(s[i]);
return ret;
}
inline void solve()
{
int n=strlen(s+1);
for(int L=1;L<=n;++L)
for(int i=1;i+L-1<=n;++i)
{
int j=i+L-1;
if(i==j) dp[i][j]=DP(1,make_string(s,i,i));//边界
else if(i+1==j&&s[i]==s[j]) dp[i][j]=DP(2,make_string(s,i,j));//边界
else if(s[i]==s[j])
{
dp[i][j].ans=dp[i+1][j-1].ans+2;
dp[i][j].s=s[i]+dp[i+1][j-1].s+s[j];
}
else
{
int k1=dp[i+1][j].ans,k2=dp[i][j-1].ans;
if(k1 > k2) dp[i][j]=dp[i+1][j];
else if(k1 < k2) dp[i][j]=dp[i][j-1];
else {
dp[i][j].ans=k1;dp[i][j].s=min(dp[i+1][j].s,dp[i][j-1].s);
}
}
}
cout<<dp[1][n].s<<endl;
}
int main()
{
while(scanf("%s",s+1)==1) solve();
return 0;
}