[bzoj1143][祭祀][ctsc2008]解题报告

传送门：ctsc2008祭祀

题面

---

在遥远的东方，有一个神秘的民族，自称Y族。他们世代居住在水面上，奉龙王为神。每逢重大庆典， Y族都会在水面上举办盛大的祭祀活动。我们可以把Y族居住地水系看成一个由岔口和河道组成的网络。每条河道连接着两个岔口，并且水在河道内按照一个固定的方向流动。显然，水系中不会有环流（下图描述一个环流的例子）。
由于人数众多的原因，Y族的祭祀活动会在多个岔口上同时举行。出于对龙王的尊重，这些祭祀地点的选择必须非常慎重。准确地说，Y族人认为，如果水流可以从一个祭祀点流到另外一个祭祀点，那么祭祀就会失去它神圣的意义。族长希望在保持祭祀神圣性的基础上，选择尽可能多的祭祀的地点。
Input
　　第一行包含两个用空格隔开的整数N、M，分别表示岔口和河道的数目，岔口从1到N编号。接下来M行，每行包含两个用空格隔开的整数u、v，描述一条连接岔口u和岔口v的河道，水流方向为自u向v。 N ≤ 100 M ≤ 1 000
Output
　　第一行包含一个整数K，表示最多能选取的祭祀点的个数。
Sample Input
4 4
1 2
3 4
3 2
4 2
Sample Output
2

---

分析

---

裸的求最长反链。
最长反链=最小链覆盖，证明 VFK orz
于是求出传递闭包后求个最小路径覆盖就吼了
（组合数学是本神书啊，感觉哪里都能遇到它

---

#include <cstdio>
#include <algorithm>
#include <iostream>
#include <vector>
#include <queue>

using namespace std;

const int N = 100 + 3;
const int inf = 0x7fffffff;

int mat[N][N];
int n, m, s, t;
int dis[N << 1];

struct edge
{
    int from, to, flow;
};

vector<edge> E;
vector<int> G[N << 1];

inline void add(int a, int b, int c)
{
    E.push_back((edge){a, b, c}); E.push_back((edge){b, a, 0});
    int k = E.size();
    G[a].push_back(k - 2); G[b].push_back(k - 1);
} 

inline void floyd()
{
    for (int i = 1; i <= n; i ++)
        for (int j = 1; j <= n; j ++)
            for (int k = 1; k <= n; k ++)
                if (mat[i][k] && mat[k][j])
                    mat[i][j] = 1;
}

bool bfs()
{
    for (int i = s; i <= t; i ++) dis[i] = -1;
    queue<int> q;
    q.push(s); dis[s] = 0;
    while (!q.empty())
    {
        int u = q.front(); q.pop();
        for (int i = 0; i < G[u].size(); i ++)
        {
            edge &e = E[G[u][i]];
            if (dis[e.to] == -1 && e.flow)
            {
                dis[e.to] = dis[u] + 1;
                q.push(e.to);
            }
        }
    }
    return dis[t] > -1;
}

int dfs(int x, int a)
{
    if (x == t || a == 0) return a;
    int flow = 0, f;
    for (int i = 0; i < G[x].size(); i ++)
    {
        edge &e = E[G[x][i]];
        if (dis[e.to] == dis[x] + 1 && (f = dfs(e.to, min(e.flow, a))))
        {
            e.flow -= f; E[G[x][i] ^ 1].flow += f;
            flow += f; a -= f;
            if (a == 0) break ;
        }
    }
    return flow;
}

int main()
{
    scanf("%d%d", &n, &m); s = 0, t = (n << 1) + 1;
    for (int i = 1; i <= m; i ++)
    {
        int a, b; scanf("%d%d", &a, &b);
        mat[a][b] = 1;
    }
    floyd();
    for (int i = 1; i <= n; i ++)
        for (int j = 1; j <= n; j ++)
        {
            if (i == j) continue ;
            else if (mat[i][j]) 
                add(i, j + n, 1);
        }
    for (int i = 1; i <= n; i ++)
    {
        add(s, i, 1);
        add(i + n, t, 1);
    }
    int ans = 0;
    while (bfs()) ans += dfs(s, inf);
    printf("%d", n - ans);
    return 0;
}