社区讨论
49求条
P4180[BJWC2010] 严格次小生成树参与者 1已保存回复 0
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 0 条
- 当前快照
- 1 份
- 快照标识符
- @mlj62q5z
- 此快照首次捕获于
- 2026/02/12 15:58 3 周前
- 此快照最后确认于
- 2026/02/12 16:11 3 周前
CPP
#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10;
const int M=3e5+10;
int n,m,f[N][20],d[N],d1[N][20],mst,d2[N][20],p[N];
struct node{
int u,v,w,flag;
}a[M];
vector<pair<int,int>> g[N];
bool cmp(struct node a,struct node b){
return a.w<b.w;
}
int find(int x){
if(p[x]!=x) p[x]=find(p[x]);
return p[x];
}
void kkk(){
sort(a+1,a+1+m,cmp);
mst=0;
int cnt=0;
for(int i=1;i<=n;i++){
p[i]=i;
}
for(int i=1;i<=m;i++){
int u=a[i].u,v=a[i].v,w=a[i].w;
u=find(u),v=find(v);
if(u!=v){
g[u].push_back({v,w});
g[v].push_back({u,w});
p[u]=v;
cnt++;
mst+=w;
a[i].flag=1;
if(cnt==n-1) break;
}
}
}
void dfs(int u,int fa,int dep){
d[u]=dep,f[u][0]=fa;
for(int i=1;i<20;i++){
f[u][i]=f[f[u][i-1]][i-1];
d1[u][i]=max(d1[u][i-1],d1[f[u][i-1]][i-1]);
if(d1[u][i-1]==d1[f[u][i-1]][i-1])
d2[u][i]=max(d2[u][i-1],d2[f[u][i-1]][i-1]);
else{
d2[u][i]=min(d1[u][i-1],d1[f[u][i-1]][i-1]);
d2[u][i]=max(d2[u][i],max(d2[u][i-1],d2[f[u][i-1]][i-1]));
}
}
for(auto t:g[u]){
int v=t.first;
int w=t.second;
if(v==fa) continue;
d1[v][0]=w,d2[v][0]=0;
dfs(v,u,dep+1);
}
}
int lca(int x,int y){
if(d[x]<d[y]) swap(x,y);
for(int i=19;i>=0;i--){
if(d[f[x][i]]>=d[y]){
x=f[x][i];
}
}
if(x==y) return x;
for(int i=19;i>=0;i--){
if(f[x][i]!=f[y][i]){
x=f[x][i],y=f[y][i];
}
}
return f[x][0];
}
void update(int u,int t,int &m1,int &m2){
for(int i=19;i>=0;i--){
if(d[f[u][i]]>=d[t]){
if(d1[u][i]>m1){
m2=max(m1,d2[u][i]);
m1=d1[u][i];
}
else if(d1[u][i]<m1&&d1[u][i]>m2){
m2=d1[u][i];
}
u=f[u][i];
}
}
}
int main(){
cin>>n>>m;
for(int i=1;i<=m;i++){
cin>>a[i].u>>a[i].v>>a[i].w;
}
kkk();
dfs(1,0,1);
int ans=2e9;
for(int i=1;i<=m;i++){
if(!a[i].flag){
int u=a[i].u,v=a[i].v,w=a[i].w;
int t=lca(u,v);
int m1=0,m2=0;
update(u,t,m1,m2);
update(v,t,m1,m2);
if(w==m1) ans=min(ans,mst+w-m2);
else ans=min(ans,mst+w-m1);
}
}
cout<<ans<<"\n";
return 0;
}
回复
共 0 条回复,欢迎继续交流。
正在加载回复...