POJ1275 Cashier Employment[差分约束系统 || 单纯形法]

Cashier Employment

Time Limit: 1000MS		Memory Limit: 10000K
Total Submissions: 7997		Accepted: 3054

Description

A supermarket in Tehran is open 24 hours a day every day and needs a number of cashiers to fit its need. The supermarket manager has hired you to help him, solve his problem. The problem is that the supermarket needs different number of cashiers at different times of each day (for example, a few cashiers after midnight, and many in the afternoon) to provide good service to its customers, and he wants to hire the least number of cashiers for this job.

The manager has provided you with the least number of cashiers needed for every one-hour slot of the day. This data is given as R(0), R(1), ..., R(23): R(0) represents the least number of cashiers needed from midnight to 1:00 A.M., R(1) shows this number for duration of 1:00 A.M. to 2:00 A.M., and so on. Note that these numbers are the same every day. There are N qualified applicants for this job. Each applicant i works non-stop once each 24 hours in a shift of exactly 8 hours starting from a specified hour, say ti (0 <= ti <= 23), exactly from the start of the hour mentioned. That is, if the ith applicant is hired, he/she will work starting from ti o'clock sharp for 8 hours. Cashiers do not replace one another and work exactly as scheduled, and there are enough cash registers and counters for those who are hired.

You are to write a program to read the R(i) 's for i=0..23 and ti 's for i=1..N that are all, non-negative integer numbers and compute the least number of cashiers needed to be employed to meet the mentioned constraints. Note that there can be more cashiers than the least number needed for a specific slot.

Input

The first line of input is the number of test cases for this problem (at most 20). Each test case starts with 24 integer numbers representing the R(0), R(1), ..., R(23) in one line (R(i) can be at most 1000). Then there is N, number of applicants in another line (0 <= N <= 1000), after which come N lines each containing one ti (0 <= ti <= 23). There are no blank lines between test cases.

Output

For each test case, the output should be written in one line, which is the least number of cashiers needed. 
If there is no solution for the test case, you should write No Solution for that case. 

Sample Input

1
1 0 1 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1
5
0
23
22
1
10

Sample Output

1

Source

Tehran 2000

---

 德黑兰的一家每天24小时营业的超市，需要一批出纳员来满足它的需求。超市经理雇佣你来帮他解决一个问题————超市在每天的不同时段需要不同数目的出纳员（例如，午夜只需一小批，而下午则需要很多）来为顾客提供优质服务，他希望雇佣最少数目的纳员。
            超市经历已经提供一天里每一小时需要出纳员的最少数量————R(),R(),...,R(）。R（）表示从午夜到凌晨1:00所需要出纳员的最少数目；R()表示凌晨1：00到2：00之间需要的；等等。每一天，这些数据都是相同的。有N人申请这项工作，每个申请者i在每天24小时当中，从一个特定的时刻开始连续工作恰好8小时。定义ti(<=ti<=)为上面提到的开始时刻，也就是说，如果第i个申请者被录用，他（或她）将从ti时刻开始连续工作8小时。
            试着编写一个程序，输入R(i),i=，...，，以及ti，i=，...，N，它们都是非负整数，计算为满足上述限制需要雇佣的最少出纳员数目、在每一时刻可以有比对应R(i)更多的出纳员在工作
输入描述:
            输入文件的第1行为一个整数T，表示输入文件中测试数据的数目(至多20个)。每个测试数据第一行为24个整数，表示R（），R（），...，R(),R(i)最大可以取到1000。接下来一行是一个整数N，表示申请者的数目，<=N<=。接下来有N行，每行为一个整数ti，<=ti<=,测试数据之间没有空行。
输出描述：
           对输入文件中的每个测试数据，输出占一行，为需要雇佣的出纳员的最少数目。如果某个测试数据没有解。则输出"No Solution"。

题意

---

看到后第一反应：这不是noi2008志愿者吗，只不过没有了每个人的花费，最小化人数，花费设为1

一个时间段一个约束 24个，一个人一个变量 1000个，x为每个人选不选，A[i][j]=1当tj<=i<tj+8(注意时间循环)

最小化 x

满足约束 Ax<R

然后有很多细节问题：

1.我靠多组数据，初始化a[][] v

2.我靠连续8小时可以循环，24之后是1，改(我的是从1到24)

3.我靠还有无解的情况，判断选所有人可不可行

4.我靠无解用的sum数组也要清空

然后终于AC了，16ms（一片0ms），然而思维难度大大降低

一点总结：这种样子的题目都碰到3个了.......

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
typedef long long ll;
const int N=,M=;
const double eps=1e-,INF=1e9;
inline int read(){
    char c=getchar();int x=,f=;
    while(c<''||c>''){if(c=='-')f=-;c=getchar();}
    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}
    return x*f;
}
int n=,m,t;
double a[M][N],b[M],c[N],v;int sum[N];
void init(){
    memset(a,,sizeof(a));
    memset(sum,,sizeof(sum));
    v=0.0;
}
void pivot(int l,int e){
    b[l]/=a[l][e];
    for(int j=;j<=n;j++) if(j!=e) a[l][j]/=a[l][e];
    a[l][e]=/a[l][e];
    for(int i=;i<=m;i++) if(i!=l&&fabs(a[i][e])>){
        b[i]-=a[i][e]*b[l];
        for(int j=;j<=n;j++) if(j!=e) a[i][j]-=a[i][e]*a[l][j];
        a[i][e]=-a[i][e]*a[l][e];
    }
    v+=c[e]*b[l];
    for(int j=;j<=n;j++) if(j!=e) c[j]-=c[e]*a[l][j];
    c[e]=-c[e]*a[l][e];
}
double simplex(){
    while(true){
        int e=,l=;
        for(e=;e<=n;e++) if(c[e]>eps) break;
        if(e==n+) return v;
        double mn=INF;
        for(int i=;i<=m;i++)
            if(a[i][e]>eps&&mn>b[i]/a[i][e]) mn=b[i]/a[i][e],l=i;
        if(mn==INF) return INF;//unbounded
        pivot(l,e);
    }
}

bool check(){
    for(int i=;i<=n;i++) if(sum[i]<c[i]) return ;
    return ;
}
int main(){
    //freopen("in.txt","r",stdin);
    int T=read();
    while(T--){
        init();
        for(int i=;i<=n;i++) c[i]=read();
        m=read();
        for(int i=;i<=m;i++){
            t=read()+;
            for(int j=;j<=;j++){
                if(t==) t=;//printf("hi %d\n",t);
                a[i][t]=; sum[t]++;
                t++;
            }
            b[i]=;
        }
        if(!check()){puts("No Solution");continue;}
        simplex();
        printf("%d\n",(int)(v+0.5));
    }
}