貪欲P1084 NOIP2012T6流行制御半分LCA + +

質問の意味：これらの力は、ブロックされた軍の一つ以上により、リーフ・ノードへのルートからのすべてのパスを見つけるように、複数のノードが、ツリーにあります。

すべての軍が同時に移動することができますので、この問題があることを単純化することができるので、すべての部隊を移動するための最も長い時間を見つけて、時間を最小限にするために

ソリューション：LCA-半分+ +貪欲

1.二部;与えられたトピックとしてデータが非常に大きいですが、また、それぞれのアプローチのための明白なテストになりますが、十分ではない、それが必要である二部

2.貪欲;限り、各ステップが上がっていくと軍隊のいずれか、その軍隊は、そう、明らかに優れたマルチノードを制御することができるようになります。アプローチ：各軍は、時間の2分以内に、できるだけ多くを行くことにしましょう、このとき、2例に分けることができます。

<1>ないルートに、それはそこに滞在しましょう

<2>軍に残っているどのくらいの時間を記録するために、ルーツになりました

軍隊がない場合は、この後者の場合には、十分な時間が、それは息子、無軍によるそれの根の息子の後にそのルーツに戻る作るために残され、それが戻ってみましょう。もしそれは戻すのに十分な時間の軍の残りの部分は、その後、去ります

次に、それはまた、すべての人の息子の根が残っているん軍はまた、点列が来ない、シーケンス行武装勢力の根の上に残ります。二つの配列を比較するために、あなたがで決定することができます軍隊ならば、すべての流行は制御下に置かすることができ、この時間

3.LCAは、ジャンプアップ後の半分の軍隊たびので、それはそれの半分を使用することが明らかにされていませんか？

コードを接続します。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
const int N=6e4;

int n,m,t,tot=0,atot=0,btot=0,ctot=0;
int d[N],query[N],f[N][20];
int ver[2*N],edge[2*N],Next[2*N],head[N];
bool ok,sta[N],need[N];
ll ans,tim[N],ned[N],dist[N][20];
pair<ll,int> h[N];
queue<int> q;

void add(int x,int y,int z){
    ver[++tot]=y,edge[tot]=z,Next[tot]=head[x],head[x]=tot;
}

void bfs(){
    q.push(1);
    d[1]=1;
    while(q.size()){
        int x=q.front();q.pop();
        for(int i=head[x];i;i=Next[i]){
            int y=ver[i];
            if(d[y]) continue;
            d[y]=d[x]+1;
            f[y][0]=x,dist[y][0]=edge[i];
            for(int j=1;j<=21&&y-(1<<j)>=1;j++){
                f[y][j]=f[f[y][j-1]][j-1];
                dist[y][j]=dist[y][j-1]+dist[f[y][j-1]][j-1];
            }
            q.push(y);
        }
    }
}

bool dfs(int x){
    bool pson=0;
    if(sta[x]) return 1;
    for(int i=head[x];i;i=Next[i]){
        int y=ver[i];
        if(d[y]<d[x]) continue;
        pson=1;
        if(!dfs(y)) return 0;
    }
    if(!pson) return 0;
    return 1;
}

bool check(ll lim)
{
    memset(sta,0,sizeof(sta));
    memset(tim,0,sizeof(tim));
    memset(ned,0,sizeof(ned));
    memset(h,0,sizeof(h));
    memset(need,0,sizeof(need));
    atot=0,btot=0,ctot=0;
    for(int i=1;i<=m;i++){
        ll x=query[i],cnt=0;
        for(int j=t;j>=0;j--){
            if(f[x][j]>1&&cnt+dist[x][j]<=lim){
                cnt+=dist[x][j];
                x=f[x][j];
            }
        }
        if(f[x][0]==1&&cnt+dist[x][0]<=lim) h[++ctot]=make_pair(lim-cnt-dist[x][0],x);
        else sta[x]=1;
    }
    for(int i=head[1];i;i=Next[i]) if(!dfs(ver[i])) need[ver[i]]=1;
    sort(h+1,h+ctot+1);
    for(int i=1;i<=ctot;i++){
        if(need[h[i].second]&&h[i].first<dist[h[i].second][0]) need[h[i].second]=0;
        else tim[++atot]=h[i].first;//军队
    }
    for(int i=head[1];i;i=Next[i]) if(need[ver[i]]) ned[++btot]=dist[ver[i]][0];//节点距离 
    if(atot<btot) return 0;
    sort(tim+1,tim+atot+1);
	sort(ned+1,ned+btot+1);
    int i=1,j=1;
    while(i<=btot&&j<=atot){
        if(tim[j]>=ned[i]) i++,j++;
        else j++;
    }
    if(i>btot) return 1;
    return 0;
}

int main()
{
	ll l=0,r=0,mid;
    cin>>n;
    t=log2(n)+1;
    for(int i=1;i<n;i++){
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        add(x,y,z),add(y,x,z);
        r+=z;
    }
    bfs();
    cin>>m;
    for(int i=1;i<=m;i++) scanf("%d",&query[i]);
    while(l<=r){
        mid=(l+r)>>1;
        if(check(mid)){
            r=mid-1;
            ans=mid;
            ok=true;
        }
        else l=mid+1;
    }
    if(!ok) cout<<-1;
    else cout<<ans;
    return 0;
}