bzoj2744 [HEOI2012]朋友圈——二分图匹配

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=2744

首先，求一个图的最大团等价于求它的补图的最大独立集，而二分图的最大独立集 = 总点数 - 最大匹配数；

所以先把图转化成补图，也就是A国奇、偶点各自成团，B国奇、偶点相互连边而其内部无边，还有些A到B的边用邻接矩阵存了；

可以发现A国最多只能选出两个点来，而A国选的那些点会影响B国的最大独立集；

发现A国点很少，不妨暴力枚举！

注意计算最大独立集时减去的数 tmp，原来写的是 tmp = 最大匹配数，其中左部点(B国奇数点)不管，右部点和A国点相连的也匹配；

想的是右部点中和A国相连的就算匹配上了也会算在减去的部分，但这样不太对呢，而且会少算左部点中和A国相连的；

所以直接再来一个 out 数组表示和A国点相连，把B国中左右部的 out 点都算上，然后匹配时不去匹配 out 的右部点；

虽然感觉挺麻烦但其实也挺好写的。

代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<vector>
using namespace std;
int const maxn=;
int A,B,m,a[],b[maxn],pre[maxn],ans,tmp,sid[][maxn];
bool vis[maxn],out[maxn];
vector<int>vb[maxn];
bool dfs(int x)
{
    for(int i=;i<vb[x].size();i++)
    {
        int u=vb[x][i];
        if(vis[u]||out[u])continue; vis[u]=;
        if(!pre[u]||dfs(pre[u])){pre[u]=x;return ;}
    }
    return ;
}
bool ck(int x,int y)
{
    int k=(x|y),ret=;
    for(;k;k-=(k&-k))ret++;
    return (ret&);
}
int main()
{
    scanf("%d%d%d",&A,&B,&m);
    for(int i=;i<=A;i++)
    {
        scanf("%d",&a[i]);
//        if(a[i]&1)a1.push_back(i);
//        else a2.push_back(i);
    }
    for(int i=;i<=B;i++)
    {
        scanf("%d",&b[i]);
        for(int j=;j<i;j++)
        {
            if(ck(b[i],b[j])||(b[i]^b[j])%==)continue;//补图上无边
            vb[i].push_back(j),vb[j].push_back(i);
        }
    }
    for(int i=,x,y;i<=m;i++){scanf("%d%d",&x,&y); sid[x][y]=;}
    for(int j=;j<=B;j++)
    {
        if(b[j]%==)continue;//不dfs右部点
        memset(vis,,sizeof vis);
        if(dfs(j))ans++;
    }
    ans=B-ans;
    for(int i=;i<=A;i++)
    {
        tmp=;
        memset(pre,,sizeof pre);
        memset(out,,sizeof out);
        for(int j=;j<=B;j++)if(!sid[i][j])out[j]=,tmp++;//补图上有边
        for(int j=;j<=B;j++)
        {
            if(b[j]%==||out[j])continue;
//            if(!sid[i][j]){tmp++; continue;}
            memset(vis,,sizeof vis);
            if(dfs(j))tmp++;
        }
        tmp=B-tmp+;
        ans=max(ans,tmp);
    }
    for(int i=;i<=A;i++)
    {
//        int i=a1[ii];
        if(a[i]%==)continue;
        for(int j=;j<=A;j++)
        {
            if(a[j]%==)continue;
//            int j=a2[jj];
            tmp=;
            memset(pre,,sizeof pre);
            memset(out,,sizeof out);
            for(int k=;k<=B;k++)if(!sid[i][j]||!sid[j][k])out[j]=,tmp++;//!
            for(int k=;k<=B;k++)
            {
                if(b[k]%==||out[k])continue;//补图上有边
                memset(vis,,sizeof vis);
                if(dfs(k))tmp++;
            }
            tmp=B-tmp+;
            ans=max(ans,tmp);
        }
    }
    printf("%d\n",ans);
    return ;
}