题目大意:
有n个人,m条有向边。对于一条边u->v表示u要从v得到一本书且仅一本(这一本数不一定是v原有的,也可以是v从别人那得到的),问在给出的边的条件下能不能使得每个人都得到一本书(非自己本人的)
题解:
目的就是让一个人可以能得到一本书,不管用什么的样的方式,只要符合条件就行,每个人都能有匹配就行了。所以只要求得最大匹配是否为n,判断一下。可以拆点建图,用二分图求或是用网络流求。
匈牙利算法(173ms)
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <queue>
using namespace std;
struct node
{
int u,v;
int nex;
}eage[40000];
int head[20000],cnt;
void ADD(int u,int v)
{
eage[cnt].u=u;
eage[cnt].v=v;
eage[cnt].nex=head[u];
head[u]=cnt++;
}
int used[20000],girl[20000];
bool dfs(int u)
{
for(int i=head[u];i!=-1;i=eage[i].nex)
{
int v=eage[i].v;
if(!used[v])
{
used[v]=1;
if(girl[v]==-1||dfs(girl[v]))
{
girl[v]=u;
return true;
}
}
}
return false;
}
int main()
{
int n,m;
while(~scanf("%d%d",&n,&m))
{
memset(head,-1,sizeof(head));
cnt=0;
int a,b;
while(m--)
{
scanf("%d%d",&a,&b);
ADD(a,b);
}
int all=0;
memset(girl,-1,sizeof(girl));
for(int i=0;i<n;i++)
{
memset(used,0,sizeof(used));
if(dfs(i))all++;
}
if(all==n)printf("YES\n");
else printf("NO\n");
}
return 0;
}
1. Hopcroft-Carp 算法(29ms)
#include <cstring>
#include <queue>
#include <cstdio>
#define INF 0x3f3f3f3f
using namespace std;
struct node
{
int u,v,nex;
}eage[30000];
int head[20000],cnt;
void ADD(int u,int v)
{
eage[cnt].u=u;
eage[cnt].v=v;
eage[cnt].nex=head[u];
head[u]=cnt++;
}
int dis,Nx,Ny;
int Mx[20000],My[20000],dx[20000],dy[20000];
int vst[20000];
bool searchP()
{
queue<int>q;
dis=INF;
memset(dx,-1,sizeof(dx));
memset(dy,-1,sizeof(dy));
for(int i=1;i<=Nx;i++)
{
if(Mx[i]==-1){q.push(i);
dx[i]=0;}
}
while(!q.empty())
{
int u=q.front();
q.pop();
if(dx[u]>dis)break;
for(int i=head[u];i!=-1;i=eage[i].nex)
{
int v=eage[i].v;
if(dy[v]==-1)
{
dy[v]=dx[u]+1;
if(My[v]==-1)dis=dy[v];
else
{
dx[My[v]]=dy[v]+1;
q.push(My[v]);
}
}
}
}
return dis!=INF;
}
bool DFS(int u)
{
for(int i=head[u];i!=-1;i=eage[i].nex)
{
int v=eage[i].v;
if(!vst[v]&&dy[v]==dx[u]+1)
{
vst[v]=1;
if(My[v]!=-1&&dy[v]==dis)continue;
if(My[v]==-1||DFS(My[v]))
{
My[v]=u;
Mx[u]=v;
return 1;
}
}
}
return 0;
}
int MaxMatch(int n)
{
int res=0;
memset(Mx,-1,sizeof(Mx));
memset(My,-1,sizeof(Mx));
Nx=Ny=n;
while(searchP())
{
memset(vst,0,sizeof(vst));
for(int i=1;i<=Nx;i++)
if(Mx[i]==-1&&DFS(i))res++;
}
return res;
}
int main()
{
int n,m;
while(~scanf("%d%d",&n,&m))
{
int a,b;
cnt=0;
memset(head,-1,sizeof(head));
while(m--)
{
scanf("%d%d",&a,&b);
ADD(a+1,b+1);
}
int all=MaxMatch(n);
if(all==n)printf("YES\n");
else printf("NO\n");
}
return 0;
}
网络流(43ms)
#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <queue>
#include <algorithm>
#include <cmath>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long LL;
const int maxn=1e5+10;
int cas=1;
int d[maxn];///标记数组
int cur[maxn];
int start,tend;///源点和汇点
struct node///建边
{
int to,cap,next;
}eage[maxn<<1];
int head[maxn];
bool vis[maxn];///标记
int cnt;
///建边函数
void ADD(int u,int v,int c)
{
eage[cnt].to=v;
eage[cnt].cap=c;
eage[cnt].next=head[u];
head[u]=cnt++;
eage[cnt].to=u;
eage[cnt].cap=0;
eage[cnt].next=head[v];
head[v]=cnt++;
}
bool BFS()
{
memset(d,-1,sizeof(d));
int Q[maxn*2];///队列
int Thead=0,Ttail=0;
Q[Ttail++]=start;
d[start]=0;///标记
while(Thead<Ttail)
{
//cout<<"sasadfasfdassgfafgdagfdd"<<endl;
int x=Q[Thead];
if(x==tend)return true;
for(int i=head[x];i!=-1;i=eage[i].next)
{
int temp=eage[i].to;
if(d[temp]==-1&&eage[i].cap>0)///没有标记,且可行流大于0
{
d[temp]=d[x]+1;
Q[Ttail++]=temp;
}
}
Thead++;
}
return false;///汇点是否成功标号,也就是说是否找到曾广路
}
int DFS(int x,int cap)
{
//cout<<"saasdadasd"<<endl;
if(x==tend)return cap;
int flow = 0 , f;
for(int i=head[x];i!=-1;i=eage[i].next)
{
int temp = eage[i].to;
if(d[temp]==d[x]+1&&eage[i].cap)
{
f=DFS(temp,min(cap-flow,eage[i].cap));
eage[i].cap -= f;
eage[i^1].cap+=f;
flow+=f;
if(flow==cap)break;
}
}
if(!flow)d[x]=-2;///防止重搜,剪枝。否则会超时
return flow;
}
int maxflow()
{
int flow=0,f;
while(BFS())
{
//cout<<"sad"<<endl;
while((f=DFS(start,INF))>0)
flow+=f;
}
return flow;
}
int main()
{
int n,m;
while(~scanf("%d%d",&n,&m))
{
cnt=0;
memset(head,-1,sizeof(head));
int u,v;
start=0;
tend=2*n+1;
for(int i=1;i<=n;i++)ADD(start,i,1);
for(int i=n+1;i<=2*n;i++)ADD(i,tend,1);
while(m--)
{
scanf("%d%d",&u,&v);///u<->v的流量是w(双向边)
ADD(u+1,v+1+n,1);
}
int ans=maxflow();
if(ans==n)printf("YES\n");
else printf("NO\n");
}
return 0;
}