bzoj 2087: [Poi2010]Sheep【凸包+极角排序+dp】

首先处理处理出来哪些边能连——能把羊分成两个偶数部分的，实现是在凸包上枚举极点，极角排序，枚举凸包上点对判断两边羊的个数的奇偶即可，设可以连边为v[i][j]=1
然后设f[i][j]为从i到j个凸包上点的方案数，初始状态是相邻点f[i][i+1]=1，转移是

\[f[i][j]=\sum_{k=i+1}^{i-1}[v[i][k]\&\&v[k][j]]f[i][k]*f[k][j] \]

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
const int N=605,M=20005;
int n,m,mod,f[N][N];
bool v[N][N];
struct dian
{
	int x,y;
	dian(int X=0,int Y=0)
	{
		x=X,y=Y;
	}
	dian operator - (const dian &a) const
	{
		return dian(x-a.x,y-a.y);
	}
	int operator * (const dian &a) const
	{
		return x*a.y-y*a.x;
	}
}a[N],b[M],d;
int read()
{
	int r=0,f=1;
	char p=getchar();
	while(p>'9'||p<'0')
	{
		if(p=='-')
			f=-1;
		p=getchar();
	}
	while(p>='0'&&p<='9')
	{
		r=r*10+p-48;
		p=getchar();
	}
	return r*f;
}
inline bool cmp(const dian &a,const dian &b)
{
	return (a-d)*(b-d)<0;
}
int main()
{
	n=read(),m=read(),mod=read();
	for(int i=1;i<=n;i++)
		a[i].x=read(),a[i].y=read();
	for(int i=1;i<=m;i++)
		b[i].x=read(),b[i].y=read();
	for(int i=1;i<=n;++i)
	{
		d=a[i];
        sort(b+1,b+m+1,cmp);
		int x=0;
		for(int j=i+1;j<=n;++j)
		{
			while(x<m&&(b[x+1]-d)*(a[j]-d)<0)
				++x;
			if((x&1)||(b[x+1]-d)*(a[j]-d)==0)
				continue;
			v[i][j]=1;
		}
    }
	for(int i=1;i<=n;i++)
		f[i][i%n+1]=1;
	for(int l=2;l<=n;l++)
		for(int i=1,j;i<=n;i++)
		{
			j=(i+l-1)%n+1;
			if(v[i][j])
				for(int k=i%n+1;k!=j;k=k%n+1)
					if(v[i][k]&&v[k][j])
						f[i][j]=(f[i][j]+f[i][k]*f[k][j])%mod;
		}
	printf("%d\n",f[1][n]);
	return 0;
}