「BZOJ1086」[SCOI2005] 王室联邦
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
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
2 1 1 3 3 3 3 2
2 1 8
题解
一开始以为要求方案数,直接吓哭了。。。
仅当n<B是无解的吧。。
发现一个省至少要有B,dfs,如果子树大小超过B,直接子树划个省,根为省会
。。。
剩余的部分小于B的话随便扔哪都是合法的吧
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 |
#include<iostream> #include<cstdio> using namespace std; int read() { int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } int n,B,cnt,top,pro; int last[1005],q[1005],size[1005],cap[1005],belong[1005]; struct data{int to,next;}e[2005]; void insert(int u,int v) { e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt; e[++cnt].to=u;e[cnt].next=last[v];last[v]=cnt; } void dfs(int x,int fa) { q[++top]=x; for(int i=last[x];i;i=e[i].next) if(e[i].to!=fa) { dfs(e[i].to,x); if(size[x]+size[e[i].to]>=B) { size[x]=0; cap[++pro]=x; while(q[top]!=x) belong[q[top--]]=pro; } else size[x]+=size[e[i].to]; } size[x]++; } void paint(int x,int fa,int c) { if(belong[x])c=belong[x]; else belong[x]=c; for(int i=last[x];i;i=e[i].next) if(e[i].to!=fa) paint(e[i].to,x,c); } int main() { n=read();B=read(); if(n<B){puts("0");return 0;} for(int i=1;i<n;i++) { int u=read(),v=read(); insert(u,v); } dfs(1,0); if(!pro)cap[++pro]=1; paint(1,0,pro); printf("%d\n",pro); for(int i=1;i<=n;i++)printf("%d ",belong[i]); printf("\n"); for(int i=1;i<=pro;i++)printf("%d ",cap[i]); printf("\n"); return 0; } |
Subscribe