【Codeforces 733F】 Drivers Dissatisfaction 【MST+bianry lifting】

题意:

      给出一个图,每条边都有权值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;
}


  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值