[ZJOI2010]贪吃的老鼠(网络流+建图)

题目描述

奶酪店里最近出现了m只老鼠！它们的目标就是把生产出来的所有奶酪都吃掉。奶酪店中一天会生产n块奶酪，其中第i块的大小为pi，会在第ri秒被生产出来，并且必须在第di秒之前将它吃掉。第j只老鼠吃奶酪的速度为sj，因此如果它单独吃完第i快奶酪所需的时间为pi/sj。老鼠们吃奶酪的习惯很独特，具体来说：

(1) 在任一时刻，一只老鼠最多可以吃一块奶酪；

(2) 在任一时刻，一块奶酪最多被一只老鼠吃。

由于奶酪的保质期常常很短，为了将它们全部吃掉，老鼠们需要使用一种神奇的魔法来延长奶酪的保质期。将奶酪的保质期延长T秒是指所有的奶酪的di变成di+T。同时，使用魔法的代价很高，因此老鼠们希望找到最小的T使得可以吃掉所有的奶酪。

题解

玄学的建图方式。

首先我们可以把时间离散化，然后把时间分成了几段，把每段看成一个类，这个类中可用的奶酪是已知的。

源点向奶酪连大小的边，奶酪向老鼠连t*v[i]大小的边。

然后考虑如何限制每种奶酪同一时刻只会被一只老鼠吃，每只老鼠同一时刻只会吃一种奶酪。

做法是把奶酪按照速度差分，这样可以保证任意一种选择方法都可以拆成几只老鼠的组合。

证明？%&%&*%&&%

代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
#include<cmath>
#define N 2002
using namespace std;
queue<int>q;
const double eps=1e-;
double p[N],d[N],v[N],b[N],D[N],r[N],sum;
int deep[N],tot,head[N],cur[N],top,n,m,id[][],t;
inline int rd(){
    int x=;char c=getchar();bool f=;
    while(!isdigit(c)){if(c=='-')f=;c=getchar();}
    while(isdigit(c)){x=(x<<)+(x<<)+(c^);c=getchar();}
    return f?-x:x;
}
struct node{int n,to;double l;}e[N*N];
inline void add(int u,int v,double l){
    e[++tot].n=head[u];e[tot].to=v;head[u]=tot;e[tot].l=l;
    e[++tot].n=head[v];e[tot].to=u;head[v]=tot;e[tot].l=;
}
inline bool bfs(int s,int t){
    memset(deep,,sizeof(deep));
    memcpy(cur,head,sizeof(head));
    q.push(s);deep[s]=;
    while(!q.empty()){
        int u=q.front();q.pop();
        for(int i=head[u];i;i=e[i].n){
            int v=e[i].to;
            if(e[i].l>eps&&!deep[v]){
                deep[v]=deep[u]+;q.push(v);
            }
        }
    }
    return deep[t];
}
double dfs(int u,int t,double l){
    if(u==t||!l)return l;
    double f,flow=;
    for(int &i=cur[u];i;i=e[i].n){
        int v=e[i].to;
        if(deep[v]==deep[u]+&&(f=dfs(v,t,min(l,e[i].l)))){
            e[i].l-=f;e[i^].l+=f;l-=f;flow+=f;
            if(l<eps)break;
        }
    }
    return flow;
}
inline bool check(double mid){
    top=;tot=;
    memset(head,,sizeof(head));
    for(int i=;i<=n;++i)D[i]=d[i]+mid,b[++top]=r[i],b[++top]=D[i];
    sort(b+,b+top+);
    top=unique(b+,b+top+)-b-;
    for(int i=;i<=n;++i)add(,i,p[i]);int num=n,T=top*m+n+;
    for(int i=;i<=top;++i)
      for(int j=;j<=m;++j)id[i][j]=++num,add(num,T,(m-j+)*(b[i]-b[i-])*v[j]);
    for(int i=;i<=n;++i){
        int x=lower_bound(b+,b+top+,r[i])-b,y=lower_bound(b+,b+top+,D[i])-b;
        for(int j=x+;j<=y;++j)
          for(int k=;k<=m;++k)add(i,id[j][k],(b[j]-b[j-])*v[k]);
    }
    double flow=;
    while(bfs(,T))flow+=dfs(,T,2e15);
    return fabs(flow-sum)<eps;
}
int main(){
    t=rd();
    while(t--){
        n=rd();m=rd();sum=;
        for(int i=;i<=n;++i)p[i]=rd(),r[i]=rd(),d[i]=rd(),sum+=p[i];
        for(int i=;i<=m;++i)v[i]=rd();sort(v+,v+m+);
        for(int i=m;i>=;--i)v[i]-=v[i-];
        double l=,r=1e9,ans=2e9;
        while(l+eps<r){
            double mid=(l+r)/;
            if(check(mid))ans=mid,r=mid;else l=mid;
        }
        printf("%.4lf\n",ans);
    }
    return ;
}