蒟蒻看到这一题只想到了模拟==看了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;
}