[BZOJ1086][SCOI2005]王室联邦（树上分块）

题目：

我是超链接

题解：

别看他给你鬼扯什么省会，实际上这道题的题目翻译一下就是：给你一棵树，要求把树分成大小在[b,3b]之内的块
差不多是树上分块的基本思想了

朴素的分块方式是贪心 能加就加 这种方法存在着严重的效率问题 可以被菊花卡成O(n)块
因此我们可以为其它的块预留位置 如果一块大小刚好>=b 就将这坨东西分成一块
首先任选一点开始深搜 维护一个栈 每个点退出递归时压栈 自下至上进行合并
如果某棵子树深搜完之后栈内元素数>=b 就把当前的栈内元素合并为一个块
但是这种方法存在一个问题 就是如果某棵子树深搜之后不到b 去深搜下一个子树 可能在下一个子树内部的某个位置超过b
这样会导致分成的块不连通
因此我们在每次进入递归时维护一个栈底，对于当前子树来说这个栈底就是整个栈的底，栈底以下的元素不能修改或弹栈
这样当一棵子树深搜过后由于子树内未分块节点不超过b，之前搜过的未分块节点数也不超过b，因此每块不超过2b
那么题目为什么给了3b呢？
深搜结束后可能会剩余一些节点，这些节点的数量不超过b，而且一定与当前分出的最后一块连通
因此我们将剩余节点分到最后一块中，可以保证最后一块的大小不超过3b
一遍深搜即可出解。—–引自PoPoQQQ

代码：

#include <cstdio>
using namespace std;
const int N=5000;
int tot,nxt[N],point[N],v[N],root[N],belong[N],stack[N],top,b,cnt;
void addline(int x,int y)
{
    ++tot; nxt[tot]=point[x]; point[x]=tot; v[tot]=y;
    ++tot; nxt[tot]=point[y]; point[y]=tot; v[tot]=x;
}
void dfs(int x,int fa)
{
    int bottom=top;
    for (int i=point[x];i;i=nxt[i])
      if (v[i]!=fa)
      {
        dfs(v[i],x);
        if(top-bottom>=b)
        {
            root[++cnt]=x; 
            while (top!=bottom) belong[stack[top--]]=cnt;
          }
      }
    stack[++top]=x;
}
int main()
{
    int n,x,y;scanf("%d%d",&n,&b);
    for (int i=1;i<n;i++)
    {
        scanf("%d%d",&x,&y);
        addline(x,y);
    }
    dfs(1,0);
    while (top) belong[stack[top--]]=cnt;
    printf("%d\n",cnt);
    for (int i=1;i<=n;i++) printf("%d ",belong[i]);
    printf("\n");
    for (int i=1;i<=cnt;i++) printf("%d ",root[i]);
}