HDU 3247 AC自动机 + 状态压缩dp

小HH：乳鸽的神题，状态很容易看出来，有50000*1024，很难保持，我用散列表超时了，用bitset刚好可以卡过，不过后来我想，只有尾结点才有效，中间的很多结点完全可以忽略，可以先用最短路吧各个尾结点之间的距离算出来，经过测试，不到50个点，马上就优化到50*1024了，本来9s多过的优化到了100多MS

#include <stdio.h>
#include <math.h>
#include <queue>
#include <string.h>
#include <algorithm>
using namespace std;
const int MAXN=60000;//结点的最大数，50000+1000*10
const int MAXL=(1<<10);
const int MAX=2;
const int MAXM=11;
typedef struct Trie_Node
{
	bool virus;//是否病毒
	int end,fail;
	int next[MAX];
}Trie;
Trie tree[MAXN];
char str[50010];
int dis[MAXN];
int pos[MAXM];
int G[MAXM][MAXM];
int dp[MAXL][MAXM];
int size,cnt;


void insert(char *s,int id)
{
	int p=0;
	int i=0;
	while(s[i]!='\0')
	{
		if(tree[p].next[s[i]-'0']==0)
		{
			tree[++size].end=0;
			tree[size].fail=0;
			tree[size].virus=false;
			for(int j=0;j<MAX;j++)
				tree[size].next[j]=0;
			tree[p].next[s[i]-'0']=size;
		}
		p=tree[p].next[s[i]-'0'];
		i++;
	}
	if(id>=0) tree[p].end=(1<<id);
	else tree[p].virus=true;
}
void BFS() //求出fail
{
	queue<int>q;
	int temp,p;
	q.push(0);
	while(!q.empty())
	{
		temp=q.front();
		q.pop();
		for(int i=0;i<MAX;i++)
		{
			if(tree[temp].next[i]!=0)
			{
				p=tree[temp].next[i];
				q.push(p);
				if(temp!=0)tree[p].fail=tree[tree[temp].fail].next[i];
				tree[p].end|=tree[tree[p].fail].end;
				tree[p].virus|=tree[tree[p].fail].virus;
			}
			else tree[temp].next[i]=tree[tree[temp].fail].next[i];
		}
	}
}
void Path(int k)
{
	queue<int>q;
	q.push(pos[k]);
	memset(dis,-1,sizeof(dis));
	dis[pos[k]]=0;
	int now,p;

	while(!q.empty())
	{
		now=q.front();
		q.pop();
		for(int i=0;i<MAX;i++)
		{
			p=tree[now].next[i];
			if(dis[p]<0&&!tree[p].virus)
			{
				dis[p]=dis[now]+1;
				q.push(p);
			}
		}
	}
	for(int i=0;i<cnt;i++)
		G[k][i]=dis[pos[i]];

}
inline int MIN(int x,int y)
{
	if(x<0||y<0)return x>y?x:y;
	else return x>y?y:x;
}
void DoIt(int n)
{
	memset(dp,-1,sizeof(dp));
	dp[0][0]=0;
	for(int i=0;i<(1<<n);i++)
		for(int j=0;j<cnt;j++)
		{
			if(dp[i][j]<0)continue;
			for(int k=0;k<cnt;k++)
			{
				if(G[j][k]<0)continue;
				int t=i|tree[pos[k]].end;
				dp[t][k]=MIN(dp[t][k],dp[i][j]+G[j][k]);
			}
		}
		int t=(1<<n)-1;
		int ans=-1;
		for(int i=0;i<cnt;i++)
			ans=MIN(ans,dp[t][i]);
		printf("%d\n",ans);
}
int main()
{
	int n,m;
	while(scanf("%d%d",&n,&m))
	{
		if(n==0&&m==0)break;
		tree[0].end=0;
		tree[0].fail=0;
		tree[0].virus=false;
		for(int i=0;i<MAX;i++)tree[0].next[i]=0;
		size=0;
		for(int i=0;i<n;i++)
		{
			scanf("%s",&str);
			insert(str,i);
		}
		for(int i=0;i<m;i++)
		{
			scanf("%s",&str);
			insert(str,-1);
		}
		BFS();
		cnt=1;
		pos[0]=0;
		for(int i=0;i<=size;i++)
			if(tree[i].end)
				pos[cnt++]=i;
		for(int i=0;i<cnt;i++)
			Path(i);
		DoIt(n);
	}
	return 0;
}