Google Code Jam 2015 Round 1A: Problem B. Haircut

蒟蒻看到这一题只想到了模拟==看了analysis后才恍然大悟竟然是二分==其实单调性也挺明显的额肯定是时间越长serve的人越多。><

首先，对于每一个时间点time，可以求出来此刻serve的人数。num=sum(time/m[i]+1)，+1是因为time和m[i]不能整除时，time时刻还在serve一个人，time和m[i]可以整除时，time时刻即将serve一个人。所以此时求出的sum是个上界，会比N大。

二分的上界是max(m[i])*N，即只有一个barber而且serve时间最长的case。注意二分的上界最大可以是1e15，所以计算sum时要用long long存，我之前脑残一看N才是1e9，然后就用int存sum了==然后WA了好久。==

求出time后，官方的analysis是说当前时刻会同时serve多个人，所以需要先相减求出该时刻serve的总人数然后遍历m[i]看哪个正好到N就是ans。

后来偷偷滴看了dreamoon男神的代码后发现，也可以直接遍历m[i]，保存恰好整除的barber，然后ans就是第倒数sum-N个。

#include<iostream>
#include<stdio.h>
#include<cstdio>
#include<stdlib.h>
#include<vector>
#include<string>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<stack>
#include<queue>
#include<ctype.h>
#include<map>
#include<time.h>
#include<bitset>
#include<set>
#include<list>
using namespace std;
//gcj R1A pro B 

const int maxn=1010;
int T;
int N;
int m[maxn];
int B;
int upbound;
int ans;
long long servednumber(long long time)//因为time是从1e15开始二分的，所以返回int会溢出
{
    if(time<0)
    {
        return 0;
    }
    long long ret=0;
    for(int i=1;i<=B;i++)
    {
        ret+=time/m[i]+1;//what if %m[i]==0?
    }
    return ret;
}
long long binarysearch()
{
    long long l=0;
    long long r=(long long)upbound*(long long)N;
    while(l<r)
    {
        long long mid=(l+r)/2;
        if(servednumber(mid)<N)
        {
            l=mid+1;
        }
        else
        {
            r=mid;
        }
    }
    //cout<<"time: "<<l<<endl;
   // cout<<l<<" "<<r<<endl;
    return min(l,r);
    //return r;
}
void solve()
{
    long long time=binarysearch();
//    long long sum=0; //it's also fine
//    vector<int>d;
//    for(int i=1;i<=B;i++)
//    {
//        sum+=time/m[i]+1;
//        if(time%m[i]==0)
//        {
//            d.push_back(i);
//        }
//    }
    //cout<<sum<<" "<<N<<endl;
//ans=d[d.size()-(sum-N)-1];
    int servedbefore=servednumber(time-1);
    int servednow=N-servedbefore;
    for(int i=1;i<=B;i++)
    {
        if(time%m[i]==0)
        {
            servednow--;
            if(servednow==0)
            {
                ans=i;
                break;
            }
            else if(servednow<0)
            {
                ans=i;
                break;
            }
        }
    }
}
int main()
{

    //freopen("input.txt","r",stdin);
    freopen("B-large-practice.in","r",stdin);
    freopen("out1.txt","w",stdout);
    scanf("%d",&T);
    for(int ca=1;ca<=T;ca++)
    {
        memset(m,0,sizeof(m));
        upbound=0;
        ans=0;
        scanf("%d %d",&B,&N);
        for(int i=1;i<=B;i++)
        {
            scanf("%d",&m[i]);
            upbound=max(upbound,m[i]);
        }
        solve();
        printf("Case #%d: %d\n",ca,ans);
    }
    return 0;
}