BZOJ 4819 新生舞会
第一句话:算出3.363636的孩子啊,你跑错流种了。
貌似上一篇我讲SDOI出原题?嘿还真是!
半个月前有一个叫WG的男人给我们搞过一场考试... ...
里面有一道题叫做保密... ...SDOI2011... ...
对于每个点求ΣS/ΣT最值然后转跑浮点数网络流... ...
是不是感觉我在讲这道题的正解... ...没错正解就是这样... ...所以我说是原题吧。
跟 保密 是一样的思路。求Σa/Σb的最值,可以用二分答案来做。
Σa/Σb的最大值怎么求呢?设一个当前答案ans;
显然如果有方案使Σa/Σb>ans则答案更优。
把式子化一下得:
Σa>Σb*ans;
即
Σa-Σb*ans>0;
即
Σ(a-b*ans)>0;
每次重新建图或者带着二分的mid跑最大费用最大流即可。
没错跑出3.363636是跑的最小费用最大流... ...别问我怎么知道的。
用KM算法算最大权匹配跑的更快奈何我不会。好在它是一个左右相等的完全二分图,用费用流做也是可以的。
出题人良心发现没有卡费用流真是资磁。
还有就是浮点数二分答案跟整数不一样,要用eps。
有哪位大佬能够告诉我 为什么我的费用流 跑的那么慢吗?
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <cstring>
#define LL long long
#define dob double
using namespace std; const int N = 110;
const int M = 100010;
const dob eps = 1e-7;
const dob dobInf = 19260817.2333;
struct Node{int to,C;double val;int next;}E[M];
int head[M],tot,n,S,T,up[M],In[M];
dob a[N][N],b[N][N],far[M],Ans; int gi()
{
int x=0,res=1;char ch=getchar();
while(ch<'0' || ch>'9'){if(ch=='-')res=-res;ch=getchar();}
while(ch>='0'&&ch<='9')x=x*10+ch-48,ch=getchar();
return x*res;
} inline void link(int u,int v,int c,double val)
{
E[++tot]=(Node){v,c,val,head[u]};
head[u]=tot;
E[++tot]=(Node){u,0,-val,head[v]};
head[v]=tot;
} inline void slink(double Eps)
{
memset(head,0,sizeof(head));tot=1;
for(int i=1;i<=n;++i)
{
link(S,i,1,0.0000000);
link(i+n,T,1,0.0000000);
for(int j=1;j<=n;++j)
link(i,n+j,1,(a[i][j])-Eps*b[i][j]);
}
} inline bool spfa()
{
for(int i=0;i<=T;++i)
far[i]=-dobInf;
memset(up,0,sizeof(up));
queue<int>Q;Q.push(S);
far[S]=0.000000;In[S]=1;
while(!Q.empty())
{
int x=Q.front();Q.pop();
for(int e=head[x];e;e=E[e].next)
if(E[e].C && far[x]+E[e].val>far[E[e].to])
{
up[E[e].to]=e;
far[E[e].to]=far[x]+E[e].val;
if(!In[E[e].to])
Q.push(E[e].to),In[E[e].to]=1;
}
In[x]=0;
}
if(!up[T])return false;
for(int i=T;i^S;i=E[up[i]^1].to)
E[up[i]].C--,E[up[i]^1].C++;
Ans-=far[T];return true;
} inline bool check()
{
Ans=0.00;
while(spfa())
if(Ans>=0.00)return true;
return Ans>=0.00;
} int main()
{
n=gi();S=2*n+1,T=2*n+2;
for(int i=1;i<=n;++i)
for(int j=1;j<=n;++j)
scanf("%lf",&a[i][j]);
for(int i=1;i<=n;++i)
for(int j=1;j<=n;++j)
scanf("%lf",&b[i][j]);
dob l=0.0,r=10000,ans=r;
while(r-l>eps)
{
double mid=(l+r)/2.0;
slink(mid);
if(check())r=mid,ans=mid;
else l=mid;
}
printf("%.6lf\n",ans);
return 0;
}
BZOJ 4819 新生舞会的更多相关文章
- 【BZOJ 4819】 4819: [Sdoi2017]新生舞会 (0-1分数规划、二分+KM)
4819: [Sdoi2017]新生舞会 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 601 Solved: 313 Description 学校 ...
- bzoj 4819: [Sdoi2017]新生舞会
Description 学校组织了一次新生舞会,Cathy作为经验丰富的老学姐,负责为同学们安排舞伴.有n个男生和n个女生参加舞会 买一个男生和一个女生一起跳舞,互为舞伴.Cathy收集了这些同学之间 ...
- [BZOJ 4819] [SDOI 2017] 新生舞会
Description 学校组织了一次新生舞会,Cathy作为经验丰富的老学姐,负责为同学们安排舞伴. 有 \(n\) 个男生和 \(n\) 个女生参加舞会买一个男生和一个女生一起跳舞,互为舞伴. C ...
- 4819: [Sdoi2017]新生舞会(分数规划)
4819: [Sdoi2017]新生舞会 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1031 Solved: 530[Submit][Statu ...
- BZOJ4819 新生舞会
4819: [Sdoi2017]新生舞会 Time Limit: 10 Sec Memory Limit: 128 MB Description 学校组织了一次新生舞会,Cathy作为经验丰富的老学 ...
- 【BZOJ4819】新生舞会(分数规划,网络流)
[BZOJ4819]新生舞会(分数规划,网络流) 题面 BZOJ Description 学校组织了一次新生舞会,Cathy作为经验丰富的老学姐,负责为同学们安排舞伴.有n个男生和n个女生参加舞会 买 ...
- [BZOJ4819][SDOI2017]新生舞会(分数规划+费用流,KM)
4819: [Sdoi2017]新生舞会 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1097 Solved: 566[Submit][Statu ...
- bzoj4819 [Sdoi2017]新生舞会
Description 学校组织了一次新生舞会,Cathy作为经验丰富的老学姐,负责为同学们安排舞伴.有n个男生和n个女生参加舞会买一个男生和一个女生一起跳舞,互为舞伴.Cathy收集了这些同学之间的 ...
- [Sdoi2017]新生舞会 [01分数规划 二分图最大权匹配]
[Sdoi2017]新生舞会 题意:沙茶01分数规划 貌似\(*10^7\)变成整数更科学 #include <iostream> #include <cstdio> #inc ...
随机推荐
- 【深度学习系列】用PaddlePaddle和Tensorflow进行图像分类
上个月发布了四篇文章,主要讲了深度学习中的"hello world"----mnist图像识别,以及卷积神经网络的原理详解,包括基本原理.自己手写CNN和paddlepaddle的 ...
- C# tostring
GUID: 即Globally Unique Identifier(全球唯一标识符) 也称作 UUID(Universally Unique IDentifier) . GUID是一个通过特定算法产生 ...
- 校验ISBN的方法
国际标准书号(International Standard Book Number,ISBN:是国际通用的图书或独立的出版物(除定期出版的期刊)代码.出版社可以通过国际标准书号清晰的辨认所有非期刊书 ...
- ShoneSharp语言(S#)的设计和使用介绍系列(1)— 开篇
ShoneSharp语言(S#)的设计和使用介绍 系列(1)- 开篇 作者:Shone 声明:原创文章欢迎转载,但请注明出处,https://www.cnblogs.com/ShoneSharp. 一 ...
- 爬起点小说 day01
先介绍下我自己爬起点小说的思路: 1.爬取所有的类型列表 a.链接存redis中 类型表:novel_list 具体每一种类型:bnovel_all_list(把novel_list和bnovel_l ...
- web兼容性测试相关知识
一.客户端兼容性 1.浏览器的兼容性测试 a.内核角度 Tridnt内核:代表作IE.腾讯.遨游.世界之窗等 Gecko内核:代表作Firefox webkit内核:代表作Safari.Chrome ...
- 深入理解php内核 编写扩展 I:介绍PHP和Zend
内容: 编写扩展I - PHP和Zend起步 原文:http://devzone.zend.com/public/view/tag/Extension Part I: Introduction to ...
- java桌面程序打包教程
首先打包成j可执行的jar文件. . 接下来找到自己生成jar文件的路径就可以看到jar文件了,我的是在桌面: 在桌面新建一个文件夹(名字随便取,一般去项目名字) 上面是我取的文件夹名字,然后把资料文 ...
- 实战-Mysql5.6.36脚本编译安装及初始化
概述 本文为centos7.3自动化编译安装mysql5.3.6的脚本及后续初始化操作,话不多少,直接上脚本. 安装脚本install.py如下: #coding=utf-8 #!/usr/bin/p ...
- 浏览器出现Cannot set property 'onclick' of null的问题
Part1: 当js文件放在head里面时,如果绑定了onclick事件,就会出现这样的错误, 是因为W3School的写法是浏览器先加载完按钮节点才执行的js,所以当浏览器自顶向下解析时,找不到on ...