luogu P2570 [ZJOI2010]贪吃的老鼠【二分+最大流】
首先考虑只满足第一个条件,二分答案,把过期时间加上mid之后的2n个时间离散,老鼠拆成每个时间的,第i个时间第j个老鼠为id[i][j],连接(s,i,p[i]),对于离散后时间(g[j-1]~g[j])在i奶酪的时间区间里的ij,连接(i,id,老鼠速度*时间段长),然后连(id,t,inf),判断合法就是dinic==sump
易证这样是满足第一个条件的,因为可以调整吃奶酪的老鼠的顺序来使其合法
然后看第二个条件,把老鼠速度从大到小排序并差分,把上个建图中的(i,id,老鼠速度*时间段长)改成(i,id,老鼠差分速度*时间段长),(id,t,inf)改成(id,t,老鼠排名*老鼠差分速度*时间段长)
这样的原理是把大的速度分成若干小段,然后把问题变成选这些小段,而小段的数量是被限制的,所以这些小段拼起来一定是若干不同的老鼠
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;
const int N=20005;
const double eps=1e-6;
int T,n,m,p[N],l[N],r[N],a[N],h[N],cnt=1,le[N],s,t,sm,id[105][105],con;
double g[N];
struct qwe
{
int ne,to;
double va;
}e[N*20];
int read()
{
int l=0,f=1;
char p=getchar();
while(p>'9'||p<'0')
{
if(p=='-')
f=-1;
p=getchar();
}
while(p>='0'&&p<='9')
{
l=l*10+p-48;
p=getchar();
}
return l*f;
}
bool cmp(const int &a,const int &b)
{
return a>b;
}
void add(int u,int v,double w)
{
cnt++;
e[cnt].ne=h[u];
e[cnt].to=v;
e[cnt].va=w;
h[u]=cnt;
}
void ins(int u,int v,double w)
{
add(u,v,w);
add(v,u,0);
}
bool bfs()
{
queue<int>q;
memset(le,0,sizeof(le));
le[s]=1;
q.push(s);
while(!q.empty())
{
int u=q.front();
q.pop();
for(int i=h[u];i;i=e[i].ne)
if(e[i].va>0&&!le[e[i].to])
{
le[e[i].to]=le[u]+1;
q.push(e[i].to);
}
}
return le[t];
}
double dfs(int u,double f)
{
if(u==t||!f)
return f;
double us=0;
for(int i=h[u];i&&us<f;i=e[i].ne)
if(e[i].va>0&&le[e[i].to]==le[u]+1)
{
double t=dfs(e[i].to,min(e[i].va,f-us));
e[i].va-=t;
e[i^1].va+=t;
us+=t;
}
if(us<eps)
le[u]=0;
return us;
}
double dinic()
{
double r=0;
while(bfs())
r+=dfs(s,1e18);
return r;
}
bool ok(double w)
{
memset(h,0,sizeof(h));
cnt=1;
for(int i=1;i<=n;i++)
ins(s,i,p[i]),g[i]=l[i],g[i+n]=r[i]+w;
sort(g+1,g+1+2*n);
for(int i=2;i<=2*n;i++)
if(g[i]-g[i-1]>eps)
for(int j=1;j<=m;j++)
{
ins(id[i][j],t,(g[i]-g[i-1])*j*a[j]);
for(int k=1;k<=n;k++)
if(g[i-1]>=l[k]&&g[i]<=r[k]+w)
ins(k,id[i][j],(g[i]-g[i-1])*a[j]);
}
// double nw=dinic();cerr<<nw<<" "<<sm<<endl;
return dinic()-sm>=0;
}
int main()
{
T=read();
while(T--)
{
con=n=read(),m=read(),sm=0;
for(int i=1;i<=n;i++)
p[i]=read(),l[i]=read(),r[i]=read(),sm+=p[i];
for(int i=1;i<=m;i++)
a[i]=read();
sort(a+1,a+1+m,cmp);
for(int i=1;i<m;i++)
a[i]-=a[i+1];
for(int i=1;i<=2*n;i++)
for(int j=1;j<=m;j++)
id[i][j]=++con;
s=0,t=con+1;
double l=0,r=1e7,ans=r;
while(r-l>eps)
{
double mid=(l+r)/2;
if(ok(mid))
r=mid,ans=mid;
else
l=mid;
}
printf("%.4f\n",ans);
}
return 0;
}
luogu P2570 [ZJOI2010]贪吃的老鼠【二分+最大流】的更多相关文章
- Luogu P2570 [ZJOI2010]贪吃的老鼠
Luogu P2570 [ZJOI2010]贪吃的老鼠 题目描述 奶酪店里最近出现了\(m\)只老鼠!它们的目标就是把生产出来的所有奶酪都吃掉.奶酪店中一天会生产\(n\)块奶酪,其中第\(i\)块的 ...
- 洛谷$P2570\ [ZJOI2010]$贪吃的老鼠 网络流+二分
正解:网络流+二分 解题报告: 传送门$QwQ$ 和上一题有点儿像,,,?$QwQ$但是比上一题要有趣很多$QwQ$ 首先把大致思路捋下?依然是.二分出每个奶酪的开始和结束时间,然后check下最大流 ...
- Luogu2570 ZJOI2010 贪吃的老鼠 二分答案+最大流
题目链接:https://www.luogu.org/problemnew/show/P2570 题意概述: 好像没什么好概述的.....很简洁? 分析: 首先想到二分时间,转化成判定性问题,在一定时 ...
- P2570 [ZJOI2010]贪吃的老鼠
传送门 →_→唯一一篇能看得懂的题解---->这里 很容易想到二分+网络流,然而并没有什么卵用--出题人的思路太神了-- 首先考虑如果一块奶酪在同一时间可以被多只老鼠吃的话,该如何建图.首先不难 ...
- Luogu2570 [ZJOI2010]贪吃的老鼠 ---- 网络流
Luogu2570 [ZJOI2010]贪吃的老鼠 题面描述 https://www.luogu.org/problemnew/show/P2570 然后题意大概就是m只老鼠,然后吃n个奶酪,已知 ...
- [ZJOI2010]贪吃的老鼠
很不错的一道网络流的题目 二分答案是显然的 首先不考虑每个饼干只能一个老鼠吃 那很显然的建图就是将时间点按照开始结束的点分成2*n-1段 然后对每一段时间建m个老鼠的点,然后s-它限流,再从它到目前可 ...
- [ZJOI2010]贪吃的老鼠 网络流
---题面--- 题解: 这是一道强题emmmm,做法非常巧妙,,,我也是看了好久大佬题解才看明白一点 首先考虑没有限制的情况,即n个老鼠可以在同一时刻吃同一块奶酪 对各个时间段拆点,连奶酪 ---& ...
- 【题解】ZJOI2010贪吃的老鼠
%%%%真的好强...看题解我都看了好久才完全明白.放一下参考的博客,谢谢神犇QAQ 1号博客 2号博客(超级赞的啦) 因为理解的过程太艰辛,所以必须记录一下这道强题:这道题目最难的两个约束就在 ...
- [ZJOI2010]贪吃的老鼠(网络流+建图)
题目描述 奶酪店里最近出现了m只老鼠!它们的目标就是把生产出来的所有奶酪都吃掉.奶酪店中一天会生产n块奶酪,其中第i块的大小为pi,会在第ri秒被生产出来,并且必须在第di秒之前将它吃掉.第j只老鼠吃 ...
随机推荐
- [2018-11-03]2018年10月28日宁波dotnet社区活动回顾及下次活动预告
离上次活动,有半年了,汗.之后尽量保证每月一次,以组织为主,多邀请嘉宾来分享. 本次活动不足之处 人手不足:由于活动组织事项受限于人手(目前就我一个,这次活动前后我又应邀给大红鹰学院应届生介绍dotn ...
- poj 1469 COURSES (二分图模板应用 【*模板】 )
COURSES Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 18454 Accepted: 7275 Descript ...
- Contiki Rtimer 模块
一.rtimer概述 The Contiki rtimer library provides scheduling and execution of real-time tasks (with pre ...
- linux 下errno各个值的意义(转) errno.h
strerror(errno):获取errno对应的错误 查看错误代码errno是调试程序的一个重要方法.当linux C api函数发生异常时,一般会将errno变量(需include errno. ...
- 理解YOLOv2训练过程中输出参数含义
原英文地址: https://timebutt.github.io/static/understanding-yolov2-training-output/ 最近有人问起在YOLOv2训练过程中输出在 ...
- Early Media and Music on Hold
Early media refers to any media that is played to the initial caller’s phone before the remote party ...
- CodeForces - 434D Nanami's Power Plant
Codeforces - 434D 题目大意: 给定一个长为n的序列,序列中的第i为上的值\(x_i\),序列第i位上的值\(x_i\in[l_i,r_i]\),价值为\(f_i(x_i)\),其中\ ...
- nginx 反向代理配置
转载一篇特别好的nginx配置博文:http://www.cnblogs.com/hunttown/p/5759959.html
- redis设置密码和redis主从复制
redis设置密码和redis主从复制 一.redis设置密码 1.Redis实用特性 安全性 主从复制(侦听器)事务处理 持久化机制 发布订阅消息 2.安全性:设置客户端连接后进行任何其他指定前需 ...
- javaCV入门指南:调用FFmpeg原生API和JavaCV是如何封装了FFmpeg的音视频操作?
通过"javaCV入门指南:序章 "大家知道了处理音视频流媒体的前置基本知识,基本知识包含了像素格式.编解码格式.封装格式.网络协议以及一些音视频专业名词,专业名词不会赘述,自行搜 ...