ZOJ3229 Shoot the Bullet

题意：一共有n天，m个妹子，一个人至少要给第i个妹子排Gi张照片，他第i天只能拍di张照片，只能给ci个人拍照，问是否能完成任务，能完成最多能拍多少张照片，第i天给按顺序输入的妹子排几张。。

上下界网络流，先要知道怎么求无源无汇网络流（不说了）。。。

这题先就简单的按照题意建个图出来，连起t->s（无源无汇），然后加上一个超级源点，超级汇点，然后按照一般的无源无汇网络流建图，跑一遍，判断一下是否可行。。如果不可行，就说明不能满足每个点流入的流量等于流出的，所以输出-1，否则，删掉超级源点，超级汇点，从原来的源点和汇点跑一遍得出答案。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int MAXN=2010;
const int MAXE=1000010;
const int INF=(1<<30)-1;
struct EDGE
{
	int v,next;
	int cap,flow;
}edge[MAXE];
int head[MAXN],size;
void init()
{
	memset(head,-1,sizeof(head));
	size=0;
}
void add_edge(int u,int v,int cap)
{
	edge[size].v=v;
	edge[size].cap=cap;
	edge[size].flow=0;
	edge[size].next=head[u];
	head[u]=size++;
	edge[size].v=u;
	edge[size].cap=0;
	edge[size].flow=0;
	edge[size].next=head[v];
	head[v]=size++;
}
int nv;
int pe[MAXN],pre[MAXN],dist[MAXN],gap[MAXN];
int ISAP(int s,int t)
{
	int cur_flow,temp,u,neck,i;
	int max_flow;
	memset(pre,-1,sizeof(pre));
	memset(gap,0,sizeof(gap));
	memset(dist,0,sizeof(dist));
	for(i=0;i<nv;i++)
		pe[i]=head[i];
	max_flow=0;
	u=s;
	while(dist[s]<nv)
	{
		if(u==t)
		{
			cur_flow=INF+1;
			for(i=s;i!=t;i=edge[pe[i]].v)
			{
				if(cur_flow>edge[pe[i]].cap)
				{
					cur_flow=edge[pe[i]].cap;
					neck=i;
				}
			}
			for(i=s;i!=t;i=edge[pe[i]].v)
			{
				temp=pe[i];
				edge[temp].cap-=cur_flow;
				edge[temp^1].cap+=cur_flow;
			}
			max_flow+=cur_flow;
			u=neck;
		}
		for(i=pe[u];i!=-1;i=edge[i].next)
		{
			if(edge[i].cap&&dist[u]==dist[edge[i].v]+1)
				break;
		}
		if(i!=-1)
		{
			pe[u]=i;
			pre[edge[i].v]=u;
			u=edge[i].v;
		}
		else
		{
			if(0==--gap[dist[u]])
				break;
			temp=nv;
			pe[u]=head[u];
			for(i=head[u];i!=-1;i=edge[i].next)
			{
				if(edge[i].cap)
					temp=min(temp,dist[edge[i].v]);
			}
			dist[u]=temp+1;
			++gap[dist[u]];
			if(u!=s)
				u=pre[u];
		}
	}
	return max_flow;
}
int n,m;
int f[MAXE],id[MAXE],down[MAXE];
void build()
{
	int i,j,val;
	memset(f,0,sizeof(f));
	init();
	int s=0;
	int t=n+m+1;
	int ss=t+1;
	int tt=ss+1;
	nv=tt+1;
	for(i=1;i<=m;i++)
	{
		scanf("%d",&val);
		f[i+n]-=val;
		f[t]+=val;
		add_edge(i+n,t,INF);
	}
	int c,d,l,r,x;
	int cnt=0;
	for(i=1;i<=n;i++)
	{
		scanf("%d%d",&c,&d);
		add_edge(s,i,d);
		while(c--)
		{
			scanf("%d%d%d",&x,&l,&r);
			x++;
			f[i]-=l;
			f[x+n]+=l;
			down[cnt]=l;
			id[cnt++]=size;
			add_edge(i,x+n,r-l);
		}
	}
	int sum=0;
	for(i=0;i<=t;i++)
	{
		if(f[i]>0)
		{
			sum+=f[i];
			add_edge(ss,i,f[i]);
		}
		else if(f[i]<0)
			add_edge(i,tt,-f[i]);
	}
	add_edge(t,s,INF);
	int ans=ISAP(ss,tt);
	if(sum!=ans)
	{
		printf("-1\n");
		printf("\n");
		return;
	}
	head[ss]=-1;
	head[tt]=-1;
	printf("%d\n",ISAP(s,t));
	for(i=0;i<cnt;i++)
		printf("%d\n",edge[id[i]^1].cap+down[i]);
	printf("\n");
}
int main()
{
	while(scanf("%d%d",&n,&m)==2)
	{
		build();
	}
	return 0;
}