BZOJ 1086: [SCOI2005]王室联邦树分块

最新推荐文章于 2020-09-17 19:00:41 发布

BlackJack_

最新推荐文章于 2020-09-17 19:00:41 发布

阅读量518

点赞数

分类专栏：块状树 —————————树

—————————树同时被 2 个专栏收录

53 篇文章 0 订阅

订阅专栏

块状树

2 篇文章 0 订阅

订阅专栏

1086: [SCOI2005]王室联邦

Time Limit: 10 Sec Memory Limit: 162 MBSec Special Judge
Submit: 1975 Solved: 1211
[Submit][Status][Discuss]

Description

　　“余”人国的国王想重新编制他的国家。他想把他的国家划分成若干个省，每个省都由他们王室联邦的一个成员来管理。他的国家有n个城市，编号为1..n。一些城市之间有道路相连，任意两个不同的城市之间有且仅有一条直接或间接的道路。为了防止管理太过分散，每个省至少要有B个城市，为了能有效的管理，每个省最多只有3B个城市。每个省必须有一个省会，这个省会可以位于省内，也可以在该省外。但是该省的任意一个城市到达省会所经过的道路上的城市（除了最后一个城市，即该省省会）都必须属于该省。一个城市可以作为多个省的省会。聪明的你快帮帮这个国王吧！

Input

　　第一行包含两个数N，B（1<=N<=1000, 1 <= B <= N）。接下来N－1行，每行描述一条边，包含两个数，即这条边连接的两个城市的编号。

Output

　　如果无法满足国王的要求，输出0。否则输出数K，表示你给出的划分方案中省的个数，编号为1..K。第二行输出N个数，第I个数表示编号为I的城市属于的省的编号，第三行输出K个数，表示这K个省的省会的城市编号，如果有多种方案，你可以输出任意一种。

Sample Input

8 2
1 2
2 3
1 8
8 7
8 6
4 6
6 5

Sample Output

3
2 1 1 3 3 3 3 2
2 1 8

转载一篇神犇(shenben) 的blog

大神们都一句话带过。

在网上扒了老长时间，才领悟。

整理一下，发给大家。

先说一个问题：

如果在当前节点，其有两个儿子。

在左儿子的子树中搜得的节点数是a-1的话，这时候我们去搜右儿子，不妨假设右儿子为根的子树是一条长为10000的链。

那么显然最下面的叶节点跟左儿子代表的子树是在同一块中的，这显然碎成渣了。

所以每一次搜到一个节点，我们都要把当前的栈的指针当做该节点的相对栈底，这样的话就能避免这个问题。

首先以任意节点为根DFS，DFS的任务是解决以当前节点为根的子树（不包括当前节点）中的节点的归属，并汇报不知道去向何方者。

具体做法是：

对于每个正在处理的节点v，定义一个等待序列，扫一遍它的孩子。在扫的时候，假设当前for循环正在处理的是孩子u，DFS(u)，然后把传回来的等待序列加到v的等待序列中。如果v的等待序列节点数超过了B，那么就让等待队列中的节点组成一个省，省会是v，但v不划入那个省中。最后，for循环结束，把v加入到等待序列中，return。在主函数中接收DFS的返回值，若等待序列为空，皆大欢喜；但事实上，等待序列不可能为空。

那么等待序列中的节点数一定不超过B + 1，怎么办？答案就是放在上一个被划分出的省中去，那么上一个被划分出的省一定<=2B - 1，而现在最后无家可归的节点数一定<= B + 1，所以放在上一个被划分出的省中去节点数一定不超过3B！而且最后被划分出来的省一定和最后无家可归的节点中的某一个相邻！那么考虑什么时候会出现无解呢？答案就是我们找不到“上一个被划分出的省”

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<set>
#include<map>
using namespace std;

inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
	return x*f;
}
void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=1010;

int last[N],ecnt;
struct EDGE{int to,nt;}e[N<<1];
inline void add(int u,int v)
{e[++ecnt]=(EDGE){v,last[u]};last[u]=ecnt;}

int n,b,size[N],prvn,bel[N],cap[N],s[N],top;

void dfs(int u,int fa)
{
	register int i;
	s[++top]=u;
	for(i=last[u];i;i=e[i].nt)
	{
		if(e[i].to==fa)continue;
		dfs(e[i].to,u);
		if(size[e[i].to]+size[u]>=b)
		{size[u]=0;cap[++prvn]=u;while(s[top]^u)bel[s[top--]]=prvn;}
		else size[u]+=size[e[i].to];
	}
	size[u]++;
}

void clr(int u,int fa,int tp)
{
	bel[u]?tp=bel[u]:bel[u]=tp;register int i;
	for(i=last[u];i;i=e[i].nt)if(e[i].to^fa)clr(e[i].to,u,tp);
}

int main()
{
	n=read();b=read();
	register int i,u,v;
	for(i=1;i<n;++i){u=read();v=read();add(u,v);add(v,u);}
	dfs(1,1);clr(1,1,prvn);
	print(prvn);puts("");
	for(i=1;i<n;++i)print(bel[i]),putchar(' ');print(bel[n]);puts("");
	for(i=1;i<prvn;++i)print(cap[i]),putchar(' ');print(cap[prvn]);puts("");
	return 0;
}
/*
8 2
1 2
2 3
1 8
8 7
8 6
4 6
6 5

3
2 1 1 3 3 3 3 2
2 1 8
*/