luogu4677山区建小学题解--区间DP

题目链接

https://www.luogu.org/problemnew/show/P4677

分析

这道题方法跟之前题不一样,我们相当于枚举一个左右端点来线性扩展,同时划分断点进行决策

\(f[i][j]\)表示在前\(i\)个村庄中建立\(j\)个小学的最小距离总和

我们将枚举到第\(i\)个村庄作为阶段,修了\(j\)所小学作为状态,通过枚举断点\(k\)来分割第\(j\)所小学与前\(j-1\)所小学

也就是说我们判断\(f[k][j-1]\)加上将新加入的第\(j\)座小学建在后面的第\(k+1\)到第\(i\)座村庄中作出的贡献（也就是新产生的距离,我们假设\(f[k][j-1]\)已经是最优的）是否更优,那么怎么这个贡献怎么求呢呢?比较显然当小学建在\([k+1,i]\)中点处产生的新距离之和最小.为了快速求我们可以先预处理出来

状态转移

memset(f,0x3f,sizeof(f));
    f[0][0]=0;
    for(ri i=1;i<=n;i++){//枚举第几座村庄 
        for(ri j=1;j<=min(i,m);j++){//枚举修了多少小学 
            for(ri k=j-1;k<i;k++){//断点,枚举前j-1所小学都建在了[1,k]这个区间村庄内 
                f[i][j]=min(f[i][j],f[k][j-1]+dis[k+1][i]);
            }
        }
    }

代码

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <cctype>
#include <cmath>
#define ll long long 
#define ri register int
using std::min;
using std::abs;
using std::max;
template <class T>inline void read(T &x){
    x=0;int ne=0;char c;
    while(!isdigit(c=getchar()))ne=c=='-';
    x=c-48;
    while(isdigit(c=getchar()))x=(x<<3)+(x<<1)+c-48;
    x=ne?-x:x;return ;
}
const int maxn=505;
const int inf=0x7fffffff;
int f[maxn][maxn],s[maxn][maxn],dis[maxn][maxn];
int n,m,pos[maxn];
/*inline int dis(int l,int r){
    int mid=(l+r)>>1,x=0;
    for(ri i=l;i<=r;i++)x+=abs(pos[i]-pos[mid]);
    return x;
}*/
int main(){
    read(n),read(m);
    for(ri i=2;i<=n;i++){
        read(pos[i]);
        pos[i]+=pos[i-1];
    }
    for(ri i=1;i<=n;i++){
    	for(ri j=i;j<=n;j++){
    		s[i][j]+=s[i][j-1]+pos[j];
        }
    }
    for(ri l=1;l<=n;l++){
        for(ri r=l;r<=n;r++){
            int mid=(l+r)>>1;
            dis[l][r]+=(mid-l)*pos[mid]-s[l][mid-1];
            dis[l][r]+=s[mid+1][r]-(r-mid)*pos[mid];
        }
    }
    memset(f,0x3f,sizeof(f));
    f[0][0]=0;
    for(ri i=1;i<=n;i++){//枚举第几座村庄 
        for(ri j=1;j<=min(i,m);j++){//枚举修了多少小学 
            for(ri k=j-1;k<i;k++){//断点,枚举前j-1所小学都建在了[1,k]这个区间村庄内 
                f[i][j]=min(f[i][j],f[k][j-1]+dis[k+1][i]);
            }
        }
    }
    printf("%d\n",f[n][m]);
    return 0;
}