POJ 1275 Cashier Employment(差分约束系统)

摘自冯威论文——《数与图的完美结合》

设num[ i ]为i时刻能够开始工作的人数，x[ i ]为实际雇佣的人数，那么x[ I ]<=num[ I ]。
设r[ i ]为i时刻至少需要工作的人数，于是有如下关系：
    x[ I-7 ]+x[ I-6 ]+x[ I-5 ]+x[ I-4 ]+x[ I-3 ]+x[ I-2 ]+x[ I-1 ]+x[ I ]>=r[ I ]
设s[ I ]=x[ 1 ]+x[ 2 ]…+x[ I ]，得到
    0<=s[ I ]-s[ I-1 ]<=num[ I ]， 0<=I<=23
    s[ I ]-s[ I-8 ]>=r[ I ]， 8<=I<=23
    s[ 23 ]+s[ I ]-s[ I+16 ]>=r[ I ]， 0<=I<=7

    对于以上的几组不等式，我们采用一种非常笨拙的办法处理这一系列的不等式(其实也是让零乱的式子变得更加整齐，易于处理)。首先我们要明白差分约束系统的应用对象（它通常针对多个二项相减的不等式的）于是我们将上面的所有式子都转化成两项未知项在左边，另外的常数项在右边，且中间用>=连接的式子，即:
    s[ I ]-s[ I-1 ]>=0            (0<=I<=23)
    s[ I-1 ]-s[ I ]>=-num[ I ]       (0<=I<=23)
    s[ I ]-s[ I-8 ]>=r[ I ]         (8<=I<=23)
    s[ I ]-s[ I+16 ]>=r[ I ]-s[ 23 ]  (0<=I<= 7)

 

论文有点小小的漏洞，也导致了它论文附带的程序是错误的，有BUG

可行方案中还差一个约束条件

S[23] - S[-1] >= sum;sum为雇佣的出纳员总数

 

将所有形如 A-B>=C 的式子 我们从节点B 引出一条有向边指向A边的权值为C  （这里注意由于左右确定，式子又是统一的>=的不等式，所以A和B是相对确定的，边是一定是指向A的） ，图就建成了 。

在程序中，我用0作为上面说的-1来构图

//差分约束系统(SPFA) //对于差分约束系统，建立的不等式是>=则在SPFA中求最长路 //反之求最短路 //注意两者对DIS初始化的不同 #include<iostream> #include<queue> #define INF 1000000000 using namespace std; int r[25],t[25],caseN,n,T; int head[25],V[700],next[700],W[700]; int cnt[25],dis[25]; bool inq[25]; int m; inline void addEdge(int u,int v,int w)//新学的数组存放邻接表的方式 { V[m] = v; W[m] = w; next[m] = head[u];//head为表头，调试后就能理解，这是倒着插的过程 head[u] = m++; } int SPFA(int ans) { memset(inq,0,sizeof(inq));//元素是否在队列中 memset(cnt,0,sizeof(cnt));//记录结点入队列的次数 for(int i = 1;i <= 24;++i) dis[i] = -INF; queue<int> q; q.push(0); inq[0] = 1; dis[0] = 0; cnt[0]++; while(!q.empty())//SPFA+邻接表模板 { int u = q.front(); q.pop(); inq[u] = 0; for(int e = head[u];e != -1;e = next[e]) { if(dis[u] + W[e] > dis[V[e]]) { dis[V[e]] = dis[u] + W[e]; if(!inq[V[e]]) { q.push(V[e]); inq[V[e]] = 1; if(++cnt[V[e]] > 24) return -1;//结点入度次数，如果超过|V|则说明存在回路 } } } } if(dis[24] == ans) return 1;//如果不存在正权回路且S[24]恰好等于ANS,则说明找到可行解，切ANS为最优解 else return 0; } void buildGraph(int ans) { m = 0;//重新建图时，边数记得回零 memset(head,-1,sizeof(head)); addEdge(0,24,ans); for(int i = 1;i <= 24;++i) { addEdge(i - 1,i,0); addEdge(i,i - 1,-t[i]); } for(int i = 1;i <= 8;++i) addEdge(i + 16,i,r[i] - ans); for(int i = 9;i <= 24;++i) addEdge(i - 8,i,r[i]); } int main() { //freopen("in.txt","r",stdin); scanf("%d",&caseN); bool ok; while(caseN--) { memset(t,0,sizeof(t)); ok = 0; for(int i = 1;i <= 24;++i) scanf("%d",&r[i]); scanf("%d",&n); for(int i = 0;i < n;++i) { scanf("%d",&T); t[T+1]++; } for(int ans = 0;ans <= n;++ans) { buildGraph(ans); if(SPFA(ans) > 0) { printf("%d/n",ans); ok = 1; break; } } if(!ok) printf("No Solution/n"); } return 0; }