【刷题】BZOJ 2744 [HEOI2012]朋友圈

Description

在很久很久以前，曾经有两个国家和睦相处，无忧无虑的生活着。一年一度的评比大会开始了，作为和平的两国，一个朋友圈数量最多的永远都是最值得他人的尊敬，所以现在就是需要你求朋友圈的最大数目。

两个国家看成是AB两国，现在是两个国家的描述：

A国：每个人都有一个友善值，当两个A国人的友善值a、b，如果a xor b mod 2=1，那么这两个人都是朋友，否则不是；
B国：每个人都有一个友善值，当两个B国人的友善值a、b，如果a xor b mod 2=0或者 (a or b)化成二进制有奇数个1，那么两个人是朋友，否则不是朋友；
A、B两国之间的人也有可能是朋友，数据中将会给出A、B之间“朋友”的情况。
在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

Solution

看A国的条件，实际上就是要求两个数奇偶性不同，那么可以二分，而且最大团为2
看B国的条件，连出的图就是奇和奇连边，偶和偶连边，一部分奇和偶连边；可以发现它的补图就是一部分奇和偶连边，可以二分。并且二分图的最大独立集等于其补图的最大团
于是
建出B国的补图
对A国分情况

一个都不选时，直接对B的补图跑最大匹配，更新答案
选一个时，枚举A国选的那个人，找出B国能够选的人，把它们从补图中拎出来，跑最大匹配，+1更新答案
选两个时，枚举A国选的两个人，接下来同上，只是+2更新答案

基本上就是这样，最后用些时间戳加速什么的

#include<bits/stdc++.h>
#define ui unsigned int
#define ll long long
#define db double
#define ld long double
#define ull unsigned long long
const int MAXA=200+10,MAXB=3000+10,inf=0x3f3f3f3f;
int a,b,m,T,e,connect[MAXB],use[MAXB],A[MAXA],B[MAXB],G[MAXA][MAXB],beg[MAXB],nex[MAXB*MAXB],to[MAXB*MAXB],ans,clk1,clk2,avail[MAXB];
template<typename T> inline void read(T &x)
{
    T data=0,w=1;
    char ch=0;
    while(ch!='-'&&(ch<'0'||ch>'9'))ch=getchar();
    if(ch=='-')w=-1,ch=getchar();
    while(ch>='0'&&ch<='9')data=((T)data<<3)+((T)data<<1)+(ch^'0'),ch=getchar();
    x=data*w;
}
template<typename T> inline void write(T x,char ch='\0')
{
    if(x<0)putchar('-'),x=-x;
    if(x>9)write(x/10);
    putchar(x%10+'0');
    if(ch!='\0')putchar(ch);
}
template<typename T> inline void chkmin(T &x,T y){x=(y<x?y:x);}
template<typename T> inline void chkmax(T &x,T y){x=(y>x?y:x);}
template<typename T> inline T min(T x,T y){return x<y?x:y;}
template<typename T> inline T max(T x,T y){return x>y?x:y;}
inline void insert(int x,int y)
{
    to[++e]=y;
    nex[e]=beg[x];
    beg[x]=e;
}
inline void Build()
{
    for(register int i=1;i<=b;++i)
        for(register int j=1;j<=b;++j)
            if(i==j||(!(B[i]&1))||(B[j]&1))continue;
            else if(((B[i]^B[j])&1)&&__builtin_popcount(B[i]|B[j])%2==0)insert(i,j);
}
inline bool dfs(int x)
{
    for(register int i=beg[x];i;i=nex[i])
        if(avail[to[i]]==clk1&&use[to[i]]!=clk2)
        {
            use[to[i]]=clk2;
            if(!connect[to[i]]||dfs(connect[to[i]]))
            {
                connect[to[i]]=x;
                return true;
            }
        }
    return false;
}
inline int solve(int a1,int a2)
{
    memset(connect,0,sizeof(connect));
    int res=0;
    ++clk1;
    for(register int i=1;i<=b;++i)
        if(G[a1][i]&&G[a2][i])res++,avail[i]=clk1;
    for(register int i=1;i<=b;++i)
        if(avail[i]==clk1&&(B[i]&1))
        {
            ++clk2;
            if(dfs(i))res--;
        }
    return res;
}
int main()
{
    ans=0;
    read(a);read(b);read(m);
    for(register int i=1;i<=a;++i)read(A[i]);
    for(register int i=1;i<=b;++i)read(B[i]);
    Build();
    for(register int i=1;i<=m;++i)
    {
        int u,v;read(u);read(v);
        G[u][v]=1;
    }
    for(register int i=1;i<=b;++i)G[0][i]=1;
    chkmax(ans,solve(0,0));
    for(register int i=1;i<=a;++i)chkmax(ans,solve(i,0)+1);
    for(register int i=1;i<=a;++i)
        for(register int j=i+1;j<=a;++j)
            if((A[i]^A[j])&1)chkmax(ans,solve(i,j)+2);
    write(ans,'\n');
    return 0;
}