hdu4694

就是求有向图的支配树啦。。。题面就不简述了传送门

上面是非常详细的论文。。。补点笔记吧

定义一个点的半支配点,是一个点的dfs祖先中,能通过四种dfs树边连到该点,且dfn最小的点

而一个点的支配点,一定在该点到其半支配点路径上的点的支配点或是半支配点

用带权并查集,按dfn逆序处理就是了= =

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<vector>
#include<queue>
#include<set>
#include<map>
#include<stack>
#include<bitset>
#include<ext/pb_ds/priority_queue.hpp>
using namespace std;

const int maxn = 5E4 + 50;
typedef long long LL;

int n,m,dfs_clock,dfn[maxn],best[maxn],semi[maxn],idom[maxn],id[maxn],fa[maxn];
LL Ans[maxn];

vector <int> v[maxn],G[maxn],h[maxn],pre[maxn],son[maxn];

int Min(int x,int y) {return dfn[semi[x]] < dfn[semi[y]]?x:y;}

int Getfa(int k)
{
	if (k == fa[k]) return k;
	int ret = Getfa(fa[k]);
	best[k] = Min(best[fa[k]],best[k]);
	return fa[k] = ret;
}

void Dfs(int x)
{
	dfn[x] = ++dfs_clock; id[dfs_clock] = x;
	for (int i = 0; i < v[x].size(); i++)
	{
		int to = v[x][i];
		if (dfn[to]) continue;
		Dfs(to); son[x].push_back(to);
	}
}

void tarjan()
{
	for (int i = dfs_clock; i; i--)
	{
		int k = id[i];
		for (int j = 0; j < pre[k].size(); j++)
		{
			int pr = pre[k][j]; if (!dfn[pr]) continue;
			if (dfn[pr] < dfn[k])
			{
				if (dfn[pr] < dfn[semi[k]]) semi[k] = pr;
			}
			else
			{
				Getfa(pr); semi[k] = semi[Min(best[pr],k)];
			}
		}
		if (k != semi[k]) G[semi[k]].push_back(k);
		for (int j = 0; j < G[k].size(); j++)
		{
			int to = G[k][j]; Getfa(to);
			if (semi[best[to]] == k) idom[to] = k;
			else idom[to] = best[to];
		}
		for (int j = 0; j < son[k].size(); j++) fa[son[k][j]] = k;
	}
	for (int i = 2; i <= dfs_clock; i++)
	{
		int k = id[i];
		if (idom[k] != semi[k]) idom[k] = idom[idom[k]];
		if (k != idom[k]) h[idom[k]].push_back(k);
	}
}

void Dfs2(int x)
{
	Ans[x] += 1LL*x;
	for (int i = 0; i < h[x].size(); i++)
	{
		int to = h[x][i];
		Ans[to] += Ans[x]; Dfs2(to);
	}
}

char s[20];
void Print(LL x)
{
	if (!x) {putchar('0'); return;}
	int len = 0;
	for (; x; x /= 10LL) s[++len] = x % 10LL;
	while (len) putchar(s[len--] + '0');
}

int getint()
{
	char ch = getchar(); int ret = 0;
	while (ch < '0' || '9' < ch) 
	{
		if (ch == EOF) return -1;
		ch = getchar();
	}
	while ('0' <= ch && ch <= '9')
		ret = ret*10 + ch - '0',ch = getchar();
	return ret;
}

void Clear()
{
	for (int i = 0; i <= n; i++)
	{
		v[i].clear(); G[i].clear();
		h[i].clear(); son[i].clear();
		dfs_clock = dfn[i] = Ans[i] = 0; pre[i].clear();
		idom[i] = semi[i] = best[i] = fa[i] = i;
	}
}

int main()
{
	#ifdef DMC
		freopen("DMC.txt","r",stdin);
	#endif
	
	for (;;)
	{
		n = getint(); if (n == -1) break; 
		m = getint(); Clear();
		while (m--)
		{
			int x = getint(),y = getint();
			v[x].push_back(y);
			pre[y].push_back(x);
		}
		Dfs(n); tarjan(); Dfs2(n);
		for (int i = 1; i < n; i++) Print(Ans[i]),putchar(' ');
		Print(Ans[n]); puts("");
	}
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值