fjwc2019 D1T3 不同的缩写（dinic+trie+dfs）

#180. 「2019冬令营提高组」不同的缩写

乍看之下没有什么好的方法鸭.......于是考虑暴力。

长度？二分似乎可行。

于是我们二分最长子串的长度（设为$len$），蓝后暴力查找。

先在每个串内练好后继边建图

for(int i=;i<=n;++i){
        int len=strlen(a[i]+);
        for(int j=;j<;++j) To[i][len][j]=-;
        for(int j=len-;j>=;--j){
            memcpy(To[i][j],To[i][j+],sizeof(To[i][j+]));
            To[i][j][a[i][j+]-'a']=j+;
        }
    }

每次用dfs查找一个串中长度不超过$len$的子串个数。

一个重要的剪枝：当长度不超过$len$的子串个数已经超过n个时，显然我们可以停止查找。因为显然无论怎么配，这个串都有解。

现在，我们筛掉了一定有解的串，那么对于剩下的串我们怎么判断答案为$len$时是否有解？

考虑暴力建出一棵trie树，保存剩下串的信息。

接下来我们实现匹配。

首先我们先新建总起点，终点$S,T$；以及$n-w$个点表示对应的子串（前面剪枝剪掉的子串（设为$w$个）就不用建点辣）

（当然你可以方便地直接开n个点）

我们遍历子串$k$，当访问到树上的某个节点$p$时，从$k$向$p$连一条流量为1的边（没错！我们等下要跑网络流），表示一种匹配。

所有串遍历完后，trie上的每个点都向$T$连一条流量为1的边。

最后，$S$向$n-w$个代表子串的点连一条流量为1的边。

于是我们就可以从$S$到$T$愉快地跑一遍dinic辣

---

找到了最短的长度，现在考虑输出其中一种方案。

直接上dfs遍历每一个串，如果找到一个还未打上结束标记的节点，就把这个节点的结束标记打上该串的编号，然后跳出dfs，遍历下个串。

最后开个字符栈，再跑遍dfs。

蓝后就是艰难(for me)的打code了

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#define stop system("pause")
using namespace std;
inline int min(int a,int b){return a<b?a:b;}
const int inf=;
#define N 305
#define M 500005
char a[N][N]; bool vis[M];
int n,ans=-,To[N][N][],tot,is[N];
int d[M],cur[M],S,T;
int tri,p[M][],fr[M];
int cnt,hd[M],nxt[M],ed[M],poi[M],cap[M];
queue <int> h;
string b,name[N];
void clear(int x){for(int i=;i<;++i)p[x][i]=;}
void add(int x,int y,int v){
    nxt[ed[x]]=++cnt; hd[x]=hd[x]?hd[x]:cnt;
    ed[x]=cnt; poi[cnt]=y; cap[cnt]=v;
}
void ins(int x,int y,int v){add(x,y,v);add(y,x,);}
#define to poi[i]
bool bfs(){
    memset(vis,,sizeof(vis));
    memset(d,-,sizeof(d));
    h.push(S); vis[S]=; d[S]=;
    while(!h.empty()){
        int x=h.front(); h.pop();
        for(int i=hd[x];i;i=nxt[i]){
            if(!vis[to]&&cap[i]>){
                vis[to]=;
                d[to]=d[x]+;
                h.push(to);
            }
        }
    }return vis[T];
}
int dfs(int x,int a){//dinic
    if(x==T||a==) return a;
    int F=,f;
    for(int &i=cur[x];i&&a>;i=nxt[i])
        if(d[to]==d[x]+&&(f=dfs(to,min(cap[i],a)))>)
            cap[i]-=f,a-=f,cap[i^]+=f,F+=f;
    return F;
}
int dinic(){
    int re=;
    while(bfs()){
        for(int i=;i<=tri+;++i) cur[i]=hd[i];
        re+=dfs(S,inf);
    }
    return re;
}
void dfs1(int d,int id,int o,int dl){
    if(d) ++tot;
    if(d==dl||tot>=n) return;
    for(int i=;i<&&tot<n;++i)
        if(To[id][o][i]!=-)
            dfs1(d+,id,To[id][o][i],dl);
}
void dfs2(int d,int id,int o,int dl,int u){
    if(d) ins(id,u,);
    if(d==dl) return;
    for(int i=;i<;++i)
        if(To[id][o][i]!=-){
            if(!p[u][i]) p[u][i]=++tri,clear(tri);
            dfs2(d+,id,To[id][o][i],dl,p[u][i]);
        }
}
bool chk(int lim){
    memset(ed,,sizeof(ed));
    memset(hd,,sizeof(hd));
    memset(nxt,,sizeof(nxt));
    cnt=; tri=n+; clear(n+);
    for(int i=;i<=n;++i){
        tot=; is[i]=;
        dfs1(,i,,lim);
        if(tot<n) is[i]=,dfs2(,i,,lim,n+);//剪枝
    }S=tri+;T=tri+;
    int tflow=;
    for(int i=;i<=n;++i)
        if(is[i]) ++tflow,ins(S,i,);
    for(int i=n+;i<=tri;++i) ins(i,T,);
    return tflow==dinic();
}
int dfs3(int d,int id,int o,int dl,int u){
    if(d&&!fr[u]){fr[u]=id; return ;}
    if(d==dl) return ;
    for(int i=;i<;++i)
        if(To[id][o][i]!=-){
            if(!p[u][i]) p[u][i]=++tri,clear(tri);
            if(dfs3(d+,id,To[id][o][i],dl,p[u][i]))
                return ;
        }
    return ;
}
void dfs4(int u){
    if(fr[u]) name[fr[u]]=b;
    for(int i=;i<;++i)
        if(p[u][i]){
            b+=(i+'a'); dfs4(p[u][i]);
            b.erase(b.size()-);
        }
}
int main(){
    freopen("diff.in","r",stdin);
    freopen("diff.out","w",stdout);
    scanf("%d",&n);
    for(int i=;i<=n;++i) scanf("%s",a[i]+);
    for(int i=;i<=n;++i){
        int len=strlen(a[i]+);
        for(int j=;j<;++j) To[i][len][j]=-;
        for(int j=len-;j>=;--j){
            memcpy(To[i][j],To[i][j+],sizeof(To[i][j+]));
            To[i][j][a[i][j+]-'a']=j+;
        }
    }
    int l=,r=n;
    while(l<r){
        int mid=(l+r)/;
        if(chk(mid)) r=mid;
        else l=mid+;
    }if(!chk(l)){printf("-1");return ;}
    ans=l;
    for(int i=;i<=n;++i)
        if(is[i])
            for(int j=hd[i];j;j=nxt[j])
                if(cap[j]==)
                    fr[poi[j]]=i;
    for(int i=;i<=n;++i)
        if(!is[i]) dfs3(,i,,ans,n+);
    b=""; dfs4(n+);
    printf("%d\n",ans);
    for(int i=;i<=n;++i) cout<<name[i]<<endl;
    return ;
}