[luogu4331]数字序列

令$a'_{i}=a_{i}+n-i$、$b'_{i}=b_{i}+n-i$，代价仍然是$\sum_{i=1}^{n}|a'_{i}-b'_{i}|$，但条件变为了$b'_{i}\le b'_{i+1}$，即不下降（以下为了方便，$a'_{i}$和$b'_{i}$仍然用$a_{i}$和$b_{i}$表示，原来的不需要考虑）

考虑暴力dp，令$f_{i,j}$表示前$i$个数且$b_{i}=j$的最小的代价，转移时先令$f_{i-1,j}=\min_{k\le j}f_{i-1,k}$，再加上绝对值，即$f_{i,j}=f_{i-1,j}+|a_{i}-j|$

由于过程是交替进行的，可以看作先加上绝对值再取min（体现在定义上就是$f_{i,j}$表示前$i$个数且$b_{i}\le j$的最小的代价），细节上由于第1次全部都是0本身就不需要取min，然后最后答案即为$f_{n,\infty}$

归纳$f_{i}$具有以下性质：其斜率单调不递增且小于等于0（即下凸但没有上升的部分）

在这一条件下，对于取绝对值再取min的过程，可以看作：对于$a_{i}$左半部分相当于斜率增加1，对于$a_{i}$右半部分斜率减少1，特别的，由于要取min，因此对于斜率为0的部分仍然不变

（做法上有一点类似[AGC049E]）

定义$f'_{i,j}=f_{i,j}-f_{i,j-1}$（$j\ge 1$），特别的$f'_{i,0}=\sum_{j=1}^{i}a_{i}$，那么$f_{i,j}=\sum_{k=0}^{j}f'_{i,k}$（其实后面两点不重要，因为求答案肯定通过构造出来的方案求更方便QAQ）

然后对于$f'_{i,j}$（$1\le j\le a_{i}$）区间减1，对于$f'_{i,j}$（$a_{i}<j$且$f'_{i,j}<0$）区间加1，直接用线段树维护即可

关于最优解的构造：先有$b_{n}=\min_{f'_{n,k}=0}k$，然后再令$b_{i}=\min(b_{i+1},\min_{f'_{i,k}=0}k-1)$即可（很明显$b_{i+1}$之前$f_{i,j}$最小的位置就是这里，同时这里的$f_{i,j}$不会是跟$f_{i,j-1}$取min的结果）

#include<bits/stdc++.h>
using namespace std;
#define N 1000005
#define L (k<<1)
#define R (L+1)
#define mid (l+r>>1) 
map<int,int>mat;
map<int,int>::iterator it;
int n,m,a[N],b[N],val[N],f[N<<2],tag[N<<2];
long long ans;
char ch[21];
int read(){
    int x=0;
    char c=getchar();
    while ((c<'0')||(c>'9'))c=getchar();
    while (('0'<=c)&&(c<='9')){
        x=x*10+c-'0';
        c=getchar();
    }
    return x;
}
void write(int k){
    int s=0;
    while (k){
        ch[++s]=k%10+'0';
        k/=10;
    }
    while (s)putchar(ch[s--]);
    putchar(' ');
}
void upd(int k,int x){
    tag[k]+=x;
    f[k]+=x;
}
void down(int k){
    upd(L,tag[k]);
    upd(R,tag[k]);
    tag[k]=0;
}
void update(int k,int l,int r,int x,int y,int z){
    if ((l>y)||(x>r))return;
    if ((x<=l)&&(r<=y)){
        upd(k,z);
        return;
    }
    down(k);
    update(L,l,mid,x,y,z);
    update(R,mid+1,r,x,y,z);
    f[k]=min(f[L],f[R]);
}
int query(int k,int l,int r){
    if (l==r)return l;
    down(k);
    if (f[R])return query(R,mid+1,r);
    return query(L,l,mid);
}
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++)a[i]=read()+n-i;
    memcpy(val,a,sizeof(val));
    sort(val+1,val+n+1);
    m=unique(val+1,val+n+1)-val-1;
    for(int i=1;i<=n;i++)a[i]=lower_bound(val+1,val+m+1,a[i])-val;
    for(int i=1;i<=n;i++){
        update(1,1,m,1,a[i],-1);
        if (i>1)update(1,1,m,a[i]+1,query(1,1,m),1);
        b[i]=query(1,1,m);
    }
    for(int i=n-1;i;i--)b[i]=min(b[i],b[i+1]);
    for(int i=1;i<=n;i++)ans+=abs(val[a[i]]-val[b[i]]);
    printf("%lld\n",ans);
    for(int i=1;i<=n;i++)write(val[b[i]]-n+i);
}