这几天为了补我弱爆了的计算几何基础,在刷刘汝佳白书上的计算几何部分。发现这道题网上的题解几乎都是n3水过去的或者有n2的做法但说是使用了容斥什么的?个人认为运用了叉积分配律的做法应该是最快且最好写的正确做法了,不过网上貌似还没看到这么写的,所以打算来给这道大水题写个题解。
计算三角形的面积其实就是计算其中任意两条边的叉积的绝对值除以2,直接枚举三角形的三个点复杂度是n3的,考虑枚举两个点为三角形左右的点,那么符合条件的第三个点就是夹在这两个点中间的点。
假设左侧点为a,右侧点为c, 中间的包含的点为b1,b2,b3…bn。那么最后符合条件的三角形的总面积(未取绝对值及除2前)为
cross(ab1, ac)+cross(ab2, ac)+…+cross(abn, ac)[cross表示叉积函数,两个参数为向量]
由于叉积符合分配律,所以该式可转化为
cross(ab1+ab2+…+abn, ac)
回到该题,应用如上性质,若对给定点先按照角度排序,在循环的时候处理一下叉积的和,即可n2的解决这道题。
#pragma GCC optimize(2)
#include <bits/stdc++.h>
#define maxn 100010
using namespace std;
typedef long long ll;
const long double eps=1e-6;
const long double PI=acos(-1.0);
int n;
double r, ans, rad[maxn];
struct Point//点&向量
{
double x, y;
Point(double x=0, double y=0): x(x),y(y) {}
}p[maxn], o(0, 0);
Point operator + (Point a, Point b) {return Point(a.x+b.x, a.y+b.y);}
Point operator - (Point a, Point b) {return Point(a.x-b.x, a.y-b.y);}
Point operator * (Point a, double b) {return Point (a.x*b, a.y*b);}
Point operator / (Point a, double b) {return Point (a.x/b, a.y/b);}
double Dot(Point a, Point b) {return a.x*b.x+a.y*b.y;}
double Cross(Point a, Point b) {return a.x*b.y-a.y*b.x;}
Point get_point(double rad)//由圆心角求点坐标
{
return Point(o.x+cos(rad)*r, o.y+sin(rad)*r);
}
int main()
{
while (cin>>n)
{
cin>>r;
ans=0;
if (n==0) break;
for (int i=1; i<=n; i++)
cin>>rad[i];
sort(rad+1, rad+1+n);
for (int i=1; i<=n; i++)
p[i]=get_point(rad[i]/360*2*PI);
for (int i=1; i<n; i++)
{
Point s=p[i+1];
for (int j=i+2; j<=n; j++)
{
ans+=fabs(Cross(p[j]-p[i], s-p[i]*(j-i-1)));
s=s+p[j];
}
}
cout<<fixed<<setprecision(0)<<ans/2<<"\n";
}
}