poj 1699 TSP(状态压缩DP + KMP)

题意:给你n个串,让你构造出一个最短的串,使得这n个串每个串都是这个串的子串,输出你构造的串的长度即可
(如果要输出构造的串中字典序最小的,就有点难了)
数据范围小,不一定要使用KMP预处理的
分两步:
1:预处理出一个串的前缀最多匹配多少长度的另一个串的后缀

2:然后就相当于旅行商问题了,每两个串之间都有一定的距离,设dp[i][j]表示i状态的字符串,最后一个字符串为j时所构造的字符串的最短长度,每次枚举一个中间点进行转移即可,具体见代码

这题用ac自动机死活过不了,蛋疼啊

#include <set>
#include <map>
#include <stack>
#include <queue>
#include <deque>
#include <cmath>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define L(i) i<<1
#define R(i) i<<1|1
#define INF  0x3f3f3f3f
#define pi acos(-1.0)
#define eps 1e-9
#define maxn 10010
#define MOD 1000000007
int n,m;
char str[16][110];
int len[20][20];
int p[110];
void getp(int m,char b[]){
    p[1]=0;
    int i,j=0;
    for(i=2;i<=m;i++){
        while(j>0&&b[j+1]!=b[i]) j=p[j];
        if(b[j+1]==b[i]) j+=1;
        p[i]=j;
    }
}
int kmp(int n,int m,char a[],char b[]){
    int i,j=0,Max=0;
    for(i=1;i<=n;i++){
        while(j>0&&b[j+1]!=a[i]) j=p[j];
        if(b[j+1]==a[i]) j+=1;
        Max=j;
        if(j==m){
            j=p[j];
        }
    }
    return Max;
}
int dp[1<<15][15];
int pre[15];
int le[20];
int main()
{

    int t,ca=1,n;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        for(int i=0;i<n;i++)     scanf("%s",str[i]+1),le[i]=strlen(str[i]+1);
        memset(len,0,sizeof(len));
        for(int i=0;i<n;i++)
        {
            for(int j=i+1;j<n;j++)
            {
                memset(p,0,sizeof(p));
                getp(le[i],str[i]);
                int tmp=kmp(le[j],le[i],str[j],str[i]);
                if(tmp>len[i][j])
                    len[i][j]=tmp;
                getp(le[j],str[j]);
                tmp=kmp(le[i],le[j],str[i],str[j]);
                if(tmp>len[j][i])
                    len[j][i]=tmp;
            }
        }
        memset(dp,INF,sizeof(dp));
        for(int i = 0; i < n; i++)
            dp[1<<i][i] = le[i];
        for(int i = 0; i < (1<<n); i++)
            for(int j = 0; j < n; j++)
                for(int k = 0; k < n; k++)
                {
                    if(!(i&(1<<k)) || !(i&(1<<j)) || j == k)
                        continue;
                    dp[i|(1<<k)][k] = min(dp[i|(1<<k)][k],dp[i][j]+le[k]-len[k][j]);
                        //dp[i][j] = min(dp[i][j],dp[i^(1<<j)][k]+le[j]-len[j][k]);

                }
        int ans = INF;
        for(int i = 0; i < n; i++)
            ans = min(ans,dp[(1<<n)-1][i]);
        printf("%d\n",ans);
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值