noip 2015 運輸計划 (lca+二分)


/*
95 最后一個點T了 qian lv ji qiong 了 
沒學過樹剖 聽chx聽xzc說的神奇的方法  Orz
首先求出每個計划的路徑長度 這里寫的倍增
然后二分答案
對於每個ans 統計>他的路徑條數 tot 並維護最大差值 dec 
並且對於每條不合法的路徑維護每個點的經過次數
然后枚舉點 如果經過次數==tot說明每一條不合法的都經過他
然后嘗試把它建成蟲洞 如果他對應邊的權值>=dec 那么我們刪掉它ans就合法了
關鍵是統計每個點在非法路徑中的經過次數 :
維護sum數組 對於每個非法的路徑起點a b LCA(a,b)==s sum[a]++ sum[b]++ sum[s]-=2
這樣網上更新的話 經過的點的sum值都變成1 祖先s的變成0 
這樣就實現了sum數組的維護 
*/
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#define maxn 300100
using namespace std;
int n,m,num,head[maxn],ans,inf;
int fa[maxn][30],dep[maxn],dis[maxn],sum[maxn],edge[maxn];
struct node
{
    int u,v,t,pre;
}e[maxn*2];
struct Ans
{
    int ai,bi,anc,di;
}lca[maxn];
int init()
{
    int x=0;char s=getchar();
    while(s<'0'||s>'9')s=getchar();
    while(s>='0'&&s<='9'){x=x*10+s-'0';s=getchar();}
    return x;
}
void Add(int from,int to,int dis)
{
    num++;
    e[num].u=from;
    e[num].v=to;
    e[num].t=dis;
    e[num].pre=head[from];
    head[from]=num;
}
void Dfs(int now,int from,int c,int Dis)
{
    fa[now][0]=from;
    dep[now]=c;dis[now]=Dis;
    for(int i=head[now];i;i=e[i].pre)
      if(e[i].v!=from)
        {
          edge[e[i].v]=i;
          Dfs(e[i].v,now,c+1,Dis+e[i].t);
        }
}
void Get_fa()
{
    for(int j=1;j<=16;j++)
      for(int i=1;i<=n;i++)
        fa[i][j]=fa[fa[i][j-1]][j-1];
}
int Get_same(int a,int t)
{
    for(int i=0;i<16;i++)
     if(t&(1<<i)) a=fa[a][i];
    return a;
}
int LCA(int a,int b)
{
    if(dep[a]<dep[b])swap(a,b);
    a=Get_same(a,dep[a]-dep[b]);
    if(a==b)return a;
    for(int i=16;i>=0;i--)
      if(fa[a][i]!=fa[b][i])
        {
          a=fa[a][i];
          b=fa[b][i];
        }
    return fa[a][0];
}
void Init()
{
    n=init();m=init();
    int u,v,t;
    for(int i=1;i<=n-1;i++)
      {
          u=init();v=init();t=init();
          Add(u,v,t);Add(v,u,t);
      }
    Dfs(1,1,0,0);
    Get_fa();
    for(int i=1;i<=m;i++)
      {
          lca[i].ai=init();lca[i].bi=init();
          lca[i].anc=LCA(lca[i].ai,lca[i].bi);
          lca[i].di=dis[lca[i].ai]+dis[lca[i].bi]-2*dis[lca[i].anc];
          inf=max(inf,lca[i].di);
      }
}
void Up_sum(int now,int from)
{
    for(int i=head[now];i;i=e[i].pre)
      if(e[i].v!=from)
        {
          Up_sum(e[i].v,now);
          sum[now]+=sum[e[i].v];
        }
}
int Judge(int x)
{
    memset(sum,0,sizeof(sum));
    int tot=0,dec=0;
    for(int i=1;i<=m;i++)
      if(lca[i].di>x)//非法路徑 
        {
          tot++;
          dec=max(dec,lca[i].di-x);//最長非法路徑與ans差值 
          sum[lca[i].ai]++;
          sum[lca[i].bi]++;
          sum[lca[i].anc]-=2;
        }
    Up_sum(1,1);//更新sum數組 
    for(int i=1;i<=n;i++)
      if(tot==sum[i]&&e[edge[i]].t>=dec)//刪掉edge[i]這條邊之后答案合法了 
        return 1;
    return 0;
}
void Solve()//二分答案 
{
    int l=0,r=inf;
    while(l<=r)
      {
          int mid=(l+r)/2;
          int tmp=Judge(mid);
          if(tmp==1)
            {
                r=mid-1;
                ans=mid;
          }
        else l=mid+1;
      }
}
void Printf()
{
    printf("%d\n",ans);
}
int main()
{
    //freopen("transport.in","r",stdin);
    //freopen("transport.out","w",stdout);
    Init();
    Solve();
    Printf();
    return 0;
}

 


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM