【校内模拟】多边形（凸包DP）（极角序扫描线）（前缀和优化DP）

简要题意：

给平面一个点集，对内空凸多边形计数。

保证没有三点共线。

$n\le 500$

---

题解：

模板题难度的凸包DP?（大概）

然而这是我第一次写凸包DP，以前都是靠口胡。

按照 $x$ 为第一关键字，$y$ 为第二关键字进行坐标排序，枚举多边形左下角的点。

设 $f[o][i][j]$ 表示起始点标号为 $o$ 按照极角序考虑的最后一条边为 $(i,j)$ 的空凸多边形有多少个，转移直接枚举下一个点看是否满足凸性和内空即可。容易发现 $o$ 那一维之间不会转移，可以丢掉。

容易发现可以用极角序扫描线+前缀和进行优化。$O(n^3)$ 绝对不算难写，但是写得轻不轻松就各凭本事了。

我这个考场忘记清空的人可能没有资格说这句话

---

代码：

#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;}