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;
}