Codeforces 939E 三分

题意：有一个可重集$multiset:s$与$m$个操作，这些操作只会是以下两个之一：

$1$.给出一个数字，现在给出的数字一定大于之前给的数字（如果在这之前给出过数字），每次给出后将这个数字加入一个$multiset$中

$2$.求出$s$中的一个子集，使得这个自己中的$max-mean$最小，其中$max$定义为这个子集中的最大元素，$mean$定义为这个可重集的元素均值

$Solution:$

最大元素仅需一个元素贡献即可，那么这个子集中一定有$s$中的最大值，那么该如何降低均值呢？在加入最大值后均值时最大值，需要加入比他小的值来拉低均值，那么首选一定是最小值，贪心的加入，均值一定减小。但是，每个时刻的最小值越加越大，均值越来越小，当某个时刻的最小值大于均值时，这个时候反而拉高均值了，所以当前的最小值大于现在的均值时停止加入即是最优方案。按顺序枚举最小值复杂度是$O(n^2)$的，要想办法优化这个过程，这个过程我们不妨看作函数$f(x)=(sum[x]+a[cnt])/(x+1)$的找最小值的过程（可重集按照升序给出，那么给到第$cnt$个，最大值自然是$a[cnt]$，那么均值就是上面的函数），按照分析，函数是先减后增的，可以用三分加速这个过程。求导二分是难实现的，因为这个函数是嵌套函数，前缀和$sum[x]$也是一个函数。

#include<bits/stdc++.h>
#define ll long long
using namespace std;

ll read()
{
	ll ret=0,base=1;
	char ch=getchar();
	while(!isdigit(ch))
	{
		if(ch=='-') base=-1;
		ch=getchar();
	}
	while(isdigit(ch))
	{
		ret=(ret<<3)+(ret<<1)+ch-48;
		ch=getchar();
	}
	return ret*base;
}

ll a[500005],sum[500005];
int cnt;

double f(int x)
{
	return (double)(sum[x]+a[cnt])/(x+1);
}

double solve()
{
//	if(cnt==1) return 0;
	ll max1=a[cnt];double ret=0;
	//第x位的函数值
	//(sum[x]+max)/(x+1) 函数先减后增
	int l=1,r=cnt-1;
	while(l<=r)
	{
		int mid1=(2*l+r)/3,mid2=(2*r+l)/3;
		double v1=f(mid1),v2=f(mid2);
		ret=max(ret,max1-v1);
		ret=max(ret,max1-v2);
		if(v1<v2) r=mid2-1;
		else l=mid1+1;
	}
	return ret;
}

int main()
{
	int q=read();
	while(q--)
	{
		if(read()==1)
		{
			a[++cnt]=read();
			sum[cnt]=sum[cnt-1]+a[cnt];
		}
		else printf("%lf\n",solve());
	}
 	return 0;
}