BZOJ4061/Gym100624F CERC2012 Farm and Factory 最短路、切比雪夫距离
设\(f_i\)表示\(i\)到\(1\)号点的最短距离,\(g_i\)表示\(i\)到\(2\)号点的最短距离,\(s_i\)表示\(n+1\)号点到\(i\)号点的最短距离,\(A=s_1,B=s_2\)
根据最短路三角形不等式,\(|f_i - A| \leq s_i \leq f_i + A , |g_i - B| \leq s_i \leq g_i + B\)
而\(s_i\)要取到最小值,所以\(s_i = \max\{|f_i - A| , |g_i - B|\}\)
所以我们要求的是\(\sum\limits_{i=1}^N \max\{|f_i - A| , |g_i - B|\}\),这相当于求一个动点\((A,B)\)到平面上\(N\)个点\((f_i,g_i)\)的最小切比雪夫距离和。
切比雪夫距离可以转为曼哈顿距离,将坐标\((x,y)\)变为\((\frac{x+y}{2} , \frac{x-y}{2})\),前者的切比雪夫距离等效于后者的曼哈顿距离。而曼哈顿距离可以直接拆开横纵坐标然后取中位数。
注意:我天真的以为2012年的题不会卡SPFA……
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<ctime>
#include<algorithm>
#include<cstring>
#include<iomanip>
#include<queue>
#define INF 0x3f3f3f3f
//This code is written by Itst
using namespace std;
inline int read(){
int a = 0;
char c = getchar();
while(!isdigit(c) && c != EOF)
c = getchar();
while(isdigit(c)){
a = a * 10 + c - 48;
c = getchar();
}
return a;
}
#define PLI pair < long long , int >
#define st first
#define nd second
const int MAXN = 1e5 + 7;
struct Edge{
int end , upEd , w;
}Ed[MAXN * 6];
int head[MAXN] , N , M , cntEd;
long long dis[2][MAXN];
priority_queue < PLI > q;
inline void addEd(int a , int b , int w){
Ed[++cntEd].end = b;
Ed[cntEd].w = w;
Ed[cntEd].upEd = head[a];
head[a] = cntEd;
}
void SPFA(int ind){
memset(dis[ind] , 0x3f , sizeof(long long) * (N + 1));
dis[ind][ind + 1] = 0;
q.push(PLI(0 , ind + 1));
while(!q.empty()){
PLI t = q.top();
q.pop();
if(-t.st != dis[ind][t.nd]) continue;
for(int i = head[t.nd] ; i ; i = Ed[i].upEd)
if(dis[ind][Ed[i].end] > dis[ind][t.nd] + Ed[i].w){
dis[ind][Ed[i].end] = dis[ind][t.nd] + Ed[i].w;
q.push(PLI(-dis[ind][Ed[i].end] , Ed[i].end));
}
}
}
inline long long abss(long long x){return x < 0 ? -x : x;}
void out(long long a , int b){
cout << a / b << '.';
a %= b;
for(int i = 1 ; i <= 8 ; ++i){
a *= 10;
cout << a / b;
a %= b;
}
putchar('\n');
}
int main(){
vector < long long > x , y;
for(int T = read() ; T ; --T){
N = read(); M = read();
memset(head , 0 , sizeof(int) * (N + 1));
cntEd = 0;
for(int i = 1 ; i <= M ; ++i){
int a = read() , b = read() , c = read();
addEd(a , b , c); addEd(b , a , c);
}
SPFA(0); SPFA(1);
x.clear(); y.clear();
long long sum = 0;
for(int i = 1 ; i <= N ; ++i){
x.push_back(dis[0][i] - dis[1][i]);
y.push_back(dis[0][i] + dis[1][i]);
}
sort(x.begin() , x.end()); sort(y.begin() , y.end());
long long mid = x[N >> 1];
for(int i = 0 ; i < N ; ++i)
sum += abss(x[i] - mid);
mid = y[N >> 1];
for(int i = 0 ; i < N ; ++i)
sum += abss(y[i] - mid);
out(sum , 2 * N);
cerr << N << ' ' << sum << endl;
}
return 0;
}
BZOJ4061/Gym100624F CERC2012 Farm and Factory 最短路、切比雪夫距离的更多相关文章
- 【BZOJ4061】[Cerc2012]Farm and factory(最短路,构造)
[BZOJ4061][Cerc2012]Farm and factory(最短路,构造) 题面 BZOJ 然而权限题QwQ. 题解 先求出所有点到达\(1,2\)的最短路,不妨记为\(d_{u,1}, ...
- [BZOJ4061][Cerc2012]Farm and factory
bzoj 鉴于是权限题,放一下题面. Description 向Byteland的国王Bitolomew致敬!国王Bitolomew认为Byteland是一个独一无二的国家.它太小了,它所有的市民(包 ...
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
- 题解-bzoj4061 CERC-2012Farm and Factory
Problem Please contact lydsy2012@163.com! 题意概要:给定\(n\)点\(m\)边无向图,设定两个起点为\(1,2\),现要求在图中增加一个点,并将这个点与其他 ...
- sgu 185 最短路建网络流
题目:给出一个图,从图中找出两条最短路,使得边不重复. 分析:既然是最短路,那么,两条路径上的所有节点的入边(s,x).出边(x,e)必定是最优的,即 dis[x] = dis[s]+edge_dis ...
- poj 2449 Remmarguts' Date(第K短路问题 Dijkstra+A*)
http://poj.org/problem?id=2449 Remmarguts' Date Time Limit: 4000MS Memory Limit: 65536K Total Subm ...
- Cogs 309. [USACO 3.2] 香甜的黄油 dijkstra,堆,最短路,floyd
题目:http://cojs.tk/cogs/problem/problem.php?pid=309 309. [USACO 3.2] 香甜的黄油 ★★ 输入文件:butter.in 输出文件 ...
- 多源最短路Floyd 算法————matlab实现
弗洛伊德(Floyd)算法是一种用于寻找给定的加权图中顶点间最短路径的算法.该算法名称以创始人之一.1978年图灵奖获得者.斯坦福大学计算机科学系教授罗伯特·弗洛伊德命名. 基本思想 通过Floyd计 ...
- 单源最短路Dijkstra算法——matlab实现
迪杰斯特拉(Dijkstra)算法是典型最短路径算法,用于计算一个节点到其他节点的最短路径. 它的主要特点是以起始点为中心向外层层扩展(广度优先搜索思想),直到扩展到终点为止. 基本思想 通过Dijk ...
随机推荐
- 如何将web项目部署到weblogic
在Eclipse中配置weblogic11g服务器: 下载并安装Eclipse:www.eclipse.org 下载并安装Weblogic Server Plugin for Eclipse:http ...
- fetch数据请求的封装
export default class HttpUtils { static get(url){ return new Promise((resolve,reject)=>{ fetch(ur ...
- 杨学明老师推出全新课程--《敏捷开发&IPD和敏捷开发结合的实践》
课时:13小时(2天) 敏捷开发&IPD和敏捷开发结合的实践 讲 师:杨学明 [课程背景] 集成产品开发(IPD).集成能力成熟度模型(CMMI).敏捷开发(Agile Developmen ...
- Nginx 出现 403 Forbidden 最终解决
Nginx 出现 403 Forbidden 最终解决 步骤一: 检查目录权限.权限不足的就加个权限吧. 例子:chmod -R 755 / var/www 步骤二: 打开nginx.conf 例子: ...
- Git&GitHub
Git是一个开源的分布式版本控制系统,可以有效.高速地处理从很小到非常大的项目版本管理. Git 是 Linus Torvalds 为了帮助管理 Linux 内核开发而开发的一个开放源码的版本控制软 ...
- ASP.NET -- WebForm -- 页面生命周期
ASP.NET -- WebForm -- 页面生命周期 ASP.NET 页运行时,此页将经历一个生命周期,在生命周期中将执行一系列处理步骤.这些步骤包括初始化.实例化控件.还原和维护状态.运行事件 ...
- Java多线程(五)线程的生命周期
点我跳过黑哥的卑鄙广告行为,进入正文. Java多线程系列更新中~ 正式篇: Java多线程(一) 什么是线程 Java多线程(二)关于多线程的CPU密集型和IO密集型这件事 Java多线程(三)如何 ...
- Fluentd初探 简介与安装
Fluentd是一个开源的数据收集器,专为处理数据流设计,有点像 syslogd ,但是使用JSON作为数据格式.它采用了插件式的架构,具有高可扩展性高可用性,同时还实现了高可靠的信息转发. 据分(Y ...
- JSP内置对象概述
JSP内置对象预先定义了九个这个的对象: request(请求) . response (响应). session (会话). application (应用程序). out . pageContex ...
- Zookeeper的一致性
转载:http://flyfoxs.iteye.com/blog/2121560 下面内容主要摘抄于<<Hadoop实战>>,红色高亮部分是本人添加的白话注释. Zookeep ...