bzoj2957 奥妙重重的线段树

https://www.lydsy.com/JudgeOnline/problem.php?id=2957

线段树的query和update竟然还可以结合起来用！

题意：小A的楼房外有一大片施工工地，工地上有N栋待建的楼房。每天，这片工地上的房子拆了又建、建了又拆。他经常无聊地看着窗外发呆，数自己能够看到多少栋房子。
　　为了简化问题，我们考虑这些事件发生在一个二维平面上。小A在平面上(0,0)点的位置，第i栋楼房可以用一条连接(i,0)和(i,Hi)的线段表示，其中Hi为第i栋楼房的高度。如果这栋楼房上任何一个高度大于0的点与(0,0)的连线没有与之前的线段相交，那么这栋楼房就被认为是可见的。
　　施工队的建造总共进行了M天。初始时，所有楼房都还没有开始建造，它们的高度均为0。在第i天，建筑队将会将横坐标为Xi的房屋的高度变为Yi(高度可以比原来大---修建，也可以比原来小---拆除，甚至可以保持不变---建筑队这天什么事也没做)。请你帮小A数数每天在建筑队完工之后，他能看到多少栋楼房？

 

开始考虑离线操作，但是后来发现这是一个三维问题，寻找x1 < x2,y1 < y2,t1 < t2的最长序列，其中x的为位置，y为斜率，t为询问的顺序，但是因为不求他的对数，觉得不能用cdq分治来做，所以考虑线段树。

当一个点发生修改的时候，对左区间是没有影响的，对右区间的影响是左区间的最大值一下的所有点都不计入sum。

考虑维护一个sum和max即可。

说不定我以后可以想到cdq分治的解法。嘿嘿嘿。

#include <map>
#include <set>
#include <ctime>
#include <cmath>
#include <queue>
#include <stack>
#include <vector>
#include <string>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <sstream>
#include <iostream>
#include <algorithm>
#include <functional>
using namespace std;
inline int read(){int now=0;register char c=getchar();for(;!isdigit(c);c=getchar());
for(;isdigit(c);now=now*10+c-'0',c=getchar());return now;}
#define For(i, x, y) for(int i=x;i<=y;i++)  
#define _For(i, x, y) for(int i=x;i>=y;i--)
#define Mem(f, x) memset(f,x,sizeof(f))  
#define Sca(x) scanf("%d", &x)
#define Sca2(x,y) scanf("%d%d",&x,&y)
#define Sca3(x,y,z) scanf("%d%d%d",&x,&y,&z)
#define Scl(x) scanf("%lld",&x);  
#define Pri(x) printf("%d\n", x)
#define Prl(x) printf("%lld\n",x);  
#define CLR(u) for(int i=0;i<=N;i++)u[i].clear();
#define LL long long
#define ULL unsigned long long  
#define mp make_pair
#define PII pair<int,int>
#define PIL pair<int,long long>
#define PLL pair<long long,long long>
#define pb push_back
#define fi first
#define se second 
typedef vector<int> VI;
const double eps = 1e-9;
const int maxn = 1e5 + 10;
const int INF = 0x3f3f3f3f;
const int mod = 1e9 + 7; 
int N,M,K;
struct Tree{int l,r;double mx;int sum;
}tree[maxn << 2];
void Build(int t,int l,int r){tree[t].l = l;tree[t].r = r;if(r == l){tree[t].mx = tree[t].sum = 0;return;}int m = (l + r) >> 1;Build(t << 1,l,m);Build(t << 1 | 1,m + 1,r);
}
int query(int t,double mx){if(tree[t].mx <= mx) return 0;if(tree[t].l == tree[t].r) return 1;int m = (tree[t].l + tree[t].r) >> 1;if(tree[t << 1].mx <=mx) return query(t << 1 | 1,mx);else return query(t << 1,mx) + tree[t].sum - tree[t << 1].sum;
}
void Pushup(int t){tree[t].mx = max(tree[t << 1].mx,tree[t << 1 | 1].mx);tree[t].sum = tree[t << 1].sum + query(t << 1 | 1,tree[t << 1].mx);
}
void update(int t,int x,double k){if(tree[t].l == tree[t].r){tree[t].sum = 1;;tree[t].mx = k;return;}int mid = (tree[t].l + tree[t].r ) >> 1;if(x <= mid) update(t << 1,x,k);else update(t << 1 | 1,x,k);Pushup(t);
}
int main()
{Sca2(N,M);Build(1,1,N);For(i,1,M){int x,y; Sca2(x,y);double t = y * 1.0 / x;update(1,x,t);Pri(tree[1].sum);}#ifdef VSCodesystem("pause");#endifreturn 0;
}