反正现在做题那么少就争取做一题写一题博客吧

看到题目发现数字种类不多,而且结合价值的要求可以容易地想到使用费用流

但是我们如果朴素地建图就会遇到一个问题,若\(i,j\)符合要求,那么给\(i,j\)连的应该是双向边,但双向边怎么跑网络流?

所以我们就要考虑怎么给边定向,我们稍加观察就会发现如果\(i,j\)合法,\(j,k\)也合法,那么\(i,k\)显然是不合法的(分四类情况讨论一下都是不合法的)

考虑那个整除出一个质数的条件,我们发现如果我们定义\(ct_i\)为\(a_i\)的所有质因数的指数之和

那么\(i,j\)满足条件当且仅当\((ct_i=ct_j+1\and a_j|a_i)\or(ct_j=ct_i+1\and a_i|a_j)\),这样一来我们就可以解释上面\(i,j,k\)三者的关系了

那么我们发现在有了\(ct_i\)之后我们就可以把数字分在两边,左边放奇数,右边放偶数,这样就构成了一个二分图,我们只从左边向右边连边即可

然后注意一下我们用EK跑最小费用最大流时每一次增广求出的最长路一定不会大于上一次增广求出的最长路,因此直接按顺序累加即可,注意最后一次的特判,详见代码

#include<cstdio>
#include<queue>
#include<iostream>
#define int long long
#define RI register int
#define CI const int&
using namespace std;
const int N=205,INF=1e18;
struct edge
{
int to,nxt,v,c;
}e[N*N<<2]; int n,head[N],s,t,cnt=1,a[N],b[N],c[N],ct[N];
inline int resolve(int x,int ret=0)
{
for (RI i=2;i*i<=x;++i) while (x%i==0) ++ret,x/=i; return ret+(x>1);
}
class Network_Flow
{
private:
queue <int> q; int dis[N],pre[N],lst[N],cap[N]; bool vis[N];
#define to e[i].to
inline bool SPFA(CI s,CI t)
{
RI i; for (i=s;i<=t;++i) dis[i]=-INF; q.push(s);
cap[s]=INF; dis[s]=0; while (!q.empty())
{
int now=q.front(); vis[now]=0; q.pop();
for (i=head[now];i;i=e[i].nxt)
if (e[i].v&&dis[now]+e[i].c>dis[to])
{
dis[to]=dis[pre[to]=now]+e[lst[to]=i].c;
if (cap[to]=min(cap[now],e[i].v),!vis[to]) q.push(to),vis[to]=1;
}
}
return dis[t]!=-INF;
}
#undef to
public:
inline void addedge(CI x,CI y,CI v,CI c)
{
e[++cnt]=(edge){y,head[x],v,c}; head[x]=cnt;
e[++cnt]=(edge){x,head[y],0,-c}; head[y]=cnt;
}
inline int solve(CI s,CI t)
{
int ret=0,ans=0; while (SPFA(s,t))
{
if (ret+dis[t]*cap[t]>=0)
{
ret+=dis[t]*cap[t]; ans+=cap[t];
for (int nw=t;nw!=s;nw=pre[nw])
e[lst[nw]].v-=cap[t],e[lst[nw]^1].v+=cap[t];
} else return ans+ret/(-dis[t]);
}
return ans;
}
}NF;
signed main()
{
RI i,j; for (scanf("%lld",&n),i=1;i<=n;++i) scanf("%lld",&a[i]),ct[i]=resolve(a[i]);
for (i=1;i<=n;++i) scanf("%lld",&b[i]); for (i=1;i<=n;++i) scanf("%lld",&c[i]);
for (t=n+1,i=1;i<=n;++i) if (ct[i]&1) NF.addedge(s,i,b[i],0); else NF.addedge(i,t,b[i],0);
for (i=1;i<=n;++i) if (ct[i]&1) for (j=1;j<=n;++j) if (!(ct[j]&1))
if ((ct[i]+1==ct[j]&&a[j]%a[i]==0)||(ct[j]+1==ct[i]&&a[i]%a[j]==0))
NF.addedge(i,j,INF,c[i]*c[j]); return printf("%lld",NF.solve(s,t)),0;
}

Luogu P4068 [SDOI2016]数字配对的更多相关文章

  1. p4068 [SDOI2016]数字配对

    传送门 分析 我们考虑对所有a[i]质因数分解,然后记cnt[i]为a[i]是由几个质数相乘得到的 然后我们建一个二分图,左面为所有cnt[i]为奇数的点,右面是为偶数的点 我们从源点向左面点连容量b ...

  2. 图论(费用流):BZOJ 4514 [Sdoi2016]数字配对

    4514: [Sdoi2016]数字配对 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 820  Solved: 345[Submit][Status ...

  3. BZOJ 4514: [Sdoi2016]数字配对 [费用流 数论]

    4514: [Sdoi2016]数字配对 题意: 有 n 种数字,第 i 种数字是 ai.有 bi 个,权值是 ci. 若两个数字 ai.aj 满足,ai 是 aj 的倍数,且 ai/aj 是一个质数 ...

  4. 【bzoj4514】: [Sdoi2016]数字配对 图论-费用流

    [bzoj4514]: [Sdoi2016]数字配对 好像正常的做法是建二分图? 我的是拆点然后 S->i cap=b[i] cost=0 i'->T cap=b[i] cost=0 然后 ...

  5. BZOJ 4514: [Sdoi2016]数字配对

    4514: [Sdoi2016]数字配对 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1606  Solved: 608[Submit][Statu ...

  6. 【BZOJ4514】[Sdoi2016]数字配对 费用流

    [BZOJ4514][Sdoi2016]数字配对 Description 有 n 种数字,第 i 种数字是 ai.有 bi 个,权值是 ci. 若两个数字 ai.aj 满足,ai 是 aj 的倍数,且 ...

  7. BZOJ4514——[Sdoi2016]数字配对

    有 n 种数字,第 i 种数字是 ai.有 bi 个,权值是 ci. 若两个数字 ai.aj 满足,ai 是 aj 的倍数,且 ai/aj 是一个质数, 那么这两个数字可以配对,并获得 ci×cj 的 ...

  8. bzoj4514 [Sdoi2016]数字配对

    Description 有 n 种数字,第 i 种数字是 ai.有 bi 个,权值是 ci. 若两个数字 ai.aj 满足,ai 是 aj 的倍数,且 ai/aj 是一个质数, 那么这两个数字可以配对 ...

  9. BZOJ4514[Sdoi2016]数字配对——最大费用最大流

    题目描述 有 n 种数字,第 i 种数字是 ai.有 bi 个,权值是 ci. 若两个数字 ai.aj 满足,ai 是 aj 的倍数,且 ai/aj 是一个质数, 那么这两个数字可以配对,并获得 ci ...

随机推荐

  1. 初学Python:Python的发展历史及下载安装

    Python作为一种计算机程序设计语言,自20世纪90年代初诞生至如今被人们逐渐悉知,经过版本更新以及功能添加,已广泛应用于各种独立的.大型项目的开发.Python 已经成为最受欢迎的程序设计语言之一 ...

  2. 《细说PHP》第四版 样章 第23章 自定义PHP接口规范 5

    23.3  接口的安全控制规范 23.2节的示例实现了一个简单接口,但是这个接口此时是在“裸奔”的.因为这个接口所有人都可以请求,不仅我们的客户端可以正常访问数据,如果有人使用如fiddler.wir ...

  3. IT兄弟连 Java语法教程 流程控制语句 分支结构语句5

    5  switch-case条件语句 Java中的第二种分支控制语句时switch语句,switch语句提供了多路支持,因此可以使程序在多个选项中进行选择.尽管一系列嵌套if语句可以执行多路测试,然而 ...

  4. pytest框架之fixture前置和后置

    一.conftest.py 定义公共的fixture,多个测试类中都可以调用 pytest提供了conftest.py文件,可以将fixture定义在此文件中 运行测试用例时,不需要去导入这个文件,会 ...

  5. 基于Node 的http转发demo,项目中请使用express+http-proxy-middleware

    var http = require("http"); var data = ''; function getData() { const options = { host: 'w ...

  6. Day02stu

    ⦁ 环境搭建之接口测试工具Jmeter搭建 1) 什么是Jmeter? Jmeter是一个接口测试工具,基于Java开发,是是Apche公司使用Java平台开发的一款测试工具. 运行环境需要配置:JD ...

  7. 多线程通信的两种方式? (可重入锁ReentrantLock和Object)

    (一)Java中线程协作的最常见的两种方式: (1)利用Object的wait().notify()和notifyAll()方法及synchronized (2)使用Condition.Reentra ...

  8. jQuery 源码分析(十八) ready事件详解

    ready事件是当DOM文档树加载完成后执行一个函数(不包含图片,css等),因此它的触发要早于load事件.用法: $(document).ready(fun) ;fun是一个函数,这样当DOM树加 ...

  9. EfficientNet: Rethinking Model Scaling for Convolutional Neural Networks

    增加模型精度的方法有增加网络的深度,特征图的通道数以及分辨率(如下图a-d所示).这篇文章研究了模型缩放,发现仔细平衡网络的深度.宽度和分辨率可以获得更好的性能(下图e).在此基础上,提出了一种新的缩 ...

  10. [转]Apache Doris资料汇总

    参考资料 https://www.toutiao.com/i6709706901770207748/?tt_from=weixin&utm_campaign=client_share& ...