从前一个和谐的班级,所有人都是搞OI的。有 nn 个是男生,有 00 个是女生。男生编号分别为 1,…,n1,…,n。
现在老师想把他们分成若干个两人小组写动态仙人掌,一个人负责搬砖另一个人负责吐槽。每个人至多属于一个小组。
有若干个这样的条件:第 vv 个男生和第 uu 个男生愿意组成小组。
请问这个班级里最多产生多少个小组?
输入格式
第一行两个正整数,n,mn,m。保证 n≥2n≥2。
接下来 mm 行,每行两个整数 v,uv,u 表示第 vv 个男生和第 uu 个男生愿意组成小组。保证 1≤v,u≤n1≤v,u≤n,保证 v≠uv≠u,保证同一个条件不会出现两次。
输出格式
第一行一个整数,表示最多产生多少个小组。
接下来一行 nn 个整数,描述一组最优方案。第 vv 个整数表示 vv 号男生所在小组的另一个男生的编号。如果 vv 号男生没有小组请输出 00。
样例一
input
10 20
9 2
7 6
10 8
3 9
1 10
7 1
10 9
8 6
8 2
8 1
3 1
7 5
4 7
5 9
7 8
10 4
9 1
4 8
6 3
2 5
output
5
9 5 6 10 2 3 8 7 1 4
样例二
input
5 4
1 5
4 2
2 1
4 3
output
2
2 1 4 3 0
限制与约定
1≤n≤5001≤n≤500,1≤m≤1247501≤m≤124750。
时间限制:1s1s
空间限制:256MB256MB
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<vector>
#include<queue>
using namespace std;
#define LL long long
#define pb push_back
#define Set(a, v) memset(a, v, sizeof(a))
#define For(i, a, b) for(int i = (a); i <= (int)(b); i++)
#define Forr(i, a, b) for(int i = (a); i >= (int)(b); i--)
#define MAXN (500+5)
struct BlossomTree{
vector<int> G[MAXN];
queue<int> q;
int n;
int pre[MAXN], link[MAXN], col[MAXN], fa[MAXN];
bool vis[MAXN];
void init(int rn){
n = rn;
Set(link, 0);
For(i, 0, n) G[i].clear();
}
void Addedges(int u, int v){
G[u].pb(v); G[v].pb(u);
}
int find(int x){
return x == fa[x]? x: (fa[x] = find(fa[x]));
}
void Combine(int u, int lca){
while(u != lca){
int f = link[u], g = pre[f];
if(find(g) != lca) pre[g] = f; //put the rever edge
if(col[f] == 1) col[f] = 2, q.push(f);
fa[find(u)] = find(f); fa[find(f)] = find(g);
u = pre[link[u]];
}
}
void Contract(int u, int v){
Set(vis, 0);
int lca = 0;
for(int i = u; i; i = pre[link[i]]) i = find(i), vis[i] = true;
for(int i = v; i; i = pre[link[i]]){
i = find(i);
if(vis[i]){lca = i; break;}
}
if(find(u) != lca) pre[u] = v;
if(find(v) != lca) pre[v] = u;
Combine(u, lca); Combine(v, lca);
}
void Bfs(int s){
while(!q.empty()) q.pop();
Set(col, 0); Set(pre, 0);
For(i, 1, n) fa[i] = i;
q.push(s);
while(!q.empty()){
int now = q.front(); q.pop();
For(i, 0, G[now].size()-1){
int v = G[now][i];
if(find(now)==find(v) || link[now]==v || col[v]==1) continue;
if(col[v] == 2) Contract(now, v);
else if(link[v]){
pre[v] = now; col[v] = 1;
col[link[v]] = 2;
q.push(link[v]);
}else{
pre[v] = now;
int g = link[now], f = now, o = v;
while(o){
link[f] = o; link[o] = f;
o = g; f = pre[o]; g = link[f];
}
return;
}
}
}
}
int MaxMatch(){
For(i, 1, n) if(!link[i]) Bfs(i);
int ret = 0;
For(i, 1, n) if(link[i]) ret++;
return ret/2;
}
void print(){
printf("%d", link[1]);
For(i, 2, n) printf(" %d", link[i]);
printf("\n");
}
}BT;
int main(){
int n, m;
scanf("%d%d", &n, &m);
BT.init(n);
For(i, 1, m){
int u, v;
scanf("%d%d", &u, &v);
BT.Addedges(u, v);
}
printf("%d\n", BT.MaxMatch());
BT.print();
return 0;
}