bailian.openjudge 4126:DNA

最新推荐文章于 2022-10-19 14:29:27 发布

zxyOI

最新推荐文章于 2022-10-19 14:29:27 发布

阅读量535

点赞数

分类专栏： # 0x22 深度优先搜索

本文链接：https://blog.csdn.net/smartzxf/article/details/102078224

版权

0x22 深度优先搜索专栏收录该内容

13 篇文章 0 订阅

订阅专栏

http://bailian.openjudge.cn/practice/4126

描述

考虑一段DNA单链，上面有N个基因片段。这里的基因片段可重叠（例如AGCTC包含AGC和CTC），不可倒置（例如AGCTC不包含TCG）。要问这样的单链最短长度是多少。

输入

输入的第一行是一个正整数T（不超过13），表示数据组数。每组数据若干行，其中第一行一个正整数N（不超过9），表示基因片段的数目，接下来N行每行一个基因片段，由AGCT四个字母组成，且长度介于1和15之间（含两端）。

输出

每组数据输出一样，表示最短的单链长度包含这N个基因片段。

样例输入

1
5
TCGG
GCAG
CCGC
GATC
ATCG

样例输出

代码

#include<algorithm>
#include<iostream>
#include<string.h>
#include<stdio.h>
using namespace std;

const long long inf=1e9+7;
long long dp[1<<10][20];
long long length[15];
string str[15];
long long n;
long long visited[15][15];

long long TMP(long long l1,long long l2)
{
    long long ans=0;
    long long i=0;
    while(i<length[l1])
    {
        long long k=i;
        long long j=0;
        bool flag=true;
        while(k<length[l1] && j<length[l2])
        {
            if(str[l1][k]!=str[l2][j])
            {
                flag=false;
                break;
            }
            k++;
            j++;
        }
        
        if(flag)
        {
            ans=j;
            break;
        }
        i++;
    }

    return ans;
}

int main()
{
    long long T;
    cin>>T;
    while(T--)
    {
        cin>>n;
        for(long long i=0; i<(1<<n); i++)
        	for(long long j=0; j<n; j++)
            	dp[i][j]=inf;
        
        for(long long i=0; i<n; i++)
        {
            cin>>str[i];
            length[i]=str[i].length();
        }
        
        for(long long i=0; i<n; i++)
            for(long long j=0; j<n; j++)
            {
                if(i==j)
                {
                    visited[i][j]=0;
                    continue;
                }
                visited[i][j]=TMP(i,j);
            }

        for(long long i=0; i<n; i++)
            dp[1<<i][i]=length[i];
        
        for(long long i=1; i<(1<<n); i++)
        	for(long long j=0; j<n; j++)
	            if(i & (1<<j)) // i 状态含有 j 
	                for(int k=0; k<n; k++)
	                {
	                    if(k==j) continue;
	                    if(i & (1<<k)) continue; // i 状态不含有 k ,用状态 dp[i][j] 更新状态 dp[i | (1<<k)][k] 
	                    if(visited[j][k]==length[k])
	                        dp[i | (1<<k)][j]=min(dp[i | (1<<k)][j],dp[i][j]);
	                    else
	                        dp[i | (1<<k)][k]=min(dp[i | (1<<k)][k],dp[i][j]+length[k]-visited[j][k]);
	                }
		        
        long long ans=inf;
        for(long long i=0; i<n; i++)
            ans=min(ans,dp[(1<<n)-1][i]);
        printf("%d\n",ans);
    }
    
    return 0;
}