• [HEOI2012][BZOJ2744] 朋友圈|匈牙利算法|最大独立集|最大团


    2744: [HEOI2012]朋友圈

    Time Limit: 30 Sec  Memory Limit: 128 MB
    Submit: 614  Solved: 175
    [Submit][Status][Discuss]

    Description

    在很久很久以前,曾经有两个国家和睦相处,无忧无虑的生活着。一年一度的评比大会开始了,作为和平的两国,一个朋友圈数量最多的永远都是最值得他人的尊敬,所以现在就是需要你求朋友圈的最大数目。
    两个国家看成是AB两国,现在是两个国家的描述:
    1.         A国:每个人都有一个友善值,当两个A国人的友善值a、b,如果a xor b mod 2=1,
    那么这两个人都是朋友,否则不是;
    2.         B国:每个人都有一个友善值,当两个B国人的友善值a、b,如果a xor b mod 2=0
    或者 (a or b)化成二进制有奇数个1,那么两个人是朋友,否则不是朋友;
    3.         A、B两国之间的人也有可能是朋友,数据中将会给出A、B之间“朋友”的情况。
    4.     在AB两国,朋友圈的定义:一个朋友圈集合S,满足

    S∈A∪ B                  ,对于所有的i,j∈  S ,i 和       j   是朋友

    由于落后的古代,没有电脑这个也就成了每年最大的难题,而你能帮他们求出最大朋 友圈的人数吗?

    Input

     

    第一行t<=6,表示输入数据总数。
    接下来t个数据:
    第一行输入三个整数A,B,M,表示A国人数、B国人数、AB两国之间是朋友的对数;第二行A个数ai,表示A国第i个人的友善值;第三行B个数bi,表示B国第j个人的友善值;
    第4——3+M行,每行两个整数(i,j),表示第i个A国人和第j个B国人是朋友。

    Output

     
    输出t行,每行,输出一个整数,表示最大朋友圈的数目。

    Sample Input

    2 4 7
    1 2
    2 6 5 4
    1 1
    1 2
    1 3
    2 1
    2 2
    2 3
    2 4

    Sample Output

    5
    【样例说明】
    最大朋友圈包含A国第1、2人和B国第1、2、3人。

    HINT

    【数据范围】

    两类数据

    第一类:|A|<=200 |B| <= 200

    第二类:|A| <= 10 |B| <= 3000

     

    Source

    纠结了一上午……

    第一次听说一个叫反图的东西。

    反图就是把原图连的边去掉,原图不连的边连上。

    我们可以看出,题目要求的是一个图的最大团,那么就等价于求反图的最大独立集。

    我们再看一遍题目。

    A国:如果a xor b mod 2=1,那么这两个人都是朋友,否则不是;

    也就是说奇数与偶数之间为朋友,而在反图中,A国部分就被划分为两部分,奇数部分是一个团,偶数部分是一个团。也就是说A国最多只能选出2个人(或者是1个,0个)。

    再看看B国。

    如果a xor b mod 2=0,或者 (a or b)化成二进制有奇数个1,那么两个人是朋友,否则不是朋友;也就是说,反图中的奇数间没有边相连,偶数间没有边相连,奇数与偶数间如果不满足“ (a or b)化成二进制有奇数个1”,那么就有边相连。于是二分图产生了。

    我们要求的是这个二分图的最大独立集,用B的总数减去二分图最大匹配(可以用匈牙利或者网络流)。

    等等……好像把A国给无视了(QAQ)。

    因为A国最多只能选出两个人,那我们就枚举就好了,枚举0个人,1个人,2个人(2个人要满足这俩人在反图中没有边相连,也就是说他俩是朋友),枚举出这些点之后,我们计算一下加入的俩点对B国中哪些点会产生影响,即不是朋友。直接去掉就好了。

    至此仅剩一个问题了:每次枚举重建图会妥妥TLE啊。

    时间戳就好了。

    膜POPOQQQ大爷。

    #include<iostream>
    #include<cstdio>
    #include<algorithm>
    #include<cstring>
    #include<cmath>
    #include<cstdlib>
    #define N 3010
    using namespace std;
    int head[N],cnt;
    int A,B,M,ans,T1,T2;
    int a[N],b[N];
    bool map[N][N];
    int res[N],sta[N],ban[N],tim[N];
    int next[N*N>>2],list[N*N>>2];
    inline int read()
    {
        int a=0,f=1; char c=getchar();
        while (c<'0'||c>'9') {if (c=='-') f=-1; c=getchar();}
        while (c>='0'&&c<='9') {a=a*10+c-'0'; c=getchar();}
        return a*f;
    }
    inline void insert(int x,int y)
    {
        next[++cnt]=head[x];
        head[x]=cnt;
        list[cnt]=y;
    }
    inline int count(int x)
    {
        int re=0;
        while (x) x-=x&-x,++re;
        return re;
    }
    bool hungary(int x)
    {
        if (ban[x]==T1) return false;
        for (int i=head[x];i;i=next[i])
            if (ban[list[i]]!=T1&&sta[list[i]]!=T2)
            {
                sta[list[i]]=T2;
                if (tim[list[i]]!=T1||!res[list[i]]||hungary(res[list[i]]))
                {
                    tim[list[i]]=T1;
                    res[list[i]]=x;
                    return 1;
                }
            }
        return 0;
    }        
    inline int maxset(int x=0,int y=0)
    {
        int re=0;
        ++T1;
        for (int i=1;i<=B;i++)
            if (map[x][i]||map[y][i])
                ban[i]=T1,++re;
        for (int i=1;i<=B;i++)
            if (b[i]&1)
            {
                ++T2;
                if (hungary(i)) ++re;
            }
        return B-re;
    }            
    int main()
    {
        int p,q;
        A=read(); B=read(); M=read();
        memset(map,1,sizeof(map));
        for (int i=1;i<=A;i++) a[i]=read();
        for (int i=1;i<=B;i++) b[i]=read();
        for (int i=1;i<=M;i++) p=read(),q=read(),map[p][q]=0;
        for (int i=1;i<=B;i++)
            if (b[i]&1)
                for (int j=1;j<=B;j++)
                    if (~b[j]&1)
                        if (~count(b[i]|b[j])&1)
                            insert(i,j);
        for (int i=1;i<=B;i++) map[0][i]=0;
        ans=maxset();
        for (int i=1;i<=A;i++)
            ans=max(maxset(i)+1,ans);
        for (int i=1;i<=A;i++)
            if (a[i]&1)
                for (int j=1;j<=A;j++)
                    if (~a[j]&1)
                        ans=max(maxset(i,j)+2,ans);
        printf("%d
    ",ans);
    }
  • 相关阅读:
    17、网卡驱动程序-DM9000举例
    16、NOR FLASH驱动框架
    15.1 linux操作系统下nand flash驱动框架2
    15、NAND FLASH驱动程序框架
    14、块设备驱动程序框架分析
    12.2 linux USB框架分析(详细注册match匹配过程)
    arm-linux-gcc: Command not found
    12、USB设备驱动程序
    POJ-2752 Seek the Name, Seek the Fame (KMP)
    POJ-2406 Power Strings (KMP)
  • 原文地址:https://www.cnblogs.com/ws-fqk/p/4759783.html
Copyright © 2020-2023  润新知