【BZOJ】1977 [BeiJing2010组队]次小生成树 Tree kruskal+LCA

题目传送门

这题好像在好久以前就讲过啊……为什么我现在才去做……

其实正解挺好想的，首先做一棵最小生成树，然后枚举一条非树边，显然这条非树边的加入可以使原本树上的唯一路径形成一个环，把这个环上严格第二小的边删掉，就是一棵严格次小生成树了。然后对所有的次小生成树求一个最小权值即可。

至于怎么查找原来的唯一路径上的严格次大值，可以用ST表维护最大值和次大值。

附上AC代码：

#include <cstdio>#include <cctype>#include <algorithm>using namespace std;typedef long long ll;const int M=3e5+10,N=1e5+10;struct note{    int x,y,w;    bool operator < (const note lyf) const {return w<lyf.w;}}a[M];struct side{    int to,w,nt;}s[N<<1];int n,m,fa[N],h[N],num,d[N],f[N][20],mx1[N][20],mx2[N][20],mn;ll ans;bool b[M];inline char nc(void){    static char ch[100010],*p1=ch,*p2=ch;    return p1==p2&&(p2=(p1=ch)+fread(ch,1,100010,stdin),p1==p2)?EOF:*p1++;}inline void read(int &a){    static char c=nc();int f=1;    for (;!isdigit(c);c=nc()) if (c=='-') f=-1;    for (a=0;isdigit(c);a=(a<<3)+(a<<1)+c-'0',c=nc());    return (void)(a*=f);}inline int gf(int x){return x==fa[x]?x:fa[x]=gf(fa[x]);}inline void add(int x,int y,int w){    s[++num]=(side){y,w,h[x]},h[x]=num;    s[++num]=(side){x,w,h[y]},h[y]=num;}inline void so(int x,int fa){    for (int i=1; i<=17; ++i){        if (d[x]<(1<<i)) break;        f[x][i]=f[f[x][i-1]][i-1];        mx1[x][i]=max(mx1[x][i-1],mx1[f[x][i-1]][i-1]);        if (mx1[x][i-1]==mx1[f[x][i-1]][i-1]) mx2[x][i]=max(mx2[x][i-1],mx2[f[x][i-1]][i-1]);        else {            mx2[x][i]=min(mx1[x][i-1],mx1[f[x][i-1]][i-1]);            mx2[x][i]=max(mx2[x][i],max(mx2[x][i-1],mx2[f[x][i-1]][i-1]));        }    }    for (int i=h[x]; i; i=s[i].nt)        if (s[i].to!=fa) d[s[i].to]=d[f[s[i].to][0]=x]+1,mx1[s[i].to][0]=s[i].w,so(s[i].to,x);}inline int lca(int x,int y){    if (d[x]<d[y]) swap(x,y);int t=d[x]-d[y];    for (int i=0; i<=17; ++i) if (t&(1<<i)) x=f[x][i];    for (int i=17; i>=0; --i) if (f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];    if (x==y) return x;    return f[x][0];}inline void calc(int x,int fa,int w){    int d1=0,d2=0,t=d[x]-d[fa];    for (int i=0; i<=17; ++i)        if (t&(1<<i)){            if (mx1[x][i]>d1) d2=d1,d1=mx1[x][i];            d2=max(d2,mx2[x][i]),x=f[x][i];        }    if (d1!=w) mn=min(mn,w-d1); else mn=min(mn,w-d2);    return;}int main(void){    read(n),read(m);    for (int i=1; i<=m; ++i) read(a[i].x),read(a[i].y),read(a[i].w);    sort(a+1,a+1+m);    for (int i=1; i<=n; ++i) fa[i]=i;    for (int i=1; i<=m; ++i){        int fx=gf(a[i].x),fy=gf(a[i].y);        if (fx!=fy) fa[fx]=fy,ans+=a[i].w,b[i]=1,add(a[i].x,a[i].y,a[i].w);    }    so(1,0),mn=2e9;    for (int i=1; i<=m; ++i)        if (!b[i]){            int t=lca(a[i].x,a[i].y);            calc(a[i].x,t,a[i].w),calc(a[i].y,t,a[i].w);        }    return printf("%lld\n",ans+mn),0;}