第一句话:算出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 新生舞会的更多相关文章

  1. 【BZOJ 4819】 4819: [Sdoi2017]新生舞会 (0-1分数规划、二分+KM)

    4819: [Sdoi2017]新生舞会 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 601  Solved: 313 Description 学校 ...

  2. bzoj 4819: [Sdoi2017]新生舞会

    Description 学校组织了一次新生舞会,Cathy作为经验丰富的老学姐,负责为同学们安排舞伴.有n个男生和n个女生参加舞会 买一个男生和一个女生一起跳舞,互为舞伴.Cathy收集了这些同学之间 ...

  3. [BZOJ 4819] [SDOI 2017] 新生舞会

    Description 学校组织了一次新生舞会,Cathy作为经验丰富的老学姐,负责为同学们安排舞伴. 有 \(n\) 个男生和 \(n\) 个女生参加舞会买一个男生和一个女生一起跳舞,互为舞伴. C ...

  4. 4819: [Sdoi2017]新生舞会(分数规划)

    4819: [Sdoi2017]新生舞会 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1031  Solved: 530[Submit][Statu ...

  5. BZOJ4819 新生舞会

    4819: [Sdoi2017]新生舞会 Time Limit: 10 Sec  Memory Limit: 128 MB Description 学校组织了一次新生舞会,Cathy作为经验丰富的老学 ...

  6. 【BZOJ4819】新生舞会(分数规划,网络流)

    [BZOJ4819]新生舞会(分数规划,网络流) 题面 BZOJ Description 学校组织了一次新生舞会,Cathy作为经验丰富的老学姐,负责为同学们安排舞伴.有n个男生和n个女生参加舞会 买 ...

  7. [BZOJ4819][SDOI2017]新生舞会(分数规划+费用流,KM)

    4819: [Sdoi2017]新生舞会 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1097  Solved: 566[Submit][Statu ...

  8. bzoj4819 [Sdoi2017]新生舞会

    Description 学校组织了一次新生舞会,Cathy作为经验丰富的老学姐,负责为同学们安排舞伴.有n个男生和n个女生参加舞会买一个男生和一个女生一起跳舞,互为舞伴.Cathy收集了这些同学之间的 ...

  9. [Sdoi2017]新生舞会 [01分数规划 二分图最大权匹配]

    [Sdoi2017]新生舞会 题意:沙茶01分数规划 貌似\(*10^7\)变成整数更科学 #include <iostream> #include <cstdio> #inc ...

随机推荐

  1. 《Linux命令行与shell脚本编程大全》第十六章 控制脚本

    一些控制脚本的方式:向脚本发送信号.修改脚本优先级,在脚本运行时切换到运行模式 16.1 处理信号 linux利用信号与运行在系统中的进程进行通信. 也可以通过对脚本进行编程,使其在收到特定信号时执行 ...

  2. thinkphp5源码解析(1)数据库

    前言 tp5的数据库操作全部通过Db类完成,比较符合国人的习惯,比如简单的Db::query().Db::execute(),还有复杂的链式操作Db::where('id=1')->select ...

  3. GIL(全局解释器锁)

    引入 现在绝大部分的Python都是CPython解释器(但不是必须使用CPython解释器),而CPython的一个特性就是有GIL,作用保证解释器级别的代码在运行时不被其他的线程进行修改,即加锁处 ...

  4. Mysql 的 IF 判断

    mysql自带很多判断逻辑,今天说一说IF的判断语句,正好今天做项目的时候也用到了 1. IF 判断 IF判断和我们代码里面写的有略微的差别,举个例子 IF('表达式','结果1','结果2')  如 ...

  5. 【luogu P1613】跑路

    https://www.luogu.org/problem/show?pid=1613 看到2k就能想到倍增.用一个数组avai[i][j][k]表示点i与点j是否存在长2k的路径,则可以递推出ava ...

  6. NGUI_Button

    十.按钮,Button 1.按钮的核心作用: 按钮能够接收单击并触发响应事件 按钮单击时能同时触发多个响应事件 按钮可以有普通.悬停.单击.禁用等多个状态的不同表现 广泛的说,按钮的核心在于接收事件 ...

  7. 从durable谈起,我是如何用搜索引擎抓住技术的关键字学习新姿势打开敏捷开发的大门

    ---又名我讨厌伸手党 我又把个人博客的子标题改为了 你可以在书和搜索引擎找到90%的问题的答案,为什么要问别人?剩下的10%或许没有答案,为什么要问别人? 这是由于最近在网上看到各种伸手,对于我这种 ...

  8. 基于C++11的线程池

    1.封装的线程对象 class task : public std::tr1::enable_shared_from_this<task> { public: task():exit_(f ...

  9. Ubuntu1204 vim中文乱码解决方法

    加入中文字符编码 sudo vi /var/lib/locales/supported.d/local #加入以下的中文字符集 zh_CN.GBK GBK zh_CN.GB2312 GB2312 zh ...

  10. Android使用XUtils框架上传照片(一张或多张)和文本,server接收照片和文字(无乱码)

    Android上传图片,这里我使用了如今比較流行的XUtils框架.该框架能够实现文件上传.文件下载.图片缓存等等,有待研究. 以下是Android端上传的代码: xUtils.jar下载 Strin ...