Gym 101606L-计算几何+dp

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
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;
}