BJ模拟:Mr. Panda and Fantastic Beats(广义SAM)

传送门

题解:
对所有串建广义SAM,我们可以知道每个串是否只在1出现。
然后在SAM上贪心即可。

#include <bits/stdc++.h>
using namespace std;
const int N=5e5+50, L=26;
int son[N][L],len[N],fail[N],id[N],mx[N],mn[N],mnlen[N],lst,tot;
inline void init() {
    for(int i=1;i<=tot;i++) {
        memset(son[i],0,sizeof(son[i]));
        len[i]=fail[i]=id[i]=mn[i]=mx[i]=0;
    } tot=lst=1;
}
inline int copy(int &x,int y) {x=++tot; memcpy(son[x],son[y],sizeof(son[x])); fail[x]=fail[y]; }
inline void extend(int c) {
    if(son[lst][c]) {
        int q=son[lst][c];
        if(len[q]==len[lst]+1) lst=q;
        else {
            int nq; copy(nq,q);
            fail[q]=nq; len[nq]=len[lst]+1;
            while(lst&&son[lst][c]==q) son[lst][c]=nq, lst=fail[lst]; 
            lst=nq;
        } 
    } else {
        int p=++tot; len[p]=len[lst]+1; 
        while(lst && !son[lst][c]) son[lst][c]=p, lst=fail[lst];
        if(!lst) fail[p]=1;
        else {
            int q=son[lst][c];
            if(len[q]==len[lst]+1) fail[p]=q;
            else {
                int np; copy(np,q);
                fail[q]=fail[p]=np; len[np]=len[lst]+1;
                while(lst && son[lst][c]==q) son[lst][c]=np, lst=fail[lst];
            }
        }
        lst=p;
    }
}
inline void W(int p) {
    if(id[p]==1) return;
    int ind=26,ml=0x3f3f3f3f;
    for(int i=0;i<L;i++) 
        if(mn[son[p][i]]==1)  
            if(ml>mnlen[son[p][i]]) 
                ind=i, ml=mnlen[son[p][i]];
    if(p==1&&ind==26) printf("Impossible");
    if(ind!=26) putchar('a'+ind), W(son[p][ind]);
}
inline void solve() {
    init(); int n; cin>>n;
    for(int i=1;i<=n;i++) {
        static char ch[N]; lst=1;
        scanf("%s",ch+1); int l=strlen(ch+1);
        for(int j=1;j<=l;j++) extend(ch[j]-'a'),id[lst]=i;
    }
    static int q[N], c[N];
    for(int i=0;i<=tot;i++) c[i]=0;
    for(int i=1;i<=tot;i++) ++c[len[i]];
    for(int i=1;i<=tot;i++) c[i]+=c[i-1];
    for(int i=1;i<=tot;i++) q[c[len[i]]--]=i;
    for(int i=tot;i>=1;i--) {
        int u=q[i]; mn[u]=id[u]; mnlen[u]=1;
        if(fail[u]) id[fail[u]]=max(id[fail[u]],id[u]);
    }
    for(int i=tot;i>=1;i--) {   
        int u=q[i];
        for(int j=0;j<L;j++)
            if(son[u][j]) {
                int v=son[u][j];
                if(mn[v]<mn[u]) mn[u]=mn[v], mnlen[u]=mnlen[v]+1;
                else if(mn[v]==mn[u]) mnlen[u]=min(mnlen[u],mnlen[v]+1);
            }
    }
    static int tl; 
    printf("Case #%d: ",++tl);
    W(1); putchar('\n');    
}
int main() {
    int T; cin>>T; while(T--) solve();  
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值