首先把相同的事件点合并,那么每个点有时间$t_i$,位置$x_i$,价值$v_i$。
考虑DP,设$f_i$表示按时间从小到大考虑每个事件,目前位于事件$i$的时间与位置时,最多能让多少个事件发生。在$t_i$秒初对方不能碰到$i$,但在这一秒结束结算时,对方可以碰到$i$。
则$f_i=\max(f_j)+v_i$。
其中,$j$需要满足一系列条件:
$1.j$要能在规定时间内到达$i$,也就是$t_i-t_j\geq|x_i-x_j|$。
将绝对值拆掉并移项,得:
$t_i-x_i\geq t_j-x_j(1)$
$t_i+x_i\geq t_j+x_j(2)$
$2.$在$t_i$秒初对方不能碰到$i$:
设$s_i$表示前$i$秒存在的事件数,那么在第$t_i$秒初,也就是第$t_i-1$秒末,对方一共前进了$s_{t_i-1}-f_j$步,所以$s_{t_i-1}-f_j<x_i$,即$f_j>s_{t_i-1}-x_i(3)$。
$3.j$在移动到$i$的过程中不能被对方追上:
最坏情况下,一定是$j$在$t_j$时刻位于$x_j$,然后一直向右移动,那么在第$p(t_j\leq p<t_i)$秒末,对方位于$s_p-f_j$,$j$位于$x_j-t_j+p$,需要满足:
$s_p-f_j<x_j-t_j+p$
移项得:
$s_p-p<x_j-t_j+f_j$
即$\max(s_p-p)<x_j-t_j+f_j$
对$s_p-p$建线段树维护区间最大值,然后在线段树上二分即可得到最大的$o$,满足$[t_j,o)$都不会被追上,设$lim_j=o$,则$lim_j\geq t_i(4)$。
对于限制$(3)$,只需要求出最大的$f_j$,然后判断是否可行即可。
对于限制$(1)$,可以通过排序解决。
对于限制$(2)$和$(4)$,可以cdq分治后扫描线+树状数组处理。
时间复杂度$O(n\log^2n)$。
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=100010,M=262150,inf=10000000;
int Case,_,n,m,cnt,i,j,ans;
int pt[N],post[N],s[N],w[N],val[N],lim[N],f[N];
int qa[N],qb[N],ca,cb;
int FLAG,bit[N],vis[N];
int mx[M],O,OFFSET;
struct P{int t,x,v;P(){}P(int _t,int _x,int _v){t=_t,x=_x,v=_v;}}a[N],b[N];
inline bool cmpa(const P&a,const P&b){return a.t==b.t?a.x<b.x:a.t<b.t;}
inline bool cmpb(const P&a,const P&b){
if(a.t-a.x!=b.t-b.x)return a.t-a.x<b.t-b.x;
return a.t+a.x<b.t+b.x;
}
inline bool cmpq(int x,int y){
return val[x]<val[y];
}
inline int lower(int x){
int l=1,r=m,mid,t;
while(l<=r)if(pt[mid=(l+r)>>1]<=x)l=(t=mid)+1;else r=mid-1;
return t;
}
inline void up(int&x,int y){x<y?(x=y):0;}
inline void add(int x,int p){for(;x<=m;x+=x&-x)if(vis[x]<FLAG)vis[x]=FLAG,bit[x]=p;else up(bit[x],p);}
inline void ask(int&t,int x){for(;x;x-=x&-x)if(vis[x]==FLAG)up(t,bit[x]);}
void build(int x,int a,int b){
if(a==b){
mx[x]=w[a];
return;
}
int mid=(a+b)>>1;
build(x<<1,a,mid);
build(x<<1|1,mid+1,b);
mx[x]=max(mx[x<<1],mx[x<<1|1]);
}
void getr(int x,int a,int b,int c){
if(c<=a){
if(mx[x]<OFFSET)return;
if(a==b){
O=a;
return;
}
}
int mid=(a+b)>>1;
if(c<=mid){
getr(x<<1,a,mid,c);
if(O)return;
}
getr(x<<1|1,mid+1,b,c);
}
inline void getlim(int x){
int offset=b[x].x+f[x]-b[x].t;
lim[x]=0;
int t=post[x];
if(w[t]>=offset)return;
O=0;
OFFSET=offset;
getr(1,1,m,t);
if(!O)O=m;
lim[x]=O;
}
void CDQ(int l,int r){
if(l==r){
if(b[l].t==0)up(f[l],0);
if(f[l]<=s[post[l]-1]-b[l].x)f[l]=-inf;
if(f[l]>=0){
f[l]+=b[l].v;
up(ans,f[l]);
getlim(l);
}
return;
}
int mid=(l+r)>>1;
CDQ(l,mid);
int i,j;
ca=cb=0;
for(i=l;i<=mid;i++)if(f[i]>=0&&lim[i])qa[++ca]=i;
for(i=r;i>mid;i--)qb[++cb]=i;
sort(qa+1,qa+ca+1,cmpq);
sort(qb+1,qb+cb+1,cmpq);
FLAG++;
for(i=j=1;i<=cb;i++){
while(j<=ca&&val[qa[j]]<=val[qb[i]]){
add(m-lim[qa[j]]+1,f[qa[j]]);
j++;
}
ask(f[qb[i]],m-post[qb[i]]+1);
}
CDQ(mid+1,r);
}
inline void solve(){
sort(a+1,a+m+1,cmpa);
for(cnt=0,i=1;i<=m;i=j){
for(j=i;j<=m&&a[i].t==a[j].t&&a[i].x==a[j].x;j++);
b[++cnt]=P(a[i].t,a[i].x,j-i);
if(b[cnt].t==0)b[cnt].v=0;
}
sort(b+1,b+cnt+1,cmpb);
for(i=1;i<=cnt;i++)pt[i]=b[i].t;
sort(pt+1,pt+cnt+1);
for(m=0,i=1;i<=cnt;i++)if(i==1||pt[i]>pt[i-1])pt[++m]=pt[i];
for(i=1;i<=cnt;i++)post[i]=lower(b[i].t);
for(i=1;i<=m;i++)s[i]=0;
for(i=1;i<=cnt;i++)s[post[i]]+=b[i].v;
for(i=1;i<=m;i++)s[i]+=s[i-1];
for(i=1;i<=m;i++)w[i]=s[i]-pt[i];
for(i=1;i<=cnt;i++)f[i]=-inf,val[i]=b[i].t+b[i].x;
build(1,1,m);
CDQ(1,cnt);
}
int main(){
scanf("%d",&Case);
for(_=1;_<=Case;_++){
scanf("%d%d",&n,&m);
for(i=1;i<=m;i++)scanf("%d%d",&a[i].t,&a[i].x);
a[++m]=P(0,1,0);
ans=0;
if(m)solve();
printf("Case #%d: %d\n",_,ans);
}
}