分析:
一共有N个旅游景点(编号1..N),由M条单向路线连接。游览线路保证是没有环路的。每一个调查团可以从任意一个景点出发,沿着计划好的游览线路依次调查。每个景点只会有一个调查团经过,不会被重复调查。求所需要的最小调查团数。题解:
将这个问题转换为二分图匹配来做:- 路径的性质:起点的入度为0,终点的出度为0,中间节点的出入度都为1。
- 转换为二分图匹配:说明除了起点之外,一条路径上的其他点入度都为1,可以转换为除了起点外,其他点都能被另一个点匹配上。没有被匹配上的点肯定是起点,而有多少起点,就有多少条路,要求最小路径覆盖,即求最少路径数->最少起点数->最少不被匹配点数->N-最大匹配数->二分图求最大匹配!
AC代码:
对比发现果真网络流快些,所以以后就只用网络流做二分图匹配,Hungary什么的就( ^_^ )/~~拜拜了(虽然它代码短QwQ)
- 网络流 ISAP做法:
- 网络流 ISAP做法:
/*************************************************************************
> File Name: test.cpp
> Author: Akira
> Mail: qaq.febr2.qaq@gmail.com
************************************************************************/
#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <cmath>
#include <vector>
#include <set>
#include <list>
#include <ctime>
typedef long long LL;
typedef unsigned long long ULL;
typedef long double LD;
#define MST(a,b) memset(a,b,sizeof(a))
#define CLR(a) MST(a,0)
#define Sqr(a) ((a)*(a))
using namespace std;
#define MaxN 100000
#define MaxM MaxN*10
#define INF 0x3f3f3f3f
#define bug cout<<88888888<<endl;
int N,M; //N个点,M条边
struct Edge{
int u,v,next;
int flow;
}edge[MaxM]; //最大边数,一般都是实际边数的2倍甚至以上
int head[MaxN];
int cont;
void init(){ //记得init
cont = 0;
MST(head, -1);
}
void add(int u, int v, int flow){
edge[cont].u = u;
edge[cont].v = v;
edge[cont].flow = flow;
edge[cont].next = head[u];
head[u] = cont++;
}
void Add(int u, int v, int flow){ //建正反两条边,反向流量为0
add(u, v, flow);
add(v, u, 0);
}
int dis[MaxN];
int num[MaxN];
int cur[MaxN];
int pre[MaxN];
void BFS(int source,int sink)
{
queue<int>q;
CLR(num);
MST(dis,-1);
q.push(sink);
dis[sink]=0;
num[0]=1;
while(!q.empty())
{
int u=q.front();
q.pop();
for(int i=head[u];i!=-1;i=edge[i].next)
{
int v = edge[i].v;
if(dis[v] == -1)
{
dis[v] = dis[u] + 1;
num[dis[v]]++;
q.push(v);
}
}
}
}
int ISAP(int source,int sink,int n) //从源点到汇点,n为总点数,返回最大流
{
memcpy(cur,head,sizeof(cur));
int flow=0, u = pre[source] = source;
BFS( source,sink);
while( dis[source] < n )
{
if(u == sink)
{
int df = INF, pos;
for(int i =source;i != sink;i = edge[cur[i]].v)
{
if(df > edge[cur[i]].flow)
{
df = edge[cur[i]].flow;
pos = i;
}
}
for(int i = source;i != sink;i = edge[cur[i]].v)
{
edge[cur[i]].flow -= df;
edge[cur[i]^1].flow += df;
}
flow += df;
//cout << flow << endl;
u = pos;
}
int st;
for(st = cur[u];st != -1;st = edge[st].next)
{
if(dis[edge[st].v] + 1 == dis[u] && edge[st].flow)
{
break;
}
}
if(st != -1)
{
cur[u] = st;
pre[edge[st].v] = u;
u = edge[st].v;
}
else
{
if( (--num[dis[u]])==0 ) break;
int mind = n;
for(int id = head[u];id != -1;id = edge[id].next)
{
if(mind > dis[edge[id].v] && edge[id].flow != 0)
{
cur[u] = id;
mind = dis[edge[id].v];
}
}
dis[u] = mind+1;
num[dis[u]]++;
if(u!=source)
u = pre[u];
}
}
return flow;
}
int main()
{
while(~scanf("%d%d", &N, &M))
{
init();
int S = 0;
int T = N*2+1;
for(int i=1;i<=N;i++) Add(S, i, 1);
for(int i=1;i<=N;i++) Add(N+i, T, 1);
int u,v;
while(M--)
{
scanf("%d%d", &u, &v);
Add(u, N+v, 1);
}
cout << N - ISAP(S,T,N*2+2) <<endl;
}
system("pause");
}
- 匈牙利算法:
/*************************************************************************
> File Name: test.cpp
> Author: Akira
> Mail: qaq.febr2.qaq@gmail.com
************************************************************************/
#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <cmath>
#include <vector>
#include <set>
#include <list>
#include <ctime>
typedef long long LL;
typedef unsigned long long ULL;
typedef long double LD;
#define MST(a,b) memset(a,b,sizeof(a))
#define CLR(a) MST(a,0)
#define Sqr(a) ((a)*(a))
using namespace std;
#define MaxN 100000
#define MaxM MaxN*10
#define INF 0x3f3f3f3f
#define bug cout<<88888888<<endl;
struct Edge
{
int from;
int to;
int weight;
Edge(int f,int t, int w):from(f),to(t), weight(w){}
};
vector<int> G[MaxN];
/* G[i] 存储顶点 i 出发的边的编号 */
vector<Edge> edges;
typedef vector<int>::iterator it;
int N,M; //总点数
int Left; //左边点
int Right; //右边点
int matching[MaxN]; // 存储求解结果
int check[MaxN];
bool DFS(int u)
{
for(it i = G[u].begin(); i!=G[u].end(); i++)//对u的每个邻接点
{
int v = edges[*i].to;
if(!check[v])//要求不在交替路中
{
check[v] = true; //放入交替路
if(matching[v] == -1 || DFS(matching[v])) //如果是未覆盖点,说明交替路为增广路,则交换路径,并返回成功
{
matching[v] = u;
matching[u] = v;
return true;
}
}
}
return false; //不存在增广路,返回失败
}
int Hungarian()
{
int ans = 0;
MST(matching, -1);
for(int u = 1;u<=Left;u++)
{
if(matching[u] == -1)
{
CLR(check);
if(DFS(u))
ans++;
}
//cout << ans << endl;
}
return ans;
}
int main()
{
while(~scanf("%d%d", &N, &M))
{
Left = N;
int u,v;
for(int i=0;i<M;i++)
{
scanf("%d%d", &u, &v);
Edge a(u,N+v,1);
edges.push_back(a);
G[u].push_back(i);
}
cout << N-Hungarian() <<endl;
}
system("pause");
}