hdu2883

hdu2883
N个客户 在si 开始 ei 结束（区间）一个机器要烤(n*t)单位时间
现在有M个烤肉架
问能不能烤完

分析
源点到每个客户n*t 限制客户烤的时间 客户到单位时间段 INF（M）为 限制 该客户可以在这个时间最多烤多少，单位时间段到汇点为M限制单位时间烤的数量。

优化
由于时间段数据很大，直接化成单位时间段会T，化为不重复时间段，该时间段到汇点为 len*M

#include <bits/stdc++.h>
//#include<cstdio>
//#include<cstring>
//#include<algorithm>
//#include<queue>
//#include<vector>
//#include<string>
//#include<iostream>
using namespace std;
typedef long long ll;
const ll INF=1e18;
const int maxn=100005;
const int maxm=maxn*10;
struct point
{
    int s,e;
    ll p;
};
point pp[maxn];
struct edge
{
    int v,next;
    ll w;
};
edge no[maxm];
set<int >se;
int temp[maxn];
int head[maxn];
int level[maxn];
int cnt;
queue <int> q;
bool bfs(int s,int t,int n)
{
    for(int i=0;i<=n;i++)
        level[i]=-1;
    while(!q.empty())
    {
        q.pop();
    }
    q.push(s);
    level[s]=0;
    while(!q.empty())
    {
        int u=q.front();
        if(u==t)
            return true;
        q.pop();
        for(int i=head[u];i;i=no[i].next)
        {
            int v=no[i].v;
            ll w=no[i].w;
            if(w>0&&level[v]==-1)
            {
                level[v]=level[u]+1;
                q.push(v);
            }
        }
    }
    return false;
}
ll dfs(int u,int t,ll max_)
{
    if(u==t)
        return max_;
    ll sum=0;
    for(int i=head[u];i;i=no[i].next)
    {
        edge & e=no[i];
        if(e.w>0&&level[e.v]==level[u]+1)
        {
            ll f=dfs(e.v,t,min(max_-sum,e.w));
            sum+=f;
            no[i].w-=f;
            no[1+((i-1)^1)].w+=f;
            if(sum==max_)
                break;
        }
    }
    if(!sum)
        level[u]=-1;
    return sum;
}
ll dinic(int s,int t,int n)
{
    ll ans=0;
    while(bfs(s,t,n))
    {
        ll f;
        while((f=dfs(s,t,INF))>0)
            ans+=f;
    }
    return ans;
}
void init()
{
    se.clear();
    cnt=1;
    for(int i=0;i<maxn;i++)
    {
        head[i]=0;
    }
    //memset(vis,0,sizeof(vis));
}
void add(int x,int y,ll w)
{
    no[cnt].v=y;
    no[cnt].w=w;
    no[cnt].next=head[x];
    head[x]=cnt;
    cnt++;
}

int main()
{
    int n,d;
    while(cin>>n>>d)
    {
        init();
        ll ans=0;
        int ss=0,t;
        for(int i=1;i<=n;i++)
        {
            int s,n,e,t;
            cin>>s>>n>>e>>t;
            pp[i].s=s;
            pp[i].e=e;
            se.insert(s);
            se.insert(e);
            pp[i].p=n*t;
            ans+=pp[i].p;
        }
        int num=0;
        for(auto it:se)
        {
            ++num;
            temp[num]=it;
        }
        t=n+num+1;
        for(int i=1;i<=n;i++)
        {
            add(ss,i,pp[i].p);
            add(i,ss,0);
            for(int j=1;j<num;j++)
            {
                if(temp[j]>=pp[i].s&&temp[j]<pp[i].e)
                {
                    add(i,j+n,INF);
                    add(j+n,i,0);
                }
            }
        }
        for(int i=1;i<num;i++)
        {
            add(i+n,t,d*(temp[i+1]-temp[i]));
            add(t,i+n,0);
        }
        ll xx=dinic(ss,t,n+num+1);
        //cout<<xx<<endl;
        if(ans==xx)
            cout<<"Yes"<<endl;
        else
            cout<<"No"<<endl;

    }
    return 0;
}