分析:
一看这道题,就有一个疑问进入我的脑海:
我们要怎么确定一个村庄到底去哪个邮局呢?
实际上邮局的存在就把村庄分成了若干区间
一个区间内的村中都去一个邮局
如果我们确定了区间,那么要找到区间最优解是不难的(枚举暴力计算)
也就是说整个序列被我们划分成了P个区间,每个区间内的最优解我们可以预处理出
所以设计状态:f[i][j]表示前i个村庄被划分成了j个区间(建了j个邮局)
转移方程:f[i][j]=min{f[k][j-1]+c[k+1][i]}
c[i][j]表示在区间(i,j)内建一个邮局的最小代价
c的预处理要注意一下(推柿子就好了)
//这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int INF=0x33333333;
int dis[305],sum[305];
int f[305][35],c[305][305];
int n,m;
void pre()
{
int i,j,k;
for (i=1;i<=n;i++)
{
c[i][i]=0;
for (j=i+1;j<=n;j++)
{
int s=sum[j]-sum[i]-dis[i]*(j-i);
c[i][j]=s;
for (k=i+1;k<=j;k++)
{
s=s+(2*k-i-j-1)*(dis[k]-dis[k-1]);
c[i][j]=min(c[i][j],s);
}
}
}
}
void doit()
{
int i,j,k;
memset(f,0,sizeof(f));
for (int i=1;i<=n;i++) f[i][1]=c[1][i];
for (i=2;i<=n;i++)
for (j=2;j<=min(i,m);j++)
{
f[i][j]=INF;
for (k=1;k<i;k++)
f[i][j]=min(f[i][j],f[k][j-1]+c[k+1][i]);
}
printf("%d\n",f[n][m]);
}
int main()
{
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++)
scanf("%d",&dis[i]);
int a=dis[1];
for (int i=1;i<=n;i++)
{
dis[i]-=a; //每个点到第一个点的距离
sum[i]=sum[i-1]+dis[i]; //dis的前缀和
}
pre();
doit();
return 0;
}
之后我又发现c数组的计算完全不用这么麻烦,
就像UVa11300中说的,取最优值得位置一定是区间的中点
所以我们实际上可以O(1)出解
但是这道题还没完,因为我们现在的效率是O(n^3)的
但是我们可以利用四边形不等式优化到O(n^2)
//这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int INF=0x33333333;
int dis[305],sum[305];
int f[305][35],c[305][305],s[305][305];
int n,m;
int get(int i,int j)
{
if ((i+j)%2==0) return sum[j]-sum[(i+j)/2]-sum[(i+j)/2-1]+sum[i-1];
else return sum[j]-sum[(i+j)/2]-sum[(i+j)/2-1]+sum[i-1]-dis[(i+j)/2];
}
void pre()
{
int i,j,k;
for (i=1;i<=n;i++)
{
c[i][i]=0;
for (j=i+1;j<=n;j++)
c[i][j]=get(i,j);
}
}
void doit()
{
int i,j,k;
memset(f,0x33,sizeof(f));
for (i=1;i<=n;i++) f[i][1]=c[1][i];
for (j=2;j<=m;j++)
{
s[n+1][j]=n-1;
for (i=n;i>=j;i--)
for (k=s[i][j-1];k<=s[i+1][j];k++)
if (f[i][j]>f[k][j-1]+c[k+1][i])
{
f[i][j]=f[k][j-1]+c[k+1][i];
s[i][j]=k;
}
}
printf("%d\n",f[n][m]);
}
int main()
{
scanf("%d%d",&n,&m);
for (int i=1;i<=n;i++)
{
scanf("%d",&dis[i]);
sum[i]=sum[i-1]+dis[i];
}
pre();
doit();
return 0;
}