bzoj1049 [HAOI2006]数字序列（dp+BIT+贪心）

这题实在是太神了orz

首先考虑第一问。使得修改元素最少，我们不妨补集转化，改为求：最多多少个点可以不修改。
对于$j<=i$，如果$a_j,a_i$不需要修改，则必须满足$a_i-a_j>=i-j$，也就是$a_i-i>=a_j-j$，我们令$b_i=a_i-i$，则就是要求$b_i>=b_j$，也即求序列b的最长不下降子序列。答案就是n-序列b的最长不下降子序列。

考虑第二问，使得修改幅度最小。不会qaq，请看神犇讲解：portal

复杂度大概最坏$O(n^3)$，但是数据随机，就水过去了qaq

细节较多，给a[0]=-inf,a[n+1]=inf会好一些。

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 35010
inline char gc(){
    static char buf[1<<16],*S,*T;
    if(S==T){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return  EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=gc();}
    while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
    return x*f;
}
int n,a[N],f[N],aa[N],m,c[N];
ll g[N],s1[N],s2[N];
vector<int>b[N];
inline void add(int x,int val){
    for(;x<=m;x+=x&-x) c[x]=max(c[x],val);
}
inline int ask(int x){
    int res=0;for(;x;x-=x&-x) res=max(res,c[x]);return res;
}
int main(){
//  freopen("a.in","r",stdin);
    n=read();
    for(int i=1;i<=n;++i) a[i]=aa[i]=read()-i;++n;a[n]=aa[n]=inf;
    sort(aa+1,aa+n+1);m=unique(aa+1,aa+n+1)-aa-1;
    for(int i=1;i<=n;++i){
        int id=lower_bound(aa+1,aa+m+1,a[i])-aa;
        f[i]=ask(id)+1;add(id,f[i]);b[f[i]].push_back(i);g[i]=1LL<<60;
    }b[0].push_back(0);a[0]=-inf;
    for(int i=1;i<=n;++i){
        for(int ii=0;ii<b[f[i]-1].size();++ii){
            int j=b[f[i]-1][ii];if(j>i) break;
            if(a[j]>a[i]) continue;s1[j]=s2[j]=0;
            for(int k=j+1;k<=i;++k)
                s1[k]=s1[k-1]+abs(a[k]-a[j]),
                s2[k]=s2[k-1]+abs(a[k]-a[i]);
            for(int k=j;k<i;++k)
                g[i]=min(g[i],g[j]+s1[k]-s1[j]+s2[i]-s2[k]);    
        }
    }printf("%d\n%lld\n",n-f[n],g[n]);
    return 0;
}