[ZJOI2010]贪吃的老鼠

很不错的一道网络流的题目

二分答案是显然的

首先不考虑每个饼干只能一个老鼠吃

那很显然的建图就是将时间点按照开始结束的点分成2*n-1段

然后对每一段时间建m个老鼠的点，然后s-它限流，再从它到目前可以运行的饼干

那么考虑加上限制每个饼干只能被一个老鼠吃

我们可以考虑一下进行差分

将s-老鼠的值变成差分后的值（当然还得乘以相应个数） 同理老鼠到饼干也类似的连（就不用乘以相应个数了）

为什么这样做是可以的呢？

因为 我们可以把每一个决策都用这里差分的值来表示出来

那还有一个问题就是有的状态实际上是看似不存在的

比如有v1,v2,v3

那么如果没有v2-v1,v1只有v3-v2这个状态实际是不存在的

但是这个一定是比v3状态吃的要小的，即用比限制时间少的v3是一定能达到这个状态的

所以这个是对的

另外：由于几乎没用过c++实数的原因 被坑的很惨

刚开始以为是eps出问题导致精度不足。。后来发现对于精度要求高这不能用cout

long double在精度上高于double

输出时要转换成double输出

#include <bits/stdc++.h>
using namespace std;
#define ll long long
const ll maxn=;
#define INF 1e9
#define exp 1e-6
ll n,m,s,t,l,sum;
ll head[maxn],d[maxn],p[maxn],r[maxn];
ll dd[maxn],ss[maxn];
bool vis[maxn];
struct re{
    ll a,b;
    long double c,flow;
}a[maxn];
void arr(ll x,ll y,long double z)
{
    a[++l].a=head[x];
    a[l].b=y;
    a[l].c=z;
    a[l].flow=;
    head[x]=l;
}
bool bfs()
{
    memset(vis,,sizeof(vis));
    queue<ll> q;
    q.push(s);
    d[s]=; vis[s]=;
    while (!q.empty())
    {
        ll x=q.front();q.pop();
        ll u=head[x];
        while (u)
        {
            ll v=a[u].b;
            if (!vis[v]&&a[u].c>a[u].flow)
            {
                vis[v]=;
                d[v]=d[x]+;
                q.push(v);
            }
            u=a[u].a;
        }
    }
    return vis[t];
}
long double dfs(ll x,long double y)
{
    if (x==t||y==) return y;
    long double flow=,f;
    ll tmp;
    ll u=head[x];
    while (u)
    {
        ll v=a[u].b;
        if (d[x]+==d[v]&&(f=dfs(v,min(y,a[u].c-a[u].flow)))>)
        {
            a[u].flow+=f;
            if (u%) tmp=u+; else tmp=u-;
            y-=f;
            a[tmp].flow-=f;
            flow+=f;
            if (y==) break;
        }
        u=a[u].a;
    }
    return flow;
}
long double maxflow()
{
    long double flow=;
    while (bfs())
    {
        flow+=dfs(s,INF);
    }
    return flow;
}
struct ree{
    long double a;
    ll b;
}b[maxn];
bool cmp(ree a,ree b)
{
    return(a.a<b.a);
}
bool cmp2(ll x,ll y)
{
    return(x<y);
}
bool check(long double x)
{
    for (ll i=;i<=n;i++)
    {
        b[i*-].a=r[i]; b[i*].a=dd[i]+x;
        b[i*].b=i*; b[i*-].b=i*-;
    }
    sort(b+,b+*n+,cmp);
    l=;memset(head,,sizeof(head));
    ll last2=;
    s=; t=(*n-)*m+n+;
    for (ll i=;i<=n;i++)
    {
        arr((*n-)*m+i,t,p[i]);
        arr(t,(*n-)*m+i,);
    }
    for (ll i=;i<=*n-;i++)
    {
        for (ll j=;j<=m;j++)
        {
            arr(,j+last2,(m-j+)*(ss[j]-ss[j-])*(b[i+].a-b[i].a));
            arr(j+last2,,);
        }
        for (ll k=;k<=n;k++)
          if (r[k]<=b[i].a&&dd[k]+x>=b[i+].a)
            for (ll j=;j<=m;j++)
            {
                arr(j+last2,(*n-)*m+k,(b[i+].a-b[i].a)*(ss[j]-ss[j-]));
                arr((*n-)*m+k,j+last2,);
            }
        last2+=m;
    }
    ll to=maxflow();
    if (to==sum) return();
    else return ();
}
int main()
{
    freopen("cheese.in","r",stdin);
    freopen("cheese.out","w",stdout);
    ll T;
    cin>>T;
    for (ll i=;i<=T;i++)
    {
        cin>>n>>m; sum=;
        for (ll i=;i<=n;i++)
        {
          cin>>p[i]>>r[i]>>dd[i];
          sum+=p[i];
        }
        for (ll i=;i<=m;i++)
          cin>>ss[i];
        sort(ss+,ss+m+,cmp2);
        long double hh=,tt=1e9;
        while ((tt-hh)>exp)
        {

            long double mid=(hh+tt)/;
            if (check(mid)) tt=mid; else hh=mid;
        }
        double tmp=tt;
        printf("%.9f\n",tmp);
    }
    return ;
}