Cake (凸包+區間dp)

題意:一個多邊形先是確定是否是凸多邊形,可以用凸包模板來判斷。然後是求最小的切割價值。

最小切割價值用區間dp來求解dp【i】【j】=min(dp【i】【j】,dp【i】【k】+dp【k】【j】

#include<bits/stdc++.h>
using namespace std;
int dp[306][306];
double add(double a,double b)
{
	if(abs(a+b)<1e-10*(abs(a)+abs(b)))
	return 0;
	return a+b;
}
struct P{
	double x,y;
	P(){}
	P(double x,double y) : x(x),y(y) {
	}
	P operator - (P p) {
		return P(add(x,-p.x),add(y,-p.y));
	}
	double det(P p)
	{
		return add(x*p.y,-y*p.x);
	}
}ps[306];
bool cmp_x(const P& p,const P& q)
{
	if(p.x!=q.x) return p.x<q.x;
	return p.y<q.y;
}
vector<P> concex_hull(P* ps,int n)
{
	sort(ps,ps+n,cmp_x);
	int k=0;vector<P> qs(n*2);
	for(int i=0;i<n;i++)
	{
		while(k>1&&(qs[k-1]-qs[k-2]).det(ps[i]-qs[k-1])<=0) k--;
		qs[k++]=ps[i];
	}
	for(int i=n-2,t=k;i>=0;i--)
	{
		while(k>t&&(qs[k-1]-qs[k-2]).det(ps[i]-qs[k-1])<=0) k--;
		qs[k++]=ps[i];
	}
	qs.resize(k-1);
	return qs;
}
int e[306][306];
int main()
{
	int n,m;
	while(~scanf("%d %d",&n,&m))
	{
		memset(dp,0x3f,sizeof(dp));
		for(int i=0;i<n;i++)
		{
			scanf("%lf %lf",&ps[i].x,&ps[i].y);
		}
		vector<P> qs=concex_hull(ps,n);
		for(int i=0;i<n;i++)
		{
			for(int j=i+2;j<n;j++)
			{
				e[i][j]=e[j][i]=(int)(abs(qs[i].x+qs[j].x)*abs(qs[i].y+qs[j].y))%m;
			}
		}
		if(n!=qs.size())
		{
			printf("I can't cut.\n");
		}
		else
		{
			for(int i=0;i<n;i++)
			dp[i][(i+1)%n]=0;
			for(int len=2;len<n;len++)
			{
				for(int i=0;i+len<n;i++)
				{
					int j=i+len;
					for(int k=i+1;k<j;k++)
					{
							dp[i][j]=min(dp[i][j],dp[i][k]+dp[k][j]+e[i][k]+e[k][j]);
					}
				}
			}
			printf("%d\n",dp[0][n-1]);
		}
	}
	return 0;
 } 

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章