bzoj1049:[HAOI2006]数字序列

传送门

第一问是sb题,求\(a[i]\)严格上升转化为求\(a[i]-i\)的不下降就好了,答案就是\(n-len\)
第二问好神仙啊,首先也没想到\(O(n^3)\)能过\(35000\)
结论也很nb,对于区间\([l,r]\),要使不降最小代价一定是存在一个分割点\(k\),使得前一半都等于\(a\),后一半为\(b\)\(a<=b\)
然后考虑到要在满足第一问的条件下,那么对于\(l,r\)的枚举就满足\(f[l]+1==f[r]且b[l]<=b[r]\)\(f[i]\)为以\(i\)结尾的最长不下降子序列的长度,\(b[i]=a[i]-i\)
然后还存在坑点
1、直接枚举会TLE(毕竟还是\(O(n^3)\)嘛),可以用vector存下每个\(f[i]\)对应的位置
2、枚举要到\(n+1\),因为枚举到\(n\)的时候,\(f[i]\)\(n\)相同的可能还没保持不降
3、记得要将\(0\)也算进去,就是保证\(f[i]==1\)可以转移
代码:

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>
using namespace std;
void read(int &x) {
    char ch; bool ok;
    for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch=='-') ok=1;
    for(x=0; isdigit(ch); x=x*10+ch-'0',ch=getchar()); if(ok) x=-x;
}
#define rg register
const int maxn=5e4+10;vector<int>s[maxn];
int n,a[maxn],b[maxn],ans1,g[maxn],len[maxn];
long long suml[maxn],sumr[maxn],f[maxn];
int main()
{
    read(n);
    for(rg int i=1;i<=n;i++)read(a[i]),b[i]=a[i]-i;
    for(rg int i=1;i<=n;i++)
    {
        if(!ans1)g[++ans1]=b[i],s[ans1].push_back(i),len[i]=1;
        else 
        {
            int t=upper_bound(g+1,g+ans1+1,b[i])-g;
            if(t>ans1)g[++ans1]=b[i];
            else g[t]=b[i];len[i]=t;
            s[t].push_back(i);
        }
    }
    printf("%d\n",n-ans1);
    memset(f,63,sizeof f);f[1]=f[0]=0,s[0].push_back(0),b[n+1]=2e9,len[n+1]=ans1+1;b[0]=-1e9;
    for(rg int i=2;i<=n+1;i++)
    {
        int t=s[len[i]-1].size();
        for(rg int j=0;j<t;j++)
        {
            int d=s[len[i]-1][j];if(d>i)break;if(b[d]>b[i])continue;
            for(rg int k=d+1;k<=i;k++)suml[k]=suml[k-1]+abs(b[k]-b[d]);
            for(rg int k=i-1;k>=d;k--)sumr[k]=sumr[k+1]+abs(b[k]-b[i]);
            for(rg int k=d;k<i;k++)f[i]=min(f[i],f[d]+suml[k]+sumr[k+1]);
            for(rg int k=d;k<=i;k++)suml[k]=sumr[k]=0;
        }
    }
    printf("%lld\n",f[n+1]);
}

猜你喜欢

转载自www.cnblogs.com/lcxer/p/10589068.html