BZOJ4326&&洛谷P2680 运输计划

倍增LCA+差分+二分答案

简称毒瘤题

我们预处理从根到每个点的dis值,然后把每条边的边权放到点上,然后处理每条航线起点终点的lca,然后我们二分一个改完权值的的路线长,去检验是否能达到,如何检验呢?
先把所有路线长大于mid的路线捡出来,然后看他和mid差多少,取一个最大值,这是要改造的道路的最小值,若比他还小,那么dis-x>mid不符合条件,然后是选出来的道路被经过次数要大于等于上面捡出来的路的条数,不然也不和题意,然后就没啦

代码

//By AcerMo
#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int M=600500;
int n,m,flag;
int to[M],nxt[M],w[M],head[M],cnt;
int dep[M],f[M][21],sum[M],pat[M],dis[M];
struct edge{int fr,to,l;}e[M];
inline int read()
{
    int x=0;char ch=getchar();
    while (ch>'9'||ch<'0') ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-'0',ch=getchar();
    return x;
}
inline void add(int x,int y,int z)
{
    to[++cnt]=y;w[cnt]=z;nxt[cnt]=head[x];head[x]=cnt;
    to[++cnt]=x;w[cnt]=z;nxt[cnt]=head[y];head[y]=cnt;
    return ;
}
inline void dfs(int x,int fa)
{
    dep[x]=dep[fa]+1;f[x][0]=fa;
    for (int i=1;i<=20;i++)
    f[x][i]=f[f[x][i-1]][i-1];
    for (int i=head[x];i;i=nxt[i])
    if (to[i]!=fa)
    {
        dis[to[i]]=dis[x]+w[i];
        dfs(to[i],x);pat[to[i]]=w[i];
    } 
    return ;
}
inline int lca(int x,int y)
{
    if (dep[x]<dep[y]) swap(x,y);
    for (int i=20;i>=0;i--)
    if (dep[f[x][i]]>=dep[y])
    x=f[x][i];
    if (x==y) return x;
    for (int i=20;i>=0;i--)
    if (f[x][i]!=f[y][i])
    x=f[x][i],y=f[y][i];
    return f[x][0]; 
} 
inline int getsum(int x,int fa,int tot,int mx)
{
    int ans=sum[x];
    for (int i=head[x];i;i=nxt[i])
    if (to[i]!=fa) ans+=getsum(to[i],x,tot,mx);
    if (ans>=tot&&mx<=pat[x]) flag=1;
    return ans;
}
inline bool che(int x)
{
    fill(sum,sum+n+1,0);
    flag=0;int mx=0,tot=0;
    for (int i=1;i<=m;i++)
    if (dis[e[i].fr]+dis[e[i].to]-2*dis[e[i].l]>x)
    {
        sum[e[i].fr]++;sum[e[i].to]++;sum[e[i].l]-=2;
        tot++;mx=max(mx,dis[e[i].fr]+dis[e[i].to]-2*dis[e[i].l]);
    }
    if (cnt==0) return 1;
    getsum(1,0,tot,mx-x);
    return flag;
}
signed main()
{
    n=read();m=read();int x,y,z;
    for (int i=1;i<n;i++) 
    x=read(),y=read(),z=read(),add(x,y,z);
    dfs(1,0);
    for (int i=1;i<=m;i++)
    e[i].fr=read(),e[i].to=read(),e[i].l=lca(e[i].fr,e[i].to);
    int l=0,r=300000000;
    while (l<r)
    {
        int mid=(l+r)>>1;
        if (che(mid)) r=mid;
        else l=mid+1;
    }
    cout<<l;
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值