[SCOI2016]妖怪

嘟嘟嘟


离NOI最后一周，把自己容易忘的知识点和板子复习一下。
（刚答完loj的笔试模拟，感觉上不了90……）
update：哦，我89……


先把式子写出来，每一个妖怪的战斗力(S(i) = A + frac{a}{b}D +D +frac{b}{a}A)。
令(k = frac{a}{b})，于是(S(i) = D *k + A *frac{1}{k} +A +D)。
这东西是一个对勾函数，在第一象限是单峰的。我们要求的是一堆对勾函数的最大值的最小值。有一个结论是单峰函数复合起来还是单峰函数（但是我不会证啊），所以就可以三分了。


刚开始我直接傻呵呵的令(L = 0, R = 1e10)，然后只得了10分。
应该把边界算准了：当(D *k = A *frac{1}{k})的时候，(k = sqrt{frac{A}{D}})，所以应该是(1e-4 sim 1e4)。
然后面向数据编程，我三分了115次才过掉这道题……

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<queue>
#include<vector>
#include<ctime>
#include<assert.h>
using namespace std;
#define enter puts("")
#define space putchar(' ')
#define Mem(a, x) memset(a, x, sizeof(a))
#define In inline
#define forE(i, x, y) for(int i = head[x], y; (y = e[i].to) && ~i; i = e[i].nxt)
typedef long long ll;
typedef double db;
const db INF = 1e18;
const db eps = 1e-10;
const int maxn = 1e6 + 5;
In ll read()
{
	ll ans = 0;
	char ch = getchar(), las = ' ';
	while(!isdigit(ch)) las = ch, ch = getchar();
	while(isdigit(ch)) ans = (ans << 1) + (ans << 3) + ch - '0', ch = getchar();
	if(las == '-') ans = -ans;
	return ans;
}
In void write(ll x)
{
	if(x < 0) putchar('-'), x = -x;
	if(x >= 10) write(x / 10);
	putchar(x % 10 + '0');
}
In void MYFILE()
{
#ifndef mrclr
	freopen("monster.in", "r", stdin);
	freopen("ha.out", "w", stdout);
#endif
}

int n;
db A[maxn], D[maxn];

In db calc(db x)
{
	db ret = 0;
	for(int i = 1; i <= n; ++i) ret = max(ret, A[i] + D[i] + x * D[i] + 1.0 / x * A[i]);
	return ret;
}

int main()
{
//	MYFILE();
	n = read();
	for(int i = 1; i <= n; ++i) A[i] = read(), D[i] = read();
	db L = 1e-4, R = 1e4;
	for(int i = 1; i <= 115; ++i)
	{
		db mid1 = (L + R) * 0.5, mid2 = (mid1 + R) * 0.5;
		if(calc(mid1) < calc(mid2)) R = mid2;
		else L = mid1;
	}
	printf("%.4lf
", min(calc(L), calc(R)));
	return 0;
}