Noip2018Day1T3軌道建設

トピックリンク

問題

私たちは、木の辺の重みを与えています。鎖の一方の鎖の右側縁は、重量値と経過として定義されます。選出する\(m個\)鎖を求めて\(m個\)を体重鎖において最小重み鎖は最大数です。

解決

まず、明らかに半分。

そして、回答のうち、半分を判断する方法を検討\(X \)それが実現可能であるかどうか。すなわち、選択されたかどうか\(m個\)チェーン、各鎖は、以上の重みである(X \)\実際には、これは貪欲です。

線形祖先は、パスのポイントツーポイントから定義されます。

各チェーンが見つかりましたかまで取得するには、直鎖、またはいくつかの点で2つの線形の合併であることができます。

貪欲な場所は各点の直鎖合成確かに優先順位が合成されますされています。次いで、最大直鎖長アップを確実にするために転送します。でも特に大きな貢献の長さを渡さので生成されただけ以上のことを行う\(1 \) これは、まず現在のサブツリー内の最大の鎖の合成を確認する必要があります。

問題はその後のサブツリー内いくつかの直鎖長を取得することになります。これらは現在、リニア二十から二量以上に結合されている\(X \)チェーン。残りの最大直線長ことを確認してください。

ここでは、半分の答えを見ることができます。また、使用することができます\(マルチセット\)プロセスを。とにかく、それは行うことができます問題があります。

コードポイントのストールの部分、BF5正のソリューションがあります。

コード

#include<set>
#include<cstdio>
#include<iostream>
#include<cstdlib>
#include<ctime>
#include<algorithm>
#include<queue>
#include<cstring>
using namespace std;
typedef long long ll;
const int N = 100010;
ll read() {
    ll x = 0,f = 1;char c = getchar();
    while(c < '0' || c > '9') {
        if(c == '-') f = -1;
        c = getchar();
    }
    while(c >= '0' && c <= '9') {
        x = x * 10 + c - '0';
        c = getchar();
    }
    return x * f;
}
int n,m;
struct node {
    int v,nxt,w;
}e[N << 1];
int head[N],ejs;
void add(int u,int v,int w) {
    e[++ejs].v = v;e[ejs].nxt = head[u];head[u] = ejs;e[ejs].w = w;
}
namespace BF1 {
    int dis[N];
    void dfs(int u,int fa) {
        for(int i = head[u];i;i = e[i].nxt) {
            int v = e[i].v;
            if(v == fa) continue;
            dis[v] = dis[u] + e[i].w;
            dfs(v,u);
        }
    }
    void main() {
        dfs(1,0);
        int x = 0;
        for(int i = 1;i <= n;++i) if(dis[i] > dis[x]) x = i;
//      cout<<x<<endl;
        memset(dis,0,sizeof(dis));
        dfs(x,0);
        int ans = 0;
        for(int i = 1;i <= n;++i) ans = max(ans,dis[i]);
        cout<<ans;
    }
}

namespace BF2 {
    int a[N],cnt;
    int check(int x) {
        int p = 1,ret = 0;
        for(int i = cnt;i > p;--i) {
            if(a[i] > x && i > p) {ret++;continue;}
            while(a[p] + a[i] < x && p < i) ++p;
            if(p < i) ret++,p++;
            else break;
        }
        return ret;
    }
    void main() {
        int l = 100000,r = 0;
        for(int i = 1;i <= ejs;i += 2) a[++cnt] = e[i].w,l = min(l,a[cnt]),r += a[cnt];
        sort(a + 1,a + cnt + 1);

        int ans = 0;
        while(l <= r) {
            int mid = (l + r) >> 1;
            if(check(mid) >= m) ans = mid,l = mid + 1;
            else r = mid - 1;
        }
        cout<<ans<<endl;
    }
}
int du[N];
namespace BF3 {
    
    int a[N],cnt;
    void dfs(int u,int fa) {
        for(int i = head[u];i;i = e[i].nxt) {
            int v = e[i].v;
            if(v == fa) continue;
            a[++cnt] = e[i].w;
            dfs(v,u);
        }
    }
    int check(int x) {
        int now = 0,ret = 0;
        for(int i = 1;i <= cnt;++i) {
            now += a[i];
            if(now >= x) now = 0,ret ++;
        }
        return ret;
    }
    void main() {
        for(int i = 1;i <= n;++i) 
            if(du[i] == 1) {dfs(i,0);break;}
        int l = 1000000,r = 0;
        for(int i = 1;i <= ejs;i += 2) {
            l = min(l,e[i].w);r += e[i].w;
        }
        
        int ans = 0;
        while(l <= r) {
            int mid = (l + r) >> 1;
            if(check(mid) >= m) ans = mid,l = mid + 1;
            else r = mid - 1;
        }
        cout<<ans<<endl;
        
    }

}
int L = 100000,R;
namespace BF5 {
    int ANS;
    int dfs(int u,int fa,int x) {
        multiset<int>s;
        int ret = 0;
//      if(!s.empty()) printf("%d\n",u);
        for(int i = head[u];i;i = e[i].nxt) {
            int v = e[i].v;
            if(v == fa) continue;
            int k = dfs(v,u,x);
            if(k + e[i].w >= x) ANS++;
            else s.insert(k + e[i].w);
        }
        while(!s.empty()) {
            multiset<int>::iterator it = s.begin();
            int k = *it;
            s.erase(it);
            multiset<int>::iterator is = s.lower_bound(x - k);
            if(is == s.end()) ret = max(ret,k);
            else ANS++,s.erase(is);
        }
//      s.clear();
        
//      printf("%d %d\n",u,ret);

        return ret;
    }
    void main() {
        int l = L,r = R,ans = 0;
        
        while(l <= r) {

            int mid =(l + r) >> 1;
            ANS = 0;dfs(1,0,mid);
            if(ANS >= m) ans = mid,l = mid + 1;
            else r = mid - 1;
        }
        cout<<ans;
    }
}
int main() {
    n = read(),m = read();
    int bz1 = 1,bz2 = 1;
    for(int i = 1;i < n;++i) {
        int u = read(),v = read(),w = read();
        L = min(L,w);R += w;    
        du[u]++;du[v]++;
    
        add(u,v,w);add(v,u,w);
        if(u != 1) bz1 = 0;
        if(v != u + 1) bz2 = 0;
    }
    
    if(m == 1) {BF1::main();return 0;}
    if(bz1) {BF2::main();return 0;}
    if(bz2) {BF3::main();return 0;}
    BF5::main();
    return 0;
}
/*
7 1
1 2 10
1 3 5
2 4 9
2 5 8
3 6 6
3 7 7
*/

おすすめ

転載: www.cnblogs.com/wxyww/p/luogu5021.html