BJ模拟：cut（Gomory–Hu 树）

传送门

题意：
给一个$n*n$的矩阵$A$，要求构造一个图$G$，使得$G$中$(i,j)$的最小割为$A_{i,j}$。 输出任意一个满足条件的图或判断无解。

题解：
任意一个满足条件的图可以转化为一颗GomoryHu树，我们只需要构造出这样的一颗树即可。
利用GomoryHu树的知识，现在的条件变为了$(i,j)$在树上路径的最小值为$A_{i,j}$ ，这类似与最大生成树，因为每个最大的边都必须先加入。 我们构造出最大生成树后判断无解即可。

利用树链剖分，时间复杂度为$O(n^2 \log n)$。

#include <bits/stdc++.h>
using namespace std;
typedef pair <int,int> pii;

const int RLEN=1<<18|1;
inline char nc() {
    static char ibuf[RLEN],*ib,*ob;
    (ib==ob) && (ob=(ib=ibuf)+fread(ibuf,1,RLEN,stdin));
    return (ib==ob) ? -1 : *ib++;
}
inline int rd() {
    char ch=nc(); int i=0,f=1;
    while(!isdigit(ch)) {if(ch=='-')f=-1; ch=nc();}
    while(isdigit(ch)) {i=(i<<1)+(i<<3)+ch-'0'; ch=nc();}
    return i*f;
}
const int N=3e2+50,INF=0x3f3f3f3f;
struct E {int x,y,c;} e[N*N];
inline bool operator <(const E &a,const E &b) {return a.c>b.c;}
int n,A[N][N],anc[N],val[N],fa[N],tot;
vector <pii> edge[N];
inline int ga(int x) {return (anc[x]==x) ? x : (anc[x]=ga(anc[x]));}
inline void dfs(int x,int f) {
    for(int e=edge[x].size()-1;e>=0;e--) {
        int v=edge[x][e].first; if(v==f) continue;
        val[v]=edge[x][e].second; dfs(v,fa[v]=x);
    }
}
inline int lca(int x,int y) {
    static int vis[N],vs; ++vs;
    while(x) {
        if(vis[x]==vs) return x;
        vis[x]=vs; x=fa[x];
        if(!x) swap(x,y);
    }
}
inline int gm(int i,int j) {
    int l=lca(i,j), mn=INF;
    for(int xx=i;xx!=l;xx=fa[xx]) mn=min(mn,val[xx]);
    for(int yy=j;yy!=l;yy=fa[yy]) mn=min(mn,val[yy]);
    return mn;
}
int main() {
    n=rd();
    for(int i=1;i<=n;i++) anc[i]=i;
    for(int i=1;i<=n;i++) 
        for(int j=1;j<=n;j++) {
            int x=rd(); A[i][j]=x; 
            if(i>j && A[i][j]!=A[j][i]) {puts("-1"); return 0;}
            if(i!=j) e[++tot]=(E){i,j,x};
        }
    sort(e+1,e+tot+1);
    for(int i=1;i<=tot;i++) {
        int x=e[i].x, y=e[i].y, c=e[i].c;
        if(ga(x)==ga(y)) continue;
        anc[ga(x)]=ga(y); 
        edge[x].push_back(pii(y,c));
        edge[y].push_back(pii(x,c));
    } dfs(1,0);
    for(int i=1;i<=n;i++) 
        for(int j=i+1;j<=n;j++) 
            if(i!=j && gm(i,j)!=A[i][j]) {puts("-1"); return 0;}
    printf("%d\n",n-1);
    for(int i=2;i<=n;i++) printf("%d %d %d\n",fa[i],i,val[i]);
}