题意:
给出一个图,每条边都有权值w和c,w为边长,c为使边长-1的钱,你拥有一定数目的钱,要你求该图的最小生成树
题解:
参考了cf官方题解的写法。
(首先我们要知道,如果我们知道哪几条边是被选中,那么肯定是找一个c最小的边,然后一直减
先是不考虑用钱换距离,单纯用kruskal裸着求一遍mst
然后我们对所有边,考虑将该条边作为要减的那条边加入到mst中,那么肯定会形成一个环,所以我们需要去掉一条边
显然,我们去掉(u-lca-v)中边长最大的边,这里用到了bianry lifting求lca的思想,nlogn时间预处理,logn时间查询
#include<set>
#include<map>
#include<cmath>
#include<stack>
#include<queue>
#include<bitset>
#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define PB push_back
#define MP make_pair
#define ll long long
#define MS(a,b) memset(a,b,sizeof(a))
#define LL (rt<<1)
#define RR (rt<<1|1)
#define lson l,mid,LL
#define rson mid+1,r,RR
#define pii pair<int,int>
#define pll pair<ll,ll>
#define lb(x) (x&(-x))
void In(){freopen("in.in","r",stdin);}
void Out(){freopen("out.out","w",stdout);}
const int N=2e5+10;
const int M=2e5+10;
const int Mbit=1e6+10;
const ll inf=1ll<<50;
const ll mod=1e9+7;
struct node
{
int a,b,c,h,id;
}edge[M];
bool cmp1(node x,node y)
{
return x.h<y.h;
}
bool cmp2(node x,node y)
{
return x.id<y.id;
}
ll sum;
int m,n,num,l;
int f[N],hei[N],V[N],dep[N],up[N][30],dp[N][30],id[N][30];
vector<node>G[N];
int Find(int x)
{
return f[x]==x?x:f[x]=Find(f[x]);
}
void kruskal()
{
sum=num=0;
int cnt=1;
sort(edge,edge+m,cmp1);
for(int i=1;i<=n;i++)f[i]=i,G[i].clear();
for(int i=0;i<m;i++){
if(cnt==n)break;
int a=edge[i].a,b=edge[i].b;
int fa=Find(a),fb=Find(b);
if(fa!=fb){
G[a].PB(edge[i]);G[b].PB(edge[i]);
f[fa]=fb;
sum+=edge[i].h;
V[num++]=edge[i].id;
cnt++;
}
}
sort(edge,edge+m,cmp2);
}
void dfs(int u,int fa)
{
int s=G[u].size(),v;
for(int i=0;i<s;i++){
node x=G[u][i];
int a=x.a,b=x.b,h=x.h;
if(a==u)v=b;else v=a;
if(v==fa)continue;
dep[v]=dep[u]+1;
up[v][0]=u;
dp[v][0]=h;
id[v][0]=x.id;
for(int i=1;i<=l;i++){
up[v][i]=up[up[v][i-1]][i-1];
if(dp[v][i-1]>dp[up[v][i-1]][i-1]){
dp[v][i]=dp[v][i-1];
id[v][i]=id[v][i-1];
}
else{
dp[v][i]=dp[up[v][i-1]][i-1];
id[v][i]=id[up[v][i-1]][i-1];
}
}
dfs(v,u);
}
}
pii query(int a,int b)
{
if(dep[a]<dep[b])swap(a,b);
int d=dep[a]-dep[b];
int ans=0,ans1;
for(int i=0;d;i++,d/=2)if(d&1){
if(dp[a][i]>ans)ans=dp[a][i],ans1=id[a][i];
a=up[a][i];
}
if(a==b)return MP(ans,ans1);
for(int i=l;i>=0;i--){
if(up[a][i]!=up[b][i]){
if(dp[a][i]>ans)ans=dp[a][i],ans1=id[a][i];
if(dp[b][i]>ans)ans=dp[b][i],ans1=id[b][i];
a=up[a][i];
b=up[b][i];
}
}
if(dp[a][0]>ans)ans=dp[a][0],ans1=id[a][0];
if(dp[b][0]>ans)ans=dp[b][0],ans1=id[b][0];
return MP(ans,ans1);
}
int main()
{
while(~scanf("%d%d",&n,&m)){
for(int i=0;i<m;i++)scanf("%d",&edge[i].h);
for(int i=0;i<m;i++)scanf("%d",&edge[i].c);
for(int i=0;i<m;i++)scanf("%d%d",&edge[i].a,&edge[i].b),edge[i].id=i;
kruskal();
l=0;while((1<<l)<=n)++l;
dfs(1,-1);
ll res=inf;int choose1=-1,choose2=-1;
int have;
scanf("%d",&have);
for(int i=0;i<m;i++){
int a=edge[i].a,b=edge[i].b,c=edge[i].c,h=edge[i].h;
pii q=query(a,b);
int x=q.first,y=q.second;
int t=have/c;
ll tmp=sum-x-t+h;
if(tmp<res)res=tmp,choose1=i,choose2=y;
}
edge[choose1].h-=have/edge[choose1].c;
printf("%lld\n",res);
printf("%d %d\n",choose1+1,edge[choose1].h);
for(int i=0;i<num;i++)
if(V[i]!=choose2)printf("%d %d\n",V[i]+1,edge[V[i]].h);
}
return 0;
}