首页 > 代码库 > BZOJ 1006 [HNOI2008]神奇的国度
BZOJ 1006 [HNOI2008]神奇的国度
Description
K国是一个热衷三角形的国度,连人的交往也只喜欢三角原则.他们认为三角关系:即AB相互认识,BC相互认识,CA
相互认识,是简洁高效的.为了巩固三角关系,K国禁止四边关系,五边关系等等的存在.所谓N边关系,是指N个人 A1A2
...An之间仅存在N对认识关系:(A1A2)(A2A3)...(AnA1),而没有其它认识关系.比如四边关系指ABCD四个人 AB,BC,C
D,DA相互认识,而AC,BD不认识.全民比赛时,为了防止做弊,规定任意一对相互认识的人不得在一队,国王相知道,
最少可以分多少支队。
Input
第一行两个整数N,M。1<=N<=10000,1<=M<=1000000.表示有N个人,M对认识关系. 接下来M行每行输入一对朋
友
Output
输出一个整数,最少可以分多少队
Sample Input
4 5
1 2
1 4
2 4
2 3
3 4
1 2
1 4
2 4
2 3
3 4
Sample Output
3
HINT
一种方案(1,3)(2)(4)
题解
我们把每个人都看做一个点,认识关系看做一条边,那么,“最少分多少支队”就是图的最小染色数。
我们知道,一般图的染色数是NPC问题,但是由于题目中的“禁止四边关系”等,可以发现其没有长度大于3的无弦环。
那么,这就是在求弦图的染色,直接套MCS算法按完美消除序列逆序染色即可。
具体细节请参照 陈丹琦 《弦图与区间图》
附代码:
#include <cctype> #include <cstdlib> #include <cstdio> #include <cstring> const int N = 10050; const int M = 1050000; struct Node{ int v, nxt; }nodes[M]; int cnt; int pre[N], to[M * 2], nxt[M * 2]; int head[N]; int f[N], P[N]; bool vis[N]; int n, m; inline void add(int h, int v) { nodes[++cnt].v = v; nodes[cnt].nxt = head[h]; head[h] = cnt; } void MCS() { memset(nodes, 0, sizeof nodes); memset(head, 0, sizeof head); memset(vis, 0, sizeof vis); memset(f, 0, sizeof f); int i, j, best = 0; for (i = 1; i <= n; ++i) add(0, i); for (j = n; j; --j) while (1) { int x; while (vis[nodes[x = head[best]].v]) head[best] = nodes[x].nxt; if (!x) { --best; continue; } x = nodes[x].v; vis[P[j] = x] = 1; for (i = pre[x]; ~i; i = nxt[i]) { if (vis[to[i]]) continue; add(++f[to[i]], to[i]); if (f[to[i]] > best) best = f[to[i]]; } break; } } int cnt2; inline void addedge(int u, int v) { nxt[cnt2] = pre[u]; to[pre[u] = cnt2++] = v; } inline int readInt() { int ans = 0; char c; do c = getchar(); while (!isdigit(c)); do { ans = ans * 10 + c - ‘0‘; c = getchar(); } while (isdigit(c)); return ans; } int mark[N], color[N]; int main() { cnt = cnt2 = 0; memset(pre, -1, sizeof pre); n = readInt(); m = readInt(); int x, y; while (m--) { x = readInt(); y = readInt(); addedge(x, y); addedge(y, x); } MCS(); memset(color, 0, sizeof color); memset(mark, 0, sizeof mark); int ans = 0; for (int j = n; j; --j) { int x = P[j]; for (int i = pre[x]; ~i; i = nxt[i]) mark[color[to[i]]] = j; for (color[x] = 1; mark[color[x]] == j; ++color[x]); if (color[x] > ans) ans = color[x]; } printf("%d\n", ans); return 0; }
BZOJ 1006 [HNOI2008]神奇的国度
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。