bzoj4953 [Wf2017]Posterize

http://www.elijahqi.win/archives/662
Description

数字图像的像素可以用三个在0到255之间的整数表示,它们分别表示红色、绿色和蓝色的强度。为了压缩图片或是
为了产生艺术效果,许多图像编辑工具收录了如下所述的”色调分离”操作。每个颜色通道会分别考虑,本题只考虑红
色通道的情况。不同于在红色通道使用0到255之间全部的整数,一张色调分离后的图片只会使用这些数字里至多 k
种整数。每个像素原来的红色强度会被替换成最相近的可用强度。图像编辑工具会选择k个整数来最小化替换过程
引起的平方误差之和。假设原图有n个像素,它们的红色取值是r1,···,rn,而 k 种可用整数为v1,···,vk ,那
么平方误差之和被定义为

你的任务是计算可以实现的最小平方误差之和,参数k和图片的红色强度会给出。
Input

第一行包含两个整数d(1≤d≤256)和k(1≤k≤d)
分别表示原图中不同的红色强度有多少种,色调分离后可以使用的红色强度有多少种。
接下来d行描述了每种红色强度在原图中占据的像素点数量。
每行包含两个整数r(0≤r≤255)和p(1≤p≤226)
这里r是一种红色强度的取值,而p是这种取值对应的像素点数量。这d行信息按照红色强度取值升序给出。
Output

输出最优的 k 种可选取值对应的平方误差之和。
Sample Input

样例1
2 1
50 20000
150 10000
样例2
2 2
50 20000
150 10000
样例3
4 2
0 30000
25 30000
50 30000
255 30000
Sample Output

样例1
66670000
样例2
0
样例3
37500000

leoly 写了斜率优化，现在排在bzoj榜首orz

原题大致的意思和曾经做过的ioi2000 postoffice有点像，和ceoi2004锯木厂也有点像
题意：
一共有256个点，现在给出其中n个点上有东西要搬到k个点中，求所需要代价最小是多少
所需要的代价就是两点间距离的平方乘以要搬运多少东西
我只写了一个非常暴力的dp
f[i][j] i表示1~i中一定有一个点安放在i点，并且已经安放了j个点
那么我们可以暴力去枚举一下放在i前面的那个点，统计答案
f[i][j]表示前i个位置，放了j个旗子，其中第i个位置一定放了旗子的最小平方误差，则f[i][j]=min{f[k][j-1]+w[k][i]|j-1<=k< i}，其中w[k][i]表示k~i这一段只有k和i位置放了旗子的平方误差。最后统计答案，枚举最后一个旗子插的位置i，加上i之后的人都要走到i位置的距离平方和，取最小值，即为答案。复杂度O(n3)

#include<cstdio>
#define N 330
#define LL long long
inline int read(){
    int x=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=getchar();}
    return x;
}
int n,k,d[N],r,p;
LL w[N][N],f[N][N];
inline LL min(LL a,LL b){
    return a<b?a:b;
}
int main(){
    freopen("bzoj4953.in","r",stdin);
    n=read();k=read();
    for (int i=1;i<=n;++i) r=read(),p=read(),d[r]=p;
    //for (int i=0;i<=255;++i) printf("%d\n",d[i]);
    for (int i=0;i<=255;++i){
        for (int j=i;j<=255;++j){
            int len=j-i-1;
            for (int z=i+1;z<j;++z){
                if ((z-i<<1)<=len) w[i][j]+=(LL)d[z]*(z-i)*(z-i);else w[i][j]+=(LL)d[z]*(j-z)*(j-z);
            }
        }
    }
    for(int i=0;i<=255;++i)
        for(int j=0;j<i;++j) 
            f[i][1]+=(LL)(i-j)*(i-j)*d[j];
    for (int j=2;j<=k;++j){
        for (int i=j-1;i<=255;++i){
            f[i][j]=f[i-1][j-1];
            for (int z=j-2;z<=i-1;++z){
                f[i][j]=min(f[z][j-1]+w[z][i],f[i][j]);
            }
        }
    }
    //printf("%lld",f[255][2]);
    LL ans=f[255][k];
    for (int i=k-1;i<=255;++i){
        LL tmp=0;
        for (int j=i+1;j<=255;++j){
            tmp+=(LL) d[j]*(j-i)*(j-i);
        }
        ans=min(ans,f[i][k]+tmp);
    }
    printf("%lld",ans);
    return 0;
}