【计算几何 + LIS】The 2019 Asia Yinchuan First Round Online Programming？- J - Nested Triangles

题目链接https://nanti.jisuanke.com/t/41294

---

题意

给出两个点$P$和$Q$，以及一些其他点$A1$，$A2$…$An$

要求出一个序列$v_1{v}_2...v_n$，使得三角形$△PQA{v}_i包含△PQA{v}_j(i<j)$，且$A{v}_j$不在线段$PA{v}_i$或者线段$QA{v}_i$上。
问序列$v_1{v}_2...v_n$最长是多少，多个最长输出字典序最小的情况。

---

题解

三角形包含$△PQA{v}_i包含△PQA{v}_j$就表明$\angle A{v}_{i}PQ>\angle A{v}_{j}PQ$且$\angle A{v}_{i}QP>\angle A{v}_{j}QP$。所以，如果对每个点都把这两个角度求出来作为两个维度，按照其中一个维度排个序，那就是对另一个维度求LIS。
但是求角度或者求cos都会造成精度的丢失，在这里可以用叉积判断，就是要手写个二分。
题目最后是求LIS反过来的字典序最小，那就反过来求一下，推一遍就行了。

---

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> piir;
const int N=2e5+7;
struct Point{
    ll x,y;
    Point operator - (const Point k)const{
        return (Point){x-k.x,y-k.y};
    }
    void input(){scanf("%lld%lld",&x,&y);}
    void output(){printf("%lld%lld",x,y);}
}p,q,a[N];
inline ll cross(const Point k1,const Point k2){
    return k1.x*k2.y-k2.x*k1.y;
}
int T,n,cs;
struct Node{
    int id;
    Point o;
}st[N],b[N],ans[N];
bool cmp1(Node k1,Node k2){
    if(cross(k1.o-p,k2.o-p)==0) return cross(k1.o-q,k2.o-q)>=0;
    return cross(k1.o-p,k2.o-p)>=0;
}
bool cmp2(Node k1,Node k2){
    if(cross(k1.o-p,k2.o-p)==0) return cross(k1.o-q,k2.o-q)<=0;
    return cross(k1.o-p,k2.o-p)<=0;
}
int tp;
int dp[N];
int tot=0;
vector<int>solve1(){
    tp=0;
    int mx=0;
    for(int i=1;i<=tot;i++){
        int lo=1,hi=tp;dp[i]=tp+1;
        while(lo<=hi){
            int mid=lo+hi>>1;
            if(cross(b[i].o-q,st[mid].o-q)<=0) hi=mid-1,dp[i]=mid;
            else lo=mid+1;
        }
        mx=max(dp[i],mx);
        if(dp[i]>tp) st[++tp]=b[i];
        else st[dp[i]]=b[i];
    }
    for(int i=1;i<=mx+1;i++){
        ans[i].id=n+1;
    }
    int now=mx;
    for(int i=tot;i>=1;i--){
        if(dp[i]==mx||cross(b[i].o-q,ans[dp[i]+1].o-q)<0){
            if(dp[i]>=now){
                if(ans[dp[i]].id>b[i].id) now=dp[i],ans[dp[i]]=b[i];
            }
            else if(dp[i]==now-1){
                ans[dp[i]]=b[i];now--;
            }
        }
    }
    vector<int>res;res.clear();
    for(int i=mx;i>=1;i--){
        res.push_back(ans[i].id);
    }
    return res;
}

vector<int>solve2(){
    tp=0;
    int mx=0;
    for(int i=1;i<=tot;i++){
        int lo=1,hi=tp;dp[i]=tp+1;
        while(lo<=hi){
            int mid=lo+hi>>1;
            if(cross(b[i].o-q,st[mid].o-q)>=0) hi=mid-1,dp[i]=mid;
            else lo=mid+1;
        }
        mx=max(dp[i],mx);
        if(dp[i]>tp) st[++tp]=b[i];
        else st[dp[i]]=b[i];
    }
    for(int i=1;i<=mx+1;i++){
        ans[i].id=n+1;
    }
    int now=mx;
    for(int i=tot;i>=1;i--){
        if(dp[i]==mx||cross(b[i].o-q,ans[dp[i]+1].o-q)>0){
            if(dp[i]>=now){
                if(ans[dp[i]].id>b[i].id) now=dp[i],ans[dp[i]]=b[i];
            }
            else if(dp[i]==now-1){
                ans[dp[i]]=b[i];now--;
            }
        }
    }
    vector<int>res;res.clear();
    for(int i=mx;i>=1;i--){
        res.push_back(ans[i].id);
    }
    return res;
}
int main()
{
    //freopen("in.in","r",stdin);
    scanf("%d",&T);
    for(int t=1;t<=T;t++){
        p.input();q.input();
        scanf("%d",&n);
        for(int i=1;i<=n;i++){
            a[i].input();
        }

        tot=0;
        for(int i=1;i<=n;i++){
            if(cross(a[i]-p,q-p)<=0){
                tot++;
                b[tot].o=a[i];
                b[tot].id=i;
            }
        }
        sort(b+1,b+1+tot,cmp1);
        vector<int>le=solve1();

        tot=0;
        for(int i=1;i<=n;i++){
            if(cross(a[i]-p,q-p)>=0){
                tot++;
                b[tot].o=a[i];
                b[tot].id=i;
            }
        }
        sort(b+1,b+1+tot,cmp2);
        vector<int>ri=solve2();

        if(le.size()<ri.size()) le=ri;
        else if(le.size()==ri.size()){
            for(int i=0;i<le.size();i++){
                if(le[i]==ri[i]) continue;
                if(le[i]<ri[i]) break;
                le=ri;break;
            }
        }
        printf("Case #%d: %d\n",++cs,le.size());
        for(int i=0;i<le.size();i++){
            printf("%d\n",le[i]);
        }
    }
}