标签：return fnd int 货车运输 fa dfn 详细 include 揭秘

一个 \(\log\) 的树剖解法！！！

好吧不滥用标题行了 .

注意到没有修改，轻重链剖分，链用 st 表维护，时间复杂度 \(O(q\log n)\) .

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <ctime>
#include <vector>
#include <queue>
#include <cmath>
#include <map>
#include <set>
using namespace std;
inline int chkmin(int& a, const int& b){if (a > b) a = b; return a;}
const int N = 12345, M = 51411, LgN = 11+4.514, INF = 0x3f3f3f3f;
typedef pair<int, int> pii;
int n, m, w[N];
int fa[N], dep[N], siz[N], son[N], top[N], dfn[N], rnk[N], cc;
vector<pii> g[N];
inline void addedge(int u, int v, int w){g[u].emplace_back(make_pair(v, w));}
inline void ade(int u, int v, int w){addedge(u, v, w); addedge(v, u, w);}
struct Edge{int u, v, w; bool operator < (const Edge& _)const{return w > _.w;}}ed[M];
struct dsu
{
	int fa[N];
	int fnd(int x){return x == fa[x] ? x : fa[x] = fnd(fa[x]);}
	inline void merge(int u, int v){fa[fnd(u)] = fnd(v);}
	inline bool same(int u, int v){return fnd(u) == fnd(v);}
	void init(int n){for (int i=1; i<=n; i++) fa[i] = i;}
}D;
struct st
{
	int f[N][LgN];
	inline void init()
	{
		for (int i=1; i<=n; i++) f[i][0] = w[rnk[i]];
		for (int j=1; (1<<j)<=n; j++)
			for (int i=1; i+(1<<j)-1<=n; i++) f[i][j] = min(f[i][j-1], f[i+(1<<(j-1))][j-1]);
	}
	inline int query(int l, int r){int _ = log2(r-l+1); return min(f[l][_], f[r-(1<<_)+1][_]);}
}T;
void dfs1(int u)
{
	siz[u] = 1;
	for (auto e : g[u])
	{
		int v = e.first, _ = e.second;
		if (v == fa[u]) continue;
		fa[v] = u; dep[v] = dep[u] + 1;
		dfs1(v);
		siz[v] += siz[u]; w[v] = _;
		if (!son[u] || (siz[son[u]] < siz[v])) son[u] = v;
	}
}
void dfs2(int u, int t)
{
	top[u] = t;
	rnk[++cc] = u; dfn[u] = cc;
	if (!son[u]) return ;
	dfs2(son[u], t);
	for (auto e : g[u])
		if ((e.first != fa[u]) && (e.first != son[u])) dfs2(e.first, e.first);
}
int query(int u, int v)
{
	int ans=INF, fu=top[u], fv=top[v];
	while (fu != fv)
	{
		if (dep[fu] > dep[fv]){chkmin(ans, T.query(dfn[fu], dfn[u])); u = fa[fu];}
		else{chkmin(ans, T.query(dfn[fv], dfn[v])); v = fa[fv];}
		fu = top[u]; fv = top[v];
	}
	if (dep[u] > dep[v]) swap(u, v);
	if (u != v) chkmin(ans, T.query(dfn[u]+1, dfn[v]));
	return ans;
}
void Kruskal()
{
	sort(ed+1, ed+1+m); D.init(n);
	for (int i=1; i<=m; i++)
	{
		int u = ed[i].u, v = ed[i].v, w = ed[i].w;
		if (D.same(u, v)) continue;
		D.merge(u, v); ade(u, v, w);
	}
}
int main()
{
#ifndef ONLINE_JUDGE
	freopen("i.in", "r", stdin);
#endif
	scanf("%d%d", &n, &m);
	for (int i=1; i<=m; i++) scanf("%d%d%d", &ed[i].u, &ed[i].v, &ed[i].w);
	Kruskal(); dfs1(1); dfs2(1, 1); T.init();
	int q, u, v; scanf("%d", &q);
	while (q--)
	{
		scanf("%d%d", &u, &v);
		if (!D.same(u, v)) puts("-1");
		else printf("%d\n", query(u, v));
	}
	return 0;
}

标签：return,fnd,int,货车运输,fa,dfn,详细,include,揭秘
来源： https://www.cnblogs.com/CDOI-24374/p/15868964.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。