AtCoder Grand Contest 015 E Mr.Aoki Incubator dp+树状数组

Description

---

给出若干个数轴上的点，每个点有起始位置x和速度v
一开始可以给若干点染色，染色点与未染色点相遇会传染颜色，问多少种染色方案使得最后所有点被染色
$$\begin{gathered} n\le 10^5 \\ \forall i\in [1,n],|x_i|,|v_i|\le 10^9 \end{gathered}$$

Solution

---

早上的模型套得很mio啊

考虑两个点会相遇的情况
$x_i<x_j,v_i<v_j$这样i永远追不上j
$x_i<x_j,v_i>v_j$这样i可以追上j
$x_i>x_j,v_i<v_j$这样j可以追上i
$x_i>x_j,v_i>v_j$这样j永远追不上i

我们可以把追击的模型转化成若干一次函数，在x=无线远处记一下投射的n个点，通过一些意识流操作可以发现某一条线能染的点一定是连续一段的，那么某一个点能染色的点也是按照v升序排序后连续的一段

于是问题变成了若干区间，问覆盖整个数轴的方案。设f[i]表示前i条线段，第i条一定选的方案，用树状数组优化转移就可以了

Code

---

#include <stdio.h>
#include <string.h>
#include <algorithm>
#define rep(i,st,ed) for (int i=st;i<=ed;++i)
#define drp(i,st,ed) for (int i=st;i>=ed;--i)
#define lowbit(x) (x&-x)
#define fi first
#define se second

typedef std:: pair <int,int> pair;
const int MOD=1000000007;
const int N=500005;

struct data {int x,v,id;} r[N];

int s[N<<1],f[N],n;

pair p[N];

int read() {
	int x=0,v=1; char ch=getchar();
	for (;ch<'0'||ch>'9';v=(ch=='-')?(-1):(v),ch=getchar());
	for (;ch<='9'&&ch>='0';x=x*10+ch-'0',ch=getchar());
	return x*v;
}

void upd(int &x,int v) {
	x+=v; (x>=MOD)?(x-=MOD):0;
}

void add(int x,int v) {
	for (;x<=n;x+=lowbit(x)) upd(s[x],v);
}

int get(int x) {
	int res=0;
	for (;x;x-=lowbit(x)) upd(res,s[x]);
	return res;
}

int ask(int l,int r) {
	int tmp=get(r)-get(l-1)+MOD;
	(tmp>=MOD)?(tmp-=MOD):0;
	return tmp;
}

bool cmp1(data a,data b) {
	return (a.x!=b.x)?(a.x<b.x):(a.v>b.v);
}

bool cmp2(data a,data b) {
	return a.v<b.v;
}

bool cmp3(pair a,pair b) {
	return (a.se!=b.se)?(a.se<b.se):(a.fi<b.fi);
}

int main(void) {
	//freopen("river.in","r",stdin);
	//freopen("river.out","w",stdout);
	n=read();
	rep(i,1,n) r[i].x=read(),r[i].v=read();
	std:: sort(r+1,r+n+1,cmp2);
	rep(i,1,n) r[i].id=i;
	std:: sort(r+1,r+n+1,cmp1);
	int mn=n+1,mx=0;
	rep(i,1,n) {
		mx=std:: max(mx,r[i].id);
		p[r[i].id].se=mx;
	}
	drp(i,n,1) {
		mn=std:: min(mn,r[i].id);
		p[r[i].id].fi=mn;
	}
	std:: sort(p+1,p+n+1,cmp3);
	rep(i,1,n) {
		if (p[i].fi==1) upd(f[i],1);
		int tmp=ask(std:: max(1,p[i].fi-1),p[i].se);
		upd(f[i],tmp);
		add(p[i].se,f[i]);
	}
	printf("%d\n", ask(n,n));
	return 0;
}