HDU2825 Wireless Password【AC自动机】

题意:求出题目所给的m个串中至少出现k个串且长度为n的字符串的个数。。

AC自动机+状压DP。。。感觉不懂状压DP还是弄懂状压DP比较好

AC自动机的核心其实就是状态转移。。

dp[i][j][k],i表示目前的长度,j表示AC自动机的状态(转移到了哪个节点),k表示已经出现了哪几个串的状态。。AC自动机在转移的时候要根据fail指针适当的合并状态,因为转移到了这个节点,这个节点的fail也应该转移到了,比如she,he,根据s->h->e,到了e之后,其实he也走过了,所以状态要合并,而s-h->e这个e的fail指针会指向h->e。。所以根据fail指针来合并就好。。然后DP感觉就是普通的状压DP。。

最后把出现哪几个串的状态中那些出现的串的个数>=k的加起来就好

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int MOD=20090717;
int size;
struct Node
{
	Node *fail;
	Node *next[26];
	int id,val;
	Node()
	{
		id=size++;
		val=0;
		for(int i=0;i<26;i++)
			next[i]=NULL;
	}
};
Node *T[110];
void insert(char *str,Node *root,int val)
{
	int i=0,index;
	Node *p=root;
	while(str[i])
	{
		index=str[i]-'a';
		if(p->next[index]==NULL)
		{
			p->next[index]=new Node();
			T[size-1]=p->next[index];
		}
		p=p->next[index];
		i++;
	}
	p->val|=(1<<val);
}
queue<Node *> q;
void build_ac_automation(Node *root)
{
	root->fail=NULL;
	q.push(root);
	while(!q.empty())
	{
		Node *temp=q.front();
		q.pop();
		if(temp!=root)
			temp->val|=temp->fail->val;
		for(int i=0;i<26;i++)
		{
			if(temp->next[i]==NULL)
			{
				if(temp==root)
					temp->next[i]=root;
				else
					temp->next[i]=temp->fail->next[i];
			}
			else
			{
				Node *p=temp->fail;
				while(p!=NULL)
				{
					if(p->next[i]!=NULL)
					{
						temp->next[i]->fail=p->next[i];
						break;
					}
					p=p->fail;
				}
				if(p==NULL)
					temp->next[i]->fail=root;
				q.push(temp->next[i]);
			}
		}
	}
}
int dp[26][110][1<<10];
int n,m,k;
int num[1100];
int solve()
{
	int i,j,p;
	memset(dp,0,sizeof(dp));
	dp[0][0][0]=1;
	for(i=0;i<n;i++)
	{
		for(j=0;j<size;j++)
		{
			for(p=0;p<(1<<m);p++)
			{
				if(dp[i][j][p]==0)
					continue;
				for(int x=0;x<26;x++)
				{
					Node *jj=T[j]->next[x];
					dp[i+1][jj->id][p|(jj->val)]+=dp[i][j][p];
					dp[i+1][jj->id][p|(jj->val)]%=MOD;
				}
			}
		}
	}
	int ans=0;
	for(p=0;p<(1<<m);p++)
	{
		if(num[p]<k)
			continue;
		for(i=0;i<size;i++)
			ans=(ans+dp[n][i][p])%MOD;
	}
	return ans;
}
char str[15];
int main()
{
	for(int i=0;i<1024;i++)
	{
		int temp=i;
		int cnt=0;
		while(temp)
		{
			cnt+=temp%2;
			temp/=2;
		}
		num[i]=cnt;
	}
	while(scanf("%d%d%d",&n,&m,&k)==3)
	{
		if(n==0&&m==0&&k==0)
			break;
		size=0;
		Node *root=new Node();
		T[size-1]=root;
		for(int i=0;i<m;i++)
		{
			scanf("%s",str);
			insert(str,root,i);
		}
		build_ac_automation(root);
		printf("%d\n",solve());
	}
	return 0;
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值