读完题后不难看出本题是个网络流模型,源点流出的总流量为\(k\),源点向每个和总部直接联系的间谍连边,每个间谍向其能传递的间谍连容量为\(m\)的边,能与德军情报部进行联系的间谍向汇点连容量为\(inf\)的边,若最大流为\(k\),则存在可行的方案。

处理可靠程度最大时,考虑用费用流解决,将每条边的安全程度看作边的费用,进行最大费用最大流即可,注意反向边的费用应该为原边费用的倒数,实现时把平时的费用流的一些部分改为乘法即可。

一些细节看代码吧。

\(code:\)

#include<bits/stdc++.h>
#define maxn 1000010
#define inf 1000000000
#define eps 1e-12
using namespace std;
typedef long double ld;
template<typename T> inline void read(T &x)
{
x=0;char c=getchar();bool flag=false;
while(!isdigit(c)){if(c=='-')flag=true;c=getchar();}
while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
if(flag)x=-x;
}
int n,k,S,s,t;
ld ans=1;
int am[maxn],d[maxn];
ld as[maxn],dis[maxn];
bool vis[maxn];
char str[30];
struct edge
{
int to,nxt,v;
ld c;
}e[maxn];
int head[maxn],edge_cnt=1;
void add(int from,int to,int val,ld cost)
{
e[++edge_cnt]=(edge){to,head[from],val,cost};
head[from]=edge_cnt;
e[++edge_cnt]=(edge){from,head[to],0,1/cost};
head[to]=edge_cnt;
}
bool spfa()
{
queue<int> q;
for(int i=s;i<=t;++i) dis[i]=0,vis[i]=d[i]=0;
q.push(s),dis[s]=1,d[s]=1,vis[s]=true;
while(!q.empty())
{
int x=q.front();
q.pop(),vis[x]=false;
for(int i=head[x];i;i=e[i].nxt)
{
int y=e[i].to,v=e[i].v;
ld c=e[i].c;
if(v&&dis[x]*c-dis[y]>eps)
{
dis[y]=dis[x]*c,d[y]=d[x]+1;
if(!vis[y]) vis[y]=true,q.push(y);
}
}
}
return dis[t];
}
int dfs(int x,int lim)
{
if(x==t) return lim;
int res=lim,flow;
for(int i=head[x];i;i=e[i].nxt)
{
int y=e[i].to,v=e[i].v;
ld c=e[i].c;
if(!v||dis[y]-dis[x]*c>eps||d[y]!=d[x]+1) continue;
if(flow=dfs(y,min(res,v)))
{
res-=flow;
e[i].v-=flow;
e[i^1].v+=flow;
if(!res) break;
}
}
return lim-res;
}
ld qp(ld x,int y)
{
ld v=1;
while(y)
{
if(y&1) v*=x;
x*=x,y>>=1;
}
return v;
}
void dinic()
{
int flow,sum=0;
while(spfa()>eps)
while(flow=dfs(s,inf))
sum+=flow,ans*=qp(dis[t],flow);
if(sum!=k) ans=0;
}
void print()
{
sprintf(str,"%.15Lf",ans);
int pos=0,cnt=0;
while(cnt<5)
{
char ch=str[pos++];
if((ch!='0'&&ch!='.')||cnt) cnt++;
}
if(str[pos]>='5') str[pos-1]+=1;
str[pos]=0;
for(int i=pos;i>=0;--i)
{
if(str[i]=='.') break;
if(str[i]>'9') str[i-1]++,str[i]='0';
}
printf("%s",str);
}
int main()
{
read(n),read(k),S=n+1,t=S+1,add(s,S,k,1);
for(int i=1;i<=n;++i) scanf("%Lf",&as[i]);
for(int i=1;i<=n;++i) read(am[i]);
for(int i=1;i<=n;++i)
if(am[i]&&as[i]>eps)
add(S,i,am[i],as[i]);
for(int i=1;i<=n;++i)
{
int v;
read(v);
if(v) add(i,t,inf,1);
}
while(1)
{
int x,y,m;
ld s;
read(x),read(y);
if(x==-1&&y==-1) break;
scanf("%Lf",&s),read(m);
add(x,y,m,s),add(y,x,m,s);
}
dinic();
if(ans>eps) print();
else puts("0");
return 0;
}

题解 洛谷 P5814 【[CTSC2001]终极情报网】的更多相关文章

  1. BZOJ2542: [Ctsc2001]终极情报网

    题解: 乘积最小只需要取对数.然后反向边就变成1/c,而不是-c了. 精度问题搞得我已经我想说什么了... 贴一份网上的pascal 代码: type ss=record x,y,c,r,next:l ...

  2. 图论(网络流):[CTSC2001]终极情报网

    [CTSC2001]终极情报网 [题目描述] 在最后的诺曼底登陆战开始之前,盟军与德军的情报部门围绕着最终的登陆地点展开了一场规模空前的情报战. 这场情报战中,盟军的战术是利用那些潜伏在敌军内部的双重 ...

  3. bzoj 2542: [Ctsc2001]终极情报网 费用流

    题目链接 2542: [Ctsc2001]终极情报网 Time Limit: 1 Sec  Memory Limit: 128 MBSubmit: 321  Solved: 125[Submit][S ...

  4. 题解 洛谷P5018【对称二叉树】(noip2018T4)

    \(noip2018\) \(T4\)题解 其实呢,我是觉得这题比\(T3\)水到不知道哪里去了 毕竟我比较菜,不大会\(dp\) 好了开始讲正事 这题其实考察的其实就是选手对D(大)F(法)S(师) ...

  5. 题解 洛谷 P3396 【哈希冲突】(根号分治)

    根号分治 前言 本题是一道讲解根号分治思想的论文题(然鹅我并没有找到论文),正 如论文中所说,根号算法--不仅是分块,根号分治利用的思想和分块像 似却又不同,某一篇洛谷日报中说过,分块算法实质上是一种 ...

  6. 题解-洛谷P5410 【模板】扩展 KMP(Z 函数)

    题面 洛谷P5410 [模板]扩展 KMP(Z 函数) 给定两个字符串 \(a,b\),要求出两个数组:\(b\) 的 \(z\) 函数数组 \(z\).\(b\) 与 \(a\) 的每一个后缀的 L ...

  7. 题解-洛谷P4229 某位歌姬的故事

    题面 洛谷P4229 某位歌姬的故事 \(T\) 组测试数据.有 \(n\) 个音节,每个音节 \(h_i\in[1,A]\),还有 \(m\) 个限制 \((l_i,r_i,g_i)\) 表示 \( ...

  8. 题解-洛谷P4724 【模板】三维凸包

    洛谷P4724 [模板]三维凸包 给出空间中 \(n\) 个点 \(p_i\),求凸包表面积. 数据范围:\(1\le n\le 2000\). 这篇题解因为是世界上最逊的人写的,所以也会有求凸包体积 ...

  9. 题解-洛谷P4859 已经没有什么好害怕的了

    洛谷P4859 已经没有什么好害怕的了 给定 \(n\) 和 \(k\),\(n\) 个糖果能量 \(a_i\) 和 \(n\) 个药片能量 \(b_i\),每个 \(a_i\) 和 \(b_i\) ...

随机推荐

  1. Linux 初始化系统 SystemV Upstart

    System V 特点 缺点: 启动时间长,init是串行启动,只有前一个进程启动完,才会启动下一个进程 启动脚本复杂,init只是执行启动脚本,不管其他事情,脚本需要自己处理各种情况,这往往使得脚本 ...

  2. IDEA SonarLint安装及使用

    SonarLint插件安装IDEA菜单栏选择File->Settings,左边栏选择Plugins 在线安装选择Browse repositories,搜索Sonar,选择SonarLint进行 ...

  3. java基础 内部类详解

    什么是内部类? 1.内部类也是一个类: 2.内部类位于其他类声明内部. 内部类的常见类型 1.成员内部类 2.局部内部类 3.匿名内部类 4.静态内部类 简单示例 /** * 外部类 * */ pub ...

  4. @Inherited 注解的作用

    @Inherited 用于放在注解上,例如 @Inherited @Documented @Target(ElementType.TYPE) public @interface InheritedAn ...

  5. layer.open弹框中的表单数据无法获取

    layer.open弹框中的表单数据无法获取 表单数据模板 layer.open() 页面效果: 当点击确定后,radio和textarea获取的值总是为空,解决办法: var setPriCustB ...

  6. JavaWEB实现qq邮箱发送验证码——qq1692700664

    // 随机验证码public String achieveCode() { String[] beforeShuffle = new String[] { "2", "3 ...

  7. 基于TCP与UDP协议的socket通信

    基于TCP与UDP协议的socket通信 C/S架构与初识socket 在开始socket介绍之前,得先知道一个Client端/服务端架构,也就是 C/S 架构,互联网中处处充满了 C/S 架构(Cl ...

  8. Kubernetes 中 搭建 EFK 日志搜索中心

    简介 Elastic 官方已经发布了Elasticsearch Operator ,简化了 elasticsearch 以及 kibana的部署与升级,结合 fluentd-kubernetes-da ...

  9. python之浅谈数据类型

    什么是数据类型 ​ 数据类型指的就是变量值的不同类型,姓名可能是一种数据类型.年龄可能是一种数据类型.爱好可能又是另一种数据类型,至于是什么数据类型我们将在下一章详细说明. 如何对数据分类 变量的是用 ...

  10. 2020阿里巴巴官方最新Redis开发规范!

    本文主要介绍在使用阿里云Redis的开发规范,从下面几个方面进行说明. 键值设计 命令使用 客户端使用 相关工具 通过本文的介绍可以减少使用Redis过程带来的问题. 一.键值设计 1.key名设计 ...