题解 luoguP2322 【[HNOI2006]最短母串问题】

传送门

看到 n &lt; = 12 n&lt;=12 n<=12的数据范围,一眼状压 d p ​ dp​ dp?随便开始敲了敲,一晚上就没了?

看到题解还有一篇状压,感觉挺麻烦的。这里给出可能稍微方便一点的。

状压状态和方程应该很好想, f [ i ] [ j ] f[i][j] f[i][j]表示状态 i i i下,最后的字符串是 j j j的最短长度,因为要输出方案,再来一个 g [ i ] [ j ] g[i][j] g[i][j]表示长度最短时的字符串。

做过几题字符串接龙状压题,应该很容易写出转移方程:
f [ t o ] [ k ] = m i n ( f [ t o ] [ k ] , f [ i ] [ j ] + c h o n g d i e ) f[to][k]=min(f[to][k],f[i][j]+chongdie) f[to][k]=min(f[to][k],f[i][j]+chongdie)
t o to to显然就是 i   ∣   ( 1 &lt; &lt; ( k − 1 ) ) i\ |\ (1&lt;&lt;(k-1)) i  (1<<(k1)),即状态中多接了一个 k k k c h o n g d i e chongdie chongdie k k k接在 j j j后面时两个字符串的重叠部分的长度,注意是有顺序的,是 k k k j j j后面。最后,当 f f f数组更新时,我们同时也更新 g g g数组。

具体分析一下转移:

1、 k k k j j j的子串,可以用 s t r i n g string string f i n d find find函数轻松解决, c h o n g d i e = 0 chongdie=0 chongdie=0

2、 k k k不是 j j j的子串,我们暴力匹配,枚举 j j j从哪一位开始与 k k k重叠,非常暴力,可以见代码。

答案显然易见,就是 m i n ( ∑ i = 1 n f [ ( 1 &lt; &lt; n ) − 1 ] [ i ] ) ) min(\sum\limits_{i=1}^{n}f[(1&lt;&lt;n)-1][i])) min(i=1nf[(1<<n)1][i])),然后取出所有长度最小的 g [ ( 1 &lt; &lt; n ) − 1 ] [ i ] g[(1&lt;&lt;n)-1][i] g[(1<<n)1][i],排遍序得到最小字典序的答案。

代码里还有部分注释,挺好懂的:

#include<bits/stdc++.h>
#define ts cout<<"ok"<<endl
#define ll long long
#define hh puts("")
using namespace std;
int n,f[5005][13],max_status,top;
string s[105],g[5005][13],ans[105];
inline int read(){
    int ret=0,ff=1;char ch=getchar();
    while(!isdigit(ch)){if(ch=='-') ff=-ff;ch=getchar();}
    while(isdigit(ch)){ret=(ret<<3)+(ret<<1)+(ch^48);ch=getchar();}
    return ret*ff;
}
inline int check(int x,int y){
    if(s[x].find(s[y])!=s[x].npos) return -1;//x包含y 
    for(int i=0;i<s[x].size();i++){//枚举j从哪一位开始与k重叠
        bool fl=1;
        for(int j=i;j<s[x].size();j++){
            if(s[y][j-i]!=s[x][j]){
                fl=0;
                break;
            }
        }
        if(fl) return s[x].size()-i;//返回重叠部分长 
    }
    return 0;
}
signed main(){
    n=read();
    for(int i=1;i<=n;i++) cin>>s[i];
    memset(f,0x3f,sizeof(f));
    f[0][0]=0;
    for(int i=1;i<=n;i++) f[1<<(i-1)][i]=s[i].length(),g[1<<(i-1)][i]=s[i];
    max_status=(1<<n)-1;
    for(int i=0;i<=max_status;i++){
        for(int j=1;j<=n;j++){//以j结尾 
            if(i&(1<<(j-1))){
                for(int k=1;k<=n;k++){//后面接k 
                    if(!(i&(1<<(k-1)))){
                        int to=i|(1<<(k-1));
                        int chong=check(j,k);//重叠部分长度 
                        if(chong==-1){//j包含k 
                            if(f[i][j]<f[to][j]||(f[i][j]==f[to][j]&&g[to][j]>g[i][j])){
                                f[to][j]=f[i][j];
                                g[to][j]=g[i][j];
                            }
                        }
                        else{
                            string t="";
                            for(int p=chong;p<s[k].size();p++) t+=s[k][p];//t是j接上k后多出来的一段字符 
                            t=g[i][j]+t;
                            if((f[i][j]+s[k].size()-chong<f[to][k])||
                            (f[i][j]+s[k].size()-chong==f[to][k]&&g[to][k]>t)){
                                f[to][k]=f[i][j]+s[k].size()-chong;
                                g[to][k]=t;
                            }
                        }
                    }
                }
            }
        }
    }
    int minn=1e9;
    for(int i=1;i<=n;i++){
        if(f[max_status][i]<minn){
            minn=f[max_status][i];
            top=1;
            ans[top]=g[max_status][i];
        }
        else if(f[max_status][i]==minn) ans[++top]=g[max_status][i];
    }
    sort(ans+1,ans+top+1);//排序找字典序最小 
    cout<<ans[1];
    return 0;
}
  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值