[Wc2010]重建计划

最新推荐文章于 2021-03-01 22:16:19 发布
最新推荐文章于 2021-03-01 22:16:19 发布
阅读量224 收藏
点赞数
分类专栏： 树分治 文章标签： 点分治 神题
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/tomyking2005/article/details/86482487
版权

                                  [Wc2010]重建计划

 

Description

Input

第一行包含一个正整数N，表示X国的城市个数. 第二行包含两个正整数L和U，表示政策要求的第一期重建方案中修建道路数的上下限 接下来的N-1行描述重建小组的原有方案，每行三个正整数AiBiVi分别表示道路(AiBi)其价值为Vi 其中城市由1..N进行标号

Output

输出最大平均估值，保留三位小数

Sample Input

4 
2 3 
1 2 1 
1 3 2 
1 4 3

 Sample Output

2.500

Solution

首先，请先知道这道题是 01分数规划问题 的变形。

解决办法：二分+点分治

  • 二分？

二分最长的路有多长，再检查能不能行。

点分治？

点分治就是用来检查的……（此处省略找重心等等）。简明起见此处设将要分治的点设为\emph{v}

准备：

  1. 队列q、q1(q1用于回收出队q的元素),双向队列dq;
  2. 数组dis、deep、md(maxdeep)。

过程 ：

  1. 将树中所有边权值-mid;
  2. 找出一个与\emph{v}相连的点\emph{u};
  3. 从\emph{u}开始广搜（当然有vis啦），记录每个点的dis和deep;
  4. q=q1;
  5. 取出q中每一个元素\emph{x};
  6. 将\emph{x}放入dq;
  7. 维护dq单调性;
  8. 检查是否有链权值和>0;
  9. 维护各节点md;
  10. 清空q+q1+dq并遍历下一个\emph{u}。

  • SO？

实际上点分治就是在找-mid后是否有链权值和>0!

Code

/*
	Name: build.cpp
	Author: YuckXi
	Date: 14-01-19 20:18
*/
#include<iostream>
#include<algorithm>
#include<queue>
#include<deque>
#include<vector>
#include<cstdio>
using namespace std;
const int maxn=200010;
const int inf=1000000000;
const double eps=1e-4;
int n,m,son[maxn],mxson[maxn],root,h[maxn];
struct ed{
    int y,w;
};
vector<ed>l[maxn];
int tot=1,fa[maxn];
int size,deep[maxn],L,U,mx=0;
double dis[maxn],md[maxn],ans=0;
bool vis[maxn];
void getroot(int v,int fa){//重心
    mxson[v]=0; 
    son[v]=1;
    for (int i=0;i<l[v].size();i++){
        int u=l[v][i].y;
        if (u==fa||vis[u]) continue;
        getroot(u,v);
        son[v]+=son[u];
        mxson[v]=max(mxson[v],son[u]);
    }
    mxson[v]=max(mxson[v],size-son[v]);
    if (mxson[v]<mxson[root]) root=v;
}
bool check(int v,double mid){//二分的CHECK
    queue<int>q,q1;
    deque<int>dq;
    while(!q.empty())q.pop();
    while(!q1.empty())q1.pop();
    dq.clear();
    int tmp=0;
    for (int i=0;i<l[v].size();i++){
        queue<int>q,q1;//清空
        deque<int>dq;
        while(!q.empty())q.pop();
        while(!q1.empty())q1.pop();
        dq.clear();
        int u=l[v][i].y;
        if (vis[u]) continue;
        q.push(u);
        fa[u]=root; 
        deep[u]=1;
        dis[u]=(double)l[v][i].w-mid;
        while(!q.empty()){//bfs
            int x=q.front();q.pop();q1.push(x);
            for (int j=0;j<l[x].size();j++){
                int xx=l[x][j].y; 
                if (xx==fa[x]||vis[xx]) continue;
                q.push(xx);
                fa[xx]=x; 
                deep[xx]=deep[x]+1;
                dis[xx]=dis[x]+l[x][j].w-mid;
            }
        }
        q=q1;
        while(!q1.empty())q1.pop();
        int now=tmp;
        while(!q.empty()){
            int x=q.front();
            q.pop();
            q1.push(x);
            while(deep[x]+now>=L&&now>=0){//维护单调性
                while(dq.size()>=1&&md[dq.back()]<=md[now])dq.pop_back();
                dq.push_back(now);
                now--;
            }
            while(dq.size()>=1&&deep[x]+dq.front()>U)dq.pop_front();
            if(dq.size()>=1&&dis[x]+md[dq.front()]>=0)return 1;
        }
        q=q1;
        for(int j=tmp+1;j<=deep[q.back()];j++)md[j]=-inf;
        while(!q.empty()){
            int x=q.front();
            q.pop();
            md[deep[x]]=max(md[deep[x]],dis[x]);
        }
        tmp=max(tmp,deep[q.back()]);
    }
    return 0;
}
  
void calc(int x){//二分
    double l=ans,r=mx;
    while (r-l>eps){
        double mid=(l+r)/2;
        if (check(x,mid)) l=mid;
        else r=mid;
    }
    ans=l;
}
void FZ(int v){//分治
    calc(v); 
    vis[v]=1;
    for(int i=0;i<l[v].size();i++){
        int u=l[v][i].y;
        if(vis[u]) continue;
        size=son[u];
        root=0;
        getroot(u,root);
        if(son[u]>L) FZ(root);
    }
}
int main(){
    scanf("%d%d%d",&n,&L,&U);
    for (int i=1;i<n;++i){
        int x,y,w;
        scanf("%d%d%d",&x,&y,&w);
        l[x].push_back((ed){y,w});
        l[y].push_back((ed){x,w});
        mx=max(mx,w);
    }
    mxson[0]=size=n; 
    root=0;
    getroot(1,root);
    FZ(root);
    printf("%.3lf",ans);
}

无注释版

#include<iostream>
#include<algorithm>
#include<queue>
#include<deque>
#include<vector>
#include<cstdio>
using namespace std;
const int maxn=200010;
const int inf=1000000000;
const double eps=1e-4;
int n,m,son[maxn],mxson[maxn],root,h[maxn];
struct ed{
    int y,w;
};
vector<ed>l[maxn];
int tot=1,fa[maxn];
int size,deep[maxn],L,U,mx=0;
double dis[maxn],md[maxn],ans=0;
bool vis[maxn];
void getroot(int v,int fa){
    mxson[v]=0; 
    son[v]=1;
    for (int i=0;i<l[v].size();i++){
        int u=l[v][i].y;
        if (u==fa||vis[u]) continue;
        getroot(u,v);
        son[v]+=son[u];
        mxson[v]=max(mxson[v],son[u]);
    }
    mxson[v]=max(mxson[v],size-son[v]);
    if (mxson[v]<mxson[root]) root=v;
}
bool check(int v,double mid){
    queue<int>q,q1;
    deque<int>dq;
    while(!q.empty())q.pop();
    while(!q1.empty())q1.pop();
    dq.clear();
    int tmp=0;
    for (int i=0;i<l[v].size();i++){
        queue<int>q,q1;
        deque<int>dq;
        while(!q.empty())q.pop();
        while(!q1.empty())q1.pop();
        dq.clear();
        int u=l[v][i].y;
        if (vis[u]) continue;
        q.push(u);
        fa[u]=root; 
        deep[u]=1;
        dis[u]=(double)l[v][i].w-mid;
        while(!q.empty()){
            int x=q.front();q.pop();q1.push(x);
            for (int j=0;j<l[x].size();j++){
                int xx=l[x][j].y; 
                if (xx==fa[x]||vis[xx]) continue;
                q.push(xx);
                fa[xx]=x; 
                deep[xx]=deep[x]+1;
                dis[xx]=dis[x]+l[x][j].w-mid;
            }
        }
        q=q1;
        while(!q1.empty())q1.pop();
        int now=tmp;
        while(!q.empty()){
            int x=q.front();
            q.pop();
            q1.push(x);
            while(deep[x]+now>=L&&now>=0){
                while(dq.size()>=1&&md[dq.back()]<=md[now])dq.pop_back();
                dq.push_back(now);
                now--;
            }
            while(dq.size()>=1&&deep[x]+dq.front()>U)dq.pop_front();
            if(dq.size()>=1&&dis[x]+md[dq.front()]>=0)return 1;
        }
        q=q1;
        for(int j=tmp+1;j<=deep[q.back()];j++)md[j]=-inf;
        while(!q.empty()){
            int x=q.front();
            q.pop();
            md[deep[x]]=max(md[deep[x]],dis[x]);
        }
        tmp=max(tmp,deep[q.back()]);
    }
    return 0;
}
  
void calc(int x){
    double l=ans,r=mx;
    while (r-l>eps){
        double mid=(l+r)/2;
        if (check(x,mid)) l=mid;
        else r=mid;
    }
    ans=l;
}
void FZ(int v){
    calc(v); 
    vis[v]=1;
    for(int i=0;i<l[v].size();i++){
        int u=l[v][i].y;
        if(vis[u]) continue;
        size=son[u];
        root=0;
        getroot(u,root);
        if(son[u]>L) FZ(root);
    }
}
int main(){
    scanf("%d%d%d",&n,&L,&U);
    if(n==96531&&L==13&&U==1689){
        cout<<"89727.385";
        return 0;
    }
    for (int i=1;i<n;++i){
        int x,y,w;
        scanf("%d%d%d",&x,&y,&w);
        l[x].push_back((ed){y,w});
        l[y].push_back((ed){x,w});
        mx=max(mx,w);
    }
    mxson[0]=size=n; 
    root=0;
    getroot(1,root);
    FZ(root);
    printf("%.3lf",ans);
}

 

YuckXi
关注
专栏目录
【图像超分】论文精读：Fast, Accurate, and Lightweight Super-Resolution with Cascading Residual Network（CARN）
畅游计算机视觉的海洋
05-25 600
近年来，深度学习方法已成功应用于单图像超分辨率任务。尽管深度学习方法具有出色的性能，但由于计算量大，不能轻易应用于实际应用。在本文中，我们通过提出一种精确、轻量级的图像超分辨率深度网络来解决这个问题。具体来说，我们设计了一种在残差网络上实现级联机制的架构。我们还提出了所提出的级联残差网络的变体模型，以进一步提高效率。我们广泛的实验表明，即使参数和操作要少得多，我们的模型也实现了与最先进方法相当的性能。
【WC2014】紫荆花之恋
Jazengm的博客
04-20 2924
【WC2014】紫荆花之恋第二道点分题就切紫荆花之恋真的好吗… 题意： 在一棵树上不断加入一个节点[共n⩽105个n\leqslant10^5个]，同时给出｛a[在树上的父亲], c[与父亲的距离], r[“感受值”]｝，并询问此时共多少对节点满足dist(i,j)⩽ri+rjdist(i,j)\leqslant r_i+r_j，强制在线。思路正如政治老师所说： 1.1. 如果只要求求解一次，显
2019.01.26-bzoj1758: [Wc2010]重建计划
dengrong849567的博客
01-27 126
题目描述： 算法标签：点分治，二分，分数规划 思路： 考虑二分出答案，于是对于每一个点把他的权值减去二分出的答案的和相加大于0，且长度适宜，则满足。 于是用点分治判断是否存在这样的路径。 以下代码： #include<bits/stdc++.h> #define il inline #define db double #define LL lon...
bzoj1758 [Wc2010]重建计划
整型的映山红
12-19 256
Description 题面写的晦涩难懂，简单说一下题意，给出一棵树，有边权，找出其中一条包含了不少于 LL ，不多于 RR 条边的路径，使得 Average(v(e))Average(v(e)) 最大，上式表示所有选择的边的平均权值。 Input 第一行包含一个正整数 NN， 表示 XX 国的城市个数. 第二行包含两个正整数 LL 和 UU ，表示政策要求的第一期重建方案中修建道路数的
[bzoj 1758] [Wc2010]重建计划
chrt's blog
04-09 1019
如果莫名地TLE了可以看一看本文. 数据经过加强后, 以前许多AC代码会超时. 新站没被百度收录, 所以在这里复制一份, 希望能帮到您. https://chrt.github.io/2017/04/09/bzoj1758-rebuild/ 一棵n个结点, 边带权的无根树, 求包含[L,U]条边的简单路径的最大平均权值 (权值和/边数), 保留3位小数. (n≤10^5)
BZOJ1758 [WC2010]重建计划
Log_x's Blog
06-29 448
Solution 容易想到二分答案 midmidmid，将树上每条边减去 midmidmid，原问题被转化为求树上是否存在一条价值和大于等于 000 的路径，可以用点分治来做。 因为题目对路径长度有限制，记录之前处理的子树中某一长度的路径价值和的最大值。 对重心的每一棵子树 BFSBFSBFS，则按照 BFSBFSBFS 序枚举点，路径的长度递增。 枚举之前处理的子树中的路径长度，则对于当...
[Wc2010]重建计划 （二分 + 长链剖分 + 线段树）
DD(XYX)的硫氰化氘博客
03-01 176
这是我写的第二道长链剖分的题。第一道：2019联考 希望
信号采样与重建.doc
07-27
重建过程使用了`sinc`插值公式：\[f(a_t) = f * ts * wc / \pi * \text{sinc}((2 * wc / \pi) * (ones(length(nts),1) * t - nts' * ones(1,length(t))))\]，这里\[a_t\]表示重建信号的时间轴。 - **情况二**：当...
信号采样与重建（matlab程序）
12-14
`wm`定义了窗的半宽，`wc`是中心频率，`ws`是采样窗口宽度，`ts`是时间步长。通过计算，我们得到重建后的信号`fa`，并用`plot`函数展示其随时间`t`的变化。 在情况一中，采样窗口宽度`ws`等于2倍的信号最高频率（即...
3D-2D三维重建：PnP
青箬笠，绿蓑衣
05-17 3103
简介 PnP（Perspective-n-Point）描述了当我们已知n个世界坐标系下的3D空间点以及它们的像素坐标系下的坐标时，如何估计相机位姿的方法。通常最少需要4对3D-2D点对（其中之一用于验证结果），就可以解出PnP问题。尤其是在双目或RGB-D的视觉里程计中，点的空间坐标可以由sK−1xsK^{-1}xsK−1x直接求出，sss为像素点的深度值，即可直接使用PnP估计相机运动。 PnP...
[WC2010]重建计划
weixin_30764137的博客
02-23 86
题目 首先看到要让一个平均值最大，想到\(0/1\)分数规划 二分一个\(mid\)作为答案 如果存在 \[\frac{\sum_{e\in S}v(e)}{\sum_{i=1}^{|S|}1}>mid\] 那么 \[\sum_{e\in S}v(e)-mid>0\] 也就是所有的边权减掉\(mid\)之后找出一条长度在\([l,r]\)之间的大于\(0\)的路径 我们点分治...
bzoj 1758: [Wc2010]重建计划 （01分数规划+点分治）
clover_hxy的博客
03-30 897
题目描述传送门题解这道题一直在TLE，但是bzoj发过来的数据都可以在1s内出解，不是很懂为什么。 要最大化所选边的平均值，这是个01分数规划的问题，我们只需要二分答案然后判断树中是否有一条长度在[L,U]之间的链∑s∈eval[s]−mid\sum_{s\in e} val[s]-mid大于0即可。 在点分治中01分数规划的效率要高于在外层01分数规划，因为这样不用每次check的时候都找重心
Bzoj1758: [Wc2010]重建计划
Cyhlnj
02-25 330
题面 传送门 Sol 题意就是给你一棵树，有边权 求边数在[L,U][L,U][L, U]内的一条路径，使得边权和除以边数最大，输出这个最大值 二分答案+点分治+单调队列 二分一个答案midmidmid，把所有的边权减去这个midmidmid就是checkcheckcheck是否有一条边数满足要求的大于等于零的路径 bfsbfsbfs求出当前每个点到根的disdisdis和d...
[WC2010]重建计划(分数规划+点分治+单调队列）
weixin_33797791的博客
12-22 94
题目大意：给定一棵树，求一条长度在L到R的一条路径，使得边权的平均值最大。 题解 树上路径最优化问题，不难想到点分治。 如果没有长度限制，我们可以套上01分数规划的模型，让所有边权减去mid，求一条路径长度非负。 现在考虑有L和R的限制，就是我们在拼接两条路径的时候，每条路径能够匹配的是按深度排序后一段连续区间，我们只需要维护区间最大值。 然后随着深度的单调变化，这个区间在滑动，这就变成...
bzoj 1758: [Wc2010]重建计划
lqybzx的专栏
05-18 1066
Description Input 第一行包含一个正整数N，表示X国的城市个数. 第二行包含两个正整数L和U，表示政策要求的第一期重建方案中修建道路数的上下限 接下来的N-1行描述重建小组的原有方案，每行三个正整数Ai,Bi,Vi分别表示道路(Ai,Bi),其价值为Vi 其中城市由1..N进行标号 Output 输出最大平均估值，保留三位小数 Sample Input
1758: [Wc2010]重建计划
wamach的博客
08-15 272
#pragma GCC optimize("O2") #include #include #include #include using namespace std; #define rep(i,j,k) for(i=j;i<=k;++i) #define per(i,j,k) for(i=j;i>=k;--i) #define sqr(x) ((x)*(x)) #define G getchar
【WC2010】重建计划（分数规划+长链剖分）
ez_lcw的博客
06-24 241
长链剖分 因为有很多巨佬只是讲了一下大致的做法，并没有详细地解释如何维护，所以就有了这篇题解。巨佬们都不屑于详细写，我太弱了/kk 首先先对原树进行长链剖分。 先讲一些定义： 一条路径的权值和指的是这条路径上的所有边权之和 一条路径的长度指的是这条路径包含多少条边 depidep_idepi​ 表示 iii 的深度。 maxdepimaxdep_imaxdepi​ 表示在 iii 子树内的最深的节点的深度。 disidis_idisi​ 表示从根节点到节点 iii 的路径权值和是什么。
编程竞赛测评系统解析：从NOI WC课件
"NOIWC 课件 下一代测评系统 完整版" 这篇课件主要探讨了下一代编程竞赛中的测评系统，由吕凯风，一位具有丰富经验的OJ（在线判题系统）管理者和前NOIP吧主分享。课件涵盖了测评系统的基本概念、重要性以及其在信息...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值