Lizard Lounge (计算几何+dp)
感受到了stl+c++17的强大
题目大意:起始点sx,sy还有n个其余的点x,y问以sx,sy为起点的射线中的点的最长上升子序列的个数
正常代码写过去了,但感觉不太美观,就整理了下思路,之后就玩了一些好玩的东西
思路:【斜率,与sx,sy的相对关系相同】同时满足时,在一条射线上,之后dis近的在前面,跑最长上升子序列即可
map的key 维护的是 斜率为二元组pii,与相对关系op
val为set集合,存放pii按照第一关键字dis排序
看起来很舒服
#include<bits/stdc++.h>
using namespace std;
#define ios ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
typedef long long ll;
typedef pair<int,int> pii;
const int maxn = 1e6+5;
int sx,sy,n,ans;
int dp[maxn];
map<tuple<pii,int>,set<pair<ll,int>>>mp;
vector<int>ve;
inline pii grad(int ex,int ey){
int x=ex-sx,y=ey-sy;
int t=__gcd(abs(x),abs(y));
return {y/t,x/t};
}
inline tuple<pii,int>get(int x,int y,int h){
int op;
if(x==sx) op=y<sy;
else op=x>sx;
return {grad(x,y),op};
}
inline ll get_dis(ll x,ll y){
return (x-sx)*(x-sx)+(y-sy)*(y-sy);
}
inline int solve(vector<int>&a){
dp[1]=*a.begin();
int ans=1;
for(int i=1;i<(int)a.size();i++){
int pos=lower_bound(dp+1,dp+ans+1,a[i])-dp;
dp[pos]=a[i];
ans=max(ans,pos);
}
return ans;
}
int main() {
ios;
cin>>sx>>sy>>n;
for(int i=1;i<=n;i++){
int x,y,h;
cin>>x>>y>>h;
mp[get(x,y,h)].insert({get_dis(x,y),h});
}
for(auto &[u,v]:mp){
ve.clear();
for(auto &[x,y]:v){
ve.push_back(y);
}
ans+=solve(ve);
}
cout<<ans<<endl;
}