題目描述

乳酪店裡最近出現了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-5;
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[102][32],t;
inline int rd(){
    int x=0;char c=getchar();bool f=0;
    while(!isdigit(c)){if(c=='-')f=1
 
;c=getchar();}
    while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);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=0;
}
inline bool bfs(int s,int t){
    memset(deep,0,sizeof(deep));
    memcpy(cur,head,sizeof(head));
    q.push(s);deep[s]=1;
    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]+1;q.push(v); 
            }
        }
    }
    return deep[t];
}
double dfs(int u,int t,double l){
    if(u==t||!l)return l;
    double f,flow=0;
    for(int &i=cur[u];i;i=e[i].n){
        int v=e[i].to;
        if(deep[v]==deep[u]+1&&(f=dfs(v,t,min(l,e[i].l)))){
            e[i].l-=f;e[i^1].l+=f;l-=f;flow+=f;
            if(l<eps)break;
        }
    }
    return flow;
}
inline bool check(double mid){
    top=0;tot=1;
    memset(head,0,sizeof(head));
    for(int i=1;i<=n;++i)D[i]=d[i]+mid,b[++top]=r[i],b[++top]=D[i];
    sort(b+1,b+top+1);
    top=unique(b+1,b+top+1)-b-1;
    for(int i=1;i<=n;++i)add(0,i,p[i]);int num=n,T=top*m+n+1; 
    for(int i=2;i<=top;++i)
      for(int j=1;j<=m;++j)id[i][j]=++num,add(num,T,(m-j+1)*(b[i]-b[i-1])*v[j]);
    for(int i=1;i<=n;++i){
        int x=lower_bound(b+1,b+top+1,r[i])-b,y=lower_bound(b+1,b+top+1,D[i])-b;
        for(int j=x+1;j<=y;++j)
          for(int k=1;k<=m;++k)add(i,id[j][k],(b[j]-b[j-1])*v[k]);
    }
    double flow=0;
    while(bfs(0,T))flow+=dfs(0,T,2e15);
    return fabs(flow-sum)<eps;
}
int main(){
    t=rd();
    while(t--){
        n=rd();m=rd();sum=0;
        for(int i=1;i<=n;++i)p[i]=rd(),r[i]=rd(),d[i]=rd(),sum+=p[i];
        for(int i=1;i<=m;++i)v[i]=rd();sort(v+1,v+m+1);
        for(int i=m;i>=1;--i)v[i]-=v[i-1];
        double l=0,r=1e9,ans=2e9;
        while(l+eps<r){
            double mid=(l+r)/2;
            if(check(mid))ans=mid,r=mid;else l=mid;
        }
        printf("%.4lf\n",ans);
    }
    return 0;
}