HDU 2457 DNA repair (AC自动机+DP)

题意:给N个串,一个大串,要求在最小的改变代价下,得到一个不含上述n个串的大串。

思路:dp,f[i][j]代表大串中第i位,AC自动机上第j位的最小代价。

 #include<algorithm>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<queue>
std::queue<int>Q;
int ch[][],end[],f[][],sz,root,fail[],n,Tcase;
char s[];
int newnode(){
for (int i=;i<=;i++)
ch[sz][i]=-;
end[sz]=;
sz++;
return sz-;
}
void clear(){
sz=;root=newnode();
}
int num(char c){
if (c=='A') return ;
else
if (c=='T') return ;
else
if (c=='G') return ;
else
if (c=='C') return ;
return ;
}
void insert(){
scanf("%s",s);
int len=strlen(s);
int now=root;
for (int i=;i<len;i++){
if (ch[now][num(s[i])]==-) ch[now][num(s[i])]=newnode();
now=ch[now][num(s[i])];
}
end[now]=;
}
void build(){
int now=root;
std::queue<int>Q;
fail[root]=root;
for (int i=;i<=;i++)
if (ch[root][i]==-) ch[root][i]=root;
else fail[ch[root][i]]=root,Q.push(ch[root][i]);
while (!Q.empty()){
int now=Q.front();Q.pop();
if (end[fail[now]]==) end[now]|=;
for (int i=;i<=;i++){
if (ch[now][i]==-) ch[now][i]=ch[fail[now]][i];
else fail[ch[now][i]]=ch[fail[now]][i],Q.push(ch[now][i]);
}
}
}
void dp(){
scanf("%s",s);
int len=strlen(s);
for (int i=;i<=len;i++)
for (int j=;j<=sz;j++)
f[i][j]=;
f[][root]=;
for (int i=;i<len;i++)
for (int j=;j<sz;j++)
if (f[i][j]<){
for (int k=;k<=;k++){
int now=ch[j][k];
if (end[now]) continue;
int tmp;
if (k==num(s[i])) tmp=f[i][j];
else tmp=f[i][j]+;
f[i+][now]=std::min(f[i+][now],tmp);
}
}
int ans=;
for (int i=;i<sz;i++)
ans=std::min(ans,f[len][i]);
if (ans==) ans=-;
printf("%d\n",ans);
}
int main(){
while (scanf("%d",&n)!=EOF){
if (n==) return ;
Tcase++;
printf("Case %d: ",Tcase);
clear();
for (int i=;i<=n;i++)
insert();
build();
dp();
}
}
上一篇:.net core 的图片处理及二维码的生成及解析


下一篇:SQL<>0查询不到NUll的值