1086: [SCOI2005]王室联邦

(大概 打了5天比赛   学习计划 树上莫队一直都没来得及学 今天做了第一道树上分块的题  写挂了5，6次  许多细节没有注意到  慢慢改才发现 orz 继续加油

1086: [SCOI2005]王室联邦

Time Limit: 10 Sec  Memory Limit: 162 MBSec  Special Judge
Submit: 2319  Solved: 1423
[Submit][Status][Discuss]

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

?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83

#include <bits/stdc++.h> #define N 1005 #define ll long long using namespace std; ll readll(){     ll x=0,f=1;char ch=getchar();     while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}     while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}     return x*f; } int readint(){     int x=0,f=1;char ch=getchar();     while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}     while(isdigit(ch))x=x*10+ch-'0',ch=getchar();     return f*x; } vector<int>vec[N]; int n,b,cnt; int num[N],pre[N]; int p[N],fp[N]; bool vis[N]; void dfs1(int v,int fa){     num[v]=1;pre[v]=fa;     for(int i=0;i<vec[v].size();i++){         int u=vec[v][i];         if(u!=fa){             dfs1(u,v);             num[v]+=num[u];         }     }     if(num[v]>=b){         if(num[v]<=2*b){             p[v]=cnt;fp[cnt]=v;num[v]=0;cnt++;vis[v]=1;         }         else{             int ans=0;             for(int i=0;i<vec[v].size();i++){                 if(vec[v][i]==fa) continue;                 ans+=num[vec[v][i]];                 if(ans>=b){                     fp[cnt]=v;p[vec[v][i]]=cnt;cnt++;num[v]-=ans;                     ans=0;                 }             }             if(!ans) p[v]=cnt-1,num[v]=0;             int t=0;             for(int i=vec[v].size()-1;i>=0;i--){                 if(vec[v][i]==fa||vis[vec[v][i]]==1) continue;                 if(p[vec[v][i]]) t=p[vec[v][i]];                 else p[vec[v][i]]=t;                 if(p[vec[v][i]]) num[vec[v][i]]=0,vis[vec[v][i]]=1;             }         }     } } void dfs2(int v){     for(int i=0;i<vec[v].size();i++){         int u=vec[v][i];         if(u!=pre[v]){             if(!p[u]) p[u]=p[v];             dfs2(u);         }     } } int main(){     ios::sync_with_stdio(false);     n=readint();b=readint();     int u,v;cnt=1;     for(int i=1;i<n;i++){         u=readint();v=readint();         vec[u].push_back(v);vec[v].push_back(u);     }     dfs1(1,-1);     cnt--;     if(!p[1]) p[1]=cnt;     dfs2(1);     printf("%d\n",cnt);     for(int i=1;i<=n;i++) printf("%d ",p[i]);     printf("\n");     for(int i=1;i<=cnt;i++) printf("%d ",fp[i]);     printf("\n");     return 0; }