HDU 4534 郑厂长系列故事——新闻净化(AC自动机+DP)

转载请注明出处，谢谢http://blog.csdn.net/ACM_cxlove?viewmode=contents by---cxlove

题目：给出一些模式串，其中有一些串必须出现在子串当中，有一些串是不可以出现在子串中。然后还有一些串有一些分值。给出母串，问最少需要删除多少个字母，能够满足条件，然后使得分值尽可能大。
http://acm.hdu.edu.cn/showproblem.php?pid=4534
妥妥的AC自动机+DP啊。
复赛第一场，妥妥的上来做了签到题之后就开始开E题了，然后就没有然后了。。。
AC自动机都能写错，妥妥地WA了几十次啊
每次建fail的时候，都要把val,canot,must进行转移啊，我SB了只转移了val
del[i][j][k],pot[i][j][k]表示前i个字符，当前在自动机中的j结点，必取串的状态为k时的最优情况。

#include<iostream>
#include<cstdio>
#include<map>
#include<cstring>
#include<cmath>
#include<vector>
#include<algorithm>
#include<set>
#include<stack>
#include<string>
#include<ctime>
#include<queue>
#include<cassert>
#define inf 0x11111111
#define maxn 210005
#define eps 1e-8
#define zero(a) fabs(a)<eps
#define Min(a,b) ((a)<(b)?(a):(b))
#define Max(a,b) ((a)>(b)?(a):(b))
#define pb(a) push_back(a)
#define mp(a,b) make_pair(a,b)
#define mem(a,b) memset(a,b,sizeof(a))
#define LL long long
#define MOD 1000000007
#define sqr(a) ((a)*(a))
#define Key_value ch[ch[root][1]][0]
#define test puts("OK");
#define pi acos(-1.0)
#define lowbit(x) ((-(x))&(x))
#pragma comment(linker, "/STACK:1024000000,1024000000")
using namespace std;
const int M=1610;
struct Trie  {
Trie *next[26];
Trie *fail;
int must,canot;
int val,word,kind;
};
Trie *que[M],s[M];
int idx,tot,pre[105];
Trie *NewNode()  {
Trie *tmp=&s[idx];
mem(tmp->next,NULL);
tmp->must=tmp->canot=tmp->val=0;
tmp->fail=NULL;
tmp->kind=idx++;
return tmp;
}
void Insert(Trie *root,char *s,int len,int k,int val)  {
Trie *p=root;
for(int i=0; i<len; i++){
if(p->next[s[i]-'a']==NULL) p->next[s[i]-'a']=NewNode();
p=p->next[s[i]-'a'];
}
if(val==999) {p->must|=1<<tot;tot++;}
else if(val==-999) p->canot|=1;
else p->val+=val;
}
void Bulid_fail(Trie *root)  {
int head=0,tail=0;
que[tail++]=root;
root->fail=NULL;
while(head<tail){
Trie *tmp=que[head++];
for(int i=0; i<26; i++){
if(tmp->next[i]){
if(tmp==root) tmp->next[i]->fail=root;
else{
Trie *p=tmp->fail;
while(p!=NULL){
if(p->next[i]){
tmp->next[i]->fail=p->next[i];
break;
}
p=p->fail;
}
if(p==NULL) tmp->next[i]->fail=root;
}
if(tmp->next[i]->fail->val) tmp->next[i]->val+=tmp->next[i]->fail->val;
if(tmp->next[i]->fail->must) tmp->next[i]->must|=tmp->next[i]->fail->must;
if(tmp->next[i]->fail->canot) tmp->next[i]->canot|=1;
que[tail++]=tmp->next[i];
}
else if(tmp==root) tmp->next[i]=root;
else tmp->next[i]=tmp->fail->next[i];
}
}
}
char str[105];
int del[2][M][1<<8];
int pot[2][M][1<<8];
void slove(){
int l=strlen(str);
mem(del[1],0x11);
mem(pot[1],0x8f);
del[1][0][0]=0;
pot[1][0][0]=0;
for(int i=0;i<l;i++){
mem(del[i&1],0x11);
mem(pot[i&1],0x8f);
for(int j=0;j<idx;j++){
for(int k=0;k<(1<<tot);k++){
if(del[(i+1)&1][j][k]>=inf) continue;
if(del[i&1][j][k]>del[(i+1)&1][j][k]+1){
del[i&1][j][k]=del[(i+1)&1][j][k]+1;
pot[i&1][j][k]=pot[(i+1)&1][j][k];
}
else if(del[i&1][j][k]==del[(i+1)&1][j][k]+1&&pot[i&1][j][k]<pot[(i+1)&1][j][k]){
del[i&1][j][k]=del[(i+1)&1][j][k]+1;
pot[i&1][j][k]=pot[(i+1)&1][j][k];
}
int id=str[i]-'a';
if(s[j].next[id]==NULL) continue;
int cur=s[j].next[id]->kind;
if(s[cur].canot) continue;
int val=s[cur].val;
int curk=k|s[cur].must;
if(del[i&1][cur][curk]>del[(i+1)&1][j][k]){
del[i&1][cur][curk]=del[(i+1)&1][j][k];
pot[i&1][cur][curk]=pot[(i+1)&1][j][k]+val;
}
else if(del[i&1][cur][curk]==del[(i+1)&1][j][k] && pot[i&1][cur][curk]<pot[(i+1)&1][j][k]+val){
del[i&1][cur][curk]=del[(i+1)&1][j][k];
pot[i&1][cur][curk]=pot[(i+1)&1][j][k]+val;
}
}
}
}
int a1=inf,a2;
int full=(1<<tot)-1,num=(l+1)&1;
for(int i=0;i<idx;i++){
if(del[num][i][full]<a1){
a1=del[num][i][full];
a2=pot[num][i][full];
}
else if(del[num][i][full]==a1&&pot[num][i][full]>a2){
a1=del[num][i][full];
a2=pot[num][i][full];
}
}
if(a1>=inf) puts("Banned");
else printf("%d %d\n",a1,a2);
}
int n,val;
int main(){
// freopen("input.txt","r",stdin);
int t,cas=0;
scanf("%d",&t);
while(t--){
idx=0;tot=0;
scanf("%d",&n);
Trie *root=NewNode();
for(int i=0;i<n;i++){
scanf("%s%d",str,&val);
Insert(root,str,strlen(str),i,val);
}
Bulid_fail(root);
scanf("%s",str);
printf("Case %d: ",++cas);
slove();
}
return 0;
}