1. 程式人生 > >poj1275 (二分+差分約束)

poj1275 (二分+差分約束)

題意:
在一家超市裡,每個時刻都需要有營業員看管,R(i) (0 <= i < 24)表示從i時刻開始到i+1時刻需要的營業員的數目,現在有N(N <= 1000)個申請人申請這項工作,並且每個申請者都有一個起始工作時間 ti,如果第i個申請者被錄用,那麼他會從ti時刻開始連續工作8小時。現在要求選擇一些申請者進行錄用,使得任何一個時刻i,營業員數目都能大於等於R(i)。求出至少需要錄用多少營業員。

分析:
每天的工作情況都是一樣的,我們只需要求出一天的即可。根據題意,令s[i]為一天內前i+1個小時錄用的人數,如果i>=7,則s[i] - s[i-8] >= R[i],如果0 <= i < 7,則可以推出s[23] - s[i+16] + s[i] >= R[i],同時每個時刻錄用的人數有個上限,假設第i時刻最多可以錄用的人數為b[i],則對於每一個i有0 <=s[i] - s[i-1] <= b[i]。
現在需要解決的一個問題是,第二個不等式中包含3個s組的變數,這該怎麼建圖呢?可以知道我們只需要列舉s[23],那麼這個量就是已知的了,因此不等式可以變為s[i] - s[i+16] >= R[i] - s[23],但是必須明白的一點是,既然s[23]是列舉的一天的錄用人數的最小數目,我們建圖之後求出的s[23]也應該為列舉的那個值,可以從0到n列舉s[23],第一個值便是答案,但是可以更高效地求解,因為問題具有單調性,直接二分.

因為s[-1] = 0,-1無法表示,只需要向右移動一個單位即可,那麼令s[0] = 0,思路很清晰了,二分s[24]搜尋最小值。

程式碼:

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <stdlib.h>
#include <algorithm>
#include <math.h>
#include <vector>
#include <queue>

using namespace std;

const
int inf =0x3f3f3f3f; int R[30],time[30],cnt,head[30],d[30],vis[30]; queue<int>q; void init() { cnt=0; memset(head,-1,sizeof(head)); memset(d,-inf,sizeof(d)); memset(vis,0,sizeof(vis)); while(!q.empty()) q.pop(); } struct Edge { int v,w,next; Edge() {}; Edge(int _v,int
_w,int _next) { v=_v,w=_w,next=_next; } } edge[100000]; void add_edge(int u,int v,int w) { edge[cnt]=Edge(v,w,head[u]); head[u]=cnt++; } int spfa(int s) { d[0]=0; q.push(0); vis[0]=1; while(!q.empty()) { int u=q.front(); q.pop(); vis[u]=0; if(u==24&&d[u] > s)//不加超時 return 0; for(int i=head[u];i!=-1; i=edge[i].next) { int v=edge[i].v; if(d[v]<d[u]+edge[i].w) { d[v]=d[u]+edge[i].w; if(!vis[v]) vis[v]=1,q.push(v); } } } return d[24]==s? 1: 0; } int main() { int T,n; scanf("%d",&T); while(T--) { memset(time,0,sizeof(time)); for(int i=0; i<=23; i++) scanf("%d",&R[i]); scanf("%d",&n); int tmp; for(int i=0; i<n; i++) { scanf("%d",&tmp); time[tmp]++; } int l =-1,r =n+1,m; while(r-l>1)//都往右移一個單位 { m=(r+l)/2; init(); for(int i=7; i<=23; i++) add_edge(i-7,i+1,R[i]); for(int i=0; i<7; i++) add_edge(17+i,i+1,R[i]-m); for(int i=0;i<=23; i++) { add_edge(i,i+1,0); add_edge(i+1,i,-time[i]); } add_edge(0,24,m); add_edge(24,0,-m); spfa(m)? r=m:l=m; } l >= n ? puts("No Solution") : printf("%d\n", r); } return 0; }