HDU2457 DNA repair AC自动机

这题真的是我的痛!一晚上调了3小时,全都WA了。。。在网上还找了一大堆数据,算的都一样抓狂抓狂原来,格式错了。。。

#include<stdio.h>
#include<string.h>
#include<queue>
using namespace std;
#define maxn 1010
int n,ch[maxn][4],val[maxn],fail[maxn],d[maxn][maxn],tot,root;
char a[maxn];
int newnode()
{
	memset(ch[tot],0,sizeof(ch[tot]));
	val[tot]=0;
	return tot++;
}
int id(char t)
{
	if(t=='A')return 0;
	else if(t=='C')return 1;
	else  if(t=='G')return 2;
	else return 3;
}
void insert()
{
	int len=strlen(a),cur=root,i;
	for(i=0;i<len;i++)
	{
		if(!ch[cur][id(a[i])])
			ch[cur][id(a[i])]=newnode();
		cur=ch[cur][id(a[i])];
	}
	val[cur]=1;
}
queue<int>q;
void getfail()
{
	fail[root]=root;
	for(int i=0;i<4;i++)
		if(ch[root][i])
		{
			fail[ch[root][i]]=root;
			q.push(ch[root][i]);
		}
	while(!q.empty())
	{
		int cur=q.front();
		q.pop();
		for(int i=0;i<4;i++)
			if(!ch[cur][i])
				ch[cur][i]=ch[fail[cur]][i];
			else
			{
				fail[ch[cur][i]]=ch[fail[cur]][i];
				val[ch[cur][i]]|=val[ch[fail[cur]][i]];
				q.push(ch[cur][i]);
			}
	}
}
int min(int a,int b)
{
	return a<b?a:b;
}

int main()
{
	int ex=1;
	while(~scanf("%d",&n)&&n)
	{
		int i,k,j;
		tot=0;root=newnode();
		for(i=0;i<n;i++)
		{
			scanf("%s",a);
			insert();
		}
		getfail();
		scanf("%s",a);
		memset(d,0x3f,sizeof(d));
		d[0][0]=0;
		int len=strlen(a);
		for(i=0;i<len;i++)
			for(j=0;j<tot;j++)
				if(d[i][j]!=0x3f3f3f3f)
					for(k=0;k<4;k++)
					{
						if(val[ch[j][k]])continue;
						d[i+1][ch[j][k]]=min(d[i+1][ch[j][k]],d[i][j]+(k==id(a[i])?0:1));
					}
		int ans=0x3f3f3f3f;
		for(j=0;j<tot;j++)
			if(ans>d[len][j])
				ans=d[len][j];
		printf("Case %d: ",ex++);
		if(ans==0x3f3f3f3f)
			printf("-1\n");
		else
			printf("%d\n",ans);
	}
	return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值