hdu3572-Task Schedule（最大流（ISAP），是否满流）

题目来源：http://howproblem.php?pid=3572

题意

给出多个任务，起始日期，终止日期，以及每个任务的任务量（一天里一台机器只能做其中一个任务的一点任务量），，问在各自的截止日期前，能不能完成。。所有任务。

思路

这道题目起初做没有头绪。。（菜。。）
但是题目中说每个任务可以由不同的机器完成，这个时候就可以想到一个图，假如，在所有任务前面加一个源点，流量分别是各自的任务量，把每一天都看成一个点，那么各个任务就可以和代表天的点相连，流量均为1，最后，把这些天与汇点相连，流量为机器的数量。。。
如图（样例）：

这样的话，就可以利用这个图是否满流进而判断能不能完成了。。。

代码

#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;
const int INF=0x3f3f3f3f;
const int maxn=10000+10;
int first[maxn],h[maxn],cur[maxn];
struct pp
{
    int v,nexx;
    int cap,flow;
}edge[1000000+10];
int m,n,full_flow,mx;
int tot;
void ADD(int u,int v,int w)
{
    edge[tot].v=v;
    edge[tot].nexx=first[u];
    edge[tot].cap=w;
    edge[tot].flow=0;
    first[u]=tot++;
}

void init(int cases)
{
    full_flow=tot=mx=0;
    scanf("%d%d",&m,&n);
    memset(first,-1,sizeof(first));
    for(int i=1;i<=m;i++)
    {
        int time,st,ed;
        scanf("%d%d%d",&time,&st,&ed);
        full_flow+=time;
        mx=max(mx,ed);
        ADD(0,i,time);
        ADD(i,0,0);
        for(int j=st;j<=ed;j++)
        {
            ADD(i,m+j,1);
            ADD(m+j,i,0);
        }
    }
    for(int i=1;i<=mx;i++)
    {
        ADD(m+i,m+mx+1,n);
        ADD(m+mx+1,m+i,0);
    }
    printf("Case %d: ",cases);
}
int num[maxn];
void height(int s,int t)
{
    queue<int> Q;
    memset(num,0,sizeof(num));
    memset(h,-1,sizeof(h));
    h[t]=0;
    ++num[0];
    while(!Q.empty()) Q.pop();
    Q.push(t);
    while(!Q.empty())
    {
        int v=Q.front();
        Q.pop();
        for(int i=first[v];~i;i=edge[i].nexx)
        {
            int u=edge[i].v;
            if(h[u]==-1)
            {
                h[u]=h[v]+1;
                ++num[h[u]];
                Q.push(u);
            }
        }
    }
}

int S[maxn];
void print_ISAP(int s,int t,int node)
{
    int top=0,ans=0,u=s;
    memcpy(cur,first,sizeof(first));
    while(h[s]<node)
    {
        if(u==t)
        {
            int d=INF,pos;
            for(int i=0;i<top;i++)
            {
                if(d>edge[S[i]].cap-edge[S[i]].flow)
                {
                    d=edge[S[i]].cap-edge[S[i]].flow;
                    pos=i;
                }
            }
            for(int i=0;i<top;i++)
            {
                edge[S[i]].flow+=d;
                edge[S[i]^1].flow-=d;
            }
            ans+=d;
            top=pos;
            u=edge[S[top]^1].v;
            continue;
        }
        int i=cur[u];
        for(;~i;i=edge[i].nexx)
        {
            int v=edge[i].v;
            if(h[u]==h[v]+1&&edge[i].cap>edge[i].flow)
            {
                cur[u]=i;
                S[top++]=i;
                u=v;
                break;
            }
        }
        if(i==-1)
        {
            int d=node;//点的个数
            for(int i=first[u];~i;i=edge[i].nexx)
            {
                if(edge[i].cap>edge[i].flow&&d>h[edge[i].v])
                {
                    d=h[edge[i].v];
                    cur[u]=i;
                }
            }
            if(--num[h[u]]==0) break;
            h[u]=d+1;
            ++num[h[u]];
            if(u!=s)
                u=edge[S[--top]^1].v;
        }
    }
    if(ans==full_flow)
        printf("Yes\n\n");
    else
        printf("No\n\n");
}

int main()
{
    int T;
    scanf("%d",&T);
    for(int i=1;i<=T;i++)
    {
        init(i);
        height(0,m+mx+1);
        print_ISAP(0,m+mx+1,m+mx+2);
    }
    return 0;
}