D. Gambler‘s Ruin

D. Gambler’s Ruin

题目链接

题目大意

需要设置两个赔率x和y。使得求出一个最佳的赔率，使得赌博公司利润最大
$$使得[s_x+s_y-max(s_x\ast x,s_y\ast y)]最大$$

题解思维

$$当某个（赌徒）p_i的赔率满足p_i\ast x>=1时，也就是x>=\frac{1}{p_i}时,赌徒会下注c_i$$

$$另一个则是y的，同理得到，y>=\frac{1}{1-p_i}时,赌徒会下注c_i$$

$$由于0<=p_i<=1，所以当pi越来越小时，\frac{1}{pi}越来越大，所以我们需要对输入的数据根据p_i进行排序【本文中使用了降序排序】$$

$$假设排序完成后的数组，中元素是以a_0,a_1,\cdot \cdot \cdot ,a_n，这样的数组$$

$$如果有p_i相同的元素，可以将它们合并$$

$$我们假设设置了\frac{1}{a_4}为赔率，那么在它【a_4】之前的元素，显然都满足下注的条件$$

$$由于y是1-p_i关联的，显然和x相反，也就是，当认定\frac{1}{1-aj}时，a_{之}后的都满足条件$$

$$问题就简化成了，找到a_i,a_j使得，使得赌博公司利润最大$$

$$对于s_x\ast x，其实可以理解为，这个排序好后的前缀和和(\frac{1}{a_i})的乘积=v1=(\frac{1}{a_i})\ast \sum _{i=1}^i{c}_i$$

$$对于s_y\ast y，其实可以理解为，这个排序好后的后缀和和(\frac{1}{1-a_j})的乘积=v2=(\frac{1}{1-a_j})\ast \sum _{i=j}^{n-j}{c}_i$$

发现，当i，j越来越大时，v1是越来越大，v2越来越小
$$由于存在max(s_x\ast x,s_y\ast y),那么显然，当固定i或j时，另外一个的变化会有一个零界点，当它越过零界点时，max的值会产生变化$$

$$假设固定了某个a_i,使得Sx\ast x固定，假设p从数组a_n开始，a_j是从a_n开始倒着的遍历,希望找到第一个a_j使得它大于或等于Sx\ast x$$

$$假设这个a_j的索引是p,那么当固定a_i（s_x\ast x）时，所能求得的最大利润显然，在p时[大的是y]，或者p+1时【大的是x】产生。$$

有了上面的思路之后，就可以进行程序的编写了了，程序如下

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

int n;
#define ll long long
const int lim  = 1e6+10;
const double eps = 1e-10;
 
pair<double, ll> a[lim];
ll pres[lim], sufs[lim];
int main()
{
	 ios::sync_with_stdio(false);
    cin.tie(nullptr), cout.tie(nullptr);
	cin>>n;	
	for(int i=1;i<=n;i++)
	{
		cin>>a[i].first>>a[i].second;
	}
	sort(a+1,a+n+1,greater<pair<double,ll>>());
	int cnt=0;
	for(int i=1;i<=n;i++)
	{
		if(i==1||fabs(a[i].first - a[i - 1].first) > eps)	
			a[++cnt] = a[i];
		else
			a[cnt].second += a[i].second;
	}
	n = cnt;
	for (int i = 1; i <= n; ++i) pres[i] = pres[i - 1] + a[i].second;
    for (int i = n; i; --i) sufs[i] = sufs[i + 1] + a[i].second;
    
    double ans = 0;
    //找到临界点，并将它和答案进行更新 
    for (int i = 1, p = n; i <= n; ++i) {
		if (i == n || fabs(a[i].first - a[i + 1].first) > eps) 
			while (p > i && pres[i] / a[i].first > sufs[p] / (1 - a[p].first)) --p;
	    ans = max(ans, pres[i] + sufs[p + 1] - pres[i] / a[i].first);
	    if (i >= p) break;
	    ans = max(ans, pres[i] + sufs[p] - sufs[p] / (1 - a[p].first));
    }
    cout << fixed << setprecision(10) << ans << '\n';
}