BZOJ 4205: 卡牌配对
来源:互联网 发布:c 语言 多线程共享数据 编辑:程序博客网 时间:2024/05/23 19:14
难以置信地1A了 挺开心。。
这道题建图太神了 一开始就觉得要在素数搞一些很哲♂学的东西然而想不到QAQ
考虑到200以内只有46个素数 我们可以用一个二元组来表示其实两种属性值的质因子
比如(a,b)表示a是A的质因子 b是B的质因子
其次题目说至多一组属性互质 也就是至少两组有不互质(有公共质因子)
其实上面的二元组应该加多一元 表示是哪两个属性 也就是AB AC BC三种可能
所以就可以挺顺理成章地建图 源点 连到X类卡牌 连到中间至多3*46*46个表示质因子的点 再连到Y类卡牌 and 汇点
(偷黄学长的分析一波)因为200<2*3*5*7,所以两侧每个点至多连出3*3*3条边。于是我们构成了一个70000个点,2000000条边的网络流。
于是乎我们用dicnic愉快地跑完这个分层图 如下图
2333
我觉得我这个代码很优美啊~也没过100行
#include<bits/stdc++.h>using namespace std;const int N=70005,M=4000005,inf=1e9+7;inline int read(){ char ch=getchar(); int x=0,f=1; while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();} while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();} return x*f;}struct edge{int y,c,next;}a[M]; int len,first[N];void ins(int x,int y,int c){ a[++len]=(edge){y,c,first[x]},first[x]=len; a[++len]=(edge){x,0,first[y]},first[y]=len;}queue<int>q; int h[N],t,p[50],pl,n,m; bool v[201];bool bfs(){ for(int i=0;i<=t;i++)h[i]=0; q.push(0); h[0]=1; while(!q.empty()){ int x=q.front(); q.pop(); for(int k=first[x];k;k=a[k].next){ int y=a[k].y; if(a[k].c && !h[y]) h[y]=h[x]+1,q.push(y); } } return h[t]>0;}int dfs(int x,int f){ if(x==t)return f; int de=0; for(int k=first[x];k;k=a[k].next){ int y=a[k].y; if(h[y]==h[x]+1){ int u=dfs(y,min(a[k].c,f-de)); de+=u,a[k].c-=u,a[k^1].c+=u; if(f==de)return f; } } if(!de)h[x]=-1; return de;}vector<int>P[201];void getp(){ for(int i=2;i<200;i++){ if(!v[i])p[++pl]=i; for(int j=1;j<=pl && i*p[j]<200;j++) { v[i*p[j]]=1; if(i%p[j]==0)break; } } for(int i=2;i<=200;i++)for(int j=1;j<=pl && p[j]<=i;j++) if(i%p[j]==0)P[i].push_back(j);}inline int u(int x1,int x2,int k){return n+m+(x1-1)*pl+x2+k*pl*pl;}int main(){ n=read(),m=read(); int i,j,k,x,y,z; getp(); t=n+m+3*pl*pl+1; len=1; for(i=1;i<=n;i++) { ins(0,i,1); x=read(),y=read(),z=read(); for(j=0;j<P[x].size();j++)for(k=0;k<P[y].size();k++) ins(i,u(P[x][j],P[y][k],0),1); for(j=0;j<P[x].size();j++)for(k=0;k<P[z].size();k++) ins(i,u(P[x][j],P[z][k],1),1); for(j=0;j<P[y].size();j++)for(k=0;k<P[z].size();k++) ins(i,u(P[y][j],P[z][k],2),1); } for(i=1;i<=m;i++) { ins(i+n,t,1); x=read(),y=read(),z=read(); for(j=0;j<P[x].size();j++)for(k=0;k<P[y].size();k++) ins(u(P[x][j],P[y][k],0),i+n,1); for(j=0;j<P[x].size();j++)for(k=0;k<P[z].size();k++) ins(u(P[x][j],P[z][k],1),i+n,1); for(j=0;j<P[y].size();j++)for(k=0;k<P[z].size();k++) ins(u(P[y][j],P[z][k],2),i+n,1); } int ans=0; while(bfs())ans+=dfs(0,inf); printf("%d\n",ans); return 0;}
2 1
- BZOJ 4205: 卡牌配对
- BZOJ P4205[FJ2015集训]卡牌配对
- BZOJ-4205 卡牌配对 最大流+线性筛+神建模
- BZOJ 1786 配对Pair
- BZOJ 1237 配对
- BZOJ 1237: [SCOI2008]配对
- bzoj 1237 [SCOI2008]配对
- BZOJ4205卡牌配对
- bzoj4205: 卡牌配对
- BZOJ4205: 卡牌配对
- Bzoj4205卡牌配对
- bzoj4205 卡牌配对
- bzoj4205 卡牌配对
- 【bzoj 1786】[Ahoi2008]Pair 配对
- 【bzoj 1237】[SCOI2008]配对 dp
- BZOJ 4514 [Sdoi2016]数字配对
- 【BZOJ】4514: [Sdoi2016]数字配对
- bzoj 1786 [Ahoi2008]Pair 配对
- C#中的DateTime如何解析两位数字的年份
- java源码解读之LinkedList------jdk 1.7
- 小米嵌入式面经(转发)
- 团队合作要服务化,拒绝代劳
- Netcat 工具用于设备间文件传输
- BZOJ 4205: 卡牌配对
- 异常处理——毕向东Java基础教程学习笔记
- IEnumerable vs IQueryable
- 内部类--毕向东Java基础教程学习笔记
- 管理类书籍读后感(散落在外的那些)
- 初学多线程练习1--电子时钟
- MySQL中删除重复数据只保留一条
- Weak & Strong Alias(强弱别名)
- js组合继承(原型链+对象冒充)