NOIP提高组2018 D1T3 【赛道修建】

颓了好几天,终于把这到题处理了一下。

话说,其实我考场上想出正解了,但是手残,算复杂度的时候多按了一个零,导致算出来是1亿多的复杂度,都不敢打。。。就把部分分都捡了一下。。。

题目描述:

C 城将要举办一系列的赛车比赛。在比赛前,需要在城内修建 m 条赛道。

C 城一共有 n 个路口,这些路口编号为 1,2,…,n1,有 n1 条适合于修建赛道的双向通行的道路,每条道路连接着两个路口。其中,第 ii 条道路连接的两个路口编号为 ai 和 bi,该道路的长度为 li。借助这 n-1n1 条道路,从任何一个路口出发都能到达其他所有的路口。

一条赛道是一组互不相同的道路 e1,e2,…,ek满足可以从某个路口出发,依次经过 道路 e1,e2,…,ek(每条道路经过一次,不允许调头)到达另一个路口。一条赛道的长度等于经过的各道路的长度之和。为保证安全,要求每条道路至多被一条赛道经过。

目前赛道修建的方案尚未确定。你的任务是设计一种赛道修建的方案,使得修建的 m 条赛道中长度最小的赛道长度最大(即 m 条赛道中最短赛道的长度尽可能大)

思路分析:

从部分分开始讲吧。

菊花图:

这个部分分还是挺好拿的,把边从大到小排序之后,直接二分答案。两个指针,i,正着来,j,倒着来,i就顺着取下去,对于每个i给它匹配一个最大的相加刚好大于二分值的j(一开始就大于二分值的要特殊处理),注意到这个决策是有单调性的,即:当i1小于i2时,那么相匹配的j1一定大于等于j2,这样的话,我们每次只要扫一遍就能得到答案(考完后有人和我说求菊花图要log2二分答案套二分,嗯,就是我们的zyk同志,其实我一直都觉得这个人的脑回路特别清奇。。。)。

分支不超过3:

这个的话其实就是正解的简化版了,考虑二分答案加树形DP。

  dp[u]:表示以u为根的子树中有多少长度大于mid的链。

  f[u]:表示以u为根的子树中,除去那些被统计到dp数组中的链,现在能传到父亲的链的长度最长的链的长度(貌似有点拗口呢。。。)。

转移应该很简单了吧。

首先,dp[u]=∑v∈u's sons 如果u的两个儿子的最长链加起来大于mid,那么再给dp[u]+1,f[u]=0,否则的话,f[u]=max(f[v1],f[v2])。

正解:

正解其实就是把两者结合一下啦。

但是特别注意的就是:我们要在满足匹配数最大的情况下,找一条最长的链传给父亲。

那么这个要求我们应该怎么来实现呢?——二分嘛!

我们可以二分我们需要传上去的那条链,然后在匹配的时候跳过这条链就行了。(其实这是成一类dalao教我的)

代码实现:(没想到吧,我转C++辣!!!)

(不得不说一句——C++压行是真的爽,啊哈哈哈!!!)

#include <bits/stdc++.h>
using namespace std;
const int maxn=50005;
int dp[maxn],f[maxn],head[maxn],nxt[2*maxn],vet[2*maxn],dist[2*maxn];
int tot,mid;
void add(int x,int y,int z){
    tot++;
    nxt[tot]=head[x];
    vet[tot]=y;
    head[x]=tot;
    dist[tot]=z;
}
bool cmp(int a,int b) {
    return a>b;
}
void dfs(int u,int father){
    vector<int> a;
    for (int i=head[u];i;i=nxt[i]){
        int v=vet[i];
        if (v==father) continue;
        dfs(v,u);
        f[v]+=dist[i];
        if (f[v]<mid) a.push_back(f[v]); else dp[u]++;
        dp[u]+=dp[v];
    }
    int l=0,r=a.size()-1;
    sort(a.begin(),a.end(),cmp);
    a.push_back(0);
    int ans=a.size()-1,Max=0,i=0,j=a.size()-2;
    while (i<j) if (a[i]+a[j]>=mid) i++,j--,Max++; else j--;
    while (l<=r){
        int now=(l+r)>>1,res=0,i=0,j=a.size()-2;
        while (i<j){
            if (i==now) {i++; continue;}
            if (j==now) {j--; continue;}
            if (a[i]+a[j]>=mid) i++,j--,res++; else j--;
        }
        if (res>=Max) Max=res,ans=now,r=now-1; else l=now+1;
    }
    dp[u]+=Max; f[u]=a[ans];
}
int main(){
    int n,m;
    scanf("%d%d",&n,&m);
    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);
    }
    int l=1,r=500000000,ans=0;
    while (l<=r){
        memset(dp,0,sizeof(dp));
        memset(f,0,sizeof(f));
        mid=(l+r)>>1;
        dfs(1,0);
        if (dp[1]>=m) ans=mid,l=mid+1; else r=mid-1;
    }
    printf("%d\n",ans);
    return 0;
}

猜你喜欢

转载自www.cnblogs.com/WR-Eternity/p/9964177.html