ZOJ 3537

这道题，就是先判断是不是凸多边形。先求出凸包，如果所有点都在凸包上，那这些点就围成了凸多边形，反之则不是。然后就是分割凸多边形了，刚开始也想到了先确定一条边所在的三角形，从而确定分割的边。但边没选好，导致后续的状态很凌乱。后来看了相关的资料。先按一个顺序将多边形的n个点编号1--n。求f[1][n]先以边n--1为基准，确定其所在三角形，然后就得到了两个多边形。得到状态转移方程，f[1][n]=min(f[1][n],f[1][k]+f[k][n])。对于f[1][k]和f[k][n]分为以分割出这两个多边形的那两条边为基准求出。这样得到的多边形的状态都是i--j（i<j,表示已凸多边形节点i连续的到j,i--j之间有条边得到的多边形）。边界的确定具体代码。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
#include <cmath>
#include <stack>
#include <vector>
#define LL long long
#define myabs(x) ((x)>0?(x):(-(x)))
using namespace std;
const int inf=0x3f3f3f3f;
const int maxn=300+10;
int f[maxn][maxn];
struct Point 
{
	int x,y;
};
Point operator - (Point a,Point b)
{
	Point tem={a.x-b.x,a.y-b.y};
	return tem;
}
Point p[maxn];	
int n,w;
int ch[maxn];
bool cmp(Point a,Point b)
{
	if(a.x==b.x) return a.y<b.y;
	else return a.x<b.x;
}
int cross(Point a,Point b)
{
	return a.x*b.y-a.y*b.x;
}
int judge_convex()
{
	sort(p,p+n,cmp);
	int i,m,k;
	ch[0]=0;
	m=0;
	for(i=1;i<n;i++)
	{
		while(m>0&&cross(p[ch[m]]-p[ch[m-1]],p[i]-p[ch[m]])<0) m--;
		ch[++m]=i;
	}
	k=m;
	for(i=n-2;i>=0;i--)
	{
		while(k>m&&cross(p[ch[k]]-p[ch[k-1]],p[i]-p[ch[k]])<0) k--;
		ch[++k]=i;
	}
	if(n>1) k--;
	if((k+1)==n) return 1;
	else return 0;
}
int cost(int x1,int y1,int x2,int y2)
{
	return (myabs(x1+x2)*myabs(y1+y2))%w;
}
int main()
{
	while(~scanf("%d%d",&n,&w))
	{
		int i;
		for(i=0;i<n;i++) scanf("%d%d",&p[i].x,&p[i].y);
		if(judge_convex())
		{
			memset(f,inf,sizeof(f));
			int j,k;
			for(i=0;i<n-1;i++)
				f[i][i+1]=-cost(p[ch[i]].x,p[ch[i]].y,p[ch[i+1]].x,p[ch[i+1]].y);
			for(i=2;i<n;i++)
			{
				for(j=0;j+i<n;j++)
				{
					for(k=j+1;k<j+i;k++)
					{
						f[j][j+i]=min(f[j][j+i],f[j][k]+f[k][j+i]+cost(p[ch[j]].x,p[ch[j]].y,p[ch[k]].x,p[ch[k]].y)+cost(p[ch[j+i]].x,p[ch[j+i]].y,p[ch[k]].x,p[ch[k]].y));
					}
				}
			}
			printf("%d\n",f[0][n-1]);
		}
		else  printf("I can't cut.\n");
	}
	return 0;				
}