BZOJ 3280: 小R的烦恼 费用流 多路增广

 

3280: 小R的烦恼

Time Limit: 10 Sec   Memory Limit: 512 MB
Submit: 486   Solved: 249
[ Submit][ Status][ Discuss]

Description

小R最近遇上了大麻烦，他的程序设计挂科了。于是他只好找程设老师求情。善良的程设老师答应不挂他，但是要求小R帮助他一起解决一个难题。

问题是这样的，程设老师最近要进行一项邪恶的实验来证明P=NP，这个实验一共持续n天，第i天需要a[i]个研究生来给他搬砖。研究生毕竟也是人，所以雇佣研究生是需要钱的，机智的程设老师已经联系好了m所大学，第j所大学共有l[j]个研究生，同时雇佣这所大学的一个研究生需要p[j]元钱。

本来程设老师满心欢喜的以为，这样捡最便宜的max{a[i]}个研究生雇来，就可以完成实验；结果没想到，由于他要求硕士生们每天工作25个小时不许吃饭睡觉上厕所喝水说话咳嗽打喷嚏呼吸空气，因此一天下来给他搬砖的所有研究生都会进入濒死状态。濒死状态的研究生，毫无疑问，就不能再进行工作了。但是机智的老师早早联系好了k家医院，第i家医院医治一个濒死的研究生需要d[i]天,并且需要q[i]元钱。

现在，程设老师想要知道，最少花多少钱，能够在这n天中满足每天的需要呢？若无法满足，则请输出”impossible”。注意，由于程设老师良心大大的坏，所以他是可以不把濒死的研究生送去医院的！

Input

       本题包含多组数据；第一行是一个数T(T<=11)，表示数据组数，以下T组数据。
对于每一组数据，第一行三个数，n,m,k;
以下一行n个数，表示a[1]…a[n]
接着一行2m个数，表示l[1],p[1]…l[n],p[n]
接着一行2k个数，表示d[1],q[1]…d[n],q[n]

Output

       对于每组数据以样例的格式输出一行，两个数分别表示第几组数据和最少钱数。

Sample Input

2
3 2 1
10 20 30
40 90 15 100
1 5
3 2 1
10 20 30
40 90 15 100
2 5

Sample Output

Case 1: 4650
Case 2: impossible

HINT

样例解释：买下90块钱的那40个研究生，另外再买10个100块钱的。这样，第一天用完的10个人全部送到医院，那么他们在第三天可以继续使用；同时，第二天和第三天都用新的研究生来弥补，这样一共需要花费40*90 + 10*100 + 5*10 = 4650元。
数据规模：
对于30%的数据中的每组数据,
满足n<=5,m,k<=2,其余数均小于等于100或者 
n<=10,m,k<=10,其余数均小于等于20.
对于100%的数据
n,m,k<=50,其余数均小于等于100.

 

建图方式和BZOJ 1221: [HNOI2001] 软件开发 费用流 多路增广几乎一模一样

因为忘记重置ans和maxflow WA一发

 

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<string>
#include<bitset>
#include<queue>
#include<map>
#include<set>
using namespace std;
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=10*x+ch-'0';ch=getchar();}
	return x*f;
}
inline void print(int x)
{if(x<0)putchar('-'),x=-x;if(x>=10)print(x/10);putchar(x%10+'0');}

const int N=20010,M=300010,inf=0X3f3f3f3f;

int ecnt=1,last[N];
struct EDGE{int to,nt,val,ct;}e[M];
inline void readd(int u,int v,int val,int ct){e[++ecnt]=(EDGE){v,last[u],val,ct};last[u]=ecnt;}
inline void add(int u,int v,int val,int ct){readd(u,v,val,ct);readd(v,u,0,-ct);}

int n,S,T=N-1;

int q[N],dis[N],ans=0;
bool inq[N],vis[N];
bool spfa()
{
	memset(dis,0X3f,sizeof(dis));
	memset(inq,0,sizeof(inq));
	memset(vis,0,sizeof(vis));
	int head=0,tail=1;q[0]=0;inq[0]=1;dis[0]=0;
	while(head<tail)
	{
		int u=q[head++];head%=N;inq[u]=0;
		for(int i=last[u];i;i=e[i].nt)
		if(e[i].val&&dis[e[i].to]>dis[u]+e[i].ct)
		{
			dis[e[i].to]=dis[u]+e[i].ct;
			if(!inq[e[i].to])q[tail++]=e[i].to,tail%=N,inq[e[i].to]=1;
		}
	}
	return dis[T]<inf;
}

int dfs(int u,int lim)
{
	if(u==T)
	{ans+=dis[u]*lim;return lim;}
	vis[u]=1;int res=0;
	for(int i=last[u];i;i=e[i].nt)
	if(e[i].val&&!vis[e[i].to]&&dis[e[i].to]==dis[u]+e[i].ct)
	{
		int tmp=dfs(e[i].to,min(lim,e[i].val));
		lim-=tmp;e[i].val-=tmp;e[i^1].val+=tmp;res+=tmp;
		if(!lim)break;
	}
	if(!res)dis[u]=-1;
	return res;
}

int flow=0;
void mincf(){while(spfa())flow+=dfs(S,inf);}

void initial()
{memset(last,0,sizeof(last));ecnt=1;flow=0;ans=0;}

int main()
{
	int t=read(),ti=0;
	while(t--)
	{
		n=read();int m=read(),k=read();int sum=0;ti++;
		for(int i=1,x;i<=n;++i)x=read(),add(S,i,x,0),add(i+n,T,x,0),sum+=x;
		for(int i=1;i<n;i++)add(i,i+1,inf,0);
		for(int i=1,x,y;i<=m;++i)
		{
			x=read(),y=read();add(S,i+n+n,x,0);
			for(int j=1;j<=n;++j)add(i+n+n,j+n,inf,y);
		}
		for(int i=1,x,y;i<=k;++i)
		{
			x=read();y=read();
			for(int j=1;j<=n;++j)if(j+x+1>n)break;else add(j,j+n+x+1,inf,y);
		}
		printf("Case %d: ",ti);
		mincf();if(flow==sum){print(ans);puts("");}else puts("impossible");
		initial(); 
	}
	return 0;
}
/*
2
3 2 1
10 20 30
40 90 15 100
1 5
3 2 1
10 20 30
40 90 15 100
2 5

Case 1: 4650
Case 2: impossible
*/