bzoj1071 [SCOI2007]组队

Description

NBA每年都有球员选秀环节。通常用速度和身高两项数据来衡量一个篮球运动员的基本素质。假如一支球队里速度最慢的球员速度为$minV$，身高最矮的球员高度为$minH$，那么这支球队的所有队员都应该满足: $A \cdot ( height – minH ) + B \cdot ( speed – minV ) \leqslant C$ 其中$A$和$B$，$C$为给定的经验值。这个式子很容易理解，如果一个球队的球员速度和身高差距太大，会造成配合的不协调。 请问作为球队管理层的你，在$N$名选秀球员中，最多能有多少名符合条件的候选球员。

Input

第一行四个数$N$、$A$、$B$、$C$ 下接$N$行每行两个数描述一个球员的$height$和$speed$

Output

最多候选球员数目。

Sample Input

4 1 2 10
5 1
3 2
2 3
2 1

Sample Output

4

HINT

数据范围： $N \leqslant 5000$ ,$heigh$t和$speed$不大于$10000$。$A$、$B$、$C$在长整型以内。

Solution

尺取好题。首先不等式原式很烦人，我们要试着转一下。移项得$A \cdot height + B \cdot speed \leqslant A \cdot minH + B \cdot minV + C$。
设$s = A \cdot height + B \cdot speed$。
暴力容易做，枚举$minV$和$minH$，看有多少个球员能加入。
考虑先枚举一维，枚举$minV$，忽略所有$speed$小于$minV$的球员。
下面枚举$minH$。如果将球员按$height$排序，$A \cdot height + B \cdot speed$并不单调，所以无法直接尺取。然而$A \cdot minH + B \cdot minV + C$是单调的，这个性质就是解题关键。我们维护两个序列，第一个序列按照$s$从小到大排序，第二个序列按照 $height$ 从小到大排序，我们从小到大枚举$minH$。然后就可以尺取了。
尺取维护两个指针，$r$表示第一个序列中$[1, r)$满足不等式的限制，$l$表示第二个序列中$[1, l)$满足$minH$的限制。
大概内容已经说完了，剩下的就看代码吧。

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

inline int read() {
    int x = 0, flag = 1; char ch = getchar();
    while (ch > '9' || ch < '0') { if (ch == '-') flag = -1; ch = getchar(); }
    while (ch <= '9' && ch >= '0') { x = x * 10 + ch - '0'; ch = getchar(); }
    return x * flag;
}
inline void write(int x) { if (x >= 10) write(x / 10); putchar(x % 10 + '0'); }

#define N 5001
#define rep(i, a, b) for (int i = a; i <= b; i++)
#define ll long long

int n;
ll A, B, C;
struct manType { int h, v; ll s; void get() { s = A * h + B * v; } }x[3][N];
bool cmpH(const manType a, const manType b) { return a.h < b.h; }
bool cmpS(const manType a, const manType b) { return a.s < b.s; }

int main() {
    n = read(), A = read(), B = read(), C = read();
    rep(i, 1, n) x[0][i].h = read(), x[0][i].v = read(), x[0][i].get(), x[1][i] = x[0][i];
    sort(x[0] + 1, x[0] + 1 + n, cmpS), sort(x[1] + 1, x[1] + 1 + n, cmpH);
    int ans = 0;
    rep(i, 1, n) {
        int l = 1, r = 1, cnt = 0;
        rep(j, 1, n) {
            int minV = x[0][i].v, minH = x[1][j].h;
            ll val = minH * A + minV * B + C;
            while (r <= n && x[0][r].s <= val)
                cnt += x[0][r].v >= minV && x[0][r].h >= x[1][l].h, r++;
                //如果x[0][r] < x[1][l],那么它在上一轮就会被删除
            while (l <= n && x[1][l].h < minH)
                cnt -= x[1][l].s <= val && x[1][l].v >= minV, l++;
            ans = max(ans, cnt);
        }
    }
    write(ans);
    return 0;
}