BZOJ2718毕业旅行

2718: [Violet 4]毕业旅行
Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 366 Solved: 204
Description

Input

Output
最多可选多少景点
Sample Input
7 6
1 2
2 3
5 4
4 3
3 6
6 7
Sample Output
2
Source
Ctsc2008 River & ural 1533. Fat Hobbits
最长反链=最小路径覆盖。。
裸题。。求最长反链。。
Dilworth定理。。
在有向无环图中，有如下的一些定义和性质：
链：一条链是一些点的集合，链上任意两个点x, y，满足要么 x 能到达 y ，要么 y 能到达 x 。
反链：一条反链是一些点的集合，链上任意两个点x, y，满足 x 不能到达 y，且 y 也不能到达 x。
一个定理：最长反链长度 = 最小链覆盖（用最少的链覆盖所有顶点）
对偶定理：最长链长度 = 最小反链覆盖
那么我们要求出的就是这个有向无环图的最小链覆盖了。最小链覆盖也就是路径可以相交的最小路径覆盖。
我们先来看路径不能相交的最小路径覆盖怎么来做：
建立一个二分图，两边都是n个点，原图的每个点 i 对应两个，在左边的叫做 i1, 在右边的叫做 i2 。
然后原图中如果存在一条边 (x, y)，那么就在二分图中建立 (x1, y2) 的边。
这样建立二分图之后，原图的点数 n - 二分图最大匹配 = 原图的最小路径覆盖（路径不能相交）。
这样为什么是对的呢？我们可以认为，开始时原图的每个点都是独立的一条路径，然后我们每次在二分图中选出一条边，就是将两条路径连接成一条路径，答案数就减少1。
而路径是不能相交的，所以我们在二分图中选出的边也是不能相交的所以就是二分图的最大匹配。
了解了路径不能相交的最小路径覆盖之后，怎么解路径可以相交的最小路径覆盖（也就是最小链覆盖）呢？
我们将原图做一次Floyd传递闭包，之后就可以知道任意两点 x, y，x 是否能到达 y。
如果两个点 x, y，满足 x 可以到达 y ，那么就在二分图中建立边 (x1, y2) 。
这样其实就是相当于将原图改造了一下，只要 x 能到达 y ，就直接连一条边 (x, y)，这样就可以“绕过”原图的一些被其他路径占用的点，直接构造新路径了。
这样就将可以相交的最小路径覆盖转化为了路径不能相交的最小路径覆盖了。
双倍经验题：BZOJ1143祭祀river
附上本蒟蒻的代码：

#include<cstdio>
#include<cstring>
using namespace std;
int n,m,h[40001],part[40001],ans=0,cnt=1,f[201][201];
bool used[40001];
struct kx
{
    int to,next;
}edge[40001];

int read()
{
    int w=0,c=1;
    char ch=getchar();
    while (ch<'0' || ch>'9')
      {
        if (ch=='-') c=-1;
        ch=getchar();
      }
    while (ch>='0' && ch<='9')
      w=w*10+ch-'0',ch=getchar();
    return w*c;
}

void add(int u,int v)
{
    cnt++,edge[cnt].next=h[u],h[u]=cnt,edge[cnt].to=v;
}

void floyed()
{
    int k,i,j;
    for (k=1;k<=n;k++)
      for (i=1;i<=n;i++)
        for (j=1;j<=n;j++)
          f[i][j]|=(f[i][k]&f[k][j]);//神奇的地方
}

bool find(int s)
{
    int k;
    for (k=h[s];k;k=edge[k].next)
      if (!used[edge[k].to])
        {
            used[edge[k].to]=true;
            if (!part[edge[k].to] || find(part[edge[k].to]))
              {
                part[edge[k].to]=s;
                return true;
              }
        }
    return false;
}

int main()
{
    int i,u,v,j;
    n=read(),m=read();
    for (i=1;i<=m;i++)
      u=read(),v=read(),f[u][v]=1;
    floyed();
    for (i=1;i<=n;i++)
      for (j=1;j<=n;j++)
        if (f[i][j])
          add(i,j+n);
    for (i=1;i<=n;i++)
      {
        memset(used,false,sizeof(used));
        if (find(i)) ans++;
      }
    printf("%d",n-ans);
    return 0;
}