zoukankan      html  css  js  c++  java
  • BZOJ-1143&&BZOJ-2718 祭祀river&&毕业旅行 最长反链(Floyed传递闭包+二分图匹配)

    蛋蛋安利的双倍经验题
    

    1143: [CTSC2008]祭祀river
    Time Limit: 10 Sec Memory Limit: 162 MB
    Submit: 1901 Solved: 951
    [Submit][Status][Discuss]

    Description
    在遥远的东方,有一个神秘的民族,自称Y族。他们世代居住在水面上,奉龙王为神。每逢重大庆典, Y族都会在水面上举办盛大的祭祀活动。我们可以把Y族居住地水系看成一个由岔口和河道组成的网络。每条河道连接着两个岔口,并且水在河道内按照一个固定的方向流动。显然,水系中不会有环流(下图描述一个环流的例子)。
    这里写图片描述
    由于人数众多的原因,Y族的祭祀活动会在多个岔口上同时举行。出于对龙王的尊重,这些祭祀地点的选择必须非常慎重。准确地说,Y族人认为,如果水流可以从一个祭祀点流到另外一个祭祀点,那么祭祀就会失去它神圣的意义。族长希望在保持祭祀神圣性的基础上,选择尽可能多的祭祀的地点。

    Input
    第一行包含两个用空格隔开的整数N、M,分别表示岔口和河道的数目,岔口从1到N编号。接下来M行,每行包含两个用空格隔开的整数u、v,描述一条连接岔口u和岔口v的河道,水流方向为自u向v。

    Output
    第一行包含一个整数K,表示最多能选取的祭祀点的个数。

    Sample Input
    4 4
    1 2
    3 4
    3 2
    4 2

    Sample Output
    2

    【样例说明】
    在样例给出的水系中,不存在一种方法能够选择三个或者三个以上的祭祀点。包含两个祭祀点的测试点的方案有两种:
    选择岔口1与岔口3(如样例输出第二行),选择岔口1与岔口4。
    水流可以从任意岔口流至岔口2。如果在岔口2建立祭祀点,那么任意其他岔口都不能建立祭祀点
    但是在最优的一种祭祀点的选取方案中我们可以建立两个祭祀点,所以岔口2不能建立祭祀点。对于其他岔口
    至少存在一个最优方案选择该岔口为祭祀点,所以输出为1011。

    HINT
    对于每个测试点:如果你仅输出了正确的被选取的祭祀点个数,那么你将得到该测试点30%的分数;如果你仅输出了正确的被选取的祭祀点个数与一个可行的方案,那么你将得到该测试点60%的分数;如果你的输出完全正确,那么你将得到该测试点100%的分数

    【数据规模】 N ≤ 100 M ≤ 1 000

    Source


    2718: [Violet 4]毕业旅行
    Time Limit: 10 Sec Memory Limit: 128 MB
    Submit: 376 Solved: 209
    [Submit][Status][Discuss]

    Description
    这里写图片描述
    Input
    这里写图片描述
    Output
    最多可选多少景点

    Sample Input
    7 6
    1 2
    2 3
    5 4
    4 3
    3 6
    6 7

    Sample Output
    2

    HINT
    这里写图片描述
    Source
    Ctsc2008 River & ural 1533. Fat Hobbits

    题意:
    给出一个有向无环图。在其中找出一个最大的点集使得点集中任意两个点之间不可达。
    思路:
    最长反链;那么何为反链? 链是一个点的集合,这个集合中任意两个元素v、u,要么v能走到u,要么u能走到v。 反链是一个点的集合,这个集合中任意两点谁也不能走到谁。
    求法:
    Dilworth定理:最长反链长度 = 最小链覆盖数
    对于最小链覆盖数,可以利用二分图匹配来求;又应用了一个结论:最小链覆盖数=总点数-最大匹配数
    实现:
    Floyed传递闭包,处理连通性问题; 匈牙利算法硬上,得结果。

    证明:(Dilworth定理)
    首先,最长反链中的每个点,一定在不同的链中,所以得到:最长反链长度 ≤ 最小链覆盖数
    下面关心的是,最长反链长度 ≥ 最小链覆盖数吗?
    用数学归纳法。
    首先空的图显然成立。
    假设对于所有点数小于图G的图都成立,现在来证G也满足条件。
    两个定义:源:没有一个点能走到它的点。(入度为0) 汇:不能走到任何一个点的点。(出度为0)
    设:
    这个图点集为V
    这个图的最长反链为A
    能走到A中某个点的点的集合为B,A中某个点能走到的点构成的集合为C。
    (1)若A中至少有一个点既不是源也不是汇
    则:
    1. B ≠ V (否则反链的点全是汇)
    2. C ≠ V (否则反链的点全是源)
    3. B ∪ C = V:反证法,假设有一个点v既不属于B也不属于C,那么说明反链上任何一点不能到v,v也不能到反链上任何一个点,反链加上v之后能成为更长的反链,与假设矛盾。
    4. B ∩ C = A:反证法,假设有一个点v既属于B也属于C但不属于A,那么说明反链上某点a能到v,v能到反链上某点b。若a = b,则与图是有向无环图矛盾。若a ≠ b,则a能先到v,在从v到达b,与a、b都在反链上矛盾。
    设B构成的子图、C构成的子图的最小链覆盖的链的集合分别为C[B]和C[C]。因为|B| < |V|,|C| < |V|,由归纳假设,C[B]、C[C]大小分别等于B构成的子图和C构成的子图的最长反链的大小。
    考虑一条链c ∈ C[B],则c中必然存在一个在最长反链上的点v,不然最小链覆盖会大于最长反链,与归纳假设矛盾。因为B ∩ C = A,所以C[B]中v不能到任何一个点。所以v必为c的汇。
    于是得到:A中的点都是C[B]中某条链的汇。
    同理:A中的点都是C[C]中某条链的源。
    又B ∪ C = V,B ∩ C = A,C[B]中某条链的汇必为C[C]中某条链的源,于是把C[B]和C[C]拼起来,就得到一个原图G的一种可行的链覆盖了!
    考虑这个链覆盖,发现这个链覆盖大小 = 最长反链长度,得到最小链覆盖大小 ≤ 最长反链长度。又由前面的最小链覆盖大小 ≥ 最长反链长度,所以G中最小链覆盖大小 = 最长反链长度
    (2)若A是所有源的集合或所有汇的集合
    若A是所有源的集合,那么从A中随便取一个元素v(这个可以随便找到吧)
    再从V中随便找一个汇u,使得v可以到u(这个也可以随便找到吧)(u可以等于v)
    然后G中去掉v、u形成图G’。则G’的最长反链大小一定小于等于|A| - 1,不然的就可以在G中找到一条不选v、u的最长反链,利用(1)的结论可以证明正确性。
    由于归纳假设,则G’中一定可以找到一个数量为|A| - 1的链覆盖。再加上{v} ∪ {u}这条链,则可以在G中找到一个数量为|A|的链覆盖,所以G的最小链覆盖 ≤ 最长反链长度。
    又由前面的最小链覆盖大小 ≥ 最长反链长度,所以G中最小链覆盖大小 = 最长反链长度。
    证完了,这就是传说中的Dilworth定理:最小链覆盖数 = 最长反链长度。
    其对偶定理:最长链长度 = 最小反链覆盖数 就不证了,网上的证明烂大街了。
    by vfk

    code:

    #include<iostream>
    #include<cstdio>
    #include<cstring>
    #include<cmath>
    #include<algorithm>
    using namespace std;
    int read()
    {
        int x=0,f=1; char ch=getchar();
        while (ch<'0' || ch>'9') {if (ch=='-') f=-1; ch=getchar();}
        while (ch>='0' && ch<='9') {x=x*10+ch-'0'; ch=getchar();}
        return x*f;
    }
    #define maxn 210
    #define maxm 1010
    int n,m,ans=0;
    int road[maxn][maxn],part[maxn],f[maxn];
    bool used[maxn];
    
    void floyed()
    {
        for (int k=1; k<=n; k++)
            for (int i=1; i<=n; i++)
                for (int j=1; j<=n; j++)            
                    if (road[i][k] && road[k][j])
                        road[i][j]=1;
    }
    
    bool find(int x)
    {
        if (f[x]) return 0;
        f[x]=1;
        for (int i=1; i<=n; i++)
            if (road[x][i])
                {
                    if (!used[i] || find(part[i]))  
                        {            
                            used[i]=1;part[i]=x;return 1;
                        }         
                }
        return 0;
    }
    
    int main()
    {
        n=read(),m=read();
        for (int i=1; i<=m; i++)
            {
                int u=read(),v=read();
                road[u][v]=1;
            }
        floyed();
        for (int i=1; i<=n; i++)
            {
                memset(f,0,sizeof(f));
                if (find(i)) ans++;
            }
        printf("%d
    ",n-ans);
        return 0;
    }
  • 相关阅读:
    【web性能优化】DNS解析与ip
    【web性能优化】雅虎军规
    【web性能优化】优化提纲及图片优化(慕课网笔记)
    【web性能优化】常用缓存方式(慕课网学习笔记)
    【前端】企业微信客户端调试
    【es6】es6使用集锦
    【前端】遇到的各种报错
    【前端】安装wampserver提示丢失MSVCR100.dll的解决方法
    【es6】将2个数组合并为一个数组
    【web】使用ionic搭建移动端项目 icon-radio 标签在ios下全部选中的问题
  • 原文地址:https://www.cnblogs.com/DaD3zZ-Beyonder/p/5346151.html
Copyright © 2011-2022 走看看