BZOJ4518：[SDOI2016]征途——题解

https://www.lydsy.com/JudgeOnline/problem.php?id=4518

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

Pine开始了从S地到T地的征途。

从S地到T地的路可以划分成n段，相邻两段路的分界点设有休息站。

Pine计划用m天到达T地。除第m天外，每一天晚上Pine都必须在休息站过夜。所以，一段路必须在同一天中走完。

Pine希望每一天走的路长度尽可能相近，所以他希望每一天走的路的长度的方差尽可能小。

帮助Pine求出最小方差是多少。

设方差是v，可以证明，v×m^2是一个整数。为了避免精度误差，输出结果时输出v×m^2。

设第i段划分路的长度为ai，则不难推导出答案为m*sigma(a^2)-(sigma(a))^2。

维护前缀和s[i]，于是设f[i]为前i条路的最少费用，就有f[i]=min(f[i],f[j]+(s[[i]-s[j])^2)。

愉快的斜率优化维护，大致操作跟BZOJ3675 & 洛谷3648 & UOJ104：[Apio2014]序列分割一样（是的这两道题几乎相同，可能式子不一样……？）复杂度O(nm)。

……然而这道题是在去apio之后看的wuvin课件写的，于是……这其实是裸的wqs二分。

二分划分代价为c，c越大划分次数越少，反之越多（因为c减小时，显然答案在原基础上变得更小，可以选择不继续划分（当然答案还是比原来小），或者可以多花点c来划分序列，于是划分一定变多。）将划分次数逼近为m则此时即为答案。

复杂度O(nlogV)快到不知道哪里去了，这题n可以为1e5。

#include<queue>
#include<cmath>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
typedef long double dl;
const int N=;
inline int read(){
    int X=,w=;char ch=;
    while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
    while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();
    return w?-X:X;
}
int n,m,q[N],g[N];
ll s[N],f[N];
inline ll sqr(ll x){return x*x;}
inline dl suan(int k,int j){
    if(s[k]==s[j])return -1e18;
    return (dl)(f[k]-f[j]+sqr(s[k])-sqr(s[j]))/(s[k]-s[j]);
}
bool pan(ll c){
    int l=,r=;
    for(int i=;i<=n;i++){
        while(l<r&&suan(q[l],q[l+])<(dl)*s[i])l++;
        f[i]=f[q[l]]+sqr(s[i]-s[q[l]])+c;
        g[i]=g[q[l]]+;
        while(l<r&&suan(q[r-],q[r])>suan(q[r],i))r--;
        q[++r]=i;
    }
    return g[n]<=m;
}
ll solve(ll l,ll r){
    ll ans;
    while(l<r){
        ll mid=(l+r)>>;
        if(!pan(mid))l=mid+;
        else{
            ans=m*(f[n]-mid*m)-sqr(s[n]);r=mid;//这里括号里面的m改为g[n]就会wa不知道为什么有谁能解答一下吗
//主要是感觉最后反正g[n]也会变成m就乘了g[n]结果就90分了……
        }
    }
    return ans;
}
int main(){
    n=read(),m=read();
    for(int i=;i<=n;i++)s[i]=s[i-]+read();
    printf("%lld\n",solve(,sqr(s[n])));
    return ;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+欢迎访问我的博客：http://www.cnblogs.com/luyouqi233/ +

+++++++++++++++++++++++++++++++++++++++++++