P5811 [IOI2019] 景点划分 题解
不妨令
先考虑原图是一棵树的时候怎么办?不妨枚举每一条边,设两边子树大小分别为
然而注意到这个事情不优美,不容易拓展到图上,我们考虑是否需要对每条边都这样考虑。事实上,我们找到树的重心,考虑将树视为以重心为根的树。考虑重心的某个儿子,如果子树大小
证明:
先证明充分性,即如果存在子树大小
接着证明如果不存在这样的子树则无解。如果不存在这样的子树,那么对于
于是树上问题就处理完了。
考虑图上问题。我们想让图和树有联系,必然是求出图的一棵 DFS 生成树,找出重心。如果重心有一个儿子子树大小
考虑 DFS 生成树的根为
证明:
考虑最终累加的子树大小必然是
于是整个题就做完了。
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <cstring>
#include <string>
#include <cassert>
#include <queue>
#include <vector>
using namespace std;
const int N = 2e5 + 5;
int n, m;
vector<int> G[N];
int A, B, C, rA, rB, rC;
int sz[N], fa[N];
int ans[N];
bool vis[N];
vector<int> NG[N];
int wc;
pair<int, int> a[5];
int id[N], idx, ra[N];
void dfs(int u, int f)
{
id[u] = ++idx;
ra[idx] = u;
fa[u] = f;
vis[u] = 1;
int maxn = 0;
sz[u] = 1;
for (auto& j : G[u])
{
if (vis[j]) continue;
NG[u].emplace_back(j);
NG[j].emplace_back(u);
dfs(j, u);
sz[u] += sz[j];
maxn = max(maxn, sz[j]);
}
maxn = max(maxn, n - sz[u]);
if (maxn <= (n / 2)) wc = u;
}
int cnt = 0;
void Get_A(int u, int fa)
{
ans[u] = a[1].second;
cnt++;
if (cnt == A) return;
for (auto& j : NG[u])
{
if (j == fa) continue;
Get_A(j, u);
if (cnt == A) return;
}
}
void Get_B(int u, int fa)
{
if (ans[u]) return;
ans[u] = a[2].second;
cnt++;
if (cnt == B) return;
for (auto& j : NG[u])
{
if (j != fa)
{
Get_B(j, u);
if (cnt == B) return;
}
}
}
class Union_Find
{
public:
int fa[N];
void Init()
{
for (int i = 0; i < N; i++) fa[i] = i;
}
int find(int u)
{
return (fa[u] == u ? u : fa[u] = find(fa[u]));
}
void merge(int u, int v)
{
fa[find(u)] = find(v);
}
}uf;
bool flag[N];
int sflag[N];
int main()
{
ios::sync_with_stdio(0), cin.tie(0);
cin >> n >> m >> A >> B >> C;
rA = A, rB = B, rC = C;
a[1] = make_pair(A, 1);
a[2] = make_pair(B, 2);
a[3] = make_pair(C, 3);
sort(a + 1, a + 4);
A = a[1].first, B = a[2].first, C = a[3].first;
for (int i = 1; i <= m; i++)
{
int u, v;
cin >> u >> v;
u++, v++;
G[u].emplace_back(v);
G[v].emplace_back(u);
}
dfs(1, 0);
for (auto& j : NG[wc])
{
int szs = 0;
if (j == fa[wc])
{
szs = n - sz[wc];
}
else szs = sz[j];
if (szs >= A)
{
Get_A(j, wc);
cnt = 0;
Get_B(wc, 0);
for (int i = 1; i <= n; i++)
{
if (!ans[i]) ans[i] = a[3].second;
}
for (int i = 1; i <= n; i++) cout << ans[i] << " ";
cout << "\n";
return 0;
}
}
int faf = fa[wc];
do
{
flag[faf] = 1;
faf = fa[faf];
} while (faf);
int szs = n - sz[wc];
vector<int> v;
for (auto& j : NG[wc])
{
if (j == fa[wc]) continue;
bool fg = 0;
for (int k = id[j]; k < id[j] + sz[j]; k++)
{
for (auto& p : G[ra[k]])
{
if (flag[p]) fg = 1;
}
}
if (!fg) continue;
szs += sz[j];
v.emplace_back(j);
if (szs >= A)
{
Get_A(fa[wc], wc);
for (auto& k : v)
{
Get_A(k, wc);
}
cnt = 0;
Get_B(wc, 0);
for (int i = 1; i <= n; i++)
{
if (!ans[i]) ans[i] = a[3].second;
}
for (int i = 1; i <= n; i++) cout << ans[i] << " ";
cout << "\n";
return 0;
}
}
for (int i = 1; i <= n; i++) cout << "0 ";
cout << "\n";
return 0;
}
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列01:轻松3步本地部署deepseek,普通电脑可用
· 25岁的心里话
· 按钮权限的设计及实现