hdu 2063(二分匹配)

<span style="font-size:18px;">#include <iostream>
#include <cstdio>
#include <vector>
#include <queue>
#include <cstring>
using namespace std;
const int maxn=1002,INF=1000000000;

//最大流模板
struct Edge
{
int from,to,cap,flow;
Edge(int u,int v,int c,int f):from(u),to(v),cap(c),flow(f){}
};

struct EdmondsKarp
{
int n,m,i;
vector<Edge> edges;
vector<int> G[maxn];
int a[maxn];
int p[maxn];

void init(int n)
{
for(i=0;i<=n;i++)G[i].clear();
edges.clear();
}

void AddEdge(int from,int to,int cap)
{
edges.push_back(Edge(from,to,cap,0));
edges.push_back(Edge(to,from,0,0));
m=edges.size();
G[from].push_back(m-2);
G[to].push_back(m-1);
}

int Maxflow(int s,int t)
{
int flow=0,u;
for(;;)
{
memset(a,0,sizeof(a));
queue<int> Q;
Q.push(s);
a[s]=INF;
while(!Q.empty())
{
int x=Q.front();
Q.pop();
for(i=0;i<G[x].size();i++)
{
Edge& e=edges[G[x][i]];
if(!a[e.to]&&e.cap>e.flow)
{
p[e.to]=G[x][i];
a[e.to]=min(a[x],e.cap-e.flow);
Q.push(e.to);
}
}
if(a[t])break;
}
if(!a[t])break;
for(u=t;u!=s;u=edges[p[u]].from)
{
edges[p[u]].flow+=a[t];
edges[p[u]^1].flow-=a[t];
}
flow+=a[t];
}
return flow;//flow是最大流。。
}
};

int main()
{
EdmondsKarp ans;
int k,M,N,s,t,i,j,x,y;
while(scanf("%d",&k)!=EOF)
{
s=0;t=1001;
if(k==0)break;
scanf("%d%d",&M,&N);
ans.init(1001);
for(i=1;i<=M;i++)ans.AddEdge(s,i,1);
for(i=1;i<=k;i++)
{
scanf("%d%d",&x,&y);
ans.AddEdge(x,y+500,1);
}
for(i=501;i<501+N;i++)ans.AddEdge(i,t,1);
printf("%d\n",ans.Maxflow(s,t));
}
return 0;
}

</span>

<span style="font-size:18px;">//摘录自紫书：针对无权图需要求出包含边数最多的匹配，即二分图的最大基数匹配。增加一个源点s，一个汇点t。从s向每个X连一条容量为1的弧，从所有的Y向t连一条容量为1的弧，最后把每条边变成一条由X指向Y的有向弧，容量为1.只要求出s到t的最大流，则原图中所有流量为1的弧对应了最大基数匹配。</span>

<span style="font-size:18px;">
</span>