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

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=4534

题意：有三种串，（1）第一种串必须不存在，（2）第二种串必须存在，（3）第三种有的话会有相应的价值得分。现在给出一个大串，删掉其中最少的字母使得（1）（2）满足时（3）最大？

思路：f1[i][j][k]、f2[i][j][k]表示长度为i到达节点j必须存在的串的状态为k的最小删掉的字母数和最多的得分。

struct node
{
    int next[26],fail,val,x,y;

    void init()
    {
        clr(next,0);
        fail=0;
        x=y=val=0;
    }
};

node a[N];
int e;


int tot;

void insert(char s[],int val)
{
    int i,k,p=0;
    for(i=0;s[i];i++)
    {
        k=s[i]-'a';
        if(a[p].next[k]==0)
        {
            a[e].init();
            a[p].next[k]=e++;
        }
        p=a[p].next[k];
    }
    if(val==999) a[p].x|=1<<tot,tot++;
    else if(val==-999) a[p].y=1;
    else a[p].val+=val;
}


queue<int> Q;

void build()
{
    int i,j,k,p,q;
    FOR0(i,26) if(a[0].next[i]) Q.push(a[0].next[i]);
    while(!Q.empty())
    {
        k=Q.front();
        Q.pop();
        for(i=0;i<26;i++)
        {
            if(a[k].next[i])
            {
                p=a[k].next[i];
                q=a[k].fail;
                Q.push(p);
                a[p].fail=a[q].next[i];
                a[p].val+=a[a[p].fail].val;
                a[p].x|=a[a[p].fail].x;
                a[p].y|=a[a[p].fail].y;
            }
            else
            {
                q=a[k].fail;
                a[k].next[i]=a[q].next[i];
            }
        }
    }
}


int f1[2][1600][260],f2[2][1600][260];
int n;
char s[N];


void init(int x)
{
    int i,j,k;
    FOR0(i,e) FOR0(j,1<<tot)
    {
        f1[x][i][j]=INF;
        f2[x][i][j]=-INF;
    }
}




void deal()
{
    int pre=0,cur=1,i,j,k,t,r;
    int L=strlen(s),x=(1<<tot);
    init(pre);
    f1[pre][0][0]=0;
    f2[pre][0][0]=0;
    FOR0(t,L)
    {
        init(cur);
        FOR0(i,e) FOR0(j,x) if(f1[pre][i][j]!=INF)
        {
            if(f1[pre][i][j]+1<f1[cur][i][j]||f1[pre][i][j]+1==f1[cur][i][j]&&f2[pre][i][j]>f2[cur][i][j])
            {
                f1[cur][i][j]=f1[pre][i][j]+1;
                f2[cur][i][j]=f2[pre][i][j];
            }
            k=s[t]-'a';
            r=a[i].next[k];
            if(a[r].y) continue;
            if(f1[pre][i][j]<f1[cur][r][j|a[r].x]||f1[pre][i][j]==f1[cur][r][j|a[r].x]&&f2[pre][i][j]+a[r].val>f2[cur][r][j|a[r].x])
            {
                f1[cur][r][j|a[r].x]=f1[pre][i][j];
                f2[cur][r][j|a[r].x]=f2[pre][i][j]+a[r].val;
            }
        }
        pre^=1;
        cur^=1;
    }
    int p=INF,q;
    FOR0(i,e) if(f1[pre][i][x-1]<p||f1[pre][i][x-1]==p&&f2[pre][i][x-1]>q)
    {
        p=f1[pre][i][x-1];
        q=f2[pre][i][x-1];
    }
    if(p==INF) puts("Banned");
    else PR(p,q);
}

int main()
{
    int C,num;
    RD(C);
    while(C--)
    {
        RD(n);
        a[0].init(); e=1; tot=0;
        int i,val;
        FOR0(i,n) RD(s),RD(val),insert(s,val);
        build();
        RD(s);
        printf("Case %d: ",++num);
        deal();
    }
    return 0;
}