[BZOJ1132][POI2008]Tro（计算几何）

题目描述

传送门

题目大意：平面上有N个点. 求出所有以这N个点为顶点的三角形的面积和。

题解

枚举每一个点，以它为原点建立平面直角坐标系，然后其他的点按照极角排序。
向量的叉积满足分配律，所以可以two pointer一下，每一次右指针移动到最后一个与左指针夹角小于180°的地方，然后每加入一个向量，计算一下这个向量与两个指针之间所有向量的叉积，利用分配律。
看了PoPoQQQ的题解才发现，其实不用这么麻烦，可以将第一、四象限和x、y轴正半轴上的点按照极角排序，然后维护一个前缀和
原来是用了愚蠢的方法死慢死慢
不过这题还卡精度，建议全程int/long long，最后先*10再手动除

代码

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
using namespace std;
#define N 6005const double pi=acos(-1.0);
struct Vector
{int x,y;double ang;Vector(int X=0,int Y=0){x=X,y=Y;ang=atan2((double)y,(double)x);}bool operator < (const Vector &a) const{return ang<a.ang;}
};
typedef Vector Point;int n,vx,vy,x[N],y[N];
long long ans;
Point q[N];long long Abs(long long x){return (x>0)?x:-x;}
int main()
{scanf("%d",&n);for (int i=1;i<=n;++i) scanf("%d%d",&x[i],&y[i]);for (int i=1;i<=n;++i){int cnt=0;for (int j=1;j<=n;++j)if (i!=j) q[++cnt]=Point(x[j]-x[i],y[j]-y[i]);sort(q+1,q+cnt+1);for (int i=1;i<=cnt;++i) q[cnt+i]=q[i],q[cnt+i].ang+=2*pi;int l=1,r=0;while (l<=cnt){if (l>r) r=l,vx=0,vy=0;while (r<=cnt<<1&&q[r].ang<q[l].ang+pi){ans+=Abs((long long)q[r].x*vy-(long long)q[r].y*vx);if (r<=cnt) vx+=q[r].x,vy+=q[r].y;++r;}vx-=q[l].x,vy-=q[l].y;++l;}}ans=ans*5/3;printf("%lld.%lld\n",ans/10,ans%10);
}

[BZOJ1132][POI2008]Tro（计算几何）相关推荐

bzoj1132 [POI2008]Tro（叉积）
bzoj1132 [POI2008]Tro 原题地址:http://www.lydsy.com/JudgeOnline/problem.php?id=1132 题意: 平面上有N个点. 求出所有以这N ...
【bzoj1132】[POI2008]Tro 计算几何
按照横坐标为第一关键字,纵坐标为第二关键字排序, 枚举每个点作为原点把再它后面的点拿出来,求每一对点与原点构成的三角形的面积 abs(a[i].x*a[j].y-a[i].y*a[j].x) 如果能 ...
bzoj1132:[POI2008]Tro
[传送门])(https://www.lydsy.com/JudgeOnline/problem.php?id=1132) 自己的计算几何还是挺渣的呢. 其实只要考虑一下叉积的式子就可以想到前缀和了, ...
1132: [POI2008]Tro 计算几何
首先O(N^3)是不科学的..我们起码要想办法优化到O(N^2logN). 考虑每次枚举一个点,以它作为一个顶点,然后维护一个前缀和来计算每个点与夹角在逆时针180°内的点(我表述的不好..边界可以二 ...
BZOJ1132: [POI2008]Tro
所以为啥用long long本地能过所有数据交上去就WA,一定要改int啊一开始写了个用点到直线距离公式的,维护两个指针扫过去 emmmmmm并不能过样例,因为有根号什么的东西,调不动-. 膜题解发 ...
bzoj1132 [POI2008]Tro
题目描述平面上有N个点. 求出所有以这N个点为顶点的三角形的面积和 N<=3000 输入第一行给出数字N,N在[3,3000] 下面N行给出N个点的坐标,其值在[0,10000] 输出保留 ...
【BZOJ1132】[POI2008]Tro 几何
[BZOJ1132][POI2008]Tro Description 平面上有N个点. 求出所有以这N个点为顶点的三角形的面积和 N<=3000 Input 第一行给出数字N,N在[3,3000 ...
[POI2008]Tro
题面 <center>2042: [POI2008]Tro<center> <center>时间限制:20秒内存限制:162MB<center> 题目 ...
bzoj 1132: [POI2008]Tro
1132: [POI2008]Tro Time Limit: 20 Sec Memory Limit: 162 MB Submit: 1663 Solved: 558 [Submit][Statu ...

[BZOJ1132][POI2008]Tro（计算几何）

题目描述

题解

代码

[BZOJ1132][POI2008]Tro（计算几何）相关推荐

最新文章

热门文章