BZOJ 2762: [JLOI2011]不等式组( 平衡树 )

对不等式变形..然后就是维护一些数, 随便找个数据结构都能写吧....用double感觉会有精度误差, 分类讨论把<改成<=了很久后弃疗了, 自己写了个分数体....然后速度就被完爆了..

------------------------------------------------------------------------------

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
 
using namespace std;
 
typedef long long ll;
 
const int maxn = 100009;
 
int N, ans, n;
char s[10];
 
struct F {
int a, b;
F(int _a = 0, int _b = 0) : a(_a), b(_b) {
if(b < 0)
a = -a, b = -b;
}
bool operator < (const F &t) const {
ll e = ll(a) * b, f = ll(t.a) * t.b;
if(!e) return f > 0;
if(!f) return e < 0;
if(e > 0 && f < 0) return false;
if(e < 0 && f > 0) return true;
return ll(a) * t.b < ll(b) * t.a;
}
bool operator == (const F &t) const {
return ll(a) * t.b == ll(b) * t.a;
}
};
 
struct O {
int t;
F v;
O(int _t = 0, F _v = F()) : t(_t), v(_v) {
}
} o[maxn];
 
struct Node {
Node* ch[2];
int sz, r;
F f;
inline void upd() {
sz = ch[0]->sz + ch[1]->sz + 1;
}
} pool[maxn], *pt, *Null, *Rt[2];
 
void Init_Treap() {
pt = pool;
Null = pt++;
Null->ch[0] = Null->ch[1] = Null;
Null->sz = 0;
Rt[0] = Rt[1] = Null;
}
 
Node* newNode(F v) {
pt->f = v;
pt->sz = 1;
pt->r = rand();
pt->ch[0] = pt->ch[1] = Null;
return pt++;
}
 
void Rotate(Node*&t, int d) {
Node* o = t->ch[d ^ 1];
t->ch[d ^ 1] = o->ch[d];
o->ch[d] = t;
t->upd(), o->upd();
t = o;
}
 
void Insert(Node*&t, F v) {
if(t == Null) {
t = newNode(v);
} else {
int d = (t->f < v);
Insert(t->ch[d], v);
if(t->ch[d]->r > t->r)
Rotate(t, d ^ 1);
}
t->upd();
}
 
void Delete(Node*&t, F v) {
int d = (v == t->f ? -1 : (t->f < v));
if(d == -1) {
if(t->ch[0] != Null && t->ch[1] != Null) {
int _d = (t->ch[0]->r > t->ch[1]->r);
Rotate(t, _d);
Delete(t->ch[_d], v);
} else {
t = (t->ch[0] != Null ? t->ch[0] : t->ch[1]);
}
} else {
Delete(t->ch[d], v);
}
if(t != Null) t->upd();
}
 
int Query(F v) {
int ret = 0;
for(Node* t = Rt[0]; t != Null; ) if(t->f < v) {
ret += t->ch[0]->sz + 1;
t = t->ch[1];
} else
t = t->ch[0];
for(Node* t = Rt[1]; t != Null; ) if(v < t->f) {
ret += t->ch[1]->sz + 1;
t = t->ch[0];
} else 
t = t->ch[1];
return ret;
}
 
int main() {
scanf("%d", &N);
ans = n = 0;
Init_Treap();
while(N--) {
scanf("%s", s);
if(s[0] == 'A') {
int a, b, c;
scanf("%d%d%d", &a, &b, &c);
c -= b;
if(!a) {
ans += (c < 0);
o[n++] = O(0, F(0, c < 0));
} else {
Insert(Rt[a < 0], F(c, a));
o[n++] = O(a < 0 ? -1 : 1, F(c, a));
}
} else if(s[0] == 'D') {
int t;
scanf("%d", &t), t--;
if(o[t].t < -1)
continue;
if(!o[t].t) {
ans -= o[t].v.b;
} else {
o[t].t > 0 ? Delete(Rt[0], o[t].v) : Delete(Rt[1], o[t].v);
}
o[t].t = -2;
} else {
int t;
scanf("%d", &t);
printf("%d\n", ans + Query(F(t ,1)));
}
}
return 0;
}

------------------------------------------------------------------------------

2762: [JLOI2011]不等式组

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 714  Solved: 205
[Submit][Status][Discuss]

Description

旺汪与旺喵最近在做一些不等式的练习。这些不等式都是形如ax+b>c 的一元不等式。当然,解这些不等式对旺汪来说太简单了,所以旺喵想挑战旺汪。旺喵给出一组一元不等式,并给出一个数值 。旺汪需要回答的是x=k 时成立的不等式的数量。聪明的旺汪每次都很快就给出了答案。你的任务是快速的验证旺汪的答案是不是正确的。
 

Input

输入第一行为一个正整数 ,代表接下来有N 行。
接下来每一行可能有3种形式:
1.“Add a b c”,表明要往不等式组添加一条不等式ax+b>c ;
2.“Del i”,代表删除第i 条添加的不等式(最先添加的是第1条)。
3.“Query k”,代表一个询问,即当x=k 时,在当前不等式组内成立的不等式的数量。
注意一开始不等式组为空,a,b,c,i,k 均为整数,且保证所有操作均合法,不会出现要求删除尚未添加的不等式的情况。
 

Output

 
对于每一个询问“Query k”,输出一行,为一个整数,代表询问的答案。

Sample Input

9
Add 1 1 1
Add -2 4 3
Query 0
Del 1
Query 0
Del 2
Query 0
Add 8 9 100
Query 10

Sample Output

1
1
0
0

HINT

1条添加到不等式组的不等式为x+1>1 ,第2条为-2x+4>3 ,所以第1个询问的时候只有第2条不等式可以成立,故输出1


然后删除第1条不等式,再询问的时候依然是只有第2条不等式可以成立,故输出1


再删除第2条不等式后,因为不等式组里面没有不等式了,所以没有不等式可以被满足,故输出0


继续加入第3条不等式8x+9>100 ,当x=k=10时有8*10+9=89<100,故也没有不等式可以被满足,依然输出0


数据范围:


20%的数据, N<=1000


40%的数据, N<=10000


100%的数据,N<=100000


a,b,c的范围为[-10^8,10^8],k的范围为[-10^6,10^6]


Source

 

posted @ 2016-01-15 21:06  JSZX11556  阅读(252)  评论(0编辑  收藏