BZOJ2151 种树（贪心+堆+链表/wqs二分+动态规划）

　　dp容易想到，但没法进一步优化了。

　　考虑贪心，每次选出价值最大的物品。但这显然是不对的因为会影响其他物品的选择。

　　于是考虑加上反悔操作。每次选出一个物品后，将其相邻两物品删除，再将原物品价值变为相邻两物品价值和-原物品价值。这样如果再次选择该物品就可以达到改为选择相邻两物品的效果。并且最优方案中相邻两物品一定要么都选要么都不选，否则不如选择原物品。

　　这种带反悔的贪心策略似乎类似地在网络流算法中出现，应该是一个比较普遍的做法，然而并不会证。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
#define N 200010
int n,m,a[N],pre[N],nxt[N],ans;
bool flag[N];
struct data
{
    int i,x;
    bool operator <(const data&a) const
    {
        return x<a.x;
    }
};
priority_queue<data> q;
int main()
{
#ifndef ONLINE_JUDGE
    freopen("bzoj2151.in","r",stdin);
    freopen("bzoj2151.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),m=read();
    if (m>n/) {cout<<"Error!";return ;}
    for (int i=;i<=n;i++) a[i]=read(),q.push((data){i,a[i]});
    for (int i=;i<n;i++) nxt[i]=i+;
    for (int i=;i<=n;i++) pre[i]=i-;
    pre[]=n,nxt[n]=;
    for (int i=;i<=m;i++)
    {
        while (flag[q.top().i]) q.pop();
        data t=q.top();q.pop();
        ans+=t.x;
        int x=pre[t.i],y=nxt[t.i];
        a[t.i]=a[x]+a[y]-a[t.i];
        q.push((data){t.i,a[t.i]});
        flag[x]=flag[y]=;
        nxt[pre[x]]=t.i,pre[t.i]=pre[x];
        pre[nxt[y]]=t.i,nxt[t.i]=nxt[y];
    }
    cout<<ans;
    return ;
}

　　不过dp真的没救了吗？看到这个恰好选k个非常容易让人联想到wqs二分。于是二分多选一个数的代价，跑正常的dp，就过掉了。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
int read()
{
    int x=,f=;char c=getchar();
    while (c<''||c>'') {if (c=='-') f=-;c=getchar();}
    while (c>=''&&c<='') x=(x<<)+(x<<)+(c^),c=getchar();
    return x*f;
}
#define N 200010
#define inf 2000000010
int n,m,a[N],f[N][][],g[N][][],ans;
void chkmax(int &x,int &y,int a,int b)
{
    if (a>x||a==x&&b<y) x=a,y=b;
}
bool check(int k)
{
    memset(f,,sizeof(f)),memset(g,,sizeof(g));
    f[][][]=a[]+k,g[][][]=;
    f[][][]=f[][][]=-inf;
    for (int i=;i<=n;i++)
    {
        f[i][][]=f[i-][][]+a[i]+k,g[i][][]=g[i-][][]+;
        if (i<n) f[i][][]=f[i-][][]+a[i]+k,g[i][][]=g[i-][][]+;
        if (f[i-][][]>f[i-][][]||f[i-][][]==f[i-][][]&&g[i-][][]<g[i-][][]) f[i][][]=f[i-][][],g[i][][]=g[i-][][];
        else f[i][][]=f[i-][][],g[i][][]=g[i-][][];
        if (f[i-][][]>f[i-][][]||f[i-][][]==f[i-][][]&&g[i-][][]<g[i-][][]) f[i][][]=f[i-][][],g[i][][]=g[i-][][];
        else f[i][][]=f[i-][][],g[i][][]=g[i-][][];
    }
    int x=f[n][][],y=g[n][][];
    chkmax(x,y,f[n][][],g[n][][]);
    chkmax(x,y,f[n][][],g[n][][]);
    chkmax(x,y,f[n][][],g[n][][]);
    if (y<=m) {ans=x-m*k;return ;}
    else return ;
}
int main()
{
#ifndef ONLINE_JUDGE
    freopen("force.in","r",stdin);
    freopen("force.out","w",stdout);
    const char LL[]="%I64d\n";
#else
    const char LL[]="%lld\n";
#endif
    n=read(),m=read();
    if (m>n/) {cout<<"Error!";return ;}
    for (int i=;i<=n;i++) a[i]=read();
    int l=-,r=;
    while (l<=r)
    {
        int mid=(l+r)/;
        if (check(mid)) l=mid+;
        else r=mid-;
    }
    cout<<ans;
    return ;
}