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+;

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] = ;

            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 ;

    }

    for(int i = ;i <= n;i ++){

        scanf("%d%d",&u,&v);

        add(u,v); add(v,u);

    }

    dfs(,);

    if(tot == ) gen[++tot] = ;

    dfs1(,,tot);

    printf("%d\n",tot);

    for(int i = ;i <= n;i ++)

        printf("%d ",blo[i]);

    printf("\n");

    for(int i = ;i <= tot;i ++)

        printf("%d ",gen[i]);

    printf("\n");

    return ;

}

巴特西