fzyjojP2963 -- [校内训练20161227]疫情控制问题

（题干中的废话已经划去）

dp显而易见

收益为负数的可以直接扔掉不管。不要一定更优

子串问题，考虑SAM

建立广义SAM

尝试匹配，匹配到的位置的parent树祖先如果有完整的串，那么可以从这个串转移

考虑祖先不好考虑

不妨考虑i对j(i<j)的贡献，就是子树了

线段树维护dfn序，区间对val取max

SAM学傻了

这个，由于是"完整出现”，所以AC自动机就可以（比广义SAM好写好调）

走一走属于第j个串的路径，对应的fail树的祖先包含i(i<j)的结束位置，那么可以转移

所以，同理，走完i之后，整个fail树子树的dp值和val取max

但是注意一点，根节点标号是0（反正我trie都这样写）

代码：

#include<bits/stdc++.h>
#define il inline
#define int long long
#define reg register int
#define mid ((l+r)>>1)
#define numb (ch^'0')
using namespace std;
typedef long long ll;
il void rd(int &x){
    char ch;bool fl=false;
    while(!isdigit(ch=getchar()))(ch=='-')&&(fl=true);
    for(x=numb;isdigit(ch=getchar());x=x*+numb);
    (fl==true)&&(x=-x);
}
namespace Miracle{
const int N=+;
int n,m;
struct node{
    int nxt,to;
}e[*N];
int hd[N],cnt;
void add(int x,int y){
    e[++cnt].nxt=hd[x];
    e[cnt].to=y;
    hd[x]=cnt;
}
char s[N],b[N];
int len[N],st[N];
ll v[N],f[N];
int dfn[N],dfn2[N];
int df;
struct seg{
    ll mx;
    ll tag;
}t[*N];
void pushup(int x){
    t[x].mx=max(t[x<<].mx,t[x<<|].mx);
}
void pushdown(int x){
    if(t[x].tag<=) return;
    t[x<<].tag=max(t[x<<].tag,t[x].tag);
    t[x<<].mx=max(t[x<<].mx,t[x].tag);

    t[x<<|].tag=max(t[x<<|].tag,t[x].tag);
    t[x<<|].mx=max(t[x<<|].mx,t[x].tag);
    t[x].tag=;
}
void upda(int x,int l,int r,int L,int R,ll c){
    if(L<=l&&r<=R){
        t[x].mx=max(t[x].mx,c);
        t[x].tag=max(t[x].tag,c);
        return ;
    }
    pushdown(x);
    if(L<=mid) upda(x<<,l,mid,L,R,c);
    if(mid<R) upda(x<<|,mid+,r,L,R,c);
    pushup(x);
}
ll query(int x,int l,int r,int p){
    if(l==r){
        return t[x].mx;
    }
    pushdown(x);
    if(p<=mid) return query(x<<,l,mid,p);
    else return query(x<<|,mid+,r,p);
}
void cle(int x,int l,int r){
    t[x].mx=;t[x].tag=;
    if(l==r) return;
    cle(x<<,l,mid);cle(x<<|,mid+,r);
}
struct AC{
    int fail[N],ch[N][];
    int cnt;
    void ins(char *s,int l){
        //int len=strlen(s+1);
        int now=;
        for(reg i=;i<=l;++i){
            int c=s[i]-'a';
            if(!ch[now][c]) ch[now][c]=++cnt;
            now=ch[now][c];
        }
    }
    void build(){
        queue<int>q;
        for(reg i=;i<;++i){
            if(ch[][i]) fail[ch[][i]]=,q.push(ch[][i]);
        }
        while(!q.empty()){
            int x=q.front();q.pop();
            for(reg i=;i<;++i){
                if(ch[x][i]){
                    fail[ch[x][i]]=ch[fail[x]][i];
                    q.push(ch[x][i]);
                }
                else ch[x][i]=ch[fail[x]][i];
            }
        }
        for(reg i=;i<=cnt;++i){
            add(fail[i],i);
        }
    }
    void clear(){
        memset(fail,,sizeof fail);
        cnt=;
        memset(ch,,sizeof ch);
    }
    void dfs(int x){
        dfn[x]=++df;
        for(reg i=hd[x];i;i=e[i].nxt){
            int y=e[i].to;
            dfs(y);
        }
        dfn2[x]=df;
    }
    ll wrk(char *s,int l,int id){
        ll ret=;
        int now=;
        for(reg i=;i<=l;++i){
            int c=s[i]-'a';
            now=ch[now][c];
    //        cout<<"now "<<now<<" i "<<i<<" dfn "<<dfn[now]<<endl;
            ret=max(ret,query(,,df,dfn[now]));
        }
    //    cout<<dfn[now]<<".. "<<dfn2[now]<<endl;
        upda(,,df,dfn[now],dfn2[now],ret+v[id]);
        return ret;
    }
}ac;
void clear(){
    ac.clear();
    if(df)cle(,,df);
    df=;cnt=;
    n=;
    memset(hd,,sizeof hd);

}
int main(){
    //return 0;
    int t;
    rd(t);
    while(t--){
        clear();
        rd(m);
        int tot=;
        for(reg i=;i<=m;++i){
            scanf("%s",b+);
            scanf("%lld",&v[i]);
            st[i]=tot+;
            len[i]=strlen(b+);
            for(reg j=tot+;j<=tot+len[i];++j){
                s[j]=b[j-tot];
            }
            tot=tot+len[i];
            ac.ins(b,len[i]);
        }
    //    cout<<" after ac.ins "<<endl;
        ac.build();
        ac.dfs();//find dfn
        ll ans=;
    //    cout<<" let's wrk "<<endl;
        for(reg i=;i<=m;++i){
    //        cout<<" ii "<<i<<"--------------------"<<v[i]<<" len "<<len[i]<<" st "<<st[i]<<endl;
            //if(v[i]<=0) continue;
            for(reg j=st[i];j<=st[i]+len[i]-;++j){
                b[j-st[i]+]=s[j];
            }
        //    cout<<b+1<<endl;
            f[i]=ac.wrk(b,len[i],i)+v[i];
        //    cout<<i<<" : "<<f[i]<<endl;
            ans=max(ans,f[i]);
        }
        printf("%lld\n",ans);
    }
    return ;
}

}
signed main(){
    freopen("2963.in","r",stdin);
    freopen("2963.out","w",stdout);
    Miracle::main();
    return ;
}

/*
   Author: *Miracle*
   Date: 2019/2/3 11:36:24
*/