CF733F Drivers Dissatisfaction

一、题目

点此看题

二、解法

因为每条边可以降低权值到负数,所以选生成树中 c c c最小的边降低权值是最赚的。

先按 w w w建出最小生成树,设生成树中最小的 c c c m c mc mc,非树边要找 c c c m c mc mc小的边,找到路径上 w w w最大的边去替换它。然后想一下细节自己写吧,实现有点麻烦。

#include <cstdio>
#include <iostream>
#include <algorithm>
using namespace std;
const int M = 200005;
#define int long long
#define pii pair<int,int>
#define inf 1e12
int read()
{
	int x=0,f=1;char c;
	while((c=getchar())<'0' || c>'9') {if(c=='-') f=-1;}
	while(c>='0' && c<='9') {x=(x<<3)+(x<<1)+(c^48);c=getchar();}
	return x*f;
}
int n,m,k,tot,sum,ans,mc,pos,pf,out,x,y,f[M],fa[M];
int used[M],dep[M],p[M][20],mx[M][20],A[M][20];
struct edge
{
	int v,c,id,next;
}e[2*M];
struct node
{
	int u,v,c,w,id;
	bool operator < (const node &b) const
	{
		return w<b.w;
	}
}a[M];
int find(int x)
{
	if(x!=fa[x]) fa[x]=find(fa[x]);
	return fa[x];
}
void dfs(int u,int pr)
{
	p[u][0]=pr;
	dep[u]=dep[pr]+1;
	for(int i=1;i<20;i++)
	{
		int tmp=p[u][i-1];
		p[u][i]=p[tmp][i-1];
		mx[u][i]=max(mx[u][i-1],mx[tmp][i-1]);
		if(mx[u][i]==mx[u][i-1]) A[u][i]=A[u][i-1];
		if(mx[u][i]==mx[tmp][i-1]) A[u][i]=A[tmp][i-1];
	}
	for(int i=f[u];i;i=e[i].next)
	{
		int v=e[i].v;
		if(v==pr) continue;
		mx[v][0]=e[i].c;
		A[v][0]=e[i].id;
		dfs(v,u);
	}
}
void upd(int u,int i)
{
	if(mx[u][i]>x)
	{
		x=mx[u][i];
		y=A[u][i];
	}
}
void lca(int u,int v)
{
	if(dep[u]<dep[v]) swap(u,v);
	x=y=0;
	for(int i=19;i>=0;i--)
		if(dep[p[u][i]]>=dep[v])
		{
			upd(u,i);
			x=max(x,mx[u][i]);
			u=p[u][i];
		}
	if(u==v) return ;
	for(int i=19;i>=0;i--)
		if(p[u][i]^p[v][i])
		{
			upd(u,i);upd(v,i);
			u=p[u][i];v=p[v][i];
		}
	upd(u,0);upd(v,0);
	return ;
}
signed main()
{
	n=read();m=read();
	for(int i=1;i<=m;i++) a[i].w=read();
	for(int i=1;i<=m;i++) a[i].c=read(),a[i].id=i;
	for(int i=1;i<=m;i++) a[i].u=read(),a[i].v=read();
	k=read();
	sort(a+1,a+1+m);
	for(int i=1;i<=n;i++) fa[i]=i;
	mc=inf;
	for(int i=1;i<=m;i++)
	{
		int x=find(a[i].u),y=find(a[i].v);
		if(x==y) continue;
		fa[x]=y;sum+=a[i].w;
		if(mc>a[i].c)
			mc=a[i].c,pos=i;
		used[i]=1;
		e[++tot]=edge{a[i].u,a[i].w,i,f[a[i].v]},f[a[i].v]=tot;
		e[++tot]=edge{a[i].v,a[i].w,i,f[a[i].u]},f[a[i].u]=tot;
	}
	ans=sum-k/mc;dfs(1,0);
	for(int i=1;i<=m;i++)
	{
		if(mc<=a[i].c) continue;
		lca(a[i].u,a[i].v);
		int tmp=sum-x+a[i].w-k/a[i].c;
		if(tmp<ans) {ans=tmp;pf=i;out=y;}
	}
	printf("%lld\n",ans);
	if(!pf)
	{
		for(int i=1;i<=m;i++)
			if(used[i])
			{
				if(i==pos) printf("%lld %lld\n",a[i].id,a[i].w-k/mc);
				else printf("%lld %lld\n",a[i].id,a[i].w);
			}
	}
	else
	{
		for(int i=1;i<=m;i++)
			if(used[i] || i==pf)
			{
				if(i==out) continue;
				if(i==pf) printf("%lld %lld\n",a[i].id,a[i].w-k/a[i].c);
				else printf("%lld %lld\n",a[i].id,a[i].w);
			}
	}
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值