CF140F 计算几何

题目大意：

给定平面上$n$个不同的点，最多删除$k$个点，在所有删除情况中，删完之后的所有点如果以坐标系上某一个点中心对称，输出所有这样子的中心点。若是有无数个这样的中心点，则输出$-1$

解题思路：

• 先将所有点按照$x$为第一关键字排序，再按照$y$为第二关键字，从小到大排序
• 可以发现未删除的点，如果有对称中心，则必然是$(p[1]+p[n])/2$，倘若$p[1]$不和$p[n]$对称，假设和$p[k]$对称，而$k<n$，则必然存在某点$p[i].x<p[1].x$使得$p[i]+p[n]=p[1]+p[k]$，但是$p[1]$本身就是$x$最小的点，所以会有矛盾
• 按照这个思路，即使删除中间的点，那么对称中心也已经由最两边的点确定，所以最终只要分别枚举两端$[l,r](1\le l\le k+1,n-k\le r\le n)$，而中心点就是$(p[l]+p[r])/2$，最后$O(n)$双指针判断$[l,r]$中哪些点必须删除即可，检查删除点个数是否合法。

AC代码：

#include <bits/stdc++.h>
#define ft first
#define sd second
#define pb push_back
#define IOS ios::sync_with_stdio(false), cin.tie(0), cout.tie(0) //不能跟puts混用
#define seteps(N) fixed << setprecision(N)
#define endl "\n"
const int maxn = 2e5 + 10;
using namespace std;
typedef long long ll;
typedef double db;
const db eps = 1e-8;
typedef pair<int, int> pii;
struct Point {
    int x, y;
    Point operator + (Point p) {return {x + p.x, y + p.y};}
    bool operator < (const Point &p) const {return make_pair(x, y) < make_pair(p.x, p.y);}
    bool operator == (const Point &p) const {
        if (x == p.x && y == p.y) return true;
        else return false;
    }
} p[maxn];
set <Point> ans;
int n, k;
int main() {
    IOS;
    cin >> n >> k;
    for (int i = 1; i <= n; i++) cin >> p[i].x >> p[i].y;
    if (n <= k) {
        cout << -1 << endl;
        return 0;
    }
    sort (p + 1, p + n + 1);
    for (int i = 1; i <= min(n, k + 1); i++) { 
        for (int j = n; j >= max(i, n - k); j--) {
            int num = i - 1 + n - j;
            if (num > k) continue;
            Point cent = (p[i] + p[j]);
            int l = i, r = j;
            while (l <= r) {
                if ((p[l] + p[r])  == cent) l++, r--;
                else if ((p[l] + p[r]) < cent) l++, num++;
                else r--, num++;
            }
            if (num > k) continue;
            ans.insert(cent);
        }
    }
    cout << ans.size() << endl;
    for (auto np : ans)
        cout << seteps(12) << np.x / 2.0 << " " << np.y / 2.0 << endl;
    return 0;
}