求出n边形有多少个三角形包含平面上某一点。
利用叉积可以判断两向量的位置关系,从而判断该点是否在n边形内部。
用总三角形数-不包含该点的三角形数即可,再次根据叉积判断该点在不在子多边形内部。
#include<iostream>
#include<string>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<iomanip>
#include<map>
#include<algorithm>
#include<queue>
#include<set>
#define inf 1000000000
#define pi acos(-1.0)
#define eps 1e-8
#define seed 131
using namespace std;
typedef pair<int,int> pii;
typedef unsigned long long ull;
typedef long long ll;
const int maxn=100005;
struct node
{
ll x,y;
}point[maxn*2];
int n,t;
ll cross(node& a,node& b,node& c)
{
ll x1=b.x-a.x;
ll y1=b.y-a.y;
ll x2=c.x-a.x;
ll y2=c.y-a.y;
return x1*y2-x2*y1;
}
int main()
{
scanf("%d",&n);
for(int i=0;i<n;i++)
{
cin>>point[i].x>>point[i].y;
// cout<<1<<endl;
point[i+n]=point[i];
}
scanf("%d",&t);
while(t--)
{
node u;
cin>>u.x>>u.y;
bool flag=true;
for(int i=0;i<n-1;i++)
{
if(cross(point[i],u,point[i+1])<=0)
{
flag=false;
break;
}
}
if(!flag)
{
printf("0\n");
continue;
}
ll ans=(ll)n*(n-1)*(n-2)/6;
int j=2;
for(int i=0;i<n;i++)
{
while(cross(point[i],u,point[j])>0)
{
j++;
}
ans-=(ll)(j-i-1)*(j-i-2)/2;
}
cout<<ans<<endl;
}
return 0;
}