博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[HEOI2012][BZOJ2744] 朋友圈|匈牙利算法|最大独立集|最大团
阅读量:4696 次
发布时间:2019-06-09

本文共 3178 字,大约阅读时间需要 10 分钟。

2744: [HEOI2012]朋友圈

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

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
#include
#include
#include
#include
#include
#define N 3010using 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\n",ans);}

 

转载于:https://www.cnblogs.com/ws-fqk/p/4759783.html

你可能感兴趣的文章
linux socket高性能服务器处理框架
查看>>
处理:Header V3 DSA/SHA1 Signature, key ID 6b8d79e6: NOKEY error: Failed dependencies:
查看>>
php.ini详解(转)
查看>>
Java 并发:Executor
查看>>
[转]基于Python的接口测试框架
查看>>
计算机网络中的性能指标
查看>>
185. Department Top Three Salaries (Hard)
查看>>
Android Retrofit使用教程
查看>>
Django之视图函数
查看>>
python 简要小结
查看>>
MyBatis---简介
查看>>
elasticsearch api
查看>>
Install guide for OpenLDAP and GOsa 2 on Ubuntu & Debian
查看>>
python初学5
查看>>
Docker容器学习与分享05
查看>>
Linux 防火墙 iptables基本操作
查看>>
Anaconda中配置Pyspark的Spark开发环境
查看>>
spring接口文档注解:@ApiOperation
查看>>
(十一)数组array
查看>>
Centos7yum安装LNMP
查看>>