HDU3572 Task Schedule 网络流最大流

先参考了别人的思路 打了一遍，自己想写更快的一种 可是 刚开始接触 还不行，先按别人思路来敲敲，明天继续刷这道题目

思路倒是看懂了

建图：把每个任务和每一天都看做一个点，添加源点和汇点。源点与每个任务之间连一条边，容量为完成该任务所需处理次数。若第i个任务可以在Si至Ei天处理，则由该任务向这些天分别连一条边，容量为1，表示此任务每天只能被处理一次。最后，从每一天连一条到汇点的边，容量为机器数M，表示每天可以处理M个任务。若求出的最大流等于所有任务需要处理的次数之和，说明能完成任务；否则，不能完成任务。

#include<iostream>#include<cstdio>#include<list>#include<algorithm>#include<cstring>#include<string>#include<queue>#include<stack>#include<map>#include<vector>#include<cmath>#include<memory.h>#include<set>#define ll long long#define LL __int64#define eps 1e-8const ll INF=9999999999999;using namespace std;#define M 400000100#define inf 0xfffffff//vector<pair<int,int> > G;//typedef pair<int,int> P;//vector<pair<int,int>> ::iterator iter;////map<ll,int>mp;//map<ll,int>::iterator p;////vector<int>G[8000];struct Node{int v,cap,next;}edge[550012];int pre[10012];int cur[10012];int head[10012];int level[10012];int gap[10012];int n,m;int vn,en;void clear(){en=0;memset(head,-1,sizeof(head));memset(pre,-1,sizeof(pre));memset(level,0,sizeof(level));memset(gap,0,sizeof(gap));}void add(int u,int v,int w){edge[en].cap=w;edge[en].v=v;edge[en].next=head[u];head[u]=en++;edge[en].cap=0;edge[en].v=u;edge[en].next=head[v];head[v]=en++;}int sap(int vs,int vt){for(int i=0;i<=vn;i++)cur[i]=head[i];int u=pre[vs]=vs,maxflow=0,aug=-1;gap[0]=vn;while(level[vs]<vn){loop:for(int &i=cur[u];i!=-1;i=edge[i].next){int v=edge[i].v;if(edge[i].cap && level[u]==level[v]+1){aug==-1?aug=edge[i].cap:(aug=min(aug,edge[i].cap));pre[v]=u;u=v;if(v==vt){maxflow+=aug;for(u=pre[u];v!=vs;v=u,u=pre[u]){edge[cur[u]].cap-=aug;edge[cur[u]^1].cap+=aug;}aug=-1;}goto loop;}}int minlevel=vn;for(int i=head[u];i!=-1;i=edge[i].next){int v=edge[i].v;if(edge[i].cap && minlevel>level[v]){cur[u]=i;minlevel=level[v];}}gap[level[u]]--;if(gap[level[u]]==0)break;level[u]=minlevel+1;gap[level[u]]++;u=pre[u];}return maxflow;}int main(void){int Case=0;int t;int vt;cin>>t;while(t--){clear();scanf("%d %d",&n,&m);int maxe=0,vs=0,sum=0;int p,s,e;for(int i=1;i<=n;i++){scanf("%d %d %d",&p,&s,&e);maxe=max(maxe,e);sum+=p;add(vs,i,p);for(int j=s;j<=e;j++)add(i,n+j,1);}vt=n+maxe+1;//汇点for(int i=1;i<=maxe;i++)add(n+i,vt,m);vn=vt+1;if(sap(vs,vt) == sum)printf("Case %d: Yes\n\n",++Case);elseprintf("Case %d: No\n\n",++Case);}}