BZOJ.1132.[POI2008]Tro(极角排序)

时间:2022-06-26 14:28:41

BZOJ

洛谷


考虑暴力,每次枚举三个点,答案就是\(\frac12\sum_{k<j<i}(i-k)\times(j-k)\)。

注意到叉积有分配率,所以固定\(k\),枚举\(i,j\),\(Ans=\frac12\sum_{k<i}(i-k)\sum_{k<j<i}(j-k)\),前缀和即可。

还有个问题是叉积是有符号的。初始时将所有点按纵坐标排序,枚举\(k\)的时候将所有向量\(i-k\)按极角排序。因为\(i>k\),\(i\)都在\(k\)的上方,向量之间不会超过\(180^{\circ}\),符号不会变,就可以直接前缀和了。

算的时候似乎会爆double,用long long存。

复杂度\(O(n^2\log n)\)。

另外可以直接求出向量的极角或斜率,不需要每次cmp算一次叉积。


//920kb	8544ms
#include <cstdio>
#include <cctype>
#include <algorithm>
#define gc() getchar()
typedef long long LL;
const int N=3005; struct Vec
{
LL x,y;
Vec(LL x=0,LL y=0):x(x),y(y) {}
Vec operator +(const Vec &a)const {return Vec(x+a.x, y+a.y);}
Vec operator -(const Vec &a)const {return Vec(x-a.x, y-a.y);}
Vec operator *(const LL a)const {return Vec(x*a, y*a);}
LL operator *(const Vec &a)const {return x*a.y-y*a.x;}
inline bool operator <(const Vec &a)const {return y<a.y||(y==a.y&&x<a.x);}
}p[N],v[N];
typedef Vec Point; inline int read()
{
int now=0;register char c=gc();
for(;!isdigit(c);c=gc());
for(;isdigit(c);now=now*10+c-48,c=gc());
return now;
}
inline bool cmp(const Vec &a,const Vec &b)
{
return a*b>0;
} int main()
{
int n=read();
for(int i=1; i<=n; ++i) p[i]=(Point){read(),read()};
std::sort(p+1,p+1+n);
LL ans=0;
for(int i=1; i+2<=n; ++i)
{
for(int j=i+1; j<=n; ++j) v[j]=p[j]-p[i];
std::sort(v+i+1,v+n+1,cmp);
Vec s=v[i+1];
for(int j=i+2; j<=n; ++j) ans+=s*v[j], s=s+v[j];
}
printf("%lld.%d\n",ans>>1,ans&1?5:0); return 0;
}