【bzoj2338】[HNOI2011]数矩形

Description

题解

通过这些点我们可以得到很多线段对吧。
两个线段如果作为矩形的对角线，这就要求它们长度相等且互相平分。
互相平分也就是中点重合。
我们把这些线段按照长度为第一关键字，中点位置为第二关键字排序，这样能够构成矩形的线段一定就在一起了，在对于每个线段向前暴力找可以与它构成矩形的线段。这样就OK了。
注意全程用long long而且中点不除以二，距离不开根号，面积用叉乘。这样才能保证精度。
千万别用double啥的。

#include<bits/stdc++.h>using namespace std;const int N = 1500 + 10;typedef long long ll;struct Point{    ll x, y;    Point(){};    Point(ll a, ll b){x = a, y = b;}    bool operator == (const Point &a) const{        return a.x == x && a.y == y;    }    ll operator * (const Point &a) const{        return x * a.y - y * a.x;    }    Point operator - (const Point &a) const{        return Point(x - a.x, y - a.y);    }    Point operator + (const Point &a) const{        return Point(x + a.x, y + a.y);    }}pt[N];struct Line{    ll len;    Point p1, p2, mdp;    bool operator < (const Line &a) const{        if(len == a.len)            return (mdp.x < a.mdp.x) || (mdp.x == a.mdp.x && mdp.y < a.mdp.y);        else return len < a.len;    }}ln[N*N>>1];int n, tot;ll sqr(ll x){return x * x;}ll dist(Point a, Point b){    return sqr(a.x - b.x) + sqr(a.y - b.y);}ll labs(ll x){return x > 0 ? x : -x;}void init(){    scanf("%d", &n);    for(int i = 1; i <= n; i++){        scanf("%lld%lld", &pt[i].x, &pt[i].y);        for(int j = 1; j < i; j++){            ln[++tot].p1 = pt[i];            ln[tot].p2 = pt[j];            ln[tot].len = dist(pt[i], pt[j]);            ln[tot].mdp = pt[i] + pt[j];        }    }}void work(){    ll ans = 0, s;    sort(ln+1, ln+tot+1);    for(int i = 1; i <= tot; i++)    for(int j = i-1; ln[i].mdp == ln[j].mdp && ln[i].len == ln[j].len && j; j--){        s = labs((ln[j].p1 - ln[i].p1) * (ln[j].p2 - ln[i].p1));        if(s > ans) ans = s;    }    printf("%lld\n", ans);}int main(){    init();    work();    return 0;}