可以发现答案就是斜率的严格上升序列长度。
单点修改斜率,考虑我们如何维护这个东西。
可以线段树来做。维护区间斜率最大值和只考虑这个区间的答案。
那么现在的问题在于如何合并两个区间。
首先左区间的答案肯定是都可以取的。如果右区间的最大值<=左区间最大值则右区间没有可取的,否则考虑右区间的左右子区间。
如果左子区间最大值<=左区间最大值,则左子区间没有可取的,直接递归处理右子区间。
否则递归处理左子区间,右子区间均可取。(指的是考虑了左子区间约束后的答案均可取,即右区间的答案-左子区间的答案)
复杂度
O(mlog2n)
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 100010
inline char gc(){
static char buf[1<<16],*S,*T;
if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
return *S++;
}
inline int read(){
int x=0,f=1;char ch=gc();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
return x*f;
}
int n,m;
struct node{
int ans;double mx;
}tr[N<<2];
inline int find(int p,int l,int r,double val){
if(tr[p].mx<=val) return 0;if(l==r) return 1;int mid=l+r>>1;
if(tr[p<<1].mx<=val) return find(p<<1|1,mid+1,r,val);
return find(p<<1,l,mid,val)+tr[p].ans-tr[p<<1].ans;
}
inline void change(int p,int l,int r,int x,int y){
if(l==r){tr[p].ans=1;tr[p].mx=y*1.0/x;return;}
int mid=l+r>>1;
if(x<=mid) change(p<<1,l,mid,x,y);
else change(p<<1|1,mid+1,r,x,y);
tr[p].mx=max(tr[p<<1].mx,tr[p<<1|1].mx);
tr[p].ans=tr[p<<1].ans+find(p<<1|1,mid+1,r,tr[p<<1].mx);
}
int main(){
// freopen("a.in","r",stdin);
n=read();m=read();
while(m--){
int x=read(),y=read();change(1,1,n,x,y);
printf("%d\n",tr[1].ans);
}return 0;
}