hdu 2457(ac自动机+dp)

时间:2023-03-09 06:51:01
hdu 2457(ac自动机+dp)

题意:容易理解...

分析:这是一道比较简单的ac自动机+dp的题了,直接上代码。

代码实现:

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
#include<iostream>
#include<queue>
using namespace std;
struct node{
int next[];
int fail;
int flag;
void init()
{
memset(next,,sizeof(next));
fail=;
flag=;
}
}a[];
char keyword[];
char S[];
int tot,len;
int dp[][]; int Min(int a,int b)
{
return a<b?a:b;
} void chushihua()
{
int i,j;
tot=;
a[].init();
for(i=;i<;i++)
for(j=;j<;j++)
dp[i][j]=;
dp[][]=;
} int hash(char x)
{
if(x=='A')
return ;
else if(x=='C')
return ;
else if(x=='G')
return ;
else
return ;
} void insert(char *str)
{
int index,p=;
for(;*str!='\0';str++)
{
index=hash(*str);
if(a[p].next[index]==)
{
a[++tot].init();
a[p].next[index]=tot;
}
p=a[p].next[index];
}
a[p].flag=;
} void build_fail()
{
queue<int>Q;
int p,cur,son,i;
Q.push();
while(!Q.empty())
{
p=Q.front();
Q.pop();
for(i=;i<;i++)
{
if(a[p].next[i]!=)
{
son=a[p].next[i];
cur=a[p].fail;
if(p==)
a[son].fail=;
else
{
while(cur&&a[cur].next[i]==)
cur=a[cur].fail;
a[son].fail=a[cur].next[i];
}
if(a[a[son].fail].flag==)
a[son].flag=;
if(a[son].flag==)
Q.push(son);
}
else
a[p].next[i]=a[a[p].fail].next[i];
}
}
} void solve(int t)
{
int i,j,k,son,res=;
char x;
for(i=;i<=len;i++)
{
for(j=;j<=tot;j++)
{
if(dp[i-][j]==||a[j].flag==)
continue;
for(k=;k<;k++)
{
son=a[j].next[k];
if(a[son].flag==)
continue;
if(k==)
x='A';
else if(k==)
x='C';
else if(k==)
x='G';
else
x='T';
if(x==S[i-])
{
if(dp[i][son]==)
dp[i][son]=dp[i-][j];
else
dp[i][son]=Min(dp[i][son],dp[i-][j]);
}
else
{
if(dp[i][son]==)
dp[i][son]=dp[i-][j]+;
else
dp[i][son]=Min(dp[i][son],dp[i-][j]+);
}
}
}
}
for(i=;i<=tot;i++)
if(dp[len][i]<res)
res=dp[len][i];
printf("Case %d: ",t);
if(res==)
printf("%d\n",-);
else
printf("%d\n",res);
} int main()
{
int n,t=;
while(scanf("%d",&n)!=EOF&&n)
{
t++;
chushihua();
getchar();
while(n--)
{
scanf("%s",keyword);
insert(keyword);
}
build_fail();
scanf("%s",S);
len=strlen(S);
solve(t);
}
return ;
}