首先利用凸包判断给出的点是否正好组成一个凸多边形,如果可以的话,定义状态dp[i][j]为从i点到j点构成的凸多边形的最优解。因为凸多边形是一个环,所以我们将节点后面再复制一遍节点,这样方便计算。然后就可以按照区间dp去做了。
转移方程: dp[i][j]=min(dp[i][j],dp[i][k]+dp[k][j]+temp);k为枚举的切割点,temp为割i到j点的代价。当i和j相邻时,显然temp为0,因为它是一条边,不需要我们割。
代码:
#include<stdio.h>
#include<iostream>
#include<algorithm>
#include<string.h>
#include<math.h>
#include<stack>
using namespace std;
#define inf 0x3f3f3f3f
#define pi acos(-1)
#define maxn 650
int xx,yy;
int dp[maxn][maxn];
struct Node
{
int x;
int y;
} node[maxn];
Node st[maxn];
bool cmp1(Node a,Node b)
{
if(a.y==b.y)return a.x<b.x;
return a.y<b.y;
}
bool cmp2(Node a,Node b)
{
if(atan2(a.y-yy,a.x-xx)==atan2(b.y-yy,b.x-xx))
return a.x<b.x;
return atan2(a.y-yy,a.x-xx)<atan2(b.y-yy,b.x-xx);
}
int chaji(Node a,Node b,Node c)
{
return (b.x-a.x)*(c.y-a.y)-(b.y-a.y)*(c.x-a.x);
}
int main()
{
int n,p;
while(~scanf("%d%d",&n,&p))
{
for(int i=1; i<=n; i++)
{
scanf("%d%d",&node[i].x,&node[i].y);
}
sort(node+1,node+n+1,cmp1);
st[1]=node[1];
xx=node[1].x;
yy=node[1].y;
sort(node+1,node+n+1,cmp2);
st[2]=node[2];
int top=2;
for(int i=3; i<=n; i++)
{
while(chaji(st[top-1],st[top],node[i])<=0)//注意此处还要判断等于
{
top--;
}
st[++top]=node[i];
}
// for(int i=1;i<=top;i++)
// {
// cout<<st[i].x<<" "<<st[i].y<<endl;
// }
if(top!=n)//全部节点不能构成一个凸包
{
printf("I can't cut.\n");
continue;
}
for(int i=1; i<=n; i++)//因为是一个环,我们在复制一遍节点,然后则这个数组上dp
{
node[i+n]=node[i];
}
memset(dp,inf,sizeof(dp));
int temp;
for(int len=3; len<=n; len++)
{
for(int i=1; i<=2*n-len+1; i++)
{
int j=i+len-1;
for(int k=i+1; k<j; k++)
{
if(k-i+1<3)dp[i][k]=0;//这部分的代价为0,因为不能构成一个三角形
if(j-k+1<3)dp[k][j]=0;//同理
if(abs(i-j)==1||abs(i-j)==n-1)temp=0;//i和j相邻,不需要割
else
{
temp=(abs(node[i].x+node[j].x)*abs(node[i].y+node[j].y))%p;
}
dp[i][j]=min(dp[i][j],dp[i][k]+dp[k][j]+temp);//维护最小值
}
}
}
int ans=inf;
for(int i=1; i<=n; i++)
{
ans=min(ans,dp[i][i+n-1]);//找出最优解
}
printf("%d\n",ans);
}
}
/*
4 6
6 1
7 9
4 12
1 9
0
*/