计算几何

1.4.5 计算几何

Description

uncle-lu对计算几何有着浓厚的兴趣。他经常对着平面直角坐标系发呆,思考一些有趣的问 题。今天,他想到了一个十分有意思的题目:

首先,uncle-lu会在x轴正半轴和y轴正半轴分别挑选nn个点。随后,他将x轴的点与y轴 的点一一连接,形成nn条线段,并保证任意两条线段不相交。uncle-lu确定这种连接方式有且仅有一种。最后,uncle-lu会给出mm个询问。对于每个询问,将会给定一个点P (x_p, y_p)P(x
p
​
,y
p
​
),问线段OPOP(O为坐标原点)与nn条线段会产生多少个交点?

Input
第1行包含一个正整数nn,表示线段的数量;

第2行包含nn个正整数,表示uncle-lu在x轴选取的点的横坐标;

第3行包含nn个正整数,表示uncle-lu在y轴选取的点的纵坐标;

第4行包含一个正整数mm,表示询问数量;

随后m行,每行包含两个正整数x_px
p
​
和y_py
p
​
,表示询问中给定的点的横、纵坐标。

Output
共mm行,每行包含一个非负整数,表示你对这条询问给出的答案。

Sample Input 1

3
4 5 3
3 5 4
2
1 1
3 3
Sample Output 1

0
3
Hint

3条线段分别为:(3, 0)−(0, 3)、(4, 0)−(0, 4)、(5, 0)−(0, 5)

(0, 0)−(1, 1)不与他们有交点,答案为0

(0, 0)−(3, 3)与三条线段均有交点,答案为3。

数据范围：

对于100%100的数据:n, m \leq 10^5,1 \leq x, y \le 231。n,m≤10
5
,1≤x,y≤231。

截距式

#include <bits/stdc++.h>
using namespace std;
int n,m;
int x[11000],y[11000];
int X,Y;
long long int cot(int t)
{
    return (long long)X*y[t]+(long long)Y*x[t]-(long long)x[t]*y[t];
}
int judge()
{
    int l=0,r=n+1,ans=0,mid;
    long long t;
    while(l+1<r)
    {
        mid=l+(r-l)/2;
        t=cot(mid);
        if(t>=0)
        {
            ans=mid;
            l=mid;
        }
        else
        {
            r=mid;
        }
    }
    return ans;
}
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        cin>>x[i];
    }
    for(int i=1;i<=n;i++)
    {
        cin>>y[i];
    }
    sort(x+1,x+n+1);
    sort(y+1,y+n+1);
    scanf("%d",&m);
    while(m--)
    {
        cin>>X>>Y;
        cout<<judge()<<endl;
    }
    return 0;
}