最小路径覆盖
题解:
假设 x , y x,y x,y之间有一条边,则把 x x x和 y y y分别拆成 x 1 , x 2 , y 1 , y 2 x_1,x_2,y_1,y_2 x1,x2,y1,y2。在 x 1 x_1 x1和 y 2 y_2 y2之间连接一条流量为 1 1 1的边,构造好二分图后答案即为 n − m a x m a t c h n-maxmatch n−maxmatch。对于方案的输出,可以利用并查集,对于 x x x和 y y y,如果流量流经 x − > y x->y x−>y,则 x − > y x->y x−>y在一条路径上,最后 O ( n 2 ) O(n^2) O(n2)输出即可。
代码:
#include<bits/stdc++.h>
using namespace std;
const int N = 100000;
const int INF = 0x3f3f3f3f;
int S,T,n,m,w[N],to[N],dep[N],head[N],nxt[N],num=1,u,v,pre[N];
void add(int u,int v,int ww){
num++;
to[num]=v;nxt[num]=head[u];w[num]=ww;head[u]=num;
num++;
to[num]=u;nxt[num]=head[v];w[num]=0;head[v]=num;
}
queue<int> q;
bool vis[N];
bool bfs(){
while(!q.empty()) q.pop();
memset(vis,0,sizeof(vis));
dep[S]=0;
vis[S]=1;
q.push(S);
while(!q.empty()){
int u=q.front();q.pop();
for(int i=head[u];i;i=nxt[i]){
int v=to[i];
if(vis[v]||w[i]<=0) continue;
vis[v]=1;
dep[v]=dep[u]+1;
q.push(v);
}
}
return vis[T];
}
int dfs(int u,int d){
if(u==T||d==0){
return d;
}
int ret=0;
for(int i=head[u];i;i=nxt[i]){
int v=to[i];
if(dep[v]!=dep[u]+1||w[i]<=0) continue;
int flow=dfs(v,min(w[i],d));//增广的流量
d-=flow;ret+=flow;
w[i]-=flow;w[i^1]+=flow;
if(d==0) break;
}
if(ret==0){//无法增广
dep[u]=-1;
}
return ret;
}
int par[N];
void init(){
for(int i=0;i<=2*n+1;i++) par[i]=i;
}
int find(int x){
return x==par[x]?x:par[x]=find(par[x]);
}
void unite(int x,int y){
x=find(x);y=find(y);
if(x==y) return ;
par[x]=y;
}
int main(){
scanf("%d %d",&n,&m);
S=0;T=2*n+1;
int ans=n;
for(int i=1;i<=m;i++){
scanf("%d %d",&u,&v);
add(u,v+n,1);
}
for(int i=1;i<=n;i++) add(S,i,1);
for(int i=n+1;i<=2*n;i++) add(i,T,1);
while(bfs()) ans-=dfs(S,INF);
init();
for(int i=1;i<=n;i++){
for(int j=head[i];j;j=nxt[j]){
if(to[j]&&!w[j]&&to[j]<=2*n){//如果i->to[j]还有流量
unite(i,to[j]-n);
}
}
}
for(int i=n;i;i--){
bool check=0;
for(int j=1;j<=n;j++){
if(find(j)==i){
check=1;
printf("%d ",j);
}
}
if(check) puts("");
}
printf("%d\n",ans);
return 0;
}