bzoj 2744: [HEOI2012]朋友圈

时间:2023-03-09 03:23:49
bzoj 2744: [HEOI2012]朋友圈
 #include<cstdio>
#include<iostream>
#define M 3010
using namespace std;
int A,B,m,a[M],b[M],map[M][M],head[M],next[M*M>>],u[M*M>>],cnt,ans,f1,f2,ban[M],f[M],k[M],pi[M];
int ji(int a1)
{
int sum=;
for(;a1;)
{
a1-=a1&-a1;
sum++;
}
return sum;
}
void jia(int a1,int a2)
{
cnt++;
next[cnt]=head[a1];
head[a1]=cnt;
u[cnt]=a2;
return;
}
bool xun(int a1)
{
for(int i=head[a1];i;i=next[i])
if(ban[u[i]]<f1&&f[u[i]]<f2)
{
f[u[i]]=f2;
if(k[u[i]]<f1||!pi[u[i]]||xun(pi[u[i]]))
{
k[u[i]]=f1;
pi[u[i]]=a1;
return ;
}
}
return ;
}
int make(int x=,int y=)
{
int re=;
f1++;
for(int i=;i<=B;i++)
if(!map[x][i]||!map[y][i])
{
ban[i]=f1;
re++;
}
for(int i=;i<=B;i++)
if(b[i]&&&ban[i]<f1)
{
f2++;
if(xun(i))
re++;
}
return B-re;
}
int main()
{
scanf("%d%d%d",&A,&B,&m);
for(int i=;i<=A;i++)
scanf("%d",&a[i]);
for(int i=;i<=B;i++)
scanf("%d",&b[i]);
for(int i=;i<=m;i++)
{
int a1,a2;
scanf("%d%d",&a1,&a2);
map[a1][a2]=;
}
for(int i=;i<=B;i++)
map[][i]=;
for(int i=;i<=B;i++)
if(b[i]&)
for(int j=;j<=B;j++)
if(~b[j]&)
if(~ji(b[i]|b[j])&)
jia(i,j);
ans=make();
for(int i=;i<=A;i++)
ans=max(make(i)+,ans);
for(int i=;i<=A;i++)
if(a[i]&)
for(int j=;j<=A;j++)
if(~a[j]&)
ans=max(make(i,j)+,ans);
printf("%d\n",ans);
return ;
}

由题A国只能取两个,一奇一偶。求最大团,最大团=补图最大独立子集。B国补图是一个二分图。枚举A国的点,在B国上跑最大独立子集=节点数-最大匹配数。