bzoj2338[HNOI2011]数矩形 计算几何

2338: [HNOI2011]数矩形

Time Limit: 20 Sec  Memory Limit: 128 MB
Submit: 1535  Solved: 693
[Submit][Status][Discuss]

Description

Input

 

Output

 

Sample Input

 

Sample Output

 

HINT

 

Source

Day2

我开始想着记每条线的斜率，然后排序来找平行线，但却不能保证构成矩形。
一种新奇的思路：记录对角线，两条对角线长度相等且中点相同时可以确定一个矩形
然后题目就变得简单起来。枚举每两点，存放它们构成的线段，排序。
每条线段向前扩展找可以形成对角线的线段，找到之后更新答案。
为了避免精度错误，一切都可以用整形来计算，即算距离时不开根，中点不除2

#include<bits/stdc++.h>
#define N 1505
#define ll long long
using namespace std;
int n,cnt;ll ans;
struct P{
    int x,y;
    bool operator < (const P &b)const{return x==b.x?y<b.y:x<b.x;}
    P operator - (const P &b)const{return (P){x-b.x,y-b.y};}
    bool operator == (const P &b)const{return x==b.x&&y==b.y;}
}p[N];
struct line{
    ll len;P a,b,mid;
    bool operator < (const line &b)const{return len==b.len?mid<b.mid:len<b.len;}
}l[N*N/2];
ll dis(P a,P b){return 1ll*(a.x-b.x)*(a.x-b.x)+1ll*(a.y-b.y)*(a.y-b.y);}
ll crs(P a,P b){return 1ll*a.x*b.y-1ll*a.y*b.x;}
void update(int i,int j){
    P a=l[i].a,c=l[j].a,d=l[j].b;
    P x=d-a,y=c-a;
    //if(!crs(y,l[i].mid))y=d-a;
    ll tmp=abs(crs(x,y));if(tmp>ans)ans=tmp;
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    scanf("%d%d",&p[i].x,&p[i].y);
    for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++){
            l[++cnt].len=dis(p[i],p[j]);
            l[cnt].a=p[i];l[cnt].b=p[j];
            l[cnt].mid=(P){p[i].x+p[j].x,p[i].y+p[j].y};
        }
    sort(l+1,l+1+cnt);int j;
    for(int i=1;i<=cnt;i++){
        j=i-1;
        while(l[j].len==l[i].len&&l[j].mid==l[i].mid)update(i,j--);
    }
    printf("%lld\n",ans);
    return 0;
}