简要题意:
给平面一个点集,对内空凸多边形计数。
保证没有三点共线。
n ≤ 500 n\leq 500 n≤500
题解:
模板题难度的凸包DP?(大概)
然而这是我第一次写凸包DP,以前都是靠口胡。
按照 x x x 为第一关键字, y y y 为第二关键字进行坐标排序,枚举多边形左下角的点。
设 f [ o ] [ i ] [ j ] f[o][i][j] f[o][i][j] 表示起始点标号为 o o o 按照极角序考虑的最后一条边为 ( i , j ) (i,j) (i,j) 的空凸多边形有多少个,转移直接枚举下一个点看是否满足凸性和内空即可。容易发现 o o o 那一维之间不会转移,可以丢掉。
容易发现可以用极角序扫描线+前缀和进行优化。 O ( n 3 ) O(n^3) O(n3) 绝对不算难写,但是写得轻不轻松就各凭本事了。
我这个考场忘记清空的人可能没有资格说这句话
代码:
#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const
using std::cerr;
using std::cout;
cs int mod=1e9+7;
inline void Inc(int &a,int b){a+=b-mod;a+=a>>31&mod;}
struct Pnt{
int x,y;Pnt(){}Pnt(int _x,int _y):x(_x),y(_y){}
friend Pnt operator+(cs Pnt &a,cs Pnt &b){return Pnt(a.x+b.x,a.y+b.y);}
friend Pnt operator-(cs Pnt &a,cs Pnt &b){return Pnt(a.x-b.x,a.y-b.y);}
friend ll operator*(cs Pnt &a,cs Pnt &b){return (ll)a.x*b.y-(ll)b.x*a.y;}
};
inline ll crs(cs Pnt &a,cs Pnt &b,cs Pnt &c){
return (ll)(b.x-a.x)*(c.y-a.y)-(ll)(c.x-a.x)*(b.y-a.y);
}
cs int N=5e2+7;
int n;
Pnt p[N];
inline ll crs(int a,int b,int c){
return crs(p[a],p[b],p[c]);
}
int id[N][N<<1|1];
int f[N][N],ans;
int ql[N],qr[N],tl,tr;
void get_lr(int o,int nw){
int ps=std::find(id[nw]+2,id[nw]+n+1,o)-id[nw];
tl=tr=0;for(int re i=ps;i<ps+n;++i)if(id[nw][i]>o){
if(crs(nw,o,id[nw][i])>0)ql[++tl]=id[nw][i];
else qr[++tr]=id[nw][i];
}
}
void Main(){scanf("%d",&n);
for(int re i=1;i<=n;++i)
scanf("%d%d",&p[i].x,&p[i].y);
std::sort(p+1,p+n+1,
[](cs Pnt &a,cs Pnt &b){return a.x<b.x||(a.x==b.x&&a.y<b.y);});
for(int i=1;i<=n;++i){
for(int re j=1;j<=n;++j)id[i][j]=j;
id[i][i]=1;
std::sort(id[i]+2,id[i]+i+1,
[&](int a,int b){return crs(i,a,b)>0;});
std::sort(id[i]+i+1,id[i]+n+1,
[&](int a,int b){return crs(i,a,b)>0;});
for(int re j=n+1;j<n+n;++j)id[i][j]=id[i][j-n+1];
}for(int re o=1;o<=n;++o){
for(int re i=o+1;i<=n;++i)
memset(f[i]+o+1,0,(n-o)<<2);
for(int re i=o+1;i<=n;++i){
int t=id[o][i];get_lr(o,t);int tmp=0;
for(int re il=1,ir=1;ir<=tr;)
if(il>tl||crs(t,ql[il],qr[ir])>0)
Inc(f[t][qr[ir]],tmp),++ir;
else
Inc(tmp,f[ql[il]][t]),++il;
for(int re j=tr,nw=0;j;--j)
if(!nw||crs(o,nw,qr[j])<0)
Inc(f[t][qr[j]],1),nw=qr[j];
else f[t][qr[j]]=0;
for(int re j=i+1;j<=n;++j)
Inc(ans,f[t][id[o][j]]);
}
}cout<<ans<<"\n";
}
inline void file(){
#ifdef zxyoi
freopen("polygon.in","r",stdin);
#endif
}signed main(){file();Main();return 0;}