ABC 289 G - Shopping in AtCoder store 数学推导+凸包

大意：
n个顾客，每个人有一个购买的欲望bi,m件物品，每一件物品有一个价值ci,每一个顾客会买商品当且仅当bi+ci>=定价.

现在要求对每一个商品定价，求出它的最大销售值（数量*定价）

n,m<=2e5

思路：

首先m件物品都是相互独立的，可以看成m个询问

另外，不妨对n个人的购买力做一个降序排序，显然它们满足单调性

不难发现，一旦我们定下了物品的价格，那么最终的销售额就由销售数量决定，也就是会有多少人买。此时在销售数量减少的情况下，我们一定会尽可能地抬高价格。从而我们得到一个结论:每一个商品i的定价一定是bj+ci(1<=j<=n).这是因为，它刚好可以使某个人会买这件商品。假设最优定价不满足这个结论，显然我们可以直接抬高它使其达到另一个bj+ci，此时我们在购买人数不变的情况下就提高了单价，这是更优的。

现在商品单价就只有n个选择了，对于商品i，我们的销售额就是max{j*(bj+ci)}(1<=j<=n),因为我们是按购买力降序排序，如果第j个人刚好买的起，那么前面的人一定都买的起（这里也不用关心购买力重复的问题，因为重复的话，后面相同购买力的的人对应的决策一定会更好）。

此时我们就转化了题意：对于一个i(1<=i<=m),找到max{j*(bj+ci)}

---

继续

 如果放到坐标系下来看的话（遇到双变量时，放到坐标系下往往是个不错的选择），我们令横坐标为ci（选择i也是可以的，因为i跟ci是一一对应的，但是考虑到我们是在对不同的价值做取舍，取ci会更加方便一点），纵坐标为对应的价值 ， 也就是 j*(bj+ci)，不同j的选择对应不同的总价值。

我们将ci表示成x，那么原本要求的值 max{j*(bj+ci)}(1<=j<=n) 转化为max{j*(bj+x)}=max{jx+j*bj}，可以发现里面是一条直线的表达式。显然我们最后是要找一个这些直线在每一个横坐标下对应的最大值，也就是转化成了一个凸包。最后的答案就是横坐标对应的凸包上的点的纵坐标了

这里借一下官方题解的图片：

求凸包的话，我们从1-n开始遍历的话，直线的斜率是单调递增的，

 不妨用单调栈来更新当前段的最大值对应的直线

假设当前栈内有两条直线L0,L1,交点为X_01,那么对于新加入的直线L',如果它与L0的交点X_1'的横坐标小于X_01的横坐标，显然就可以把L1淘汰掉了，因为它后面也不会有比L‘更大的机会了。

关于这个判断，就只要计算一下交点横坐标就可以了。

最后我们得到了一个凸包，那么对于一个ci，我们去二分找到它在那一段线段上就可以了

时间复杂度O(n+mlogn)

code

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define endl '\n'
#define mk make_pair
const ll N=2e5+10;
ll n,m;
ll b[N];
ll c[N];
struct P
{
	double x,y;
};
vector<pair<double,P>> vt;
double cross(P p1,P p2,P p3) {
    return (p2.x-p1.x)*(p3.y-p1.y)-(p3.x-p1.x)*(p2.y-p1.y);
}
bool judge(ll x,ll tar)
{
	return vt[x].first<=tar;
}
void solve()
{
	cin>>n>>m;
	for(int i=1;i<=n;++i) cin>>b[i];
	sort(b+1,b+1+n,greater<ll>());
	for(int i=1;i<=m;++i) cin>>c[i];
	
	for(int i=1;i<=n;++i)
	{
		P p={(double)i,(double)i*b[i]};//y=ix+i*b[i]
        //上文说过直线可以用jx+j*b[j],所以我们就直接用两个参数来表示对应的直线，可以说p就代表一条直线
		while(vt.size()>=2&&cross(p,vt.rbegin()->second,(next(vt.rbegin()))->second)<0) vt.pop_back();
//cross函数是在计算新加入的直线是否要把里面这一条直线弹出，也就是在计算交点横坐标是否满足条件。这个可以自己推一下
		//弹出无用的直线
		double x=0;
		if(vt.size())
		{
			P pp=vt.back().second;
			x=(pp.y-p.y)/(p.x-pp.x);
		} 
		vt.push_back(make_pair(x,p));
	}	
	ll len=vt.size();
//	for(auto i:vt)
//	{
//		cout<<i.second.x<<" "<<i.second.y<<endl;
//	}
	for(int i=1;i<=m;++i)
	{
		ll l=0,r=len-1;
		while(l<=r)
		{
			ll mid=l+r>>1;
			if(judge(mid,c[i])) l=mid+1;
			else r=mid-1;
		}
		P op=vt[l-1].second;
		cout<<(ll)(op.x*c[i]+op.y)<<" ";
	}
}
int main()
{
	ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
//	ll t;cin>>t;while(t--)
	solve();
	return 0;
}