Codeforces 678F Lena and Queries

题意：

你有一个点集，有三种操作：

• 往集合里插入一个点\((x, y)\)
• 从集合中删除第\(i\)次操作插入的点
• 对于给出的\(q\)，询问点集中\(x \cdot q + y\)的最大值

分析：

先不考虑插入删除操作，对于一个给定的点集，如何寻找\(x \cdot q + y\)最大值

这是一个线性规划的问题，只是可行域变成了离散的点。

设\(x \cdot q + y = z\)，其中\(z\)是优化目标。

\(y = -q \cdot x + z\)，使得经过点\((x, y)\)斜率为\(-q\)的直线的截距最大。

那么作为最优解的点一定在点集的凸包上，所以可以用单调栈求出凸包，然后三分求最大值。

因为每次操作都可能导致点集发生变化，不可能每次都求一遍凸包。

每个节点对应一个生存期\([L, R]\)，即在第\(L\)次操作到第\(R\)次操作中点集中有该点。

然后把它插入到一个线段树的区间中，这样线段树的每个点都对应一段操作区间的一个点集。

用这个点集的凸包更新所有这个区间的询问。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
using namespace std;
typedef long long LL;
const int maxn = 300000 + 10;
const LL INF = 1LL << 61;
struct Point
{
	LL x, y;
	Point(LL x = 0, LL y = 0): x(x), y(y) {}
	void read() { scanf("%lld%lld", &x, &y); }
	bool operator < (const Point& t) const {
		return x < t.x || (x == t.x && y < t.y);
	}
	Point operator + (const Point& t) const {
		return Point(x + t.x, y + t.y);
	}
	Point operator - (const Point& t) const {
		return Point(x - t.x, y - t.y);
	}
};
LL Cross(const Point& A, const Point& B) {
	return A.x * B.y - A.y * B.x;
}
LL Dot(const Point& A, const Point& B) {
	return A.x * B.x + A.y * B.y;
}
int type[maxn], top;
Point p[maxn], S[maxn];
vector<Point> v[maxn * 4];
bool del[maxn], empty[maxn];
LL ans[maxn];
void insert(int o, int L, int R, int qL, int qR, int v) {
	if(qL <= L && R <= qR) {
		::v[o].push_back(p[v]);
		return;
	}
	int M = (L + R) / 2;
	if(qL <= M) insert(o<<1, L, M, qL, qR, v);
	if(qR > M) insert(o<<1|1, M+1, R, qL, qR, v);
}
void query(int x) {
	int L = 1, R = top;
	while(R - L >= 3) {
		int mid1 = (L * 2 + R) / 3;
		int mid2 = (L + R * 2) / 3;
		if(Dot(p[x], S[mid1]) < Dot(p[x], S[mid2])) L = mid1;
		else R = mid2;
	}
	for(int i = L; i <= R; i++)
		ans[x] = max(ans[x], Dot(p[x], S[i]));
}
void solve(int o, int L, int R) {
	if(L < R) {
		int M = (L + R) / 2;
		solve(o<<1, L, M);
		solve(o<<1|1, M+1, R);
	}
	sort(v[o].begin(), v[o].end());
	top = 0;
	for(int i = 0; i < v[o].size(); i++) {
		while(top > 1 && Cross(S[top]-S[top-1], v[o][i]-S[top]) >= 0) top--;
		S[++top] = v[o][i];
	}
	for(int i = L; i <= R; i++) if(type[i] == 3 && !empty[i])
		query(i);
}
int main()
{
	int n; scanf("%d", &n);
	int cnt = 0;
	for(int i = 1; i <= n; i++) {
		scanf("%d", type + i);
		if(type[i] == 1) {
			p[i].read();
			cnt++;
		} else if(type[i] == 2) {
			int x; scanf("%d", &x);
			del[x] = true;
			cnt--;
			insert(1, 1, n, x, i, x);
		} else {
			scanf("%lld", &p[i].x);
			p[i].y = 1LL;
			if(!cnt) empty[i] = true;
		}
	}
	for(int i = 1; i <= n; i++)
		if(type[i] == 1 && !del[i])
			insert(1, 1, n, i, n, i);
	for(int i = 1; i <= n; i++) ans[i] = -INF;
	solve(1, 1, n);
	for(int i = 1; i <= n; i++) if(type[i] == 3) {
		if(empty[i]) printf("EMPTY SET\n");
		else printf("%lld\n", ans[i]);
	}
	return 0;
}