【BZOJ 2744】【HEOI2012】朋友圈

题目链接：

　　　　　TP

题解：

　　对于A国，我们发现，最大团一定不大于2。对于B国，发现同奇偶性点之间都有边，不同奇偶性之间可能有边，也就是说对于B国是一个二分图最大团，也就是求B国补图的二分图最大独立集。然后，我们枚举使用A国的人员，将其与B国连接的点做一个补图，跑跑匈牙利即可。

　　【注】大视野上测试点和题面不一样啊！MMP没有t读入，只有一组数据，日哦！！

代码：

　　

#include <cstdio>
#include <iostream>
#include <cstring>

using namespace std;

inline int read(){
    int s=0,k=1;char ch=getchar();
    while(ch<'0'|ch>'9')    ch=='-'?k=-1:0,ch=getchar();
    while(ch>47&ch<='9')    s=s*10+(ch^48),ch=getchar();
    return s*k;
}

const int N=3e3+3;

int A,B,M;

struct edges{
    int v;edges *last;
}edge[N*N],*head[N<<2];int cnt;

inline void push(int u,int v){
    edge[++cnt]=(edges){v,head[u]};head[u]=edge+cnt;
}

struct Country{
    int sgl[N<<2],cpl[N<<2];
    int cnt_sgl,cnt_cpl;
    int val[N<<2],re[N<<2];
    inline void add(int x,int pos){
        if(x&1){
            sgl[++cnt_sgl]=pos;
            val[pos]=x;re[pos]=cnt_sgl;
        }else{
            cpl[++cnt_cpl]=pos;
            val[pos]=x;re[pos]=cnt_cpl;
        }
    }
    inline void clear(){
        cnt_sgl=cnt_cpl=0;
    }
}a,b;

int f[N<<2];bool vis[N<<2];

inline bool find(int x){

    for(edges *i=head[x];i;i=i->last){
        if(!vis[i->v]){
            vis[i->v]=true;
            if(f[i->v]==-1||find(f[i->v])){
                f[i->v]=x;
                return true;
            }
        }
    }
    return false;
}

int v[500][N<<2],size[N<<2],tt[N<<2];

inline void build(){
    for(int i=1;i<=b.cnt_sgl;i++)
        if(tt[b.sgl[i]]==2){
        for(int j=1;j<=b.cnt_cpl;j++)
            if(tt[b.cpl[j]]==2){
                int x=b.val[b.sgl[i]]|b.val[b.cpl[j]];
                int t(0);
                while(x) t+=x&1,x>>=1;
                if((~t)&1)
                    push(i,j);
            }
        }
}

inline int solve(){
    memset(f,-1,sizeof(f));
    int ret=0;
    for(int i=1;i<=B;i++)
        ret+=tt[i]==2;
    for(int i=1;i<=b.cnt_sgl;i++){
        if(tt[b.sgl[i]]==2){
        memset(vis,0,sizeof(vis));
            if(find(i))
                ret--;
        }
    }
    return ret;
}

inline void clear(){
    memset(head,0,sizeof(head));
    cnt=0;
}


int main(){
    // freopen(".in","r",stdin);
    // freopen(".out","w",stdout);
        clear();
        a.clear();b.clear();
        memset(v,0,sizeof(v));
        memset(size,0,sizeof(size));
        A=read(),B=read(),M=read();
        for(int i=1;i<=A;i++)    a.add(read(),i);
        for(int i=1;i<=B;i++)    b.add(read(),i);
        for(int i=1,u,vv;i<=M;i++){
            u=read(),vv=read();
            v[u][++size[u]]=vv;
        }
        for(int i=1;i<=B;i++)
            tt[i]=2;        
        build();
        int ans=(a.cnt_sgl>0)+(a.cnt_cpl>0);
        ans=max(ans,solve());
        memset(tt,0,sizeof(tt));
        for(int i=1;i<=a.cnt_sgl;i++){
            for(int j=1;j<=size[a.sgl[i]];j++)
                tt[v[a.sgl[i]][j]]=2;
            clear();
            build();
            ans=max(ans,solve()+1);
            for(int j=1;j<=size[a.sgl[i]];j++)
                tt[v[a.sgl[i]][j]]=0;
        }
        for(int j=1;j<=a.cnt_cpl;j++){
            for(int k=1;k<=size[a.cpl[j]];k++)
                tt[v[a.cpl[j]][k]]=2;
            clear();
            build();
            ans=max(ans,solve()+1);
            for(int k=1;k<=size[a.cpl[j]];k++)
                tt[v[a.cpl[j]][k]]=0;
        }
        for(int i=1;i<=a.cnt_sgl;i++){
            for(int j=1;j<=size[a.sgl[i]];j++)
                tt[v[a.sgl[i]][j]]++;
            for(int j=1;j<=a.cnt_cpl;j++){
                for(int k=1;k<=size[a.cpl[j]];k++)
                    tt[v[a.cpl[j]][k]]++;
                clear();
                build();
                ans=max(ans,solve()+2);
                for(int k=1;k<=size[a.cpl[j]];k++)
                    tt[v[a.cpl[j]][k]]--;
            }
            for(int j=1;j<=size[a.sgl[i]];j++)
                tt[v[a.sgl[i]][j]]--;
        }
        printf("%d
",ans);
    
}