[CF 249D]Donkey and Start解题报告

最新推荐文章于 2024-05-29 20:34:53 发布

原创最新推荐文章于 2024-05-29 20:34:53 发布 · 885 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#codeforces #解题报告

本文介绍了一种利用最长上升子序列（LIS）算法解决特定几何问题的方法。通过坐标变换，将原问题转化为求解二维平面上点的最长严格递增子序列问题。文章详细解释了如何使用树状数组来高效实现这一算法。

题意

第一象限上有n<=10^5个点。以分数的形式给定两个斜率k1和k2（这样可以表示tan90°的情况）。保证对应倾斜角在[0°,90°]内。从原点开始，从我们所在的点作两条斜率分别为k1和k2的射线，并在射线所夹区域内部（不含射线上）选一个点，并继续这样做下去，直到无点可选。问途中最多经过多少个点。

像这样：

它的坐标是x轴往左，y轴往上……我估计这么做的原因是为了迷惑你……

分析

这道题的关键是很明显的：假设过原点，斜率为k1和k2的单位向量分别为v1和v2，那我们就以(v1,v2)为基底建立一组新的坐标系。而在新的坐标系中，规则就变成了：下一步只能选择x和y坐标均严格大于当前点的点。

换言之，将所有点按x坐标为第一关键字，y坐标为第二关键字排序后，我们要求的也就是从原点出发的最长上升序列。当然需要注意，该序列中的x坐标也必须单调上升。

所以按照求LIS的经典方法，用一个树状数组维护即可。

代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;
typedef long long LL;
const int SIZEN=100010;
class BIT{
public:
	int n;
	int s[SIZEN];
	#define lowbit(x) ((x)&(-x))
	void clear(int _n){
		n=_n;
		memset(s,0,sizeof(s));
	}
	void modify(int x,int t){
		for(;x<=n;x+=lowbit(x)) s[x]=max(s[x],t);
	}
	int pref_max(int x){
		int ans=0;
		for(;x;x-=lowbit(x)) ans=max(ans,s[x]);
		return ans;
	}
};
class Point{
public:
	LL x,y;
	int id;
};
void print(const Point &p){
	cout<<"("<<p.x<<","<<p.y<<" "<<p.id<<")";
}
bool operator < (const Point &a,const Point &b){
	if(a.x==b.x) return a.y<b.y;
	return a.x<b.x;
}
LL operator * (const Point &a,const Point &b){
	return a.x*b.x+a.y*b.y;
}
int N;
Point P[SIZEN];
Point alpha1,alpha2;
LL Y[SIZEN]={0};int tot=0;
void coor_trans(void){//坐标变换
	for(int i=0;i<=N;i++){
		LL a=P[i].x*alpha2.y-alpha2.x*P[i].y;
		LL b=alpha1.x*P[i].y-P[i].x*alpha1.y;
		a*=-1,b*=-1;//强行改规则,只选x,y坐标都小于它的
		P[i].x=a,P[i].y=b;
		Y[tot++]=b;
	}
	sort(Y,Y+tot);
	tot=unique(Y,Y+tot)-Y;
	for(int i=0;i<=N;i++) P[i].y=lower_bound(Y,Y+tot,P[i].y)-Y+1;
}
int dp[SIZEN]={0};
BIT T;
void work(void){
	sort(P,P+1+N);
	T.clear(N+1);
	for(int i=0,j=0;i<=N;i++){
		while(P[j].x<P[i].x){
			T.modify(P[j].y,dp[P[j].id]);
			j++;
		}
		dp[P[i].id]=T.pref_max(P[i].y-1)+1;
	}
	printf("%d\n",dp[0]-1);
}
void read(void){
	P[0].x=P[0].y=0;P[0].id=0;
	scanf("%d",&N);
	scanf("%I64d/%I64d",&alpha1.y,&alpha1.x);
	scanf("%I64d/%I64d",&alpha2.y,&alpha2.x);
	for(int i=1;i<=N;i++){
		scanf("%I64d%I64d",&P[i].x,&P[i].y);
		P[i].id=i;
	}
}
int main(){
	//freopen("t1.in","r",stdin);
	read();
	coor_trans();
	work();
	return 0;
}