HDU - 2457 DNA repair(AC自动机+dp)

题目链接：点击查看

题目大意：给出 n 个匹配串，再给出一个模式串，问最少修改模式串中多少个字母可以使得模式串中不含有任意一个匹配串

题目分析：因为又是模式串与匹配串的题目，虽然与一般意义上的匹配不太一样，但鉴于数据都非常小，所以可以视为动态规划，先用AC自动机跑出每个可行状态，再在trie树上dp就好了，dp[ i ][ j ]代表到第 i 个字母为止，状态为 j 时最小需要改动的字符个数，然后就是很简单的dp了

代码：
 

#include<iostream>
#include<cstdio>
#include<string>
#include<ctime>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<stack>
#include<climits>
#include<queue>
#include<map>
#include<set>
#include<sstream>
using namespace std;
   
typedef long long LL;
  
typedef unsigned long long ull;
   
const int inf=0x3f3f3f3f;
   
const int N=1e5+100;
 
char s[1010];
 
int fail[1010],trie[1010][26],mp[150],cnt;

bool vis[1010];

int dp[1100][1100];//dp[i][j]:长度为 i 的字符串中状态为 j 时满足条件所需要删除的最少次数 
 
void insert_word()
{
	int len=strlen(s);
	int pos=0;
	for(int i=0;i<len;i++)
	{
		int to=mp[s[i]];
		if(!trie[pos][to])
			trie[pos][to]=++cnt;
		pos=trie[pos][to];
	}
	vis[pos]=true;
}
 
void getfail()
{
	queue<int>q;
	for(int i=0;i<26;i++)
	{
		if(trie[0][i])
		{
			fail[trie[0][i]]=0;
			q.push(trie[0][i]);
		}
	}
	while(!q.empty())
	{
		int cur=q.front();
		q.pop();
		if(vis[fail[cur]])
			vis[cur]=true;
		for(int i=0;i<26;i++)
		{
			if(trie[cur][i])
			{
				fail[trie[cur][i]]=trie[fail[cur]][i];
				q.push(trie[cur][i]);
			}
			else
				trie[cur][i]=trie[fail[cur]][i];
		}
	}
}
 
void init()
{
	cnt=0;
	memset(vis,false,sizeof(vis));
	memset(trie,0,sizeof(trie));
	mp['A']=0,mp['C']=1,mp['G']=2,mp['T']=3;
}

int main()
{
//#ifndef ONLINE_JUDGE
//  freopen("input.txt","r",stdin);
//    freopen("output.txt","w",stdout);
//#endif
//  ios::sync_with_stdio(false);
	int n;
	int kase=0;
	while(scanf("%d",&n)!=EOF&&n)
	{
		init();
		for(int i=1;i<=n;i++)
		{
			scanf("%s",s);
			insert_word();
		}
		getfail();
		scanf("%s",s+1);
		int len=strlen(s+1);
		for(int i=0;i<=len;i++)
			for(int j=0;j<=cnt;j++)
				dp[i][j]=inf;
		dp[0][0]=0;
		for(int i=0;i<=len;i++)
			for(int j=0;j<=cnt;j++)
				if(dp[i][j]<inf)
					for(int k=0;k<4;k++)
					{
						int nj=trie[j][k];
						if(vis[nj])
							continue;
						if(k==mp[s[i+1]])
							dp[i+1][nj]=min(dp[i+1][nj],dp[i][j]);
						else
							dp[i+1][nj]=min(dp[i+1][nj],dp[i][j]+1);
					}
		int ans=inf;
		for(int i=0;i<=cnt;i++)
			ans=min(ans,dp[len][i]);
		if(ans==inf)
			ans=-1;
		printf("Case %d: %d\n",++kase,ans);
	}
	
	
 
      
      
      
      
      
      
      
      
      
    return 0;
}