題意:一個多邊形先是確定是否是凸多邊形,可以用凸包模板來判斷。然後是求最小的切割價值。
最小切割價值用區間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;
}