P2570 [ZJOI2010]贪吃的老鼠
在Ta的博客查看
显然二分，最大流判定
要满足两个条件：
(1) 在任一时刻，一只老鼠最多可以吃一块奶酪；
(2) 在任一时刻，一块奶酪最多被一只老鼠吃。
先按照奶酪的边界进行离散化， 变成num个块，就可以知道每个时间有哪些奶酪了
把每个老鼠拆成num个点，
初步：
每个老鼠的每个时间的奶酪连接O(len*speed)
首先每个老鼠每个时间段吃的是有限的，显然保证（1）
但是不能保证(2)
改进：
考虑让每一个奶酪在时间段只能
证明不会咕了
（总感觉不能一定能满足存在一种方案使得总共时间不会超出）
卡精度啊，，，，，
inf设太大了
并且为了防止被inf卡，可以直接记录ret表示流出流量，直接返回ret即可
#include<bits/stdc++.h> #define reg register int #define il inline #define fi first #define se second #define mk(a,b) make_pair(a,b) #define numb (ch^'0') #define pb push_back #define solid const auto & #define enter cout<<endl #define pii pair<int,int> #define double long double using namespace std; typedef long long ll; template<class T>il void rd(T &x){char ch;x=0;bool fl=false;while(!isdigit(ch=getchar()))(ch=='-')&&(fl=true);for(x=numb;isdigit(ch=getchar());x=x*10+numb);(fl==true)&&(x=-x); } template<class T>il void output(T x){if(x/10)output(x/10);putchar(x%10+'0');} template<class T>il void ot(T x){if(x<0) putchar('-'),x=-x;output(x);putchar(' ');} template<class T>il void prt(T a[],int st,int nd){for(reg i=st;i<=nd;++i) ot(a[i]);putchar('\n');}namespace Miracle{ const int N=33; const int P=33*33*2+44; const double inf=1e9; const double eps=1e-8; int n,m; struct node{int nxt,to;double w; }e[2*(P*N+N+P)]; int hd[P],cnt=1; int p[N],st[N],nd[N]; int sp[N]; double mem[2*N]; int num; int sum; void add(int x,int y,double z){e[++cnt].nxt=hd[x];e[cnt].to=y;e[cnt].w=z;hd[x]=cnt;e[++cnt].nxt=hd[y];e[cnt].to=x;e[cnt].w=0;hd[y]=cnt; } int d[P]; int s,t; double dfs(int x,double flow){if(x==t) return flow;double res=flow;for(reg i=hd[x];i&&res>eps;i=e[i].nxt){int y=e[i].to;if(e[i].w>eps&&d[y]==d[x]+1){double k=dfs(y,min(res,e[i].w));if(!k) d[y]=0;res-=k;e[i].w-=k;e[i^1].w+=k;}}return flow-res; } int q[P],l,r; bool bfs(){memset(d,0,sizeof d);l=1;r=0;q[++r]=s;d[s]=1;while(l<=r){int x=q[l++];for(reg i=hd[x];i;i=e[i].nxt){int y=e[i].to;if(e[i].w>eps&&!d[y]){d[y]=d[x]+1;q[++r]=y;if(y==t) return true;}}}return false; } int id(int x,int y){return m+(x-1)*num+y; } bool che(double mid){memset(hd,0,sizeof hd);cnt=1;num=0;for(reg i=1;i<=m;++i){mem[++num]=st[i];mem[++num]=nd[i]+mid;}sort(mem+1,mem+num+1);num=unique(mem+1,mem+num+1)-mem-1;--num;//warning!!! s=0,t=id(n,num)+1;for(reg i=1;i<=m;++i){add(i,t,p[i]);}for(reg i=1;i<=n;++i){for(reg j=1;j<=num;++j){for(reg k=1;k<=m;++k){if(st[k]<=mem[j]&&mem[j+1]<=nd[k]+mid){add(id(i,j),k,(double)sp[i]*(mem[j+1]-mem[j]));}}add(s,id(i,j),(double)sp[i]*i*(mem[j+1]-mem[j]));}}double flow=0,ret=0;while(bfs()){while(1){flow=dfs(s,inf);if(flow<eps) break;ret+=flow;}}if(ret+eps>sum) return true;return false; } void clear(){sum=0;s=0;t=0; } bool cmp(int x,int y){return x>y; } int main(){int T;rd(T);while(T--){clear();rd(m);rd(n); for(reg i=1;i<=m;++i){rd(p[i]);rd(st[i]);rd(nd[i]);sum+=p[i];}for(reg i=1;i<=n;++i){rd(sp[i]);}sort(sp+1,sp+n+1,cmp);sp[n+1]=0;for(reg i=1;i<=n;++i){sp[i]-=sp[i+1];}double L=0.0,R=1e7+3;for(reg i=1;i<=60;++i){double mid=(R+L)/2;if(che(mid)) R=mid;else L=mid;}printf("%.10Lf\n",L);}return 0; }} signed main(){Miracle::main();return 0; }/*Author: *Miracle* */
 Example)
)