社区讨论

90分求助(帮我调出来送10元红包)

P4180[BJWC2010] 严格次小生成树参与者 13已保存回复 12

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
12 条
当前快照
1 份
快照标识符
@mi6ze2m4
此快照首次捕获于
2025/11/20 13:18
4 个月前
此快照最后确认于
2025/11/20 15:45
4 个月前
查看原帖
#pragma GCC optimize(2) #include<bits/stdc++.h> #define N 155555 #define M 355555 #define inf 0x7fffffff #define ll long long using namespace std; int n,m,tot,cnt,mn=inf; ll ans; int f[N],head[N],deep[N],fa[N][19],d1[N][19],d2[N][19]; struct data{int x,y,v;bool sel;}a[M]; struct edge{int to,next,v;}e[N*2]; bool cmp(data a,data b) {return a.v<b.v;} void ins(int u,int v,int w) {e[++cnt].to=v;e[cnt].next=head[u];e[cnt].v=w;head[u]=cnt;} void insert(int u,int v,int w) {ins(u,v,w);ins(v,u,w);} int find(int x){return x==f[x]?x:find(f[x]);} void dfs(int x,int f) { for(int i=1;i<=18;i++) { if(deep[x]<(1<<i))break; fa[x][i]=fa[fa[x][i-1]][i-1]; d1[x][i]=max(d1[x][i-1],d1[fa[x][i-1]][i-1]); if(d1[x][i-1]==d1[fa[x][i-1]][i-1]) d2[x][i]=max(d2[x][i-1],d2[fa[x][i-1]][i-1]); else { d2[x][i]=min(d1[x][i-1],d1[fa[x][i-1]][i-1]); d2[x][i]=max(d2[x][i-1],d2[x][i]); d2[x][i]=max(d2[x][i],d2[fa[x][i-1]][i-1]); } } for(int i=head[x];i;i=e[i].next) if(e[i].to!=f) { fa[e[i].to][0]=x; d1[e[i].to][0]=e[i].v; deep[e[i].to]=deep[x]+1; dfs(e[i].to,x); } } int lca(int x,int y) { if(deep[x]<deep[y])swap(x,y); int t=deep[x]-deep[y]; for(int i=0;i<=18;i++) if((1<<i)&t)x=fa[x][i]; for(int i=18;i>=0;i--) { if(fa[x][i]!=fa[y][i]) {x=fa[x][i];y=fa[y][i];} } if(x==y)return x; return fa[x][0]; } void cal(int x,int f,int v) { int mx1=0,mx2=0; int t=deep[x]-deep[f]; for(int i=0;i<=18;i++) { if(t&(1<<i)) { if(d1[x][i]>mx1) { mx2=mx1; mx1=d1[x][i]; } mx2=max(mx2,d2[x][i]); x=fa[x][i]; } } if(mx1!=v)mn=min(mn,v-mx1); else mn=min(mn,v-mx2); } void solve(int t,int v) { int x=a[t].x,y=a[t].y,f=lca(x,y); cal(x,f,v);cal(y,f,v); } int main() { scanf("%d%d",&n,&m); for(int i=1;i<=n;i++) f[i]=i; for(int i=1;i<=m;i++) scanf("%d%d%d",&a[i].x,&a[i].y,&a[i].v); sort(a+1,a+m+1,cmp); for(int i=1;i<=m;i++) { int p=find(a[i].x),q=find(a[i].y); if(p!=q) { f[p]=q; ans+=a[i].v; a[i].sel=1; insert(a[i].x,a[i].y,a[i].v); tot++; if(tot==n-1)break; } } dfs(1,0); for(int i=1;i<=m;i++) if(!a[i].sel)solve(i,a[i].v); printf("%lld",ans+mn); return 0; }

回复

12 条回复,欢迎继续交流。

正在加载回复...