【XSY3920】简单的几何题（几何，凸包）

题面

简单的几何题

题解

易知 $v\ne 0$，那么直接考虑条件：

$$\begin{array}{rl}\left(x_i{u}^2+x_i-uv-y_{i}v\right)b^2& \le \left(x_i{a}^2+x_i-ab-y_{i}b\right)v^2\\ \frac{x_i{u}^2+x_i-uv-y_{i}v}{v^2}& \le \frac{x_i{a}^2+x_i-ab-y_{i}b}{b^2}\end{array}$$

设 $u^{\prime }=\frac{u}{v}$，$v^{\prime }=\frac{1}{v}$，$a^{\prime }=\frac{a}{b}$，$b^{\prime }=\frac{1}{b}$，则原不等式变为：

$$x_i{u}^{\prime 2}+x_i{v}^{\prime 2}-u^{\prime }-y_i{v}^{\prime }\le x_i{a}^{\prime 2}+x_i{b}^{\prime 2}-a^{\prime }-y_i{b}^{\prime }$$

考虑配方：

$$\begin{array}{rl}{x}_i{u}^{\prime 2}+x_i{v}^{\prime 2}-u^{\prime }-y_i{v}^{\prime }& \le x_i{a}^{\prime 2}+x_i{b}^{\prime 2}-a^{\prime }-y_i{b}^{\prime }\\ u^{\prime 2}+v^{\prime 2}-\frac{1}{x_i}{u}^{\prime }-\frac{y_i}{x_i}{v}^{\prime }& \le a^{\prime 2}+b^{\prime 2}-\frac{1}{x_i}{a}^{\prime }-\frac{y_i}{x_i}{b}^{\prime }\\ \left(u^{\prime 2}-\frac{1}{x_i}{u}^{\prime }+\frac{1}{4x_i^2}\right)+\left(v^{\prime 2}-\frac{y_i}{x_i}{v}^{\prime }+\frac{y_i^2}{4x_i^2}\right)& \le \left(a^{\prime 2}-\frac{1}{x_i}{a}^{\prime }+\frac{1}{4x_i^2}\right)+\left(b^{\prime 2}-\frac{y_i}{x_i}{b}^{\prime }+\frac{y_i^2}{4x_i^2}\right)\\ {\left(u^{\prime }-\frac{1}{2x_i}\right)}^2+{\left(v^{\prime }-\frac{y_i}{2x_i}\right)}^2& \le {\left(a^{\prime }-\frac{1}{2x_i}\right)}^2+{\left(b^{\prime }-\frac{y_i}{2x_i}\right)}^2\end{array}$$

其中左右分别可以看成是 $(u^{\prime },v^{\prime })$ 和 $(a^{\prime },b^{\prime })$ 到同一个点 $\left(\frac{1}{2x_i},\frac{y_i}{2x_i}\right)$ 的距离，其中 $(a^{\prime },b^{\prime })$ 到 $\left(\frac{1}{2x_i},\frac{y_i}{2x_i}\right)$ 的距离是已知的，记为 $R_i$。

那么这个不等式就是要求 $(u^{\prime },v^{\prime })$ 这个点在以 $\left(\frac{1}{2x_i},\frac{y_i}{2x_i}\right)$ 为圆心、$R_i$ 为半径的圆内或圆上，记这个圆为圆 $O_i$（$O_i=\left(\frac{1}{2x_i},\frac{y_i}{2x_i}\right)$），那么 $\odot O_i$ 过点 $(a^{\prime },b^{\prime })$。

那么如果 $i$ 是答案的一种，当且仅当存在点 $(u^{\prime },v^{\prime })$ 满足它被 $\odot O_i$ 包含（“包含” 指含边界），且 $\forall j\ne i$ 有 $(u^{\prime },v^{\prime })$ 在 $\odot O_j$ 外。也就是存在点 $(u^{\prime },v^{\prime })$ 恰好只被圆 $i$ 包含。

如何找到这些 $i$ 呢？

我们考虑求出点集 { O , O 1 , O 2 , ⋯   , O n } \{O,O_1,O_2,\cdots,O_n\} {O,O1​,O2​,⋯,On​} 的凸包，那么 $i$ 是答案的一种当且仅当 $O_i$ 在凸包上。

感性证明：（以下证明都十分感性，真正的证明应该是用反演变换化曲为直来证）

分情况讨论：

• 若 $O$ 不在凸包上。

  那么考虑一个凸包内（不含边界）的点 $A$，显然一定可以找到两个在凸包上的相邻点 $B,C$ 使得 $A$ 在 $△OBC$ 内，如图：

  

  我们把这个三角形抽出来单独看，作出 $\odot A,\odot B,\odot C$：

  

  设 $\odot A$ 的过点 $O$ 的直径交 $\odot A$ 于另一点 $P_A$，同理设出 $P_B$ 和 $P_C$。

  设 $\odot B$ 与 $\odot C$ 交于点 $Q$，那么容易得到 $\angle P_{B}QO=\angle P_{C}QO=90°$，那么 $P_B,Q,P_C$ 三点共线。

  容易看出 $△P_{B}O{P}_C$ 是由 $△BOC$ 位似变换得到，所以 $A$ 在 $△BOC$ 内运动可以看作是 $P_A$ 在 $△P_{B}O{P}_C$ 内运动：（颜色可能有点变化，请不要在意）

  

  那么容易看出 $\odot B$ 与 $\odot C$ 的并集一定包含 $\odot A$。

  所以 $\odot A$ 所对应的 $i$ 肯定不符合要求。

• 若 $O$ 在凸包上。

  与第一种情况是类似的，对于凸包内（不含边界）的某一个点 $A$，同样一定可以找到两个在凸包上的相邻点 $B,C$ 使得 $A$ 在 $△OBC$ 内，如图：

  

  接下来的证明过程和第一种情况类似。

感性地证毕。

整道题的大致流程就是这样。

还有一个小细节，你发现如果用 double存点的话精度会爆炸。

然后又发现每一个点都是 $\left(\frac{x}{z},\frac{y}{z}\right)$ 的形式，所以我们用三元组 $(x,y,z)$ 存储一个点即可。

然后叉积判断等用 __int128即可，具体详见代码。

代码如下：（细节很多）

#include<bits/stdc++.h>
 
#define N 100010
#define ll long long
#define lll __int128
 
using namespace std;
 
inline int read()
{
    int x=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-') f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=(x<<1)+(x<<3)+(ch^'0');
        ch=getchar();
    }
    return x*f;
}
 
struct Point
{
    ll x,y,z;
    int id;
}p[N],st[N];
 
inline int check(Point a,Point b,Point c)
{
    lll x=(lll)(a.x*b.z-b.x*a.z)*(lll)(a.y*c.z-a.z*c.y);
    lll y=(lll)(a.x*c.z-c.x*a.z)*(lll)(a.y*b.z-a.z*b.y);
    if(x==y) return 0;
    return x<y?-1:1;
}
 
inline bool inside(Point a,Point b,Point c)
{
    if(a.x*b.z==b.x*a.z||b.x*c.z==c.x*b.z) return 1;
    return !((a.x*b.z<b.x*a.z)^(b.x*c.z<c.x*b.z));
}
 
inline bool operator < (Point a,Point b)
{
    int tmp=check(p[1],a,b);
    if(!tmp)
    {
        if(p[1].x*a.z<=a.x*p[1].z&&p[1].x*b.z<=b.x*p[1].z) return a.x*b.z<b.x*a.z;
        else return a.x*b.z>b.x*a.z;
    }
    return tmp>0;
}
 
inline bool operator == (Point a,Point b)
{
    return a.x*b.z==b.x*a.z&&a.y*b.z==b.y*a.z;
}
 
int n;
int top;
bool vis[N];
bool ban[N];
 
void Graham()
{
    for(int i=2;i<=n;i++)
        if(p[i].y*p[1].z<p[1].y*p[i].z||
        (p[i].y*p[1].z==p[1].y*p[i].z&&p[i].x*p[1].z<p[1].x*p[i].z))
            swap(p[1],p[i]);
    sort(p+2,p+n+1);
    st[++top]=p[1];
    for(register int i=2;i<=n;i++)
    {
        while(top>1)
        {
            int tmp=check(st[top-1],st[top],p[i]);
            if(!tmp&&st[top]==p[i]) ban[st[top].id]=ban[p[i].id]=1;
            if(tmp<0||(!tmp&&inside(st[top-1],st[top],p[i]))) top--;
            else break;
        }
        st[++top]=p[i];
    }
    for(register int i=1;i<=top;i++) vis[st[i].id]=1;
}
 
int main()
{
    n=read()+1;
    p[1].x=read(),p[1].y=1,p[1].z=read();
    for(register int i=2;i<=n;i++)
        p[i].x=1,p[i].z=2ll*read(),p[i].y=read(),p[i].id=i-1;
    Graham();
    for(register int i=1;i<n;i++)
        if(vis[i]&&!ban[i]) printf("%d ",i);
    return 0;
}
/*
1 2 3
4 5
*/
/*
2 2 3
4 5
6 7
*/
/*
2 1 666666666
333333333 1
455943374 506861437
*/
/*
2 236701642 822463349
213494995 807598793
728662831 347575608
*/