题目:http://acm.hdu.edu.cn/showproblem.php?pid=5594

完全不会啊TAT。。

其实官方题解已经说的很清楚了。。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <queue>
#include <cmath>
#define rep(i,l,r) for (int i=l;i<=r;i++)
#define down(i,l,r) for (int i=l;i>=r;i--)
#define clr(x,y) memset(x,y,sizeof(x))
#define maxn 405
#define ll long long
#define inf int(1e9)
using namespace std;
int a[maxn],head[maxn],cur[maxn],uu[maxn],odd[maxn],even[maxn];
int tot,t,s;
struct data{int obj,pre,c;
}e[];
int read(){
int x=,f=; char ch=getchar();
while (!isdigit(ch)){ if (ch=='-') f=-; ch=getchar();}
while (isdigit(ch)) {x=x*+ch-''; ch=getchar();}
return x*f;
}
void insert(int x,int y,int z){
e[++tot].obj=y; e[tot].c=z; e[tot].pre=head[x]; head[x]=tot;
e[++tot].obj=x; e[tot].c=; e[tot].pre=head[y]; head[y]=tot;
}
bool bfs(){
queue<int >q; clr(uu,-); uu[s]=; q.push(s);
while (!q.empty()){
int u=q.front(); q.pop();
for (int j=head[u];j;j=e[j].pre){
int v=e[j].obj;
if (uu[v]==-&&e[j].c) uu[v]=uu[u]+,q.push(v);
}
}
if (uu[t]==-) return ;
return ;
}
int dfs(int x,int mx){
if (x==t) return mx;
int used=;
for (int j=cur[x];j;j=e[j].pre){
int v=e[j].obj;
if (uu[v]==uu[x]+){
int w=dfs(v,min(e[j].c,mx-used));
used+=w; e[j].c-=w; e[j^].c+=w;
if (e[j].c) cur[x]=j;
if (used==mx) return mx;
}
}
if (!used) uu[x]=-;
return used;
}
int dinic(){
int ans=;
while (bfs()){
rep(i,s,t) cur[i]=head[i];
ans+=dfs(s,inf);
}
return ans;
}
bool jud(int x){
int s=int(sqrt(x))+;
rep(i,,s) if (x%i==) return ;
return ;
}
bool work(){
tot=; clr(head,); clr(odd,); clr(even,); //clr(one,0);
int odd_cnt=,even_cnt=,one_in=,one_out=,one_cnt=;
int n=read();
rep(i,,n) {
a[i]=read();
if (a[i]>&&(a[i]&)) odd[++odd_cnt]=a[i];
else if (a[i]>&&((a[i]&)==)) even[++even_cnt]=a[i];
else one_cnt++;
}
int need=even_cnt-odd_cnt;
rep(i,,one_cnt){
if (need) {
one_in++;
odd[++odd_cnt]=;
need--;
}
else one_out++;
}
if ((!one_in && (one_out==||one_out==)) || need!=) return ;
s=; t=*odd_cnt+;
int cnt=odd_cnt;
rep(i,,cnt) insert(s,i,);
rep(i,,cnt) insert(i+cnt,t,);
rep(i,,cnt) rep(j,,cnt){
if (jud(odd[i]+even[j])) {
if (odd[i]==&&one_out>) insert(i,j+cnt,);
else insert(i,j+cnt,);
}
}
if (dinic()!=*cnt) return ;
return ;
}
int main(){
int T=read();
while (T--){
if (work()) puts("YES"); else puts("NO");
}
return ;
}

BC#65 T5 ZYB's Prime的更多相关文章

  1. hdu 5594 ZYB's Prime 最大流

    ZYB's Prime Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://acm.hdu.edu.cn/showproblem.php?pid=5 ...

  2. BC 65 ZYB's Premutation (线段树+二分搜索)

    题目简述:有一个全排列,一直每个前缀区间的逆序对数,还原这个排列. fi记录逆序对数,pi记录该位置数值,则k=fi-f(i-1)表示前i-1个数比pi大的数的个数,那么只要在剩余元素求出按大小顺序第 ...

  3. BestCoder Round #65 (ZYB's Premutation)

    ZYB's Premutation Accepts: 220 Submissions: 983 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: ...

  4. BestCoder Round #65 (ZYB's Game)

    ZYB's Game Accepts: 672 Submissions: 1207 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536 ...

  5. BestCoder Round #65 (ZYB's Biology)

    ZYB's Biology Accepts: 848 Submissions: 1199 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 13 ...

  6. BestCoder Round #65 HDOJ5592 ZYB&#39;s Premutation(树状数组+二分)

    ZYB's Premutation Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 131072/131072 K (Java/Othe ...

  7. BC 65 game

    主持人一直某个数字在1到n范围,假设甲乙已经知道,甲先猜乙后,都采用最优策略,主持人说偏大还是偏小,不断缩小范围,问最后乙能会获胜的X的取值的个数. 如果n为奇数,那么仅当x=n/2乙必然获胜,若为奇 ...

  8. OpenSSL CVE-2022-0778漏洞问题复现与非法证书构造

    本文介绍CVE-2022 0778漏洞及其复现方法,并精心构造了具有一个非法椭圆曲线参数的证书可以触发该漏洞. 本博客已迁移至CatBro's Blog,那是我自己搭建的个人博客,欢迎关注.本文链接 ...

  9. FF D8 FF FE 00 24 47 00转图片

    String[] img = "FF D8 FF FE 00 24 47 00 9D 0C 00 00 00 00 00 00 00 00 00 00 00 00 00 00 00 F0 0 ...

随机推荐

  1. iOS UIAlertController在Tableview中显示缓慢,迟钝,延迟

    在UITableViewCell中弹窗Alert延迟.在cellForRow中:cell.selectionStyle = UITableViewCellSelectionStyleNone; 或者在 ...

  2. 开发中关于Git那些事(续:Git变基)

    其实上一篇写的内容仅仅是Git的冰山一角,如果你认为Git就是简简单单的几行命令,那只能说明你还没有真正了解Git这个强大的内容寻址文件系统.这篇文章,还是接着介绍一些实用但是很少有人知晓的一些命令, ...

  3. 高度-宽度关系,同一div、不同div高度与宽度关系控制函数

    //对象1的高度等于对象2的高度n倍,调用方法:Ht1DivideHt2('#div2','#div1',3)//div2的高度是div1高度的3倍function Ht1DivideHt2(obj1 ...

  4. 三十天学不会TCP,UDP/IP网络编程-TraceRoute的哲学

    新年快乐,继续来部分粘贴复制我的这一系列文章啦,如果对和程序员有关的计算机网络知识,和对计算机网络方面的编程有兴趣,欢迎去gitbook(https://www.gitbook.com/@rogerz ...

  5. springBoot系列教程07:异常捕获

    发生异常是很正常的事,异常种类也是千奇百怪,发生异常并不可怕,只要正确的处理,并正确的返回错误信息并无大碍,如果不进行捕获或者处理,分分钟服务器宕机是很正常的事 所以处理异常时,最基本的要求就是发生异 ...

  6. requireJS教程

    目录[-] 使用 RequireJS 优化 Web 应用前端 AMD 简介 传统 JavaScript 代码的问题 AMD 的引入 清单 1. AMD 规范 RequireJS 简介 实战 Requi ...

  7. promise入门demo

    <!DOCTYPE html> <html> <head> <title></title> </head> <body&g ...

  8. 未能解析此远程名称: 'www.***.com' 解决办法

    今天发布网站的时候在本机IIS上调试完全没问题,但发布到远程服务器上却出现了如下错误: 原因:远程服务器的DNS解析失败,不能对类似于'www.***.com' 的网站解析,所以我们要进行手动解析. ...

  9. Linux下查找文件的方法

    在Linux环境下查找一个文件的方法:find 路径 -name 'filename',filename不清楚全名的话可以用*号进行匹配,如“tomcat.*”.如果不清楚路径的话可以用"/ ...

  10. Netty之多用户的聊天室(三)

    Netty之多用户的聊天室(三) 一.简单说明 笔者有意将Netty做成一个系列的文章,因为笔者并不是一个善于写文章的人,而且笔者学习很多技术一贯的习惯就是敲代码,很多东西敲着敲着就就熟了,然后再进行 ...