BZOJ4326&&洛谷P2680 运输计划

倍增LCA+差分+二分答案

简称毒瘤题

我们预处理从根到每个点的dis值,然后把每条边的边权放到点上,然后处理每条航线起点终点的lca,然后我们二分一个改完权值的的路线长,去检验是否能达到,如何检验呢?

先把所有路线长大于mid的路线捡出来,然后看他和mid差多少,取一个最大值,这是要改造的道路的最小值,若比他还小,那么dis-x>mid不符合条件,然后是选出来的道路被经过次数要大于等于上面捡出来的路的条数,不然也不和题意,然后就没啦

代码

//By AcerMo
#include
#include
#include
#include
#include
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]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;i1,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 (lint mid=(l+r)>>1;
        if (che(mid)) r=mid;
        else l=mid+1;
    }
    cout<return 0;
}

你可能感兴趣的:(差分,倍增LCA)