NOIP2015-D2T3运输计划

最新推荐文章于 2019-11-04 21:06:38 发布
最新推荐文章于 2019-11-04 21:06:38 发布
阅读量90 收藏
点赞数
原文链接：http://www.cnblogs.com/w-h-h/p/7580666.html
版权

题目背景

公元 2044 年,人类进入了宇宙纪元。

题目描述

L 国有 n 个星球,还有 n-1 条双向航道,每条航道建立在两个星球之间,这 n-1 条航道连通了 L 国的所有星球。

小 P 掌管一家物流公司,该公司有很多个运输计划,每个运输计划形如:有一艘物

流飞船需要从 ui 号星球沿最快的宇航路径飞行到 vi 号星球去。显然,飞船驶过一条航道 是需要时间的,对于航道 j,任意飞船驶过它所花费的时间为 tj,并且任意两艘飞船之 间不会产生任何干扰。

为了鼓励科技创新,L 国国王同意小 P 的物流公司参与 L 国的航道建设,即允许小 P 把某一条航道改造成虫洞,飞船驶过虫洞不消耗时间。

在虫洞的建设完成前小 P 的物流公司就预接了 m 个运输计划。在虫洞建设完成后, 这 m 个运输计划会同时开始,所有飞船一起出发。当这 m 个运输计划都完成时,小 P 的 物流公司的阶段性工作就完成了。

如果小 P 可以自由选择将哪一条航道改造成虫洞,试求出小 P 的物流公司完成阶段 性工作所需要的最短时间是多少?

输入输出格式

输入格式：

输入文件名为 transport.in。

第一行包括两个正整数 n、m,表示 L 国中星球的数量及小 P 公司预接的运输计划的数量,星球从 1 到 n 编号。

接下来 n-1 行描述航道的建设情况,其中第 i 行包含三个整数 ai, bi 和 ti,表示第

i 条双向航道修建在 ai 与 bi 两个星球之间,任意飞船驶过它所花费的时间为 ti。

接下来 m 行描述运输计划的情况,其中第 j 行包含两个正整数 uj 和 vj,表示第 j个 运输计划是从 uj 号星球飞往 vj 号星球。

输出格式：

输出 共1行,包含1个整数,表示小P的物流公司完成阶段性工作所需要的最短时间。

输入输出样例

输入样例#1： 
6 3 
1 2 3 
1 6 4 
3 1 7 
4 3 6 
3 5 5 
3 6 
2 5 
4 5
输出样例#1： 
11

说明

所有测试数据的范围和特点如下表所示

请注意常数因子带来的程序效率上的影响。

这题首先用的是二分+暴力判断

二分答案然后把大于m的树链一一枚举，找可行的边取交集

80分

  1 #include<cstdio>
  2 #include<cstdlib>
  3 #include<algorithm>
  4 #include<cstring>
  5 #define MAXN 300005
  6 #define LOG 20
  7 #define pii pair<int,int>
  8 using namespace std;
  9 int first[MAXN],Next[MAXN*2],to[MAXN*2],W[MAXN*2],cnt;
 10 int n,m;
 11 int fa[LOG][MAXN],d[LOG][MAXN],dep[MAXN];
 12 int mid;
 13 struct Lian{
 14     int x,y,lca;
 15     int Val;
 16     friend bool operator < (const Lian &p1,const Lian &p2){
 17         return (p1.Val<p2.Val);
 18     }
 19     friend bool operator > (const Lian &p1,const Lian &p2){
 20         return !(p1.Val<p2.Val);
 21     }
 22 } s[MAXN];
 23 void Add(int x,int y,int w){
 24     Next[++cnt]=first[x]; first[x]=cnt; to[cnt]=y; W[cnt]=w;
 25     Next[++cnt]=first[y]; first[y]=cnt; to[cnt]=x; W[cnt]=w;
 26     // double edge
 27 }
 28 void dfs(int x){
 29     for(int e=first[x];e;e=Next[e]){
 30         int y=to[e],w=W[e];
 31         if(y==fa[0][x]){
 32             continue;
 33         }
 34         fa[0][y]=x;
 35         d[0][y]=w;
 36         dep[y]=dep[x]+1;
 37         dfs(y);
 38     }
 39 }
 40 pii LCA(int x,int y){
 41     if(dep[x]<dep[y]){
 42         swap(x,y);
 43     }
 44     int L=0;
 45     for(int k=dep[x]-dep[y],p=0;k;k>>=1,p++){
 46         if(k&1){
 47             L+=d[p][x];
 48             x=fa[p][x];
 49         }
 50     }
 51     if(x==y){
 52         return make_pair(L,x);
 53     }
 54     for(int k=LOG-1;k>=0;k--){
 55         if(fa[k][x]!=fa[k][y]){
 56             L+=d[k][x];
 57             L+=d[k][y];
 58             x=fa[k][x];
 59             y=fa[k][y];
 60         }
 61     }
 62     return make_pair(L+d[0][x]+d[0][y],fa[0][x]);
 63 }
 64 int b[MAXN];
 65 bool check(){
 66     memset(b,0,sizeof(b));
 67     Lian t;t.Val=mid;
 68     int Pos=upper_bound(s+1,s+m+1,t)-s;
 69     if(Pos>m){
 70         return 1;
 71     }
 72     int num=m-Pos+1;
 73     for(int i=Pos;i<=m;i++){
 74         int lc=s[i].x,rc=s[i].y;
 75         int lca=s[i].lca,L=s[i].Val;
 76         while(lc!=lca){
 77             if(L-d[0][lc]<=mid){
 78                 b[lc]++;
 79                 if(b[lc]>=num){
 80                     return 1;
 81                 }
 82             }
 83             lc=fa[0][lc];
 84         }
 85         while(rc!=lca){
 86             if(L-d[0][rc]<=mid){
 87                 b[rc]++;
 88                 if(b[rc]>=num){
 89                     return 1;
 90                 }
 91             }
 92             rc=fa[0][rc];
 93         }
 94     }
 95     return 0;
 96 }
 97 int main()
 98 {
 99 //    freopen("T1.in","r",stdin);
100 //    freopen("my.out","w",stdout);
101     scanf("%d%d",&n,&m);
102     for(int i=1;i<n;i++){
103         int x,y,w;
104         scanf("%d%d%d",&x,&y,&w);
105         Add(x,y,w);
106     }
107     dep[1]=1;
108     dfs(1);
109     for(int k=1;k<LOG;k++){
110         for(int i=1;i<=n;i++){
111             fa[k][i]=fa[k-1][fa[k-1][i]];
112             d[k][i]=d[k-1][i]+d[k-1][fa[k-1][i]];
113         }
114     }
115     for(int i=1;i<=m;i++){
116         scanf("%d%d",&s[i].x,&s[i].y);
117         pii t=LCA(s[i].x,s[i].y);
118         s[i].lca=t.second;
119         s[i].Val=t.first;
120     }
121     sort(s+1,s+m+1);
122 //    for(int i=1;i<=n;i++){
123 //        mid=s[i].Val;
124 //        if(check()){
125 //            printf("1\n");
126 //        }
127 //        else{
128 //            printf("0\n");
129 //        }
130 //    }
131     int L=1,R=s[m].Val;
132     while(L<R-1){
133         mid=(L+R)/2;
134         if(check()){
135             R=mid;
136         }
137         else{
138             L=mid;
139         }
140     }
141     mid=L; if(check()){
142         printf("%d\n",L);
143     }
144     else{
145         printf("%d\n",R);
146     }
147     return 0;
148 }
View Code

实际上判断用树上差分提高效率即可AC

  1 #include<cstdio>
  2 #include<cstdlib>
  3 #include<algorithm>
  4 #include<cstring>
  5 #define MAXN 300005
  6 #define LOG 20
  7 #define pii pair<int,int>
  8 using namespace std;
  9 int first[MAXN],Next[MAXN*2],to[MAXN*2],W[MAXN*2],cnt;
 10 int n,m;
 11 int fa[LOG][MAXN],d[LOG][MAXN],dep[MAXN];
 12 int mid;
 13 struct Lian{
 14     int x,y,lca;
 15     int Val;
 16     friend bool operator < (const Lian &p1,const Lian &p2){
 17         return (p1.Val<p2.Val);
 18     }
 19     friend bool operator > (const Lian &p1,const Lian &p2){
 20         return !(p1.Val<p2.Val);
 21     }
 22 } s[MAXN];
 23 void Add(int x,int y,int w){
 24     Next[++cnt]=first[x]; first[x]=cnt; to[cnt]=y; W[cnt]=w;
 25     Next[++cnt]=first[y]; first[y]=cnt; to[cnt]=x; W[cnt]=w;
 26     // double edge
 27 }
 28 void dfs(int x){
 29     for(int e=first[x];e;e=Next[e]){
 30         int y=to[e],w=W[e];
 31         if(y==fa[0][x]){
 32             continue;
 33         }
 34         fa[0][y]=x;
 35         d[0][y]=w;
 36         dep[y]=dep[x]+1;
 37         dfs(y);
 38     }
 39 }
 40 pii LCA(int x,int y){
 41     if(dep[x]<dep[y]){
 42         swap(x,y);
 43     }
 44     int L=0;
 45     for(int k=dep[x]-dep[y],p=0;k;k>>=1,p++){
 46         if(k&1){
 47             L+=d[p][x];
 48             x=fa[p][x];
 49         }
 50     }
 51     if(x==y){
 52         return make_pair(L,x);
 53     }
 54     for(int k=LOG-1;k>=0;k--){
 55         if(fa[k][x]!=fa[k][y]){
 56             L+=d[k][x];
 57             L+=d[k][y];
 58             x=fa[k][x];
 59             y=fa[k][y];
 60         }
 61     }
 62     return make_pair(L+d[0][x]+d[0][y],fa[0][x]);
 63 }
 64 int b[MAXN];
 65 int num;
 66 int find(int x){
 67     int ret=0;
 68     for(int e=first[x];e;e=Next[e]){
 69         int y=to[e];
 70         if(y==fa[0][x]) continue;
 71         ret=max(ret,find(y));
 72         b[x]+=b[y];
 73     }
 74     if(b[x]==num){
 75         ret=max(ret,d[0][x]);
 76     }
 77     return ret;
 78 }
 79 bool check(){
 80     memset(b,0,sizeof(b));
 81     Lian t;t.Val=mid;
 82     int Pos=upper_bound(s+1,s+m+1,t)-s;
 83     if(Pos>m){
 84         return 1;
 85     }
 86     num=m-Pos+1;
 87     for(int i=Pos;i<=m;i++){
 88         int lc=s[i].x,rc=s[i].y;
 89         int lca=s[i].lca,L=s[i].Val;
 90         b[lc]++,b[rc]++;
 91         b[lca]-=2;
 92     }
 93     int temp=find(1);
 94     return (s[m].Val-temp<=mid);
 95 }
 96 int main()
 97 {
 98  //   freopen("T1.in","r",stdin);
 99 //    freopen("my.out","w",stdout);
100     scanf("%d%d",&n,&m);
101     for(int i=1;i<n;i++){
102         int x,y,w;
103         scanf("%d%d%d",&x,&y,&w);
104         Add(x,y,w);
105     }
106     dep[1]=1;
107     dfs(1);
108     for(int k=1;k<LOG;k++){
109         for(int i=1;i<=n;i++){
110             fa[k][i]=fa[k-1][fa[k-1][i]];
111             d[k][i]=d[k-1][i]+d[k-1][fa[k-1][i]];
112         }
113     }
114     for(int i=1;i<=m;i++){
115         scanf("%d%d",&s[i].x,&s[i].y);
116         pii t=LCA(s[i].x,s[i].y);
117         s[i].lca=t.second;
118         s[i].Val=t.first;
119     }
120     sort(s+1,s+m+1);
121 //    for(int i=1;i<=n;i++){
122 //        mid=s[i].Val;
123 //        if(check()){
124 //            printf("1\n");
125 //        }
126 //        else{
127 //            printf("0\n");
128 //        }
129 //    }
130     int L=0,R=s[m].Val;
131     while(L<R-1){
132         mid=(L+R)/2;
133         if(check()){
134             R=mid;
135         }
136         else{
137             L=mid;
138         }
139     }
140     mid=L; if(check()){
141         printf("%d\n",L);
142     }
143     else{
144         printf("%d\n",R);
145     }
146     return 0;
147 }
View Code

 

 

 

 

转载于:https://www.cnblogs.com/w-h-h/p/7580666.html

weixin_30607659
关注
[二分+LCA+差分]【NOIP2015D2T3】运输计划 题解
Wflower的Blog //What are you waiting for?
10-25 438
传送门： 洛谷：https://daniu.luogu.org/problemnew/show/2680#sub （链接给的是大牛分站的，主站上T了一个点，大牛分站A了，“高性能”，卡常数） UOJ：http://uoj.ac/problem/150题目分析给出一棵nn个节点且有边权的数和m个从sis_i到tit_i的运输计划，完成所有计划所花费的时间为运输计划中的最大时间。现在可以选中树上的
#NOIP2015#运输计划（树剖 + 二分 + 差分）
ElaineZXY的博客
08-12 286
#150. 【NOIP2015】运输计划  统计  描述 提交 自定义测试 公元 2044 年，人类进入了宇宙纪元。 L 国有 nn 个星球，还有 n−1n−1 条双向航道，每条航道建立在两个星球之间，这 n−1n−1 条航道连通了 LL 国的所有星球。 小 P 掌管一家物流公司， 该公司有很多个运输计划，每个运输计划形如：有一艘物流飞船需要从 uiui 号星球沿最快的宇
NOIP2015 D2T3运输计划
啊
10-20 963
常数是检验人傻的唯一标准…… 但是今天我要推翻这句话。 我这么机智，怎么还有巨大常数……#include<iostream> #include<cstdio> #include<algorithm> #include<cstring> #define Max(x,y) x>y?x:y; #define maxn 300005 using namespace std; void read(int
NOIP2015D2T3-运输计划
LiuKairui的博客
07-31 514
问题描述 公元 2044 年，人类进入了宇宙纪元。 L 国有 nn 个星球，还有 n-1n−1 条双向航道，每条航道建立在两个星球之间，这 n-1n−1 条航道连通了 LL 国的所有星球。 小 P 掌管一家物流公司， 该公司有很多个运输计划，每个运输计划形如：有一艘物流飞船需要从 ui​号星球沿最快的宇航路径飞行到 vi号星球去。显然，飞船驶过一条航道是需要时间的，对于航道 jj ，任意飞船...
NOIP 2015 [D2 T3]运输计划
SuperGate的博客
09-09 2010
【NOIP2015 Day2】运输计划 Time Limit:20000MS  Memory Limit:262144K Total Submit:10 Accepted:7 Case Time Limit:1000MS Description         公园2044年，人类进入了宇宙纪元。         L过有n个星球，还有n-1条双向航道，每条航道建立在两个星球之间，这n-
【转载】NOIP2015 TG D2T3 transport
MEET YOU FIRST TIME
08-08 1219
本题是2015年提高组复赛第二天的最后一题，也就是说，是整体而言最难的题目。而其思路也证明了这一点。但是，在江老师一句“LCA”的启发后，我历尽了半节课，终于想了出来一个做法。【题意】 给定一棵无向的连通的一共有n个节点的树，并给出m个“运输计划”，每个运输计划要从一个点一路将货物运输到另一个点，费用为沿途所有边的权值之和。 现在，可以允许你将某一条边的权值设为0，问最小的总费用和是多少？【分析】
「NOIP 2015 D2T3」运输计划
forever_dreams的博客
11-04 186
算法标签：二分，树上差分，dfs序
【BZOJ 4326】【NOIP 2015 d2t3】运输计划
Ripped的博客
02-20 449
题意给你一棵树和边权，有若干条航线，从树上的一个点到另一个点。令运输时间为所有航线所需要的时间中的最长时间。你可以将某一条边的权值变为0，求这个最长时间的最小值。样例输入6 3 1 2 3 1 6 4 3 1 7 4 3 6 3 5 5 3 6 2 5 4 5样例输出11SOL先讲一下部分分毕竟这是一道联赛题。 1、首先想到就是预处理lca，然后枚举变为0的边计算每条航线的权值，时
NOIP2015 D2T3 运输计划
weixin_30856725的博客
05-21 121
拿到题目的第一眼 首先这是一棵n个节点的树（别说你看不出来） 然后对于树上的m条链我们可以选取树上的唯一一条边使它的边权变为0 求处理后最长链的长度 20分 m=1好啦，好像可做一眼望去全是水 只需求出一条链上的所有边并计算边权和及最大边权（暴力往上跳并记录即可）边权和减去最大边权即为答案那么我们就可以O(n)过掉这道题了（不嫌麻烦的话也可以O(log n)搞树上路径） 60分？ ...
NOIP2013D2T3-华容道
LiuKairui的博客
07-30 496
问题描述 小 B 最近迷上了华容道，可是他总是要花很长的时间才能完成一次。于是，他想到用编程来完成华容道：给定一种局面， 华容道是否根本就无法完成，如果能完成， 最少需要多少时间。 小 B 玩的华容道与经典的华容道游戏略有不同，游戏规则是这样的： 在一个 n×mn×mn \times m 棋盘上有n×mn×m n \times m 个格子，其中有且只有一个格子是空白的，其余 n×m−1n×m...
NOIP2012D2T3-疫情控制
LiuKairui的博客
10-17 1126
问题描述 H 国有 n 个城市，这 n 个城市用 n−1 条双向道路相互连通构成一棵树， 1 号城市是首都，也是树中的根节点。 H 国的首都爆发了一种危害性极高的传染病。当局为了控制疫情，不让疫情扩散到边境城市（叶子节点所表示的城市），决定动用军队在一些城市建立检查点，使得从首都到边境城市的每一条路径上都至少有一个检查点，边境城市也可以建立检查点。但特别要注意的是，首都是不能建立检查点的。 ...
【NOIP2018提高组D2T3】保卫王国
Luckfort
11-18 1593
倍增 换根
NOIP 2014 D2T3 解方程 Hash大法好
世界
11-10 3729
题目大意：给定高次方程an*x^n+...+a1*x^1+a0*x^0=0 求[1,m]区间内有多少个整数根 ai 懒得高精，考场上写的long double乱搞……30分打底50分顶天QAQ 当我终于搞定了各种非官方数据之后，我只能长跪大地，手捧鲜花，仰望上苍高喊：哈希大法好！ 首先阿贝尔在200年前告诉我们 五次以上方程没有求根公式 于是我们只能枚举1~m 这个是100W 然后100
2025年软考高级 - 信息系统项目管理师考试备考全攻略
11-09
2025年软考高级 - 信息系统项目管理师考试备考全攻略
MySQL 5.7从入门到精通 第23章 新闻发布系统数据库设计 共6页.pptx
最新发布
11-09
【课程大纲】 第1章 初始MySQL 共19页.pptx 第2章 MySQL的安装与配置 共14页.pptx 第3章 数据库的基本操作 共11页.pptx 第4章 数据表的基本操作 共26页.pptx 第5章 数据类型和运算符 共17页.pptx 第6章 MySQL函数 共76页.pptx 第7章 查询数据 共48页.pptx 第8章 插入、更新与删除数据 共10页.pptx 第9章 索引 共11页.pptx 第10章 存储过程和函数 共19页.pptx 第11章 视图 共20页.pptx 第12章 触发器 共11页.pptx 第13章 用户管理 共25页.pptx 第14章 数据备份与还原 共21页.pptx 第15章 MySQL日志 共22页.pptx 第16章 性能优化 共18页.pptx 第17章 MySQL Workbench5.2 的使用 共15页.pptx 第18章 MySQL Replication 共27页.pptx 第19章 MySQL Cluster 共49页.pptx 第20章 MySQL管理利器——MySQL Utilities 共5页.pptx 第21章 读写分离的利器——MySQL Proxy 共5页.pptx 第22章 PHP操作MySQL数据库 共7页.pptx 第23章 新闻发布系统数据库设计 共6页.pptx 第24章 论坛管理系统数据库设计 共6页.pptx
高分springboot毕设+vue的游戏创意工坊与推广平台的设计与实现_orv论文-Java源码.zip
11-09
本项目是一个基于Spring Boot和Vue的游戏创意工坊与推广平台的设计与实现。该项目旨在为游戏开发者和玩家提供一个集中的平台，使他们能够分享创意、展示作品并获取反馈。平台的核心功能包括游戏创意的提交与管理、游戏作品的展示与评价、用户间的互动交流以及推广活动的组织与管理。 在技术实现上，后端采用Spring Boot框架，利用其快速开发和部署的特点，确保系统的稳定性和高效性。前端则使用Vue.js，以其灵活的数据绑定和组件化开发方式，为用户提供流畅的交互体验。数据库设计充分考虑了数据的安全性和扩展性，以支持大量用户和作品的存储需求。 此外，项目还集成了多种实用工具和插件，如用户认证、权限管理、文件存储等，以提升平台的整体功能和用户体验。通过这个项目，用户不仅能够锻炼自己的编程技能，还能深入了解游戏开发和运营的全过程。
考研助手--论文.zip
11-09
基于SSM的毕业设计源码
C# AutoCAD DWG打印成PDF 代码
11-09
自己写的转PDF代码，调试，可用。
PowerQuery的使用与M语言介绍.docx
11-09
PowerQuery的使用与M语言介绍.docx
2023年NOIP CSP-J CSP-S 初赛学习资料大全
NOIP CSP-J CSP-S 初赛学习资料集（N）-2023.09.02.pdf 这份学习资料集主要针对 NOIP CSP-J CSP-S 初赛的学习和备考，涵盖了计算机结构与组成、计算机科学速成课、CSP 初赛高频考点解析、CSP-J 初赛练习、备战 CSP...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值