BZOJ 1086：[SCOI2005]王室联邦（DFS树分块）

http://www.lydsy.com/JudgeOnline/problem.php?id=1086

题意：给出n个点的树，让你对树进行分块，每块的大小范围在[b, 3b]之间。

思路：一开始想着维护一个sz[u]代表以u为根的子树（不包括u本身）的大小，如果在范围之内就分成一块，但是这样写感觉一些细节上的东西没考虑清楚，WA了很久。

看了别人的做法：http://blog.****.net/popoqqq/article/details/42772237，很简洁。

大概是维护一个栈，每次dfs结束的时候入栈，在每次遍历子树之前标记一下当前的栈顶，当遍历子树的时候，如果当前的栈顶 - 之前的栈顶 >= b的话，就可以分成一块了。最后剩下的元素和最后一块连通。

结果的正确性：因为之前的子树大小加起来小于b，当前的子树大小小于b，所以当前的块大小小于2*b。最后剩余的元素小于b，最后的块小于2*b，因此小于3*b。

因为b小于n的，所以可以保证一定有答案。

 #include <cstdio>
 #include <cstring>
 #include <vector>
 #include <algorithm>
 using namespace std;
 #define N 1010
 struct Edge {
     int v, nxt;
 } edge[N*];
 int n, b, head[N], tot, belong[N], stk[N], top;
 vector<int> ans;

 void Add(int u, int v) {
     edge[tot] = (Edge) {v, head[u]}; head[u] = tot++;
     edge[tot] = (Edge) {u, head[v]}; head[v] = tot++;
 }

 void dfs(int u, int fa) {
     int base = top; // 之前的栈顶
     for(int i = head[u]; ~i; i = edge[i].nxt) {
         int v = edge[i].v;
         if(v == fa) continue;
         dfs(v, u);
         if(top - base >= b) {
             ans.push_back(u);
             while(base < top) belong[stk[top--]] = ans.size();
         }
     }
     stk[++top] = u;
 }

 int main() {
     while(~scanf("%d%d", &n, &b)) {
         memset(head, -, sizeof(head)); tot = top = ;
         memset(belong, -, sizeof(belong));
         for(int i = ; i < n; i++) {
             int u, v; scanf("%d%d", &u, &v);
             Add(u, v);
         }
         dfs(, -);
         while(top) belong[stk[top--]] = ans.size(); // 最后元素
         printf("%d\n", ans.size());
         for(int i = ; i <= n; i++) printf("%d%c", belong[i], i == n ? '\n' : ' ');
         for(int i = ; i < ans.size(); i++) printf("%d%c", ans[i], ans.size() -  == i ? '\n' : ' ');
     }
     return ;
 }