ans=n-二分图最大匹配
对还要输出方案
#include<queue>
#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
const int maxn=2010;
const int inf=0x7fffffff/3;
int n,m,h[maxn],tot=1,vis[maxn],S,T;
struct edge{int to,next,w;}G[maxn*100];
int xi[maxn],cur[maxn];
bool mark[maxn];
void add(int x,int y,int z){
G[++tot].to=y;G[tot].next=h[x];h[x]=tot;G[tot].w=z;
G[++tot].to=x;G[tot].next=h[y];h[y]=tot;G[tot].w=0;
}
bool bfs(){
for (int i=S;i<=T;++i) vis[i]=-1;
queue<int>q; q.push(S); vis[S]=0;
while (!q.empty()){
int u=q.front(); q.pop();
for (int i=h[u];i;i=G[i].next){
int v=G[i].to;
if (vis[v]==-1&&G[i].w){
vis[v]=vis[u]+1;
q.push(v);
}
}
}return vis[T]!=-1;
}
int dfs(int x,int f){
if (x==T||!f) return f;
int w,used=0;
for (int i=cur[x];i;i=G[i].next)
if (vis[G[i].to]==vis[x]+1){
w=f-used;
w=dfs(G[i].to,min(G[i].w,w));
if(w){xi[x]=G[i].to;if(G[i].to>n)mark[G[i].to-n]=1;}
G[i].w-=w; G[i^1].w+=w;
used+=w; if (G[i].w) cur[x]=i;
if (used==f) return f;
}
if (!used) vis[x]=-1;
return used;
}
int dinic(){
int ret=0;
while (bfs()){
for (int i=S;i<=T;++i) cur[i]=h[i];
ret+=dfs(S,inf);
}return ret;
}
int main(){
scanf("%d%d",&n,&m);
S=0; T=2*n+1;
for (int i=1;i<=m;++i){
int x,y; scanf("%d%d",&x,&y);
add(x,y+n,inf);
}
for (int i=1;i<=n;++i)
add(S,i,1),add(i+n,T,1);
int ret=dinic();
for (int i=1;i<=n;++i){
if (mark[i]) continue;
printf("%d ",i);
int tmp=i;
while(xi[tmp])printf("%d ",tmp=xi[tmp]-n);
printf("\n");
}
printf("%d\n",n-ret);
}