一开始觉得是网络流..仔细一看应该是最短路,再看数据范围..呵呵不会写...这道题是最大生成树+最近公共祖先。第一次写..表示各种乱..

因为要求运输货物质量最大,所以路径一定是在最大生成树上的。然后就用LCA求两点之间的能运输的最大重量。预处理O(nlogn),查询O(logn).

----------------------------------------------------------------------------------------

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<vector>
#include<cctype>
#define rep(i,n) for(int i=0;i<n;i++)
#define clr(x,c) memset(x,c,sizeof(x))
#define Rep(i,l,r) for(int i=l;i<r;i++)
#define addEdge(u,v,d) MST.edges.push_back((Edge){u,v,d})
#define jud(a,b) MST.find(a)-MST.find(b)
using namespace std;
const int maxn=10000+5,maxs=20;
const int inf=0x7fffffff;
struct Edge {
int u,v,d;
Edge(int _u,int _v,int _d):u(_u),v(_v),d(_d) {}
bool operator < (const Edge &x) const {
return d>x.d;
}
};
struct LCA {
int p[maxn][maxs];
int depth[maxn];
int d[maxn][maxs];
bool vis[maxn];
int n;
vector<int> g[maxn];
vector<Edge> edges;
void init(int _n) {
n=_n;
edges.clear();
rep(i,n) { g[i].clear(); vis[i]=false; }
}
void add(int u,int v,int d) {
edges.push_back( (Edge) {u,v,d} );
edges.push_back( (Edge) {v,u,d} );
int m=edges.size();
g[u].push_back(m-2);
g[v].push_back(m-1);
}
void dfs(int x) {
int t=1;
vis[x]=1;
while(depth[x]>=(1<<t)) {
p[x][t]=p[p[x][t-1]][t-1];
d[x][t]=min(d[x][t-1],d[p[x][t-1]][t-1]);
t++;
}
rep(i,g[x].size()) {
Edge &e=edges[g[x][i]];
if(vis[e.v]) continue;
p[e.v][0]=x;
d[e.v][0]=e.d;
depth[e.v]=depth[x]+1;
dfs(e.v);
}
}
void DFS() { rep(i,n) if(!vis[i]) { depth[i]=0; dfs(i); } }
int query(int a,int b) {
int tmp,log=1;
if(depth[a]<depth[b]) swap(a,b);
while((1<<(log+1))<=depth[a]) log++;
int ans=inf;
for(int i=log;i>=0;--i) if(depth[a]-(1<<i)>=depth[b]) {
ans=min(ans,d[a][i]);
a=p[a][i];
}
if(a==b) return ans;
for(int i=log;i>=0;--i) if(p[a][i]!=-1 && p[a][i]!=p[b][i]) {
ans=min(ans,d[a][i]); a=p[a][i];
ans=min(ans,d[b][i]); b=p[b][i];
}
return ans=min(ans,min(d[a][0],d[b][0]));
}
};
LCA lca;
  
struct KRUSKAL {
int n;
int p[maxn];
vector<Edge> edges;
void init(int _n) {
n=_n;
edges.clear();
}
int find(int x) { return x==p[x] ? x : p[x]=find(p[x]); }
void kruskal() {
rep(i,n) p[i]=i;
sort(edges.begin(),edges.end());
rep(i,edges.size()) {
Edge &e=edges[i];
int x=find(e.u),y=find(e.v);
if(x!=y) {
p[x]=y;
   lca.add(e.u,e.v,e.d);
}
}
}
};
  
KRUSKAL MST;
int read() {
char c=getchar();
int ans=0,f=1;
while(!isdigit(c)) { if(c=='-') f=-1; c=getchar(); }
while(isdigit(c)) { (ans*=10)+=c-'0'; c=getchar(); }
return f*ans;
}
   
int main()
{
freopen("test.in","r",stdin);
freopen("test.out","w",stdout);
int n=read(),m=read();
MST.init(n);
lca.init(n);
rep(i,m) {
int u=read(),v=read(),d=read();
addEdge(--u,--v,d);
}
MST.kruskal();
lca.DFS();
n=read();
rep(i,n) {
int a=read(),b=read();
--a; --b;
if(jud(a,b)) printf("-1\n");
else {
printf("%d\n",lca.query(a,b));
}
}
return 0;
}

----------------------------------------------------------------------------------------

3287 货车运输

2013年NOIP全国联赛提高组

 时间限制: 1 s
 空间限制: 128000 KB
 题目等级 : 钻石 Diamond

[NOIP2013提高组] CODEVS 3287 火车运输(MST+LCA)的更多相关文章

  1. JZOJ 3534. 【NOIP2013提高组day1】货车运输

    Description A 国有 n 座城市,编号从 1 到 n,城市之间有 m 条双向道路.每一条道路对车辆都有重量限制,简称限重.现在有 q 辆货车在运输货物,司机们想知道每辆车在不超过车辆限重的 ...

  2. [NOIp2013提高组]积木大赛/[NOIp2018提高组]铺设道路

    [NOIp2013提高组]积木大赛/[NOIp2018提高组]铺设道路 题目大意: 对于长度为\(n(n\le10^5)\)的非负数列\(A\),每次可以选取一个区间\(-1\).问将数列清零至少需要 ...

  3. 1154 能量项链 2006年NOIP全国联赛提高组 codevs

    1154 能量项链  2006年NOIP全国联赛提高组 codevs 题目描述 Description 在Mars星球上,每个Mars人都随身佩带着一串能量项链.在项链上有N颗能量珠.能量珠是一颗有头 ...

  4. [NOIP2013 提高组] 华容道 P1979 洛谷

    [NOIP2013 提高组] 华容道 P1979 洛谷 强烈推荐,更好的阅读体验 经典题目:spfa+bfs+转化 题目大意: 给出一个01网格图,和点坐标x,y空格坐标a,b,目标位置tx,ty要求 ...

  5. Codevs 3287 货车运输 2013年NOIP全国联赛提高组(带权LCA+并查集+最大生成树)

    3287 货车运输 2013年NOIP全国联赛提高组 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 传送门 题目描述 Description A 国有 n 座 ...

  6. Codevs 3287 货车运输 == 洛谷P1967

    3287 货车运输 2013年NOIP全国联赛提高组 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 钻石 Diamond 题目描述 Description A 国有 n 座城市,编 ...

  7. cogs 2109. [NOIP 2015] 运输计划 提高组Day2T3 树链剖分求LCA 二分答案 差分

    2109. [NOIP 2015] 运输计划 ★★★☆   输入文件:transport.in   输出文件:transport.out   简单对比时间限制:3 s   内存限制:256 MB [题 ...

  8. [NOIP2012提高组] CODEVS 1200 同余方程(扩展欧几里德算法)

    数论题..所有数论对我来说都很恶心..不想再说什么了.. ------------------------------------------------ #include<iostream&g ...

  9. NOIP2010提高组] CODEVS 1069 关押罪犯(并查集)

    这道这么简单的题目还写了这么久.. 将每个会发生冲突的两人的怒气进行排序,然后从怒气大到小,将两个人放到不同监狱中.假如两人都已经被放置且在同一监狱,这就是答案. ------------------ ...

随机推荐

  1. C语言入门(11)——switch分支语句

    C语言提供了一种用于多分支选择的switch语句, 其一般形式为: switch(表达式) { case 常量表达式1:语句1; break; case 常量表达式2:语句2; break; .... ...

  2. 2014第16周三CSS布局再学习摘录

    今天尝试写了下前端页面,费了不少时间,做出的结果仍然惨不忍睹,感觉很简单的几个页面,在现有框架多个样式混杂下就是感觉很不自在随意,晚上回来又看了些div+css方面的基础知识. 1.CSS的class ...

  3. Make Hadoop 1.2.1 run, my first try

    经历两天努力,8月25日下午2点40分,终于让hadoop1.2.1跑起来. 用的是<Hadoop实战第2版>(陆嘉恒)里面的WordCount例子,虽然书是2013年出的,但用的例子还是 ...

  4. 2014 Web开发趋势

    本文翻译自:http://www.pixelstech.net/article/1401629232-Web-design-trends-for-2014 如今,已然到了Web横行的时代.越来越多的资 ...

  5. asp.net MVC Razor 语法(1)

    Razor 不是编程语言.它是服务器端标记语言. 什么是 Razor ? Razor 是一种允许您向网页中嵌入基于服务器的代码(Visual Basic 和 C#)的标记语法. 当网页被写入浏览器时, ...

  6. DLL导出函数和类 之 __declspec(dllexport)

    可利用__declspec(dllexport)导出函数或类. 若要指定C类型约定导出,则需在前面加extern “C”. 若要导出函数,__declspec(dllexport) 关键字必须出现在调 ...

  7. tigervnc*

    yum install -y tigervnc* tigervnc-server vncserver www.webmin.cn http://blog.sina.com.cn/s/blog_4ba5 ...

  8. matlab GUI之常用对话框(二)-- 进度条的使用方法

    常用对话框(二) 进度条   waitbar 调用格式: h = waitbar(x,'message')  waitbar(x,'message','CreateCancelBtn','button ...

  9. Oralce Exp 与 Imp 的使用方法

    1.完全:EXP  SYSTEM/SYSTEM@ORCL  FILE=C:\FULL.DMP  LOG=C:\FULL.DMP.LOG  FULL=Y  BUFFER=819200如果要执行完全导出, ...

  10. Java DecimalFormat 格式化数字

    我们经常要将数字进行格式化,比如取2位小数,这是最常见的.Java 提供 DecimalFormat 类,帮你用最快的速度将数字格式化为你需要的样子.下面是一个例子: importjava.text. ...