HDU 3572 Task Schedule

http://acm.hdu.edu.cn/showproblem.php?pid=3572

题意：给N个任务，M台机器。每个任务有最早才能开始做的时间S，deadline E，和持续工作的时间P。问存不存在可行的工作时间。

题解：最大流。主要问题在于建模……

　　　建模：源点与每一个任务建边，容量为P；每一天与汇点建边，容量为M；每一个任务和天数建边，容量为1。

　　　判断：判断最后求得最大流是否和所有任务处理的次数之和相等，如果相等就输出Yes。

第一次写ISP，以后做模板用了……

#include <iostream>
using namespace std;

const int INF=0xfffffff;
const int maxn=1010;
const int maxm=500010;

struct edge{
    int v,next,val;
} G[maxm];
int pre[maxn],idx, sum;
int N,M;
int level[maxn];
int gap[maxn];
int s,t;

void add_edge(int from,int to,int val)
{
    G[idx].v=to;
    G[idx].val=val;
    G[idx].next=pre[from];
    pre[from]=idx++;
    
    G[idx].v=from;
    G[idx].val=0;
    G[idx].next=pre[to];
    pre[to]=idx++;
}

int dfs(int pos,int cost, int cnt)
{
    if (pos==t)
    {
        return cost;
    }
    
    int j,minh=cnt-1,lv=cost,d;
    
    for (j=pre[pos];j!=-1;j=G[j].next)
    {
        int v=G[j].v,val=G[j].val;
        if(val>0)
        {
            if (level[v]+1==level[pos])
            {
                if (lv<G[j].val) d=lv;
                else d=G[j].val;
                
                d=dfs(v,d,cnt);
                G[j].val-=d;
                G[j^1].val+=d;
                lv-=d;
                if (level[s]>=cnt) return cost-lv;
                if (lv==0) break;
            }
            
            if (level[v]<minh)    minh=level[v];
        }
    }
    
    if (lv==cost)
    {
        --gap[level[pos]];
        if (gap[level[pos]]==0) level[s]=cnt;
        level[pos]=minh+1;
        ++gap[level[pos]];
    }
    
    return cost-lv;
    
}

int sap(int s,int t, int cnt)
{
    int flow=0;
    gap[s]=cnt;
    while (level[s]<cnt)
    {
        flow+=dfs(s,INF, cnt);
    }
    
    return flow;
}

void init()
{
    memset(pre, -1, sizeof(pre));
    memset(gap,0,sizeof(gap));
    memset(level,0,sizeof(level));
    sum = idx = 0;
}

int main()
{
    //freopen("/Users/apple/Desktop/暑假/24.1/24.1/in","r",stdin);
    int T, ca = 1;
    scanf("%d", &T);
    while (T--)
    {
        init();
        scanf("%d %d", &N, &M);
        int  P, S, E;
        int mint=INF;
        int maxt=-1;
        for (int i = 1; i <= N; ++i)
        {
            scanf("%d %d %d", &P, &S, &E);
            sum += P;
            if(mint>S) mint=S;
            if(maxt<E) maxt=E;
            add_edge(0, i, P);
            for (int j = S; j <= E; ++j)
            {
                add_edge(i, N+j, 1);
            }
        }
        int days=maxt-mint+1;
        s=0;
        t=N+days+1;
        for (int i = 1; i <= days; ++i)
        {
            add_edge(i+N, N+days+1, M);
        }

        printf("Case %d: ",ca++);
        int flow1=sap(s,t,t+1);
        if(flow1==sum){
            puts("Yes");
            puts("");
        }
        else{
            puts("No");
            puts("");
        }
    }
    return 0;
}