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

#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;
}

转载于:https://www.cnblogs.com/wang9897/p/8454622.html