Codeforces 1284E New Year and Castle Building (计算几何)
作者:互联网
题目链接
https://codeforces.com/contest/1284/problem/E
题解
我们计算选出 \(3\) 个点构成三角形覆盖的点数之和,这个值乘以 \(\frac{(n-4)}{2}\) 就是答案。这是因为对于任意一个(凸或凹,根据题意凹的多种凹法只算一次)四边形,从 \(4\) 个顶点中选出 \(3\) 个构成三角形的 \(4\) 种方案中每个被四边形覆盖的点恰好被 \(2\) 种方案覆盖。
然后就比较套路了,补集转化,枚举一个点求有多少个三角形不包含它,然后其余的点以它为中心按极角排序,枚举最顺时针的那个,算从它的一侧选出两个点的方案数。
不过似乎没发现三角形那个性质大力扫描线也能做,好题->垃圾题?
时间复杂度 \(O(n^2\log n)\).
代码
#include<bits/stdc++.h>
#define llong long long
#define mkpr make_pair
#define riterator reverse_iterator
#define y1 Lorem_ipsum_dolor
using namespace std;
inline int read()
{
int x = 0,f = 1; char ch = getchar();
for(;!isdigit(ch);ch=getchar()) {if(ch=='-') f = -1;}
for(; isdigit(ch);ch=getchar()) {x = x*10+ch-48;}
return x*f;
}
const int mxN = 2500;
struct Point
{
llong x,y;
Point() {}
Point(llong _x,llong _y):x(_x),y(_y) {}
int quadrant() {return x>0?(y<0):(y<=0);}
} a[mxN+3],b[(mxN<<1)+3];
typedef Point Vector;
Point operator +(Point x,Point y) {return Point(x.x+y.x,x.y+y.y);}
Point operator -(Point x,Point y) {return Point(x.x-y.x,x.y-y.y);}
llong Cross(Vector x,Vector y) {return x.x*y.y-x.y*y.x;}
int n;
bool cmp_ang(Point x,Point y) {return x.quadrant()^y.quadrant()?x.quadrant()<y.quadrant():Cross(x,y)>0;}
int main()
{
scanf("%d",&n);
for(int i=1; i<=n; i++) scanf("%I64d%I64d",&a[i].x,&a[i].y);
llong ans = 0ll;
for(int i=1; i<=n; i++)
{
for(int j=1; j<=n; j++) b[j] = a[j]-a[i]; swap(b[i],b[1]);
sort(b+2,b+n+1,cmp_ang);
// printf("i=%d\n",i);
// for(int j=2; j<=n; j++) printf("(%lld,%lld)\n",b[j].x,b[j].y);
for(int j=2; j<=n; j++) b[n-1+j] = b[j];
int j,k; llong cur = 0ll;
for(j=2,k=3; j<=n; j++)
{
while(k<=j||Cross(b[j],b[k])>0) {k++;}
int cnt = j+n-1-k;
// printf("j=%d k=%d cnt=%d\n",j,k,cnt);
cur += 1ll*cnt*(cnt-1ll)/2ll;
}
ans += cur;
}
ans = (1ll*n*(n-1ll)*(n-2ll)*(n-3ll)/6ll-ans)*(n-4ll)/2ll;
printf("%I64d\n",ans);
return 0;
}
标签:Building,cnt,ch,int,Codeforces,1284E,1ll,llong,define 来源: https://www.cnblogs.com/suncongbo/p/12575009.html