首页 > 代码库 > BZOJ 1086:[SCOI2005]王室联邦(DFS树分块)

BZOJ 1086:[SCOI2005]王室联邦(DFS树分块)

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

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

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

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

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

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

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

 1 #include <cstdio> 2 #include <cstring> 3 #include <vector> 4 #include <algorithm> 5 using namespace std; 6 #define N 1010 7 struct Edge { 8     int v, nxt; 9 } edge[N*2];10 int n, b, head[N], tot, belong[N], stk[N], top;11 vector<int> ans;12 13 void Add(int u, int v) {14     edge[tot] = (Edge) {v, head[u]}; head[u] = tot++;15     edge[tot] = (Edge) {u, head[v]}; head[v] = tot++;16 }17 18 void dfs(int u, int fa) {19     int base = top; // 之前的栈顶20     for(int i = head[u]; ~i; i = edge[i].nxt) {21         int v = edge[i].v;22         if(v == fa) continue;23         dfs(v, u);24         if(top - base >= b) {25             ans.push_back(u);26             while(base < top) belong[stk[top--]] = ans.size();27         }28     }29     stk[++top] = u;30 }31 32 int main() {33     while(~scanf("%d%d", &n, &b)) {34         memset(head, -1, sizeof(head)); tot = top = 0;35         memset(belong, -1, sizeof(belong));36         for(int i = 1; i < n; i++) {37             int u, v; scanf("%d%d", &u, &v);38             Add(u, v);39         }40         dfs(1, -1);41         while(top) belong[stk[top--]] = ans.size(); // 最后元素42         printf("%d\n", ans.size());43         for(int i = 1; i <= n; i++) printf("%d%c", belong[i], i == n ? \n :  );44         for(int i = 0; i < ans.size(); i++) printf("%d%c", ans[i], ans.size() - 1 == i ? \n :  );45     }46     return 0;47 }

 

BZOJ 1086:[SCOI2005]王室联邦(DFS树分块)