codeforces 733F LCA+最小生成树

点击打开链接

题意：给一个图，及每条边的边权，然后每条边都可以使自己的边权减少1，但是需要费用ci，然后你有S的费用，问你如何分配这些费用使得最小生成树的值最小

思路：首先要明确对于这所有的费用，我是全部用在一条边上的，也就是说只用修改一条边的权值，那么我们可以先将原图的最小生成树找出来，然后建个LCA 的图，之后枚举每条边，如果是最小生成树上的边，那么直接减然后判断即可， 但是如果不是的话，那么仔LCA的那个树上就会形成一个环，在这个环中找出最大边权的那个边即可   PS：第一发超时，然后加了两个剪枝，第一个是如果边权大于你拥有的费用，那么这条边是没用的，第二个是如果你减边权的这个边减完之后，价值还是比之前最优的大，那么这个边也没用了，因为你这个边的权值之前没在最小生成树中，所以它的边权肯定是比环上的任意一个大的   加了这两个之后快了好多

#include <vector>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int inf=0x3f3f3f3f;
const ll INF=0x3f3f3f3f3f3f3f3f;
const int maxn=200010;
bool vis[maxn];
int n,m,k;
int L[maxn*2],E[maxn*2],H[maxn],dis[maxn],dp[2*maxn][20],num[maxn],sum[maxn],head[maxn],kkk,pre[maxn];
struct node{
    int to,next,cost;
}EE[maxn*4];
void add_edge(int u,int v,int cost){
    EE[kkk].to=v;EE[kkk].cost=cost;EE[kkk].next=head[u];head[u]=kkk++;
}
void dfs(int t,int deep){
    k++;E[k]=t;L[k]=deep;H[t]=k;
    for(int i=head[t];i!=-1;i=EE[i].next){
        int tt=EE[i].to;
        if(!vis[tt]){
            vis[tt]=1;pre[tt]=t;
            dis[tt]=dis[t]+EE[i].cost;
            dfs(tt,deep+1);
            k++;E[k]=t;L[k]=deep;
        }
    }
}
void RMQ_init(){
    for(int i=1;i<=2*n-1;i++) dp[i][0]=i;
    for(int i=1;(1<<i)<=2*n-1;i++){
        for(int j=1;j+(1<<i)-1<=2*n-1;j++){
            if(L[dp[j][i-1]]<L[dp[j+(1<<(i-1))][i-1]]) dp[j][i]=dp[j][i-1];
            else dp[j][i]=dp[j+(1<<(i-1))][i-1];
        }
    }
}
int RMQ(int le,int ri){
    le=H[le];ri=H[ri];
    if(le>ri) swap(le,ri);
    int kk=0;
    while((1<<(kk+1))<=ri-le+1) kk++;
    if(L[dp[le][kk]]<L[dp[ri-(1<<kk)+1][kk]]) return E[dp[le][kk]];
    else return E[dp[ri-(1<<kk)+1][kk]];
}
struct edge1{
    int u,v,cost,id;
};
edge1 es[maxn];
bool cmp(const edge1 &a,const edge1 &b){
    return a.cost<b.cost;
}
int V,EEE,f[maxn],fvis[maxn];
int find1(int x){
    if(x!=f[x]) f[x]=find1(f[x]);
    return f[x];
}
void unite(int a,int b){
    int aa=find1(a);
    int bb=find1(b);
    if(aa==bb) return ;
    f[aa]=bb;
}
ll kruskal(){
    sort(es,es+EEE,cmp);
    ll ans=0;
    for(int i=0;i<=V;i++) f[i]=i;
    for(int i=0;i<EEE;i++){
        edge1 e=es[i];
        if(find1(e.u)!=find1(e.v)){
            ans+=(ll)(es[i].cost);
            fvis[es[i].id]=1;
            unite(e.u,e.v);
            add_edge(e.u,e.v,e.cost);
            add_edge(e.v,e.u,e.cost);
        }
    }
    return ans;
}
inline int getint(){
    int res=0;
    char c=getchar();
    bool mi=false;
    while(c<'0' || c>'9') mi=(c=='-'),c=getchar();
    while('0'<=c && c<='9') res=res*10+c-'0',c=getchar();
    return mi ? -res : res;
}
int C[maxn],val[maxn];
int main(){
    int mon;
    while(scanf("%d%d",&n,&m)!=-1){
        EEE=m;V=n;
        for(int i=0;i<=n;i++) vis[i]=0,dis[i]=0;
        for(int i=0;i<=n;i++) head[i]=-1;
        for(int i=0;i<m;i++) es[i].cost=getint(),fvis[i]=0;
        for(int i=0;i<m;i++) val[i]=getint();
        for(int i=0;i<m;i++) es[i].u=getint(),es[i].v=getint(),es[i].id=i;
        mon=getint();
        ll fans=kruskal();
        ll fffans=fans;
        int fpos1=inf,fpos2=inf,fpos3=inf;
        kkk=0;k=0;vis[1]=1;dfs(1,1);RMQ_init();
        for(int i=0;i<m;i++){
            ll ftmp=fans;
            if(val[es[i].id]>mon) continue;
            if(fvis[es[i].id]){
                ftmp-=mon/val[es[i].id];
                if(ftmp<fffans){
                    fffans=ftmp;
                    fpos1=es[i].id;fpos2=inf;fpos3=inf;
                }
            }else{
                if(ftmp-mon/val[es[i].id]>fffans) continue;
                ftmp=ftmp+es[i].cost-mon/val[es[i].id];
                int ffa=RMQ(es[i].u,es[i].v);
                int min1=inf;
                int uuu=es[i].u,vvv=es[i].v;
                while(uuu!=ffa){
                    int tt1=dis[uuu]-dis[pre[uuu]];
                    if(ftmp-tt1<fffans){
                        fffans=ftmp-tt1;
                        fpos1=es[i].id;
                        fpos2=uuu;fpos3=pre[uuu];
                    }
                    uuu=pre[uuu];
                }
                while(vvv!=ffa){
                    int tt1=dis[vvv]-dis[pre[vvv]];
                    if(ftmp-tt1<fffans){
                        fffans=ftmp-tt1;
                        fpos1=es[i].id;
                        fpos2=vvv;fpos3=pre[vvv];
                    }
                    vvv=pre[vvv];
                }
            }
        }
        printf("%I64d\n",fffans);
        if(fpos2==inf){
            if(fpos1==inf){
                for(int i=0;i<m;i++){
                    if(fvis[es[i].id]) printf("%d %d\n",es[i].id+1,es[i].cost);
                }
            }else{
                for(int i=0;i<m;i++){
                    if(fvis[es[i].id]){
                        if(fpos1==es[i].id) printf("%d %d\n",es[i].id+1,es[i].cost-mon/val[es[i].id]);
                        else printf("%d %d\n",es[i].id+1,es[i].cost);
                    }
                }
            }
        }else{
            fvis[fpos1]=1;
            for(int i=0;i<m;i++){
                if((fpos2==es[i].u&&fpos3==es[i].v)||(fpos2==es[i].v&&fpos3==es[i].u)) fvis[es[i].id]=0;
            }
            for(int i=0;i<m;i++){
                if(fvis[es[i].id]){
                    if(fpos1==es[i].id) printf("%d %d\n",es[i].id+1,es[i].cost-mon/val[es[i].id]);
                    else printf("%d %d\n",es[i].id+1,es[i].cost);
                }
            }
        }
    }
    return 0;
}