bzoj 1086: [SCOI2005]王室联邦 (分块+dfs)

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
 
思路:
在树上dfs从根节点往下遍历,如果遍历到某结点的几棵子树加起来大于B,那么就把这几棵子树里的点扔到一个块里面也就是归为一个省(因为计算子树大小是从下到上的所以不用担心大于3B),这个块的省会就是当前点,这样处理完后还会剩下一些大小小于B的子树,我们再进行一遍dfs,把这些树里的点扔到其他能扔的块里面去。
 
实现代码;
#include<bits/stdc++.h>
using namespace std;
const int M = 1e4+10;
int n,B,tot,cnt,top;
int head[M],siz[M],gen[M],blo[M],q[M];

struct node{
    int to,next;
}e[M];

void add(int u,int v){
    e[++cnt].to = v;e[cnt].next = head[u];head[u] = cnt;
}

void dfs(int u,int fa){
    q[++top] = u;
    for(int i = head[u];i;i=e[i].next){
        int v = e[i].to;
        if(v == fa) continue;
        dfs(v,u);
        if(siz[u] + siz[v] >= B){
            siz[u] = 0;
            gen[++tot] = u;
            while(q[top]!=u)
                blo[q[top--]] = tot;
        }
        else
            siz[u] += siz[v];
    }
    siz[u] ++;
}

void dfs1(int u,int fa,int c){
    if(blo[u]) c = blo[u];
    else blo[u] = c;
    for(int i = head[u];i;i=e[i].next){
        int v = e[i].to;
        if(v == fa) continue;
        dfs1(v,u,c);
    }
}

int main()
{
    int u,v;
    scanf("%d%d",&n,&B);
    if(n < B){
        printf("0\n"); return 0;
    }
    for(int i = 2;i <= n;i ++){
        scanf("%d%d",&u,&v);
        add(u,v); add(v,u);
    }
    dfs(1,0);
    if(tot == 0) gen[++tot] = 1;
    dfs1(1,0,tot);
    printf("%d\n",tot);
    for(int i = 1;i <= n;i ++)
        printf("%d ",blo[i]);
    printf("\n");
    for(int i = 1;i <= tot;i ++)
        printf("%d ",gen[i]);
    printf("\n");
    return 0;
}

 

posted @ 2018-11-11 17:32  冥想选手  阅读(...)  评论(... 编辑 收藏