1049: [HAOI2006]数字序列
Time Limit: 10 Sec Memory Limit: 162 MBSubmit: 397 Solved: 151
[ Submit][ Status]
Description
现在我们有一个长度为n的整数序列A。但是它太不好看了,于是我们希望把它变成一个单调严格上升的序列。但是不希望改变过多的数,也不希望改变的幅度太大。
Input
第一行包含一个数n,接下来n个整数按顺序描述每一项的键值。
Output
第一行一个整数表示最少需要改变多少个数。 第二行一个整数,表示在改变的数最少的情况下,每个数改变的绝对值之和的最小值。
Sample Input
4
5 2 3 5
5 2 3 5
Sample Output
1
4
【数据范围】
90%的数据n<=6000。
100%的数据n<=35000。
保证所有数列是随机的。
4
【数据范围】
90%的数据n<=6000。
100%的数据n<=35000。
保证所有数列是随机的。
HINT
Source
首先,我们先把原序列-a[i]
原题转变为求单调不降序列。
原因:假象原序列的答案连成一条斜率为1的线,那么-完刚好变成平行x轴直线,同理反推,定理得证。
数字序列定律:
[i,j]中间修改的数要么为a[i],要么为a[j],否则将答案改成它不会更劣
于是就改。。。
从昨天做到今天。。
树链剖分2h,LIS 2 day。。。
再一次证明思路的重要性。。。
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<functional>
#include<iostream>
#include<cmath>
#include<cctype>
#include<ctime>
using namespace std;
#define For(i,n) for(int i=1;i<=n;i++)
#define Fork(i,k,n) for(int i=k;i<=n;i++)
#define Rep(i,n) for(int i=0;i<n;i++)
#define ForD(i,n) for(int i=n;i;i--)
#define RepD(i,n) for(int i=n;i>=0;i--)
#define Forp(x) for(int p=pre[x];p;p=next[p])
#define Lson (x<<1)
#define Rson ((x<<1)+1)
#define MEM(a) memset(a,0,sizeof(a));
#define MEMI(a) memset(a,127,sizeof(a));
#define MEMi(a) memset(a,128,sizeof(a));
#define INF (2139062143)
#define F (100000007)
#define MAXN (35000+10)
long long mul(long long a,long long b){return (a*b)%F;}
long long add(long long a,long long b){return (a+b)%F;}
long long sub(long long a,long long b){return (a-b+(a-b)/F*F+F)%F;}
typedef long long ll;
ll n,a[MAXN],f[MAXN],d[MAXN],len;
void LIS()
{
memset(d,127,sizeof(d));
d[len=0]=-INF;
For(i,n)
{
int l=0,r=len,ans=0;
while (l<=r)
{
int m=l+r>>1;
if (d[m]<=a[i]) ans=m,l=m+1;
else r=m-1;
}
f[i]=ans+1;len=max(len,f[i]);
d[f[i]]=min(d[f[i]],a[i]);
}
}
ll g[MAXN],y[MAXN];
int next[MAXN]={0};
int main()
{
// freopen("bzoj1049.in","r",stdin);
scanf("%d",&n);
a[1]=-INF,a[n+2]=INF,n+=2;
Fork(i,2,n-1) scanf("%d",&a[i]),a[i]=a[i]-i;
LIS();
printf("%d\n",n-len);
memset(d,0,sizeof(d));
For(i,n) if (!d[f[i]]) d[f[i]]=i;else next[i]=d[f[i]],d[f[i]]=i;
memset(g,127,sizeof(g));
g[0]=g[1]=0;//Fork(i,2,n) g[i]=g[i-1]+abs(a[1]-a[i]);
For(i,n)
{
if (f[i]==1) continue;
int j=d[f[i]-1];
y[i]=0;
ForD(k,i-1) y[k]=y[k+1]+abs(a[i]-a[k]);
for(;j;j=next[j])
if (a[j]<=a[i])
{
if (j+1==i) {g[i]=min(g[i],g[j]);continue;}
ll tot=0;
Fork(k,j,i-1)
{
tot+=abs(a[j]-a[k]);
g[i]=min(g[i],g[j]+tot+y[k+1]);
}
}
}
cout<<g[n]<<endl;
return 0;
}