bzoj 1195: [HNOI2006]最短母串【状压dp】

我有病吧……明明直接枚举是否匹配就可以非要写hash，然后果然冲突了（……我个非酋居然还敢用hash
设f[s][i]为已选串状态为s并且最后一个串是i，还有预处理出g[i][j]表示最长有长为g[i][j]的i串后缀等于j串前缀这里，直接暴力匹配即可……
然后注意到比较麻烦的事要求字典序最小，但是因为空间限制我们又不能给每个f存一个串，所以我们设t[s][i][j]为状态使f[s][i]长度最小且字典序最小的选串顺序的第j个选的是那个串，这个更新f[s][i]的时候直接更新，注意如果有等于当前f[s][i]的情况，要分别把旧的t[s][i]和当前情况对应的两个串都求出来比一下字典序
长度dp部分随便转移一下就好

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=15;
int n,f[5005][N],t[5005][N][N],g[N][N],len[N],ans=1e9,b[N];
bool v[N];
char c[N][55],s1[1005],s2[1005],a[1005];
struct qwe
{
    int l;
    char c[55];
};
bool cmp(const qwe &a,const qwe &b)
{
    for(int i=1;i<=min(a.l,b.l);i++)
        if(a.c[i]!=b.c[i])
            return a.c[i]<b.c[i];
    return a.l<b.l;
}
bool ok(int i,int j,int l)
{
    for(int k=1;k<=l;k++)
        if(c[i][len[i]-l+k]!=c[j][k])
            return 0;
    return 1;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%s",c[i]+1);
        len[i]=strlen(c[i]+1);
        b[i]=1<<(i-1);
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            for(int l=min(len[i],len[j]);l>=0;l--)
                if(ok(i,j,l))
                {
                    g[i][j]=l;
                    break;
                }
    // dfs(1,0);
    memset(f,0x3f,sizeof(f));
    int inf=f[0][0];
    for(int i=1;i<=n;i++)
        f[b[i]][i]=len[i],t[b[i]][i][1]=i;
    for(int s=1,le=(1<<n)-1;s<=le;s++)
    {
        int con=0;
        for(int i=1;i<=n;i++)
            if(s&b[i])
                con++;
        for(int i=1;i<=n;i++)
            if(s&b[i])
            {
                for(int j=1;j<=n;j++)
                    if(!(s&b[j]))
                    {
                        int nw=f[s][i]+len[j]-g[i][j],w=(s|b[j]);
                        if(nw<f[w][j])
                        {
                            f[w][j]=nw;
                            for(int k=1;k<=con;k++)
                                t[w][j][k]=t[s][i][k];
                            t[w][j][con+1]=j;
                        }
                        else if(nw==f[w][j])
                        {
                            int t1=0,t2=0,fl=0;
                            for(int k=1;k<=con+1;k++)
                                for(int l=1+g[t[w][j][k-1]][t[w][j][k]];l<=len[t[w][j][k]];l++)
                                    s1[++t1]=c[t[w][j][k]][l];
                            for(int k=1;k<=con;k++)
                                for(int l=1+g[t[s][i][k-1]][t[s][i][k]];l<=len[t[s][i][k]];l++)
                                    s2[++t2]=c[t[s][i][k]][l];
                            for(int l=1+g[t[s][i][con]][j];l<=len[j];l++)
                                s2[++t2]=c[j][l];
                            for(int k=1;k<=nw;k++)
                            {
                                if(s1[k]<s2[k])
                                    break;
                                if(s1[k]>s2[k])
                                {
                                    fl=1;
                                    break;
                                }
                            }
                            if(fl)
                            {
                                for(int k=1;k<=con;k++)
                                    t[w][j][k]=t[s][i][k];
                                t[w][j][con+1]=j;
                            }
                        }
                    }
            }
    }
    // printf("%d\n",ans);
    int t1=0,w=(1<<n)-1,ans=1e9;
    for(int j=1;j<=n;j++)
    {
        if(f[w][j]<ans)
        {
            ans=f[w][j];
            int t1=0;
            for(int k=1;k<=n;k++)
                for(int l=1+g[t[w][j][k-1]][t[w][j][k]];l<=len[t[w][j][k]];l++)
                    a[++t1]=c[t[w][j][k]][l];
        }
        else if(f[w][j]==ans)
        {
            int t1=0,fl=0;
            for(int k=1;k<=n;k++)
                for(int l=1+g[t[w][j][k-1]][t[w][j][k]];l<=len[t[w][j][k]];l++)
                    s1[++t1]=c[t[w][j][k]][l];
            for(int k=1;k<=ans;k++)
            {
                if(a[k]<s1[k])
                    break;
                if(a[k]>s1[k])
                {
                    fl=1;
                    break;
                }
            }
            if(fl)
                for(int k=1;k<=ans;k++)
                    a[k]=s1[k];
        }
    }
    for(int i=1;i<=ans;i++)
        printf("%c",a[i]);
    return 0;
}

转载于:https://www.cnblogs.com/lokiii/p/9644795.html