计算几何（广州大学第十八届ACM大学生程序设计竞赛）

题目描述

2023年赛季中，污渍与小夨相约，区域赛上一定要先看SUA的计算几何题，并且成功偷鸡；遗憾的是，赛季结束后，两人只能举起可乐向着一轮残月：“****，退钱！！！”；为了弥补遗憾，小夨决定出一道简单的计算几何题，并且期待赛场上的朋友们能够将其通过。

以上为题目背景；

给定 n 个点（编号 1∼n），你可以进行若干次操作：选取任意两个不同的点连线；

保证线段两两不相交的前提下，至多可以有多少条线段？同时给出具体的连线方案。

输入描述:

第一行输入一个整数 n (1≤n≤100)，代表一共有 nnn 个点；

接下来输入 n 行，每行输入两个整数 xi,yi (0≤xi,yi≤106)，表示点 i 的坐标为 (xi,yi)。 

输出描述:

第一行输出一个整数 tot，代表至多可以有 tot条线段；

接下来输出 tot行；

每行输出两个正整数 id1,id2 分别表示每条线段的两个顶点对应的编号。

示例1

输入

3
0 0
1 1
3 0

输出

1
1 2

备注:

保证数据不存在相同坐标的点。

思路

只要保证所有还没找过的点都在这条线段的同一侧，这样就可以保证这条线段，和其他还没找过的点能组成的线段，一定不相交。

然后把这两个点删了，从接下来的集合里面去找，不断迭代就好了。

AC代码：

#include <iostream>
 
using namespace std;
 
using pii = pair<int, int>;
using ll = long long;
 
const int N = 110, M = 10010;
 
int x[N], y[N], n, m, used[N];
 
pii edges[M];
 
bool judge(int p1, int p2)
{
    int less = 0, greater = 0;
    int xa = x[p1], ya = y[p1];
    int xb = x[p2], yb = y[p2];
    int d1 = 0, d2 = 0;
    for (int i = 0; i < n; ++i)
    {
        if (i == p1 || i == p2)
            continue;
        if (used[i])
            continue;
        int xc = x[i], yc = y[i];
        ll left = (ll)(yc - ya) * (xb - xa);
        ll right = (ll)(yb - ya) * (xc - xa);
        if (left == right) {
            int xmin = min(xa, xb), xmax = max(xa, xb);
            int ymin = min(ya, yb), ymax = max(ya, yb);
            if (xc >= xmin && xc <= xmax && yc >= ymin && yc <= ymax)
                return false;
            ll dd1 = (ll)(xc - xa) * (xc - xa) + (ll)(yc - ya) * (yc - ya);
            ll dd2 = (ll)(xc - xb) * (xc - xb) + (ll)(yc - yb) * (yc - yb);
            if (dd1 < dd2) d1 = 1;
            if (dd1 > dd2) d2 = 1;
        }
        if (d1 == 1 && d2 == 1) return false;
        if (left < right)
            ++less;
        if (left > right)
            ++greater;
    }
//     printf("%d-%d, less: %d, greater: %d\n", p1, p2, less, greater);
    if (less == 0 || greater == 0) return true;
    return false;
 
}
 
int main()
{
    scanf("%d", &n);
    for (int i = 0; i < n; ++i) {
        scanf("%d%d", x + i, y + i);
    }
    while (true)
    {
        bool flag = false;
        for (int i = 0; i < n; ++i)
        {
            if (used[i])
                continue;
            for (int j = i + 1; j < n; ++j)
            {
                if (used[j])
                    continue;
                if (judge(i, j))
                {
                    flag = true;
                    used[i] = used[j] = true;
                    edges[m++] = {i, j};
                    break;
                }
            }
            if (flag) break;
        }
        if (!flag)
            break;
    }
    printf("%d\n", m);
    for (int i = 0; i < m; ++i)
    {
        printf("%d %d\n", edges[i].first + 1, edges[i].second + 1);
    }
}

还有一种方法

n=int(input())
a=[]
for in range(n):
    x,y=map(int,input().split())
    a.append([x,y])
b=a.copy()
b.sort()
print(n//2)
for i in range(1,n,2):
    d=a.index(b[i-1])
    f=a.index(b[i])
    print(d+1,f+1)