zoj3545

/*
AC自动机
相当暴力的 解法:
mark[i][j][k]表示长度为i,走到节点j上,目前还有的基因组合为k的状态是否可达,true表示可达,
false表示不可达。
*/
#include <cstdio>
#include <cstring>
#include <map>
using namespace std;
const int N=105*10; //节点个数的最大值
const int S=4; //不同的字符 个数
const int K=10;
struct node{
	node *sons[S], *fail;
	int num;
}nodes[N], *root;
int cnt;  //cnt是树中节点的 个数
node *que[N];
map<char, int> mci;
bool mark[2][N][1<<K];
int vs[1<<K], tv[K];
int n, l;

void init(){
	mci['A']=0;
	mci['G']=1;
	mci['T']=2;
	mci['C']=3;
}
void clear(){
	cnt=0;
	root=NULL;
}
node* newNode(){
	node* ans=&nodes[cnt++];
	memset(ans->sons, 0, S*sizeof(node*));
	ans->fail=NULL;
	ans->num=0;
	return ans;
}
int hash(char ch){ //字符的哈希函数,根据不同的需要而定
	return mci[ch];
}
int id(node*& t){ //计算该节点的id
	return t-nodes;
}
//j表示当前关键字的标号
void insert(node*& root, char* str, int j){
	node* t=root;
	int i, k;
	for(i=0; str[i]; i++){
		if(t->sons[k=hash(str[i])]==NULL){
			t->sons[k] = newNode();
		}
		t=t->sons[k];
	}
	t->num|=(1<<j);
}
void getFail(node*& root){
	int l, r, i;
	node *t;
	l=r=0;
	root->fail=root; //这样可以保证每个节点的fail指针都是非空的
	for(que[r++]=root; l!=r; ){
		t=que[l++];
		for(i=0; i<S; i++){
			if(t->sons[i]){
				que[r++]=t->sons[i];
				if(t==root){
					t->sons[i]->fail=t;
				}else{
					t->sons[i]->fail=t->fail->sons[i];
				}
				t->sons[i]->num|=t->sons[i]->fail->num;
			}else{ //增设虚拟节点
				if(t==root) t->sons[i]=t;
				else t->sons[i]=t->fail->sons[i];
			}
		}
	}
}

bool input(){
	if(scanf("%d%d", &n, &l) == EOF) return false;
	int i;
	char str[105];
	clear();
	root=newNode();
	for(i=0; i<n; i++){
		scanf("%s%d", str, tv+i);
		insert(root, str, i);
	}
	return true;
}
void solve(){
	int i, sum, k, u, v, j, g, s;
	sum=1<<n;
	for(vs[0]=0, i=1; i<sum; i++){
		for(k=0; k<n; k++){
			if(i&(1<<k)){
				vs[i]=vs[i^(1<<k)]+tv[k];
				break;
			}
		}
	}
	getFail(root);
	for(i=0; i<cnt; i++){
		for(k=0; k<sum; k++){
			mark[0][i][k]=false;
		}
	}
	mark[0][0][0]=true;
	for(u=i=0; i<l; i++){
		v=u^1;
		for(j=0; j<cnt; j++){
			for(k=0; k<sum; k++){
				mark[v][j][k]=false;
			}
		}
		for(j=0; j<cnt; j++){
			for(k=0; k<sum; k++){
				if(!mark[u][j][k]) continue;
				for(g=0; g<S; g++){
					s=id(nodes[j].sons[g]);
					mark[v][s][k|nodes[j].sons[g]->num]=true;
				}
			}
		}
		u=v;
	}
	int ans=-1;
	for(j=0; j<cnt; j++){
		for(k=0; k<sum; k++){
			if(!mark[u][j][k]) continue;
			if(ans<vs[k]){
				ans=vs[k];
			}
		}
	}
	if(ans < 0) printf("No Rabbit after 2012!\n");
	else printf("%d\n", ans);
}
int main(){
//	freopen("in.txt", "r", stdin);
	init();
	while(input()) solve();
	return 0;
}
 
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值