poj1275 （二分+差分约束）

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

借鉴了一位大佬的博客： https://blog.csdn.net/consciousman/article/details/53812818

题意：
在一家超市里，每个时刻都需要有营业员看管，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;
}