Poj1275 差分约束

 

题目链接：http://poj.org/problem?id=1275

 

题意：一商店二十四小时营业，但每个时间段需求的雇员数不同，现有n个人申请这份工作，其可以从固定时间t连续工作八个小时，问在满足需求的情况下最小需要多少个雇

员。

 

Discuss from aoxboxcox：

为避免负数，时间计数1～24。

令：

R[i] i时间需要的人数 （1<=i<=24）

T[i] i时间应聘的人数 （1<=i<=24）

x[i] i时间录用的人数 （0<=i<=24），其中令x[0]=0

再设s[i]=x[0]+x[1]+……+x[i] （0<=i<=24），

由题意，可得如下方程组：

(1) s[i]-s[i-8]>=R[i] (8<=i<=24)

(2) s[i]-s[16+i]>=R[i]-s[24] (1<=i<=7)

(3) s[i]-s[i-1]>=0 (1<=i<=24)

(4) s[i-1]-s[i]>=-T[i] (1<=i<=24)

 

这个差分约束有个特殊的地方，(2)的右边有未知数s[24]。这时可以通过枚举s[24]=ans来判断是否有可行解。

即(2)变形为(2') s[i]-s[16+i]>=R[i]-ans (1<=i<=7)再通过SPFA求解(1)(2')(3)(4)。

不过最后有可能出现这种情况：(1)(2')(3)(4)虽然有解，但求出的s[24]小于代入(2')里的ans！

这时，显然得到的s[]不满足原来的(2)了（请仔细比较(2)与(2')）。不过虽然得到的解不满足原方程组，

但这并不代表(1)(2)(3)(4)在s[24]=ans时没有可行解！此外，值得注意的是，当得到的s[24]>ans时，

虽然s[24]不一定是最优解，但把ans置成s[24]后，确实是可行解。所以，简单把(2)置换成(2')是有问题的！

为了等价原命题，必须再加上条件：s[24]>=ans这就是所谓加出来的那条边（5） s[24]－s[0]>=ans最后说一下，

SPFA后判dis[24]==ans其实是没有必要的。

 

Code :

 

可以用二分答案来优化，这里直接从小到大枚举了：

#include<stdio.h>
#include<string.h>
#include<stdlib.h>

typedef struct{
	int v,c,nxt;
}Edge;

#define M 32
#define E 10024

Edge edg[E];
int hed[M];
int head[M];
int dis[M];
int ins[M],inn[M];
int stk[M],top;
int r[M],t[M],e,te;

void AddEdge(int *hed,int u,int v,int c)
{
	edg[e].v=v;edg[e].c=c;
	edg[e].nxt=hed[u];hed[u]=e++;
}

int SPFA(int r)
{
	int u,v,tmp;

	memset(dis,-11,sizeof(dis));
	memset(ins,0,sizeof(ins));
	memset(inn,0,sizeof(inn));
	top=0;dis[0]=0;
	ins[0]=inn[0]=1;
	stk[top++]=0;

	while(top){
		u=stk[--top];ins[u]=0;
		for(tmp=head[u];tmp!=-1;tmp=edg[tmp].nxt){
			v=edg[tmp].v;
			if(dis[v]<dis[u]+edg[tmp].c){
				dis[v]=edg[tmp].c+dis[u];
				if(!ins[v]){
					ins[v]=1;
					inn[v]++;
					stk[top++]=v;
					if(inn[v]>25) 
						return false;
				}
			}
		}
	}
	return true;
}

int main()
{
	int tt,ts;
	int i,x,n,tmp;

	scanf("%d",&tt);
	while(tt--){
		memset(t,0,sizeof(t));
		for(i=1;i<=24;i++)
			scanf("%d",r+i);
		scanf("%d",&ts);n=ts;
		while(ts--){
			scanf("%d",&x);
			t[x+1]++;
		}
		e=0;
		memset(hed,-1,sizeof(hed));

		for(i=1;i<=24;i++){
			if(i>7) AddEdge(hed,i-8,i,r[i]);
			AddEdge(hed,i-1,i,0);
			AddEdge(hed,i,i-1,-t[i]);
		}
		te=e;
		for(tmp=0;tmp<=n;tmp++){
			e=te;
			memmove(head,hed,sizeof(hed));
			for(i=1;i<8;i++)
				AddEdge(head,i+16,i,r[i]-tmp);
			AddEdge(head,0,24,tmp);
			if(SPFA(tmp)) break;
		}
		printf(tmp<=n?"%d\n":"No Solution\n",tmp);
	}
	return 0;
}