BZOJ4016【最短路】【点分治】

原创 2015年07月08日 18:31:02
/* I will wait for you */

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#include <ctime>
#include <algorithm>
#include <iostream>
#include <fstream>
#include <vector>
#include <queue>
#include <deque>
#include <set>
#include <map>
#include <string>
#define make(a,b) make_pair(a,b)
#define fi first
#define se second

using namespace std;

typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> pii;
typedef map<int, int> mii;

const int maxn = 100010;
const int maxm = 1010;
const int maxs = 26;
const int inf = 0x3f3f3f3f;
const int P = 1000000007;
const double error = 1e-9;

inline ll read()
{
	ll x = 0, f = 1; char ch = getchar();
	while (ch > '9' || ch < '0' )
		f = (ch == '-' ? -1 : 1), ch = getchar();
	while (ch <= '9' && ch >= '0')
		x = x * 10 + ch - '0', ch = getchar();
	return f * x;
}

struct edge
{
	int v, w, next;
} e[maxn * 2];

int n, m, cnt, sum, root, head[maxn], deep[maxn], del[maxn],
    size[maxn], val[maxn], K, _ans[2], g[maxn][2], 
    f[maxn][2], dis[maxn], vis[maxn], _max[maxn];

priority_queue<pii, vector<pii>, greater<pii> > q, _next[maxn];

void insert(int u, int v, int w)
{
	e[cnt] = (edge) {v, w, head[u]}, head[u] = cnt++;
}

void dijkstra()
{
	memset(dis, 0x3f, sizeof dis), dis[1] = 0;

	for (int i = 1; i <= n; i++)
		q.push(make(dis[i], i));

	for (int i = 1; i <= n; i++) {
		while (q.top().fi != dis[q.top().se])
			q.pop();
		int u = q.top().se; q.pop();

		for (int i = head[u]; ~i; i = e[i].next) {
			int v = e[i].v;
			if (dis[v] > dis[u] + e[i].w) {
				dis[v] = dis[u] + e[i].w;
				q.push(make(dis[v], v));
			}
		}
	}
}

void _dfs(int u)
{
	vis[u] = 1;
	while (!_next[u].empty()) {
		pii v = _next[u].top(); _next[u].pop();
		if (!vis[v.fi]) {
			insert(u, v.fi, v.se);
		       	insert(v.fi, u, v.se);
			_dfs(v.fi);
		}
	}
}

void find_root(int u, int p)
{
	size[u] = 1, _max[u] = 0;
	for (int i = head[u]; ~i; i = e[i].next) {
		int v = e[i].v;
		if (v != p && !del[v]) {
			find_root(v, u), size[u] += size[v];
			_max[u] = max(_max[u], size[v]);
		}
	}
	_max[u] = max(_max[u], sum - size[u]);
	if (_max[u] < _max[root]) root = u;
}

void __dfs(int u, int p)
{
	if (val[u] > g[deep[u]][0])
		g[deep[u]][0] = val[u], g[deep[u]][1] = 0;
	if (val[u] == g[deep[u]][0])
		g[deep[u]][1] += 1;

	for (int i = head[u]; ~i; i = e[i].next) {
		int v = e[i].v;
		if (v != p && !del[v]) {
			deep[v] = deep[u] + 1;
			val[v] = val[u] + e[i].w;
			__dfs(v, u);
		}
	}
}

void get_size(int u, int p)
{
	size[u] = 1;
	for (int i = head[u]; ~i; i = e[i].next) {
		int v = e[i].v;
		if (v != p && !del[v])
			get_size(v, u), size[u] += size[v];
	}
}

void solve(int u)
{
	del[u] = 1; get_size(u, 0);

	for (int i = 0; i <= size[u]; i++)
		f[i][0] = f[i][1] = 0;
	f[0][0] = 0, f[0][1] = 1;

	for (int i = head[u]; ~i; i = e[i].next) {
		int v = e[i].v;
		if (!del[v]) {
			for (int j = 0; j <= size[v]; j++)
				g[j][0] = g[j][1] = 0;
			deep[v] = 1, val[v] = e[i].w, __dfs(v, u);

			for (int j = 1; j <= K && j <= size[v]; j++) {
				if (g[j][0] + f[K - j][0] > _ans[0])
					_ans[0] = g[j][0] + f[K - j][0], _ans[1] = 0;
				if (g[j][0] + f[K - j][0] == _ans[0])
					_ans[1] += g[j][1] * f[K - j][1];
			}
			
			for (int j = 1; j <= size[v]; j++) {
				if (g[j][0] > f[j][0])
					f[j][0] = g[j][0], f[j][1] = 0;
				if (g[j][0] == f[j][0])
					f[j][1] += g[j][1];
			}
		}
	}

	for (int i = head[u]; ~i; i = e[i].next) {
		int v = e[i].v;
		if (!del[v]) {
			sum = size[u], root = 0;
			find_root(v, u), solve(root);
		}
	}
}	

int main()
{
	n = read(), m = read(), K = read() - 1;

	memset(head, -1, sizeof head);	
	for (int i = 1; i <= m; i++) {
		int u =read(), v = read(), w = read();
		insert(u, v, w), insert(v, u, w);
	}

	dijkstra();
	
	for (int i = 0; i < cnt; i++)
		if (dis[e[i].v] == dis[e[i ^ 1].v] + e[i].w)
			_next[e[i ^ 1].v].push(make(e[i].v, e[i].w));
	memset(head, -1, sizeof head);
	cnt = 0, _dfs(1);
	
	sum = _max[0] = n, root = 0;	
	find_root(1, 0), solve(root);

	printf("%d %d\n", _ans[0], _ans[1]);
	

	return 0;
}

【BZOJ4016】最短路径树问题,最短路+点分治

聪哥说要去有宫保鸡丁的地方(???)
  • xym_CSDN
  • xym_CSDN
  • 2016年08月21日 17:36
  • 509

【jzoj5088】【GDOI2017第四轮模拟day2】【最小边权和 】【图论】

题目大意有一张n个点m条边的有向图,每条边有一个互不相同的边权w,有q个询问,要求你从点a经过不超过c条边到点b,要求经过的边权递增并和尽量小,求出最小的边权和,如果没有合法方案则输出-1。解题思路设...
  • chunkitlau
  • chunkitlau
  • 2017年04月25日 15:07
  • 164

【点分治】的学习笔记和众多例题

包含点分治多道例题: 【POJ 1741 & BZOJ 1468 & BZOJ 3365】 【BZOJ 2152】 【BZOJ 2599】 【BZOJ 1316】 很详细的一篇文章,20分钟教会你点分...
  • nixinyis
  • nixinyis
  • 2017年03月23日 16:45
  • 1204

点分治时间复杂度

前言: 半年前学习的时候没有管这个问题,现在才搞懂。 结论: 点分治的时间复杂度为O(nlogn)O(nlogn)。 大致证明: 由于每次都是找数的重心,所以处理完一个大小为nn的树后,它的...
  • baidu_36797646
  • baidu_36797646
  • 2018年01月16日 15:47
  • 54

BZOJ4016 最短路径树问题(点分治)

题目链接:BZOJ4016题解:先跑一个最短路图,然后按照结点编号从小到大dfs一遍,dfs树即为题目所要求。然后就是点分治的经典做法,路径分为经过根节点的和不经过根节点的,不经过根节点的路径一定属于...
  • XHRlyb
  • XHRlyb
  • 2017年08月24日 16:09
  • 131

【点分治总结】

点分治教程:例题 给定一棵带权树,显然共有N*(N-1)/2条边,问:第k小的边边长多长? N...
  • ALPS233
  • ALPS233
  • 2016年05月13日 21:22
  • 2836

【蒟蒻的点分治专题训练】----5道题题解

点分治
  • DERITt
  • DERITt
  • 2016年05月13日 22:04
  • 638

洛谷P3806 【模板】点分治1 【点分治】

题目背景 感谢hzwer的点分治互测。 题目描述 给定一棵有n个点的树 询问树上距离为k的点对是否存在。 输入输出格式 输入格式: n,m 接下来n-...
  • qq_38678604
  • qq_38678604
  • 2017年12月17日 10:24
  • 89

bzoj2006 [ NOI2010 ] && bzoj3784 --点分治+线段树+堆

bzoj2006: 定义一个四元组{x,l,r,w},表示左端点在x,右端点在[l,r]的超级和弦的最大美妙度在将w作为右端点时取到,w可以用前缀和+线段树/ST表求出。 对于每个i,我们将{i,i+...
  • gjghfd
  • gjghfd
  • 2017年05月26日 16:42
  • 134

树分治(点分治模板)poj-1741 Tree

首先讲解一下树分治,以下的内容转自:http://blog.sina.com.cn/s/blog_6d5aa19a0100o73m.html 对于一棵有根树, 树中满足要求的一个数对所对应的一条...
  • u010660276
  • u010660276
  • 2015年04月07日 15:22
  • 5198
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:BZOJ4016【最短路】【点分治】
举报原因:
原因补充:

(最多只允许输入30个字)