-
7 7 1 2 1 3 2 4 3 4 4 5 4 6 5 7
样例输出
-
3
H市一共有N个旅游景点(编号1..N),由M条单向游览路线连接。在一个景点游览完后,可以顺着游览线路前往下一个景点。
为了避免游客重复游览同一个景点,游览线路保证是没有环路的。
每一个调查团可以从任意一个景点出发,沿着计划好的游览线路依次调查,到达终点后再返回。每个景点只会有一个调查团经过,不会重复调查。
举个例子:
上图中一共派出了3个调查团:
1. 蓝色:调查景点;2
2. 橙色:调查景点;1->3->4->6
3. 绿色:调查景点;5->7
当然对于这个图还有其他的规划方式,但是最少也需要3个调查团。
由于小组内的人数有限,所以大家希望调查团的数量尽可能少,同时也要将所有的景点都进行调查。
输入
第1行:2个整数N,M。1≤N≤500,0≤M≤20,000。
第2..M+1行:2个数字u,v,表示一条有向边(u,v)。保证不会出现重复的边,且不存在环。
输出
第1行:1个整数,表示最少需要的调查团数量。
对于一条路径,起点的入度为0,终点的出度为0,中间节点的出入度都为1。
每一个点最多只能有1个后继,同时每一个点最多只能有1个前驱。
假如我们选择了一条边(u,v),也就等价于把前驱u和后继v匹配上了。这样前驱u和后继v就不能和其他节点匹配。
利用这个我们可以这样来构图:
将每一个点拆分成2个,分别表示它作为前驱节点和后继节点。将所有的前驱节点作为A部,所有后继节点作为B部。
接下来进行连边,若原图中存在一条边(u,v),则连接A部的u和B部的v。
在这个上面做一个最大二分匹配:
其中实线表示被选中的匹配,虚线表示未被选中的。
有没有发现,和原图刚好有着对应的关系。
这样在匹配结束的时候,我们就可以直接通过匹配的情况来确定选中的路径。
但是如何保证这样就能得到最小的路径覆盖呢?
如果一个点是路径起点的话,它在B部的节点一定是没有匹配上的。
经过最大匹配算法后,B部剩下没有被匹配的点一定是最少的,也就对应了最小需要的路径数。
所以最小路径覆盖的结果才是N-最大匹配数。
正是这样,这样问题也就解决了。接下来第二个问题,怎么用网络流来解决二分匹配呢?
上一次我们讲了二分多重匹配,二分匹配不就是它的简化版么。
只需要把源点s到A部的边和B部到汇点t的边容量限定为1就可以了!
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#define MAX 1010
#define MAXCF 101
#define min(a,b) (a)>(b)?(b):(a)
using namespace std;
int cf[MAX][MAX];//存储图
int queue[MAX];//搜索队列
int path[MAX];//保存路径
int capacity[MAX];//流量数组,保存经过该点的最小流量
bool visited[MAX];//记录访问数组
int findAugmentPath(int T)
{
int i = 0, tail = 0;
memset(visited, 0, sizeof(visited));
queue[tail] = 1;//将源点加入队列
capacity[1] = MAXCF;
visited[1] = true;
while (i <= tail)
{
int u = queue[i];
if (u == T)
return capacity[T];//找到一条增广路径,返回该路径最小流量
for (int v = 2; v <= T; v++)
{
if (!visited[v] && cf[u][v] > 0)
{
path[v] = u;
capacity[v] = min(cf[u][v], capacity[u]);//记录路径上的最小残余流量
visited[v] = true;
tail++;
queue[tail] = v;
}
}
i++;
}
return 0;
}
void modifyGraph(int T)
{
int flow = capacity[T];
int now = T;
while (now != 1)
{
int fa = path[now];
cf[fa][now] -= flow;
cf[now][fa] += flow;
now = fa;
}
}
int main()
{
int N, M;
int i, j;
int u, v;
while (cin >> N >> M)
{
int S = 1;
int E = S + N + N + 1;
int result = 0;
memset(cf, 0, sizeof(cf));
for (i = 0; i < N; ++i)
{
cf[S][S + i + 1] = 1;
cf[S + N + i + 1][E] = 1;
}
for (i = 0; i < M; ++i)
{
cin >> u >> v;
cf[S + u][S + N + v] = 1;
}
int maxFlow = 0;
int delta = 0;
while (delta = findAugmentPath(E))
{
maxFlow += delta;
modifyGraph(E);
}
cout << N - maxFlow << endl;
}
return 0;
}