POJ 1275 Cashier Employment

又是一道以区间和作为定点的差分约束，看了冯威的论文才会的，不过他的论文里漏了一个条件

 

设s[i]为[0,i]雇佣的人数和，num[i]为i时刻可雇佣人数，need[i]为i时刻需要人数，则

0<=s[i]-s[i-1]<=num[i]

 

s[i]-s[i-8]>=need[i] (8<=i<=23)

 

s[i]+(s[23]-s[i+16])>=need[i] (0<=i<=8)

 

另外,s[23]-s[0]>=sum(总雇佣人数)

 

第三个等式化为s[i]-s[i+16]>=need[i]-s[23],只要二分s[23]，即最后答案就可以解决问题

 

代码：

#include<iostream> #include<queue> using namespace std; const int MAX=25; const int inf=1<<30; struct node { int v,w,next; }g[MAX*100]; int num[MAX],need[MAX],adj[MAX],dis[MAX],cnt[MAX],vis[MAX]; int e,n; void add(int u,int v,int c) { g[e].v=v; g[e].w=c; g[e].next=adj[u]; adj[u]=e++; } int spfa() { int i,u,v; queue<int>q; for(i=1;i<=24;i++) dis[i]=-inf; memset(vis,0,sizeof(vis)); memset(cnt,0,sizeof(cnt)); dis[0]=0; q.push(0); vis[0]=1; cnt[0]++; while(!q.empty()) { u=q.front(); q.pop(); vis[u]=0; for(i=adj[u];i!=-1;i=g[i].next) { v=g[i].v; if(dis[v]<dis[u]+g[i].w) { dis[v]=dis[u]+g[i].w; if(!vis[v]) { q.push(v); vis[v]=1; cnt[v]++; if(cnt[v]>24) return 0; } } } } return 1; } int ok(int x) { int i; e=0; //cout<<"yes"<<endl; memset(adj,-1,sizeof(adj)); add(0,24,x); for(i=1;i<=24;i++) { add(i-1,i,0); add(i,i-1,-num[i]); } for(i=9;i<=24;i++) add(i-8,i,need[i]); for(i=1;i<=8;i++) add(i+16,i,need[i]-x); if(spfa()) return 1; else return 0; } int main() { int i,j,T; scanf("%d",&T); while(T--) { memset(num,0,sizeof(num)); e=0; for(i=1;i<=24;i++) scanf("%d",&need[i]); scanf("%d",&n); for(i=1;i<=n;i++) { scanf("%d",&j); num[j+1]++; } int l=1,r=n,mid,ans=-1; while(l<=r) { mid=(l+r)/2; if(ok(mid)) { r=mid-1; ans=mid; } else l=mid+1; } if(ans!=-1) printf("%d/n",ans); else puts("No Solution"); } return 0; }