HDU 3572 Task Schedule ([最大流]任务分配,判断满流)

题意：有M个机器，有N个任务。每个任务必须在Si 或者以后开始做，在Ei 或者之前完成，完成任务必须处理Pi 个时间单位。其中，每个任务可以在任意（空闲）机器上工作，每个机器的同一时刻只能工作一个任务，每个任务在同一时刻只能被一个机器工作，而且任务做到一半可以打断，拿去其他机器做。能否在规定时间内把任务做完。
思路：关键是抽象出模型，建立网络图。以0为源点，与每个任务之间连一条容量是p的边，任务与之对应的时间段的每一天连一个容量是1的边，最后所有的天数，与汇点（最大的点加1）连一条容量是m的边，判断sum(pi)是不是等于最大流。
**注意一点的是，模板的数组开小会TLE

#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
#define cl(a,b) memset(a,b,sizeof(a))
const int maxn=20006;
const int inf=1<<28;
const int nv=1100;
const int ne=501000;
struct isap{
    int n,size;
    int head[nv];
    int dis[nv],gap[nv],cur[nv],pre[nv];
    int maxflow;

    struct edge{
        int v,w,next;
        edge(){}
        edge(int _v,int _w,int _next):v(_v),w(_w),next(_next){}
    }E[ne];

    void init(int x){
        n=x,size=0;
        //for(int i=0;i<=n;i++)head[i]=-1;
        cl(head,-1);
    }

    void insert(int u,int v,int w){
        E[size]=edge(v,w,head[u]);
        head[u]=size++;
        E[size]=edge(u,0,head[v]);
        head[v]=size++;
    }
    int maxFlow(int src,int des){
        maxflow=0;
        for(int i=0;i<=n;i++){
            dis[i]=gap[i]=0;
            cur[i]=head[i];
        }
        int u=pre[src]=src;
        int aug=-1;
        while(dis[src]<n){
           loop:for(int &i=cur[u];i!=-1;i=E[i].next){
                int v=E[i].v;
                if(E[i].w&&dis[u]==dis[v]+1){
                    aug=min(aug,E[i].w);
                    pre[v]=u;
                    u=v;
                    if(v==des){
                        maxflow+=aug;
                        for(u=pre[u];v!=src;v=u,u=pre[u]){
                            E[cur[u]].w-=aug;
                            E[cur[u]^1].w+=aug;
                        }
                        aug=inf;
                    }
                    goto loop;
                }
            }
            int mdis=n;
            for(int i=head[u];i!=-1;i=E[i].next){
                int v=E[i].v;
                if(E[i].w&&mdis>dis[v]){
                    cur[u]=i;
                    mdis=dis[v];
                }
            }
            if(--gap[dis[u]]==0)break;
            gap[dis[u]=mdis+1]++;
            u=pre[u];
        }
        return maxflow;
    }
}G;
int main(){
    int cas=1;
    int T;
    scanf("%d",&T);
    while(T--){
        int n,m;
        scanf("%d%d",&n,&m);
        G.init(n);
        int mx=0;
        int sum=0;
        for(int i=1;i<=n;i++){
            int pi,si,ei;
            scanf("%d%d%d",&pi,&si,&ei);
            mx=max(mx,ei);
            sum+=pi;
            G.insert(0,i,pi);
            for(int j=si;j<=ei;j++){
                G.insert(i,n+j,1);
            }
        }
        for(int i=1;i<=mx;i++){
            G.insert(n+i,mx+n+1,m);
        }
        G.n=mx+n+1;
        if(G.maxFlow(0,mx+n+1)==sum){
            printf("Case %d: Yes\n",cas++);
        }
        else {
            printf("Case %d: No\n",cas++);
        }
        printf("\n");
    }
    return 0;
}