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

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

但是我们如果朴素地建图就会遇到一个问题,若\(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. docker搭建zookeeper集群

    1.在官网拉取镜像 docker pull zookeeper 2.根据镜像启动zookeeper容器 docker run -itd --name zookeeper1 -h zookeeper1 ...

  2. 抓包工具之fiddler实战1-基本设置

    Fiddler概述 百度搜索fiddler能找到官网网站,百度软件中心也提供了下载,本人去下载了基本和官网的版本一致,但还是建议大家下载软件一定去官网进行下载. Fiddler是干什么的 在百度百科里 ...

  3. 《细说PHP》第四版 样章 第18章 数据库抽象层PDO 8-1

    18.6.5  获取数据 PDO的数据获取方法与其他数据库扩展非常类似,只要成功执行SELECT查询,都会有结果集对象生成.不管使用PDO对象中的query()方法,还是使用prepare()和exe ...

  4. mongodb简单监控

    一.mongodb提供了自己的监控工具 mongostat mongostat是mongodb自带的状态检测工具,在命令行下使用.它会间隔固定时间获取mongodb的当前运行状态,并输出.如果你发现数 ...

  5. docker registry 删除镜像 垃圾回收

    操作步骤 通过环境变量修改默认配置,允许删除 获取image的sha值 进入registry容器中,执行垃圾回收 删除残留目录 #环境变量 REGISTRY_STORAGE_DELETE_ENABLE ...

  6. Web前端——Html常用标签及属性

    html 常用的标题等标签就不记录了,只记录一下比较少见的标签以及属性 表格 table td 单元格 tr 表的行 th 表头 td或th可以下面的两个属性达到跨行或跨列 表格跨行 rowspan ...

  7. String与StringBuilder性能比对

    //String与StringBuilder性能比对package seday01;/** * String修改字符串带来的性能开销 * @author xingsir * */public clas ...

  8. C#中巧用Lambda表达式实现对象list进行截取

    场景 有一个对象的list,每个对象有唯一的属性Id,并且是从1递增,现在要根据此Id属性进行截取. 其中DataTreeNode 实现 Global.Instance.PrepareCompareD ...

  9. layui 学习笔记一:layui table 查询、新增、编辑、删除

    一.table数据的呈现(对应查询) 页面代码: @{ ViewBag.Title = "TableGrid"; } @section styles{ <link href= ...

  10. JS的with关键字到底是什么?

    with关键字 with在JS中通常被当做重复引用同一个对象多个属性的快捷方式. var obj = { a: 1, b: 2, c: 3 }; // 重复引用obj进行属性赋值 obj.a = 3; ...