YbtOJ NOIP模拟 Sequence (dp+问题转移）

题意

$给⼀个⻓度为n的序列a，你可以删去任意多个数$
$最⼤化剩余序列中，数值等于位置标号的个数。$

input

5
1 1 2 5 4

output（可变为1 2 5 4）

3

subtask

对于 50% 的数据，$n\le 1e3$ ；
​对于 100% 的数据，$n\le 5e5$;
其中 $a_i\le 1e9$ ;

solution

首先发现贪心不可行，锁定dp，很容易想到一个线性dp的做法：

令状态 $f_i$ 表示处理 $1$ 到$i$ 这些数，并且处理后的序列以 $a_i$ 结尾且数字$a_i$ 移动到位置 $a_i$ 上时的最大答案。 其中由于任意$a_i$的位置只能往前移动，所以只有满足$a_i\le i$的状态才合法。

转移易得 $f_i=\max (f_j+1)$

其中转移条件为:

$$\{\begin{array}{r}j<i\\ a_j<a_i\\ j-a_j\le i-a_i\end{array}$$

条件很好理解：$f_i$由于下标$i$的递增只能从小于 $a_i$ 的$a_j$的状态转移来；而让 数字 $a_j$ 移动到位置$a_j$所需删数的数量不能超过$a_i$的，即$j$位置前删除的数不能超过$i$位置前删除的数。

初态为$d{p}_i=1(a_i\le i)$，因为对于任意一个满足$a_i\le i$的数都一定可以从位置 $i$ 移动到位置 $a_i$

那么就有了一个朴素转移的$O(n^2)$的dp,可以过掉subtask1

//50pts subtask1
#include<bits/stdc++.h>
using namespace std;
const int maxn=5e5+10;
int n,a[maxn],f[maxn],ans;
signed main()
{
    cin>>n;
    for(int i=1;i<=n;i++)
    {
    	cin>>a[i];
    	if(a[i]<=i) f[i]=1,ans=1;
    }
    for(int i=1;i<=n;i++)
        if(a[i]<=i)
            for(int j=1;j<i;j++)
                if(a[j]<a[i]&&a[i]-a[j]<=i-j)
                    f[i]=max(f[i],f[j]+1),ans=max(ans,f[i]);
    cout<<ans<<endl;
    return 0;
}

那么怎么优化呢，可以从转移条件入手

$$\{\begin{array}{r}j<i\\ a_j<a_i\\ j-a_j\le i-a_i\end{array}$$

对于这种类似三维偏序的转移，可以用二分数据结构或者CDQ分治优化，可以做到 $O(n{\log }^{2}n)$，虽然本题数据不能AC，这种优化方式还是很值得学习的。

进一步的，凭借 小学 学习的不等式，很容易得到后两个不等式是第一个的充分不必要条件，所以只需要考虑后两个不等式即可

$$\{\begin{array}{r}{a}_j<a_i\\ j-a_j\le i-a_i\end{array}$$

采用二维偏序的思想，按照 $a_i$ 为关键字排序，这样一定满足在新序列从前向后转移时，$a_j<a_i且j<i$。又由于每次转移贡献为1，所以只要在新序列的顺序中以 $(i-a_i)$ 为元素值求最长不降序列，其长度即为所求答案。

人尽皆知的，最长不降序列问题可以通过树状数组或者单调栈做到$O(n\log n)$，就能AC此题。

//100pts subtask2
#include <bits/stdc++.h>
using namespace std;
const int maxn=5e5+10;
const int INF=0x3f3f3f3f;
struct node{int p, a,val;}a[maxn];
int n,m,h[maxn];
inline bool cmp(const node &A,const node &B)
{
    if(A.a==B.a) return A.p>B.p;
    else return A.a<B.a;
}
signed main()
{
    memset(h, 0x3f, sizeof(h));
    cin>>n;int x;
    for(int i=1;i<=n;i++)
        cin>>x,a[i]={i,x,i-x};
    sort(a+1,a+n+1,cmp);
    for(int i=1;i<=n;i++)
        if(a[i].val>=0)
            h[(upper_bound(h+1,h+n+1,a[i].val)-h)]=a[i].val;
    for (int i=n;i>=0;i--)
        if(h[i]<INF)
            return cout<<i<<endl,0;
    cout<<0<<endl;
    return 0;
}