[HDU3045]Picnic Cows

题面描述

传送门

思路

这道题看一看，貌似就是

小的跟小的混，大的跟大的混嘛

简单来说就是从小到大排序。

状态转移方程

$s$为前缀和

$$F_i=\min (F_j-(i-j)\ast a_{j+1}+s_i-s_j)$$

决策单调性

设

$$F_j-(i-j)\ast a_{j+1}+s_i-s_j\ge F_k-(i-k)\ast a_{k+1}+s_i-s_k(k>j)$$

在未来状态$t$下，

证明:

$$F_j-(t-j)\ast a_{j+1}-s_t+s_j\ge F_k-(t-k)\ast a_{k+1}+s_t-s_k$$

$$F_j-(t-j)\ast a_{j+1}-s_j\ge F_k-(t-k)\ast a_{k+1}-s_k$$

$$t=i+val(val>0)$$

仅需证明:

$$-val\ast a_{j+1}\ge -val\ast a_{k+1}$$

$$a_{j+1}\le a_{k+1}$$

证毕。

踢队头

$$F_j-(i-j)\ast a_{j+1}-s_j\ge F_k-(i-k)\ast a_{k+1}-s_k(k>j)$$

$$F_j-F_k+j\ast a_{j+1}-k\ast a_{k+1}-s_j+s_k\ge i\ast a_{j+1}-i\ast a_{k+1}$$

$$calc(j,k)=\frac{F_j-F_k+j\ast a_{j+1}-k\ast a_{k+1}-s_j+s_k}{a_{j+1}-a_{k+1}}\le i$$

$k$优于$j$

即

$$calc(q_{head},q_{head+1})\le i$$

$q_{head+1}$优于$q_{head}$，踢去$head$

维护一个$calc(q_{head},q_{head+1})\ge i$队列

同时我们可以发现，$calc(q_{head},q_{head+1})$随$i$的增大，$head$不断增大，而不断增大。

因此，队列维护的应该是一个下凸壳。

踢队尾

由于维护一个下凸壳，

所以

$$calc(q_{tail-1},q_{tail})\le calc(q_{tail},i-m)$$

反之，

当

$$calc(q_{tail-1},q_{tail})\ge calc(q_{tail},i-m)$$

踢去队尾。

AC code

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#define gc getchar()
#define ll long long 
using namespace std;
const int N=4e5+10;
inline void qr(ll &x)
{
    x=0;char c=gc;int f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=gc;}
    while(c>='0'&&c<='9'){x=x*10+(c^48);c=gc;}
    x*=f;
}
void qw(ll x)
{
    if(x<0)x=-x,putchar('-');
    if(x/10)qw(x/10);
    putchar(x%10+48);
}
ll f[N],s[N],a[N];int q[N],l,r,n,m;
inline ll calc(int j,int k)
{
    return f[k]-f[j]-s[k]+s[j]+k*a[k+1]-j*a[j+1];
}
inline ll sum(int j,int k)
{
    return a[k+1]-a[j+1];
}
bool check(int i,int j,int k)
{
	return calc(j,i)*sum(k,j)<=calc(k,j)*sum(j,i);
}
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        f[0]=0;
        for(int i=1;i<=n;i++)qr(a[i]);
        sort(a+1,a+n+1);s[0]=0;
        for(int i=1;i<=n;i++)s[i]=s[i-1]+a[i];
        for(int i=min((m<<1)-1,n);i>=m;i--)f[i]=s[i]-a[1]*i;
        l=1;r=1;q[1]=m;
        for(int i=m<<1;i<=n;i++)
        {
            while(l<r&&calc(q[l],q[l+1])<i*sum(q[l],q[l+1]))++l;
            f[i]=f[q[l]]+(s[i]-s[q[l]])-(i-q[l])*a[q[l]+1];
            while(l<r&&check(i-m+1,q[r],q[r-1]))--r;
            q[++r]=i-m+1;
        }
        qw(f[n]);puts("");
    }
    return 0;
}