您的位置:首页 > 其它

[2016ACM多校] HDU5784 极角序 正难则反

2016-08-05 18:16 267 查看

题意

给N个点,问组成多少个锐角三角形。N<=2000

思路

一开始想的是枚举线段,然后构成锐角三角形的另一个点在线段两端点的垂线内,还在以线段为直径的圆外。垂线好搞,圆卡到比赛结束。赛后发现直接用三角形总数减去非锐角的个数就可以了,没用正难则反这个思维就去坑第五题T到爆。


在算非锐角个数是,平角是个特殊情况,所以不能算钝角直角相加,而要算非锐角再除以2。最后代码只写了40行,感觉挺不错的。

AC代码 C++

#include <stdio.h>
#include <math.h>
#include <algorithm>

using namespace std;

#define MAXN 2002
#define PI 3.1415926535
#define EPS (1e-10)
#define EDIT_BORD(ang) ((ang) > PI ? (ang) - 2.0 * PI : (ang))

pair<int, int> point[MAXN];
double ang[MAXN];

int main()
{
int n, i, j, k;
long long sum;
double l, r;
while(scanf("%d", &n) > 0)
{
for(i=0; i<n; i++)
scanf("%d%d", &point[i].first, &point[i].second);
for(i=sum=0; i<n; i++)
{
for(j=k=0; j<n; j++)
if(i != j)
ang[k++] = atan2(point[j].second-point[i].second, point[j].first-point[i].first);
sort(ang, ang+k);
for(j=0; j<k; j++)
{
l = EDIT_BORD(ang[j] + PI * 0.5 - EPS);
r = EDIT_BORD(ang[j] + PI * 1.5 + EPS);
sum += max(l < r ? upper_bound(ang, ang+k, r) - upper_bound(ang, ang+k, l) : k - (upper_bound(ang, ang+k, l) - upper_bound(ang, ang+k, r)), 0LL);
}
}
printf("%I64d\n", (long long)n * (n - 1) * (n - 2) / 6 - sum / 2);
}
return 0;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签:  acm