F. Drivers Dissatisfaction+最小生成樹+lca求樹上兩點的最大值
阿新 • • 發佈:2018-07-13
continue tmp rst spa 題解 去掉 題意 drive .com
題目鏈接:F. Drivers Dissatisfaction
題意:n個點,m條邊,每條邊有一個w,代表這條路的不滿意度,每一條邊我們可以花費c來使不滿意讀-1;然後問你有s,找到一棵生成樹是滿意度最小
題解:對於s,我們可以知道花費在c最小的邊上價值最優,我們可以先求一顆最小生成樹,然後枚舉沒有用到的邊,把連接這兩點的最長邊去掉判段能否更新最小值
這裏求樹上兩點的最短路經過的最大值,我們可以有lca,或者樹鏈拋分都可以
#include<bits/stdc++.h> #include<set> #include<cstdio> #include<iomanip> #include<iostream> #include<string> #include<cstring> #include<algorithm> #define pb push_back #define ll long long #define fi first #define se second #define PI 3.14159265 #define ls l,m,rt<<1 #define rs m+1,r,rt<<1|1 #define eps 1e-7 #definepii pair<int,int> typedef unsigned long long ull; const int mod=1e3+5; const ll inf=0x3f3f3f3f3f3f3f; const int maxn=2e5+5; using namespace std; int n,m,w[maxn],c[maxn],f[maxn],dep[maxn],s,head[maxn],cnt; bool vis[maxn]; int fa[20][maxn],mx[20][maxn];//mx[i][x]保存x往上走2^i最長邊的id struct data{int u,v,w,c,id; bool operator<(const data b) { return w<b.w; } }da[maxn]; struct edge { int to,nxt,w,id; }e[maxn<<1]; int find(int x) { if(f[x]==x)return x; return f[x]=find(f[x]); } void link(int x,int y) { f[find(y)]=find(x); } void add_edge(int u,int v,int id) { e[++cnt].to=v;e[cnt].nxt=head[u];head[u]=cnt;e[cnt].id=da[id].id;e[cnt].w=da[id].w; e[++cnt].to=u;e[cnt].nxt=head[v];head[v]=cnt;e[cnt].id=da[id].id;e[cnt].w=da[id].w; } void dfs(int v) { for(int i=1;i<20;i++) { fa[i][v]=fa[i-1][fa[i-1][v]]; if(w[mx[i-1][fa[i-1][v]]]<w[mx[i-1][v]])mx[i][v]=mx[i-1][v]; else mx[i][v]=mx[i-1][fa[i-1][v]]; } for(int i=head[v];i;i=e[i].nxt) { int to=e[i].to; if(to!=fa[0][v]) { dep[to]=dep[v]+1; fa[0][to]=v; mx[0][to]=e[i].id; dfs(to); } } } int lca(int x,int y) { if(dep[x]<dep[y])swap(x,y); int t=dep[x]-dep[y],temp=0,ans=0; for(int i=0;i<20;i++) { if(t>>i&1) { if(w[mx[i][x]]>temp)ans=mx[i][x],temp=w[mx[i][x]]; x=fa[i][x]; } } if(x==y)return ans; for(int i=19;i>=0;i--) { if(fa[i][x]!=fa[i][y]) { if(w[mx[i][x]]>temp)ans=mx[i][x],temp=w[ans]; if(w[mx[i][y]]>temp)ans=mx[i][y],temp=w[ans]; x=fa[i][x]; y=fa[i][y]; } } if(w[mx[0][x]]>temp)ans=mx[0][x],temp=w[ans]; if(w[mx[0][y]]>temp)ans=mx[0][y],temp=w[ans]; return ans; } int main() { ios::sync_with_stdio(false); cin.tie(0);cout.tie(0); cin>>n>>m; for(int i=1;i<=m;i++)cin>>w[i],da[i].w=w[i]; for(int i=1;i<=m;i++)cin>>c[i],da[i].c=c[i]; for(int i=1;i<=m;i++) cin>>da[i].u>>da[i].v,da[i].id=i; cin>>s; sort(da+1,da+1+m); for(int i=0;i<=n;i++)f[i]=i; ll ans=0,num=-1;int cn=n-1; for(int i=1;i<=m;i++) { if(find(da[i].v)!=find(da[i].u)) { vis[da[i].id]=true; link(da[i].v,da[i].u); add_edge(da[i].v,da[i].u,i); if(num==-1)num=da[i].id; else if(c[num]>da[i].c) num=da[i].id; ans+=da[i].w; if(--cn==0)break; } } ll sum=ans; ans=sum-s/c[num]; dfs(1); int del=0; for(int i=1;i<=m;i++) { if(!vis[da[i].id]) { int id=da[i].id; int t=lca(da[i].u,da[i].v); if(t==0)continue; ll tmp=sum-w[t]+w[id]-s/c[id]; if(tmp<ans) { del=t;num=da[i].id; ans=tmp; } } } if(del)vis[del]=false,vis[num]=true; cout<<ans<<endl; w[num]-=s/c[num]; for(int i=1;i<=m;i++) { if(vis[i]) { cout<<i<<" "<<w[i]<<endl; } } return 0; }
F. Drivers Dissatisfaction+最小生成樹+lca求樹上兩點的最大值