HDU -- 2883 kebab（最大流判满流）

题目大意：

有n个人来买烤串，烤串机一次能烤m个烤串，给出每个人到来的时间si，离开的时间ei，每个人买ni个烤串，每个烤串所序的时间为ti，问能否满足所有顾客的需求；

思路分析：

这道题和hdu的3572很像，只不过这道题的si和ei范围太大，所以不能采用时间点，而是用时间区间；

①：设立一个源点，连向每一个人，权值为ni*ti，表示一个人所有烤串所需要的时间；

②：把每个人的时间si，ei都存下来，进行排序，然后一个人的时间区间完全包括一个相邻时间点的时间段，就从这个人向这个时间段连一条边，权值为INF，

③：每一个时间段连向汇点，权值为（（time[i]-time[i-1）*m）,表示烤串机一次能烤m个串，这一段时间能烤的烤串的总时间；

跑一遍最大流，当最大流等于所有的总时间，输出Yes，否则是No；

代码实现：

#include<cstdio>
#include<cstring>
#include<queue>
#include<iostream>
#include<algorithm>
#define Max(a,b) ((a)>(b)?(a):(b))
using namespace std;
const int N=610;
const int M=161200;
const int INF=0x3f3f3f3f;
int n,m,s,t,top,head[N],gap[N],cur[N],pre[N],dis[N],time[N],si[N],ei[N];

struct Edge{
    int to,next,flow;
}edge[M];

void Addedge(int from,int to,int val){
    edge[top].to=to,edge[top].next=head[from],edge[top].flow=val,head[from]=top++;
    edge[top].to=from,edge[top].next=head[to],edge[top].flow=0,head[to]=top++;
}

void Bfs(){
    queue<int> q;
    memset(gap,0,sizeof(gap));
    memset(dis,-1,sizeof(dis));
    gap[0]=1,dis[t]=0;
    q.push(t);
    while(!q.empty()){
        int u=q.front();
        q.pop();
        for(int i=head[u];i+1;i=edge[i].next){
            if(dis[edge[i].to]==-1){
                dis[edge[i].to]=dis[u]+1;
                gap[dis[edge[i].to]]++;
                q.push(edge[i].to);
            }
        }
    }
}

int Sap(){
    Bfs();
    memset(pre,-1,sizeof(pre));
    for(int i=0;i<=t;++i) cur[i]=head[i];
    int u=s,i,cur_flow,max_flow=0,neck,tmp;
    while(dis[s]<=t){
        if(u==t){
            cur_flow=INF;
            for(int i=s;i!=t;i=edge[cur[i]].to){
                if(cur_flow>edge[cur[i]].flow){
                    neck=i;
                    cur_flow=edge[cur[i]].flow;
                }
            }
            for(int i=s;i!=t;i=edge[cur[i]].to){
                tmp=cur[i];
                edge[tmp].flow-=cur_flow;
                edge[tmp^1].flow+=cur_flow;
            }
            max_flow+=cur_flow;
            u=neck;
        }
        int i;
        for(i=cur[u];i!=-1;i=edge[i].next)
            if(edge[i].flow&&dis[u]==dis[edge[i].to]+1) break;
        if(i!=-1){
            cur[u]=i;
            pre[edge[i].to]=u;
            u=edge[i].to;
        }else{
            if(--gap[dis[u]]==0) break;
            cur[u]=head[u];
            int mindis=t;
            for(i=head[u];i!=-1;i=edge[i].next){
                if(edge[i].flow&&mindis>dis[edge[i].to])
                    mindis=dis[edge[i].to];
            }
            dis[u]=mindis+1;
            gap[dis[u]]++;
            if(u!=s) u=pre[u];
        }
    }
    return max_flow;
}

int main(){
    int ni,ti;
    while(~scanf("%d%d",&n,&m)){
        memset(head,-1,sizeof(head));
        top=0,s=0;
        int sum=0,cnt=0,c=0;
        for(int i=1;i<=n;++i){
            scanf("%d%d%d%d",&si[i],&ni,&ei[i],&ti);
            sum+=ni*ti;
            Addedge(s,i,ni*ti);
            time[cnt++]=si[i],time[cnt++]=ei[i];
        }
        sort(time,time+cnt);
        for(int i=1;i<cnt;++i)
            if(time[c]!=time[i]) time[++c]=time[i];
        t=n+c+1;
        for(int i=1;i<=c;++i)
            Addedge(i+n,t,(time[i]-time[i-1])*m);
        for(int i=1;i<=n;++i)
            for(int j=1;j<=c;++j)
                if(si[i]<=time[j-1]&&ei[i]>=time[j]) Addedge(i,j+n,INF);
        int res=Sap();
        if(res>=sum) printf("Yes\n");
        else printf("No\n");
    }
}