HDU3572 Task Schedule 网络流最大流

最新推荐文章于 2019-10-29 13:31:52 发布

韦我独尊-德天独厚

最新推荐文章于 2019-10-29 13:31:52 发布

阅读量1k

点赞数

分类专栏：网络流文章标签：最大流

本文链接：https://blog.csdn.net/yitiaodacaidog/article/details/13159157

版权

网络流专栏收录该内容

5 篇文章 0 订阅

订阅专栏

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

思路倒是看懂了

建图：把每个任务和每一天都看做一个点，添加源点和汇点。源点与每个任务之间连一条边，容量为完成该任务所需处理次数。若第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-8
const 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);
		else
			printf("Case %d: No\n\n",++Case);
	}
}