【BZOJ4205】卡牌配对
Description
现在有一种卡牌游戏,每张卡牌上有三个属性值:A,B,C。把卡牌分为X,Y两类,分别有n1,n2张。
两张卡牌能够配对,当且仅当,存在至多一项属性值使得两张卡牌该项属性值互质,且两张卡牌类别不同。
比如一张X类卡牌属性值分别是225,233,101,一张Y类卡牌属性值分别为115,466,99。那么这两张牌是可以配对的,因为只有101和99一组属性互质。
游戏的目的是最大化匹配上的卡牌组数,当然每张卡牌只能用一次。
Input
数据第一行两个数n1,n2,空格分割。
接下来n1行,每行3个数,依次表示每张X类卡牌的3项属性值。
接下来n2行,每行3个数,依次表示每张Y类卡牌的3项属性值。
Output
输出一个整数:最多能够匹配的数目。
Sample Input
2 2
2 2 2
2 5 5
2 2 5
5 5 5
Sample Output
2
【提示】
样例中第一张X类卡牌和第一张Y类卡牌能配对,第二张X类卡牌和两张Y类卡牌都能配对。所以最佳方案是第一张X和第一张Y配对,第二张X和第二张Y配对。
另外,请大胆使用渐进复杂度较高的算法!
HINT
对于100%的数据,n1,n2≤ 30000,属性值为不超过200的正整数
Solution
又一道场上想不出来的题。
首先考虑普通网络流匹配:左边是一类卡牌,右边是另一类卡牌,两两之间看能否配对连一条容量为1的边,做一次Dinic出解。
但是完整数据中两边点数都是3万,肯定无法两两连边。这时候我们想到优化建图:我们肯定是要构造一些特殊信息点,将符合信息的卡牌分别连向它们,以此减少边数、完成匹配。
照例从源点向每一个X卡牌连一条容量为1的边,每一个Y卡牌向汇点连一条容量为1的边,限制匹配次数。
考虑中间辅助匹配的信息点是什么:至多一组数互质等价于至少两组数不互质。不互质,意味着有公约数;而至少两组数的话,我们可以分拆成三种情况:一对卡牌的AB满足条件、AC满足条件、BC满足条件。
公约数怎么处理?我们可以把所有公约数二元组看成一个信息点。以12为例,我们在中间建一列信息点\((x,y)\),所有满足\(x|A\)且\(y|B\)的X卡牌向这个点连一条容量为1的边,相应地从这个点向满足条件的Y卡牌连一条容量为1的边。我们发现,这样恰好可以使得满足AB都不互质的卡牌们进行配对过程。
同理建出其他两组信息点。三组信息点同时存在时,我们发现不论如何,一张X卡牌最多通过一种方式也就是一组信息点和另一边的某一张Y卡牌配对,同样地Y卡牌不论如何也只会通过一种方式配对。即使是ABC都满足的一对卡牌,也最多只会选择从某一组信息点流过。
建模完成,过程优美。
但是\((x,y)\)的点数足足有\(3*200*200\)个,太多了。其实有公约数不就是有公质因子吗?200以内的素数仅仅有46个,所以现在我们只需要把信息点的\((x,y)\)变更成质数的组合\((p_i,p_j)\),按上述连边即可。信息点点数仅仅有\(3*46*46=6348\)个。
Dinic是跑得过的,有极大的速度优势。
Code
#include <cstdio>
#include <queue>
using namespace std;
const int N=30005,INF=1000000000,SZ=68000,M=6000000;
int p[N],pcnt,minp[N],id[N];
bool vis[N];
int n1,n2,a[N*2][3];
int h[SZ],tot=1;
int S,T,dis[SZ],cur[SZ];
queue<int> q;
struct Edge{int v,f,next;}e[M*2];
inline void addEdge(int u,int v,int f){
e[++tot]=(Edge){v,f,h[u]}; h[u]=tot;
e[++tot]=(Edge){u,0,h[v]}; h[v]=tot;
}
void sieve(){
minp[1]=1;
for(int i=2;i<=200;i++){
if(!vis[i]){
p[++pcnt]=i;
minp[i]=i;
id[i]=pcnt;
}
for(int j=1;j<=pcnt&&i*p[j]<=200;j++){
int x=i*p[j];
vis[x]=true;
if(i%p[j]==0){
minp[x]=p[j];
break;
}
minp[x]=p[j];
}
}
}
bool bfs(){
for(int i=1;i<=T;i++) dis[i]=-1,cur[i]=h[i];
while(!q.empty()) q.pop();
q.push(S);
dis[S]=0;
while(!q.empty()){
int u=q.front(); q.pop();
for(int i=h[u],v;i;i=e[i].next)
if(e[i].f&&dis[v=e[i].v]==-1){
dis[v]=dis[u]+1;
if(v==T) return true;
q.push(v);
}
}
return dis[T]!=-1;
}
int dfs(int u,int flow){
if(u==T) return flow;
int res=0,get;
for(int i=cur[u],v;i&&flow;i=e[i].next)
if(e[i].f&&dis[v=e[i].v]==dis[u]+1){
get=dfs(v,min(flow,e[i].f));
e[i].f-=get;
e[i^1].f+=get;
if(e[i].f) cur[u]=i;
flow-=get;
res+=get;
}
if(!res) dis[u]=-1;
return res;
}
int dinic(){
int res=0;
while(bfs())
res+=dfs(S,INF);
return res;
}
void deal(int x,int *l,int &cnt){
cnt=0;
while(x!=1){
l[++cnt]=id[minp[x]];
int p=minp[x];
while(x%p==0) x/=p;
}
}
int main(){
sieve();
scanf("%d%d",&n1,&n2);
for(int i=1;i<=n1+n2;i++) scanf("%d%d%d",&a[i][0],&a[i][1],&a[i][2]);
S=n1+n2+46*46*3+1; T=S+1;
for(int i=1;i<=n1;i++) addEdge(S,i,1);
for(int i=1;i<=n2;i++) addEdge(n1+i,T,1);
int l1[5],l2[5],l3[5],t1,t2,t3;
for(int u=1;u<=n1+n2;u++){
deal(a[u][0],l1,t1);
deal(a[u][1],l2,t2);
deal(a[u][2],l3,t3);
for(int i=1;i<=t1;i++)
for(int j=1;j<=t2;j++){
int id=46*(l1[i]-1)+(l2[j]-1)+1;
u<=n1?addEdge(u,n1+n2+id,1):addEdge(n1+n2+id,u,1);
}
for(int i=1;i<=t1;i++)
for(int j=1;j<=t3;j++){
int id=46*(l1[i]-1)+(l3[j]-1)+1;
u<=n1?addEdge(u,n1+n2+46*46+id,1):addEdge(n1+n2+46*46+id,u,1);
}
for(int i=1;i<=t2;i++)
for(int j=1;j<=t3;j++){
int id=46*(l2[i]-1)+(l3[j]-1)+1;
u<=n1?addEdge(u,n1+n2+2*46*46+id,1):addEdge(n1+n2+2*46*46+id,u,1);
}
}
printf("%d\n",dinic());
return 0;
}
【BZOJ4205】卡牌配对的更多相关文章
- BZOJ4205卡牌配对——最大流+建图优化
题目描述 现在有一种卡牌游戏,每张卡牌上有三个属性值:A,B,C.把卡牌分为X,Y两类,分别有n1,n2张. 两张卡牌能够配对,当且仅当,存在至多一项属性值使得两张卡牌该项属性值互质,且两张卡牌类别不 ...
- BZOJ4205 : 卡牌配对
对于两张卡牌,如果存在两种属性值不互质,则可以匹配. 只考虑200以内的质数,一共有46个,可以新建3*46*46个点来表示一类属性值中有这两种质数的卡牌. 然后对于每张卡牌,枚举它的质因子,最多只有 ...
- 【BZOJ4205】卡牌配对 最大流
[BZOJ4205]卡牌配对 Description 现在有一种卡牌游戏,每张卡牌上有三个属性值:A,B,C.把卡牌分为X,Y两类,分别有n1,n2张. 两张卡牌能够配对,当且仅当,存在至多一项属性值 ...
- [BZOJ4205][FJ2015集训]卡牌配对
题目:卡牌配对 传送门:None 题目大意:有$n_1$张$X$类牌和$n_2$张$Y$类类牌,每张卡牌上有三个属性值:$A,B,C$.两张卡牌能够配对,当且仅当,存在至多一项属性值使得两张卡牌该项属 ...
- BZOJ 4205: 卡牌配对
4205: 卡牌配对 Time Limit: 20 Sec Memory Limit: 512 MBSubmit: 173 Solved: 76[Submit][Status][Discuss] ...
- 刷题总结——卡牌配对(bzoj4205网络流)
题目: Description 现在有一种卡牌游戏,每张卡牌上有三个属性值:A,B,C.把卡牌分为X,Y两类,分别有n1,n2张. 两张卡牌能够配对,当且仅当,存在至多一项属性值使得两张卡牌该项属性值 ...
- [BZOJ4205][FJ2015集训] 卡牌配对 [建图+最大流]
题面 这是bzoj权限题,题面可以去下面的离线题库找 离线4205,只有题面,不能提交 思路 二分图匹配 这道题模型显然就是个二分图匹配嘛 那我们两两判断一下然后连边匹配.....就只有30分了 因为 ...
- Unity3D_(游戏)卡牌01_启动屏界面
卡牌2D游戏展示 (游戏代码放到 卡牌04_游戏界面 文章最后面~) 游戏项目已托管到github上(里面有个32bit可执行文件) 传送门 规则 开始游戏每张卡牌初始翻开展示 展示几秒后卡牌 ...
- 使用UIKit制作卡牌游戏(三)ios游戏篇
译者: Lao Jiang | 原文作者: Matthijs Hollemans写于2012/07/13 转自朋友Tommy 的翻译,自己只翻译了这第三篇教程. 原文地址: http://www.ra ...
随机推荐
- MATLAB复制图片时边框大的问题
当使用MATLAB画图时,需要将图片复制到word中,会发现图片有一个白色的边框,在论文的排版中是一个影响美观的问题 例如: >> x = 0:10; >> y = sin(x ...
- maven项目的标准目录结构
maven项目的标准目录结构如下:
- 织梦调用多个栏目typeid="1,2,3"不支持的解决方法
织梦arclist调用副栏目不显示的解决办法: 打开/include/taglib/arclist.lib.php,代码约位于295-296行,查找以下两行代码: if($CrossID=='') $ ...
- [Github] Github使用教程
前言 Github是一个面向开源及私有软件项目的托管平台.它可以免费使用,并且速度快速,拥有超多的用户.是目前管理软件开发和发现已有代码的首选平台.下面将向Github新手介绍相关操作. 正文 注册 ...
- GO/GOLANG程序员笔记大全
---------------------------------------- go 并发 // 注解:go 语言天生为程序并发所设计,可以说go的强项就是在cpu并发上的处理. // go 语言层 ...
- [mysql] 归档工具pt-archiver,binlog格式由mixed变成row
pt-archiver官方地址:https://www.percona.com/doc/percona-toolkit/3.0/pt-archiver.html 介绍:归档数据,比如将一年前的数据备份 ...
- iOS中使用RNCryptor对资源文件加密(先加密后拖进项目中)
概述:IPA 在发布时,业务相关的敏感资源文件以明文的形式存储,由于没有加密保护,这些文件在应用发布后 可能被其他人获取,并结合其他漏洞和手段产生真实攻击.所以我们要 1.在设计.开发阶段,集合业务确 ...
- final发布视频展示博客
Part One [探路者]选题展示视频链接: http://v.youku.com/v_show/id_XMzIxMDM2MTQ1Ng==.html?spm=a2h3j.8428770.341605 ...
- JavaScript实现弹出层(以layer.open为例)
首先,引用layer,自行下载. 添加如下两行 <script src=" ../layer/jquery.min.js"></script> <sc ...
- quartusII13.0使用教程
1.新建工程项目,填写项目存储路径和工程名,不要出现中文路径 2.添加已存在文件(可选),在[File name]下选择已经存在的工程项目,利用[Add]或[Add all]命令添加文件到新工程,点击 ...