二分图匹配一道~用的网络流解。
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
using namespace std;
const int inf=1<<29;
const int maxn=40010;
int n,k,e,st,des,pnt[maxn],nxt[maxn],head[maxn],flow[maxn],level[maxn];
queue<int> q;
void addEdge(int u,int v,int f)
{
pnt[e]=v;nxt[e]=head[u];flow[e]=f;head[u]=e++;
pnt[e]=u;nxt[e]=head[v];flow[e]=0;head[v]=e++;
}
bool BFS()
{
memset(level,0,sizeof(level));
while(!q.empty())
q.pop();
q.push(st);
level[st]=1;
while(!q.empty())
{
int u=q.front();
q.pop();
if(u==des)
return true;
for(int i=head[u];i!=-1;i=nxt[i])
if(flow[i]&&!level[pnt[i]])
{
level[pnt[i]]=level[u]+1;
q.push(pnt[i]);
}
}
return level[des];
}
int DFS(int u,int sum)
{
if(u==des)
return sum;
for(int i=head[u],t;i!=-1;i=nxt[i])
{
if(flow[i]&&level[pnt[i]]==level[u]+1&&(t=DFS(pnt[i],min(sum,flow[i]))))
{
flow[i]-=t;
flow[i^1]+=t;
return t;
}
}
return level[u]=0;
}
int maxflow()
{
int ans=0;
while(BFS())
while(1)
{
int val=DFS(st,inf);
if(!val)
break;
ans+=val;
}
return ans;
}
int main()
{
while(scanf("%d%d",&n,&k)!=EOF)
{
e=st=0;
des=n*2+1;
memset(head,-1,sizeof(head));
for(int i=0;i<k;i++)
{
int u,v;
scanf("%d%d",&u,&v);
addEdge(u,n+v,1);
}
for(int i=1;i<=n;i++)
addEdge(0,i,1);
for(int i=1;i<=n;i++)
addEdge(n+i,des,1);
printf("%d\n",maxflow());
}
return 0;
}