ZJOI 2010 贪吃的老鼠最大流

题目链接：点我点我:-) —只有洛谷有此题

题目描述：
奶酪店里最近出现了m只老鼠！它们的目标就是把生产出来的所有奶酪都吃掉。奶酪店中一天会生产n块奶酪，其中第i块的大小为pi，会在第ri秒被生产出来，并且必须在第di秒之前将它吃掉。第j只老鼠吃奶酪的速度为sj，因此如果它单独吃完第i快奶酪所需的时间为pi/sj。老鼠们吃奶酪的习惯很独特，具体来说：
(1) 在任一时刻，一只老鼠最多可以吃一块奶酪；
(2) 在任一时刻，一块奶酪最多被一只老鼠吃。
由于奶酪的保质期常常很短，为了将它们全部吃掉，老鼠们需要使用一种神奇的魔法来延长奶酪的保质期。将奶酪的保质期延长T秒是指所有的奶酪的di变成di+T。同时，使用魔法的代价很高，因此老鼠们希望找到最小的T使得可以吃掉所有的奶酪。

输入格式：
输入文件的第一行包含一个整数K，表示输入文件中数据的组数。
每组数据的第一行包含两个整数n和m，分别表示奶酪和老鼠的数量。接下来的n行每行包含三个整数pi,ri,di。最后m行每行包含一个整数，表示sj。pi,ri,di,sj的含义如上文所述。

输出格式：
包含K行，每行包含一个实数，表示你找到的最小的T。你的答案和标准答案的绝对误差不应超过10−410^{-4}。

数据规模：
30%的数据中，1≤n,m≤5；
100%的数据中，1≤K≤5,1≤n,m≤30,1≤pi≤10510^5, 0 ≤ri < di≤10710^7,1≤sj≤10510^5。

思路：
网络流好题
应该容易想到二分答案和离散化奶酪的时间点

先考虑朴素的建图，也就是不考虑每个时刻可以有多只老鼠在不同的小时刻吃同一个奶酪
1.我们将源点每个奶酪链流量为p[i]的边
2.然后把每只老鼠的每个时间段拆为一个点，再将其与对应的奶酪连边，流量为这只老鼠可以在这个时间段吃多少奶酪
3.所有的老鼠向汇点连边，无限流量

然而题目有限制，每个时刻可以有多只老鼠在不同的小时刻吃同一个奶酪，做法很巧妙，我们可以对老鼠的速度再进行拆点，把s[]从大到小排序然后与下一个作差，如”9 6 2 1” -> “3 4 1 1”，然后再将其与汇点连边的流量改为（rat为老鼠编号，dur为时间段大小，s为作差后的数组）rat*s[rat]*dur即可

为什么呢？
可以发现速度进行过如此处理之后，
9 = 3+4+1+1
6 = 4+1+1
2 = 1+1
1 = 1
9+6+2+1 = 3*1 + 4*2 + 1*3 + 1*4
(下面所有的s数组均为作差后的数组)
1.所以，当我们把老鼠点到汇点的流量限制设为rat*s[rat]*dur后，对于每一只老鼠，它‘吃’的量不超过实际可以‘吃’的量。（上述式子解释了，实际包括编号比它小的老鼠的吃的量）
2.每个老鼠点与奶酪的连边为s[rat]*dur，这样可以保证每只老鼠实际不可以超过自己这段时间可以吃的奶酪。
3.（对于第2点的解释）实际这段时间奶酪被吃的量不超过最快的老鼠的速度在这个时间段可以吃的量，如果每只老鼠在2点所述的边满流，那么实际相当于吃的最快的老鼠吃了这一整个时间段。而任何其它的情况均可以用别的老鼠凑出来，且由于第一点不会不合法

最后，每次二分答案，重新构图跑最大流，检查是否满流即可。

感想：
很好的网络流题目，然而连网络流都没想到是什么鬼？？？
原来网络流构图还可以这么神奇，把每个点使劲拆成这个鬼样子！！！
网络流的代码应该没什么问题了，但是构图能力还要再提高啦！
吐槽一下我调了3个小时的错误：
1.精度eps，写错居然可以导致最大流超过了我的p数组值的和！！！
2.我的天哪：(f=dfs(e.to, fmin(Mint, e.cap-e.flow)))>eps 这句话我不小心把 > eps写到最后一个括号里面了。。

代码：

//miaomiao 2017.2.3
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>using namespace std;#define pb push_back
#define Set(a, v) memset(a, v, sizeof(a))
#define For(i, a, b) for(int i = (a); i <= (int)(b); i++)#define MN (30+5)
#define N (3000+5)
#define eps (1e-6)struct Edge{int from, to;double cap, flow;
};struct Dinic{vector<int> G[N];vector<Edge> edges;int m, s, t;int cur[N], d[N]; void init(){For(i, 0, N-1) G[i].clear();edges.clear(); }void AddEdge(int u, int v, double cap){edges.pb((Edge){u, v, cap, 0});    edges.pb((Edge){v, u, 0, 0});m = edges.size();G[u].pb(m-2); G[v].pb(m-1);}bool Bfs(){queue<int> q;q.push(s); Set(d, 0);d[s] = 1; int now;while(!q.empty()){now = q.front(); q.pop();For(i, 0, G[now].size()-1){Edge &e = edges[G[now][i]];if(!d[e.to] && e.cap-e.flow > eps){d[e.to] = d[now]+1;q.push(e.to);}}}return d[t];}double dfs(int now, double Mint){if(Mint < eps || now == t) return Mint;double f;double ret = 0.0;for(int &i = cur[now]; i < G[now].size(); i++){Edge &e = edges[G[now][i]];if(d[e.to]==d[now]+1 && (f=dfs(e.to, fmin(Mint, e.cap-e.flow)))>eps){e.flow += f; Mint -= f;edges[G[now][i]^1].flow -= f;ret += f;if(Mint < eps) return ret;}}return ret;}double Maxflow(){double ret = 0.0;while(Bfs()){Set(cur, 0);ret += dfs(s, (1LL<<60)*1.0);}return ret;}
}Din;int n, m;
double L, R, rp[MN], rr[MN], rd[MN], rs[MN], Ti[MN*2];void solve(){double mid, la, sum = 0;int pn;For(i, 1, n) sum += rp[i];while(R-L > eps){mid = (L+R)/2.0;Din.init();For(i, 1, n) Din.AddEdge(0, i, rp[i]);For(i, 1, n){Ti[2*i-1] = rr[i]; Ti[2*i] = rd[i]+mid;            }sort(Ti+1, Ti+2*n+1);Din.s = 0; Din.t = m*2*n+n+1; pn = n;For(rat, 1, m){la = Ti[1];For(i, 2, 2*n){if(Ti[i]-la < eps) continue;++pn; Din.AddEdge(pn, Din.t, rat*rs[rat]*(Ti[i]-la));For(j, 1, n)if(rr[j]-la < eps && (rd[j]+mid)-Ti[i] > -eps)Din.AddEdge(j, pn, rs[rat]*(Ti[i]-la));la = Ti[i];}}if(sum-Din.Maxflow() < eps) R = mid;else L = mid;}printf("%lf\n", L);
}int main(){int T;double tot;scanf("%d", &T);while(T--){tot = 0.0;scanf("%d%d", &n, &m);For(i, 1, n) scanf("%lf%lf%lf", &rp[i], &rr[i], &rd[i]), tot += rp[i];For(i, 1, m) scanf("%lf", &rs[i]); sort(rs+1, rs+m+1, greater<int>());R = tot/rs[1] + 1.0; L = 0;For(i, 1, m-1) rs[i] -= rs[i+1];solve();}return 0;
}

ZJOI 2010 贪吃的老鼠最大流相关推荐

[ZJOI2010]贪吃的老鼠
P2570 [ZJOI2010]贪吃的老鼠在Ta的博客查看显然二分,最大流判定要满足两个条件: (1) 在任一时刻,一只老鼠最多可以吃一块奶酪: (2) 在任一时刻,一块奶酪最多被一只老鼠吃. ...
Luogu P2570 [ZJOI2010]贪吃的老鼠
Luogu P2570 [ZJOI2010]贪吃的老鼠题目描述奶酪店里最近出现了$m$只老鼠!它们的目标就是把生产出来的所有奶酪都吃掉.奶酪店中一天会生产$n$块奶酪,其中第$i$块的 ...
洛谷 P2570 [ZJOI2010]贪吃的老鼠
洛谷 P2570 [ZJOI2010]贪吃的老鼠题目题目描述奶酪店里最近出现了m只老鼠!它们的目标就是把生产出来的所有奶酪都吃掉.奶酪店中一天会生产n块奶酪,其中第i块的大小为pi,会在第ri秒 ...
[ZJOI2010] 贪吃的老鼠（二分+差分+神仙建图网络流）
problem luogu-P2570 solution 卧槽网络流尼玛神题首先这个最小延长时间 TTT ,套路地考虑二分,将问题转化为判定性问题. 其次 n,mn,mn,m 和奶酪存在时间 [l, ...
22Exchange Server 2010跨站点部署-邮件流测试及重定向
11.5.2 OWA登陆发送测试(重定向) 使用https://mail.robin.com/owa登陆下上海分支机构的用户可以看到,登陆成功,这里就使用到了重定向 11.5.3 SH-GZ邮件流 ...
洛谷$P2570\ [ZJOI2010]$贪吃的老鼠网络流+二分
正解:网络流+二分解题报告: 传送门$QwQ$ 和上一题有点儿像,,,?$QwQ$但是比上一题要有趣很多$QwQ$ 首先把大致思路捋下?依然是.二分出每个奶酪的开始和结束时间,然后check下最大流 ...
【网络流】LGP2570 [ZJOI2010]贪吃的老鼠
[题目] 原题地址原题说的很清楚了,在这里也贴出来吧. 奶酪店里最近出现了m只老鼠!它们的目标就是把生产出来的所有奶酪都吃掉.奶酪店中一天会生产n块奶酪,其中第i块的大小为pi,会在第ri秒被生产出 ...
P2570 [ZJOI2010]贪吃的老鼠
传送门 →_→唯一一篇能看得懂的题解---->这里很容易想到二分+网络流,然而并没有什么卵用--出题人的思路太神了-- 首先考虑如果一块奶酪在同一时间可以被多只老鼠吃的话,该如何建图.首先不难 ...
[ZJOI 2010]count 数字计数
Description 题库链接问你 $[l,r]$ 区间内所有整数中各个数码出现了多少次. $1\leq a\leq b\leq 10^{12}$ Solution 数位 $DP$ . ...
DP专练2 （大理石 + [ZJOI 2010]数字计数）
你肯定以为DP专练会有很多题, 但是请考虑一下本仙女的DP码力,一次性能更几个题... 来吧,别害怕呀~~ 文章目录大理石题目题解代码实现数字计数题目题解代码实现大理石题目林老师 ...

ZJOI 2010 贪吃的老鼠最大流

ZJOI 2010 贪吃的老鼠最大流相关推荐

最新文章

热门文章

ZJOI 2010 贪吃的老鼠 最大流

ZJOI 2010 贪吃的老鼠 最大流相关推荐

最新文章

热门文章

ZJOI 2010 贪吃的老鼠最大流

ZJOI 2010 贪吃的老鼠最大流相关推荐