Luogu5286

\(2019.4.14\)，新生第一题，改了\(3\)个小时

题解-租酥雨，和出题人给的正解一模一样

枚举\(AD\)，分别考虑鱼身\(BC\)和鱼尾\(EF\)
到\(E\)，\(F\)距离相等的判断要避免精度问题，可以离散化其平方以后用一种莫队思想搞
关键是怎么枚举垂直，还要控制中点也在\(AD\)上

判垂直和平行的关键就是
与向量\((a,b)\)点积相同的点在一条垂直向量\((a,b)\)的直线上，与向量\((a,b)\)叉积相同的点在一条平行向量\((a,b)\)的直线上

然后加一些转化来限制鱼身的范围：
中点在\(AD\)上；
端点相对于直线\(AD\)在线段\(AD\)上(这样的取值是一个区间)
手推(逆推)就比较简单了

细节非常多，等以后学了计算几何，有时间再打一遍

\(0.\)考场上没看清题
\(1.\)有些(可能是少数)精度问题(不要先除)可以利用一些离散的方法避免；相应地，统计的方法也可以稍加变换

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long LL;
const int INF=1e9+7;
inline LL read(){register LL x=0,f=1;register char c=getchar();while(c<48||c>57){if(c=='-')f=-1;c=getchar();}while(c>=48&&c<=57)x=(x<<3)+(x<<1)+(c&15),c=getchar();return f*x;
}const int N=1005;
const double eps=1e-10;//这题坐标1e9,所以精度至少要1e-10
const double Pi=acos(-1);struct Point{LL x,y;double k;//极角inline void cal(){k=atan2(y,x);}inline bool operator< (const Point &b) const {return k<b.k;}
}p[N],q[N<<1];struct Node{LL a,b,c,d;inline bool operator< (const Node &z) const{ //这个const不能少if(a!=z.a) return a<z.a;if(b!=z.b) return b<z.b;if(c!=z.c) return c<z.c;return d<z.d;}
}S[N*N];LL len[N],tmpl[N];int cnt[N];
int n,m,Scnt,Lcnt;LL ans,sum;inline void add(int x){if(x>m) x-=m;sum+=cnt[len[x]];cnt[len[x]]++;}//这样实现非常优秀,避免了精度问题
inline void del(int x){if(x>m) x-=m;cnt[len[x]]--;sum-=cnt[len[x]];}inline LL sqr(int x){return 1ll*x*x;}int main(){n=read();for(int i=1;i<=n;i++) p[i]=(Point){read(),read()};for(int i=1;i<=n;i++)for(int j=i+1;j<=n;j++){LL a=p[i].x-p[j].x,b=p[i].y-p[j].y,d=__gcd(abs(a),abs(b));a/=d,b/=d;if(a<0) a=-a,b=-b;if(a==0) b=abs(b);LL P=1ll*a*(p[i].x+p[j].x)+1ll*b*(p[i].y+p[j].y);//相应的要预处理出P和Q的值LL Q=1ll*a*p[i].y-1ll*b*p[i].x;S[++Scnt]=(Node){a,b,P,Q};}sort(S+1,S+Scnt+1);for(int i=1;i<=n;i++){m=0;for(int j=1;j<=n;j++){if(i!=j) q[++m]=(Point){p[j].x-p[i].x,p[j].y-p[i].y},q[m].cal();}sort(q+1,q+m+1);for(int j=m+1;j<=m+m;j++) q[j]=q[j-m],q[j].k+=2*Pi;//开两倍避免讨论Lcnt=0;for(int j=1;j<=m;j++) len[j]=tmpl[++Lcnt]=sqr(q[j].x)+sqr(q[j].y);sort(tmpl+1,tmpl+Lcnt+1);Lcnt=unique(tmpl+1,tmpl+Lcnt+1)-tmpl-1;for(int j=1;j<=m;j++) len[j]=lower_bound(tmpl+1,tmpl+Lcnt+1,len[j])-tmpl;//把相同长度的编号设为相同sum=0;memset(cnt,0,sizeof cnt);for(int j=1,p1=0,p2=0;j<=m;j++){//确定ADwhile(q[p1+1].k+eps<q[j].k+1.5*Pi) add(++p1);while(q[p2+1].k-eps<q[j].k+0.5*Pi) del(++p2);LL a=-q[j].y,b=q[j].x,d=__gcd(abs(a),abs(b));//枚举BC:与AD垂直a/=d,b/=d;if(a<0) a=-a,b=-b;if(a==0) b=abs(b);//这些都只是一种规定排序的方法LL P=1ll*a*(p[i].x+(p[i].x+q[j].x))+1ll*b*(p[i].y+(p[i].y+q[j].y));//中点在直线AD上LL L=1ll*a*p[i].y-1ll*b*p[i].x, R=1ll*a*(p[i].y+q[j].y)-1ll*b*(p[i].x+q[j].x);//两端点在AD范围内if(L>R) swap(L,R);ans+=sum*(lower_bound(S+1,S+Scnt+1,(Node){a,b,P,R}) - upper_bound(S+1,S+Scnt+1,(Node){a,b,P,L}));}}printf("%lld\n",ans<<2);
}