(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++){//确定AD
while(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
",ans<<2);
}