poj1275 (二分+差分約束)
阿新 • • 發佈:2019-01-08
題意:
在一家超市裡,每個時刻都需要有營業員看管,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;
}