4514: [Sdoi2016]数字配对 费用流
链接
https://www.lydsy.com/JudgeOnline/problem.php?id=4514
思路
EK直接贪心做
<0的时候加上剩余返回
二分图a->b的时候
把b->a也连接上
最后除2
整除和贪心可只知道它是对的
代码
#include <bits/stdc++.h>
#define ll long long
#define iter vector<int>::iterator
using namespace std;
const ll N=5e5+7;
ll read() {
ll x=0,f=1;char s=getchar();
for(;s>'9'||s<'0';s=getchar()) if(s=='-') f=-1;
for(;s>='0'&&s<='9';s=getchar()) x=x*10+s-'0';
return x*f;
}
ll n,S,T,a[2007],b[2007],c[2007];
struct node {
ll u,v,nxt,cap,cost;
}e[N<<1];
ll head[N<<1],tot=1;
void add_edge(ll u,ll v,ll cap,ll cost) {
e[++tot].v=v;
e[tot].u=u;
e[tot].cap=cap;
e[tot].cost=cost;
e[tot].nxt=head[u];
head[u]=tot;
}
void Add(ll u,ll v,ll cap,ll cost) {
// cout<<u<<" "<<v<<"\n";
add_edge(u,v,cap,cost);
add_edge(v,u,0,-cost);
}
ll dis[N];
ll vis[N],frm[N];
bool spfa() {
memset(dis,-0x3f,sizeof(dis));
memset(vis,0,sizeof(dis));
queue<int> q;
q.push(S);
dis[S]=0;
while(!q.empty()) {
ll u=q.front();
q.pop();
vis[u]=0;
for(ll i=head[u];i;i=e[i].nxt) {
ll v=e[i].v;
if(e[i].cap&&dis[v]<dis[u]+e[i].cost) {
dis[v]=dis[u]+e[i].cost;
frm[v]=i;
if(!vis[v]) {
vis[v]=1;
q.push(v);
}
}
}
}
return dis[T]>=-0x3f3f3f3f3f3f3f;
}
ll EK() {
ll flow=0,ans=0;
while(spfa()) {
ll now_flow=0x3f3f3f3f;
for(ll i=frm[T];i;i=frm[e[i].u])
now_flow=min(now_flow,(ll)e[i].cap);
for(ll i=frm[T];i;i=frm[e[i].u]) {
e[i].cap-=now_flow;
e[i^1].cap+=now_flow;
}
if(ans+1LL*now_flow*dis[T]<0) {
ll ok=ans/-dis[T];
return flow+ok;
} else
ans+=now_flow*dis[T],flow+=now_flow;
}
return flow;
}
set<int> dsr[250];
int main() {
freopen("a.in","r",stdin);
n=read();
for(ll i=1;i<=n;++i) a[i]=read();
for(ll i=1;i<=n;++i) b[i]=read();
for(ll i=1;i<=n;++i) c[i]=read();
S=1001,T=1002;
for(ll i=1;i<=n;++i) {
ll tmp=a[i];
for(ll j=2;j*j<=tmp;++j) {
if(tmp%j==0) dsr[i].insert(j);
while(tmp%j==0) tmp/=j;
}
if(tmp!=1) dsr[i].insert(tmp);
}
for(ll i=1;i<=n;++i) Add(S,i,b[i],0);
for(ll i=1;i<=n;++i) Add(i+n,T,b[i],0);
for(ll i=1;i<=n;++i) {
for(ll j=i+1;j<=n;++j) {
ll x=i,y=j;
if(a[x]<a[y]) swap(x,y);
if(a[y]==0) continue;
if(a[x]%a[y]==0) {
if(dsr[x].count(a[x]/a[y])) {
Add(x,y+n,0x3f3f3f3f3f3fLL,c[x]*c[y]);
Add(y,x+n,0x3f3f3f3f3f3fLL,c[x]*c[y]);
}
}
}
}
printf("%lld\n",EK()/2LL);
return 0;
}
4514: [Sdoi2016]数字配对 费用流的更多相关文章
- BZOJ 4514: [Sdoi2016]数字配对 [费用流 数论]
4514: [Sdoi2016]数字配对 题意: 有 n 种数字,第 i 种数字是 ai.有 bi 个,权值是 ci. 若两个数字 ai.aj 满足,ai 是 aj 的倍数,且 ai/aj 是一个质数 ...
- BZOJ.4514.[SDOI2016]数字配对(费用流SPFA 二分图)
BZOJ 洛谷 \(Solution\) 很显然的建二分图后跑最大费用流,但有个问题是一个数是只能用一次的,这样二分图两部分都有这个数. 那么就用两倍的.如果\(i\)可以向\(j'\)连边,\(j\ ...
- 【BZOJ4514】[Sdoi2016]数字配对 费用流
[BZOJ4514][Sdoi2016]数字配对 Description 有 n 种数字,第 i 种数字是 ai.有 bi 个,权值是 ci. 若两个数字 ai.aj 满足,ai 是 aj 的倍数,且 ...
- bzoj4514: [Sdoi2016]数字配对--费用流
看了一眼题目&数据范围,觉得应该是带下界的费用流 原来想拆点变成二分图,能配对的连边,跑二分图,可行性未知 后来看到另外一种解法.. 符合匹配要求的数要满足:质因子的个数相差为1,且两者可整除 ...
- 【BZOJ 4514】[Sdoi2016]数字配对 费用流
利用spfa流的性质,我直接拆两半,正解分奇偶(妙),而且判断是否整除且质数我用的是暴力根号,整洁判断质数个数差一(其他非spfa流怎么做?) #include <cstdio> #inc ...
- 图论(费用流):BZOJ 4514 [Sdoi2016]数字配对
4514: [Sdoi2016]数字配对 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 820 Solved: 345[Submit][Status ...
- BZOJ 4514: [Sdoi2016]数字配对
4514: [Sdoi2016]数字配对 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1606 Solved: 608[Submit][Statu ...
- 【BZOJ4514】【SDOI2016】数字配对 [费用流]
数字配对 Time Limit: 10 Sec Memory Limit: 128 MB[Submit][Status][Discuss] Description 有 n 种数字,第 i 种数字是 ...
- [bzoj4514]数字配对[费用流]
今年SDOI的题,看到他们在做,看到过了一百多个人,然后就被虐惨啦... 果然考试的时候还是打不了高端算法,调了...几天 默默地yy了一个费用流构图: 源连所有点,配对的点连啊,所有点连汇... 后 ...
随机推荐
- ReentrantLock源码(一)
一.简介. ReentrantLock 是一个互斥锁,在基本行为和机制上与synchonized一样,只不过synchonized用方法和声明访问了隐式的锁监视器,但是ReentrantLock 做了 ...
- SpringBoot+SpringCloud实现登录用户信息在微服务之间的传递
实现思路: 1:准备一个ThreadLocal变量,供线程之间共享. 2:每个微服务对所有过来的Feign调用进行过滤,然后从请求头中获取User用户信息,并存在ThreadLocal变量中. 3:每 ...
- 使用淘宝npm镜像
我们都知道淘宝大量采用了Nodej技术,所以它的镜像还是值得肯定的!更多相关可查看:淘宝NPM镜像 npm的官方镜像在国内访问是比较慢的,所以替代方法是使用淘宝npm镜像! 淘宝 NPM 镜像是一个完 ...
- linux环境下安装qt过程
linux(虚拟机fedora9)环境下安装qt的过程主要是按照下面几网页上的教程完成的. http://mobile.51cto.com/symbian-272869.htm http://www. ...
- 项目方说性能达到百万TPS,如何测试它的可信度?
项目方说性能达到百万TPS,如何测试它的可信度? 应用系统性能提升的关键在于运维端的接入管理模型(AAA,认证 Authentication.授权 Authorization.计费 Accountin ...
- LCA 最近公共祖先 (模板)
#include <iostream> #include <stdio.h> #include <cstring> #include <vector> ...
- JXNU暑期选拔赛
最小的数 Time Limit : 3000/1000ms (Java/Other) Memory Limit : 65535/32768K (Java/Other) Total Submissi ...
- Integer诡异特性
package 代码测试; public class ceshi { public static void main(String[] args) { Integer i1=100; Integer ...
- OGG 12.3中支持系统procedure复制的几点说明
如果需要同步系统级别的过程和package,则需要满足以下条件: 要求使用OGG12.3及以后的版本 需要使用oracle db12.2及以后的版本 需要使用集成抽取和集成投递 在DBA_GG_SUP ...
- kali linux下 hachcat安装
网上关于hachcat的简单使用方法介绍很多,然而却很少有在kali linux上的安装教程,找了好长时间,终于安装成功了,特此将中间借鉴的内容记录如下: #首先安装p7z,用于解压下载的p7z包 # ...