BZOJ 1109 [POI2007]堆积木Klo

题目链接

https://www.lydsy.com/JudgeOnline/problem.php?id=1109

思路

神奇的dp转化：

设$f_i$表示第$i$块正确的时候，第$i$块前面有多少块摆放正确。

朴素dp方程：

$$f_i=f_{j}+1[i>j(1),a_i>a_j(2),a_i-i\le a_j-j(3)]$$
很容易看出，若 $(2)(3)$两条件满足， $(1)$条件一定满足。

若按照$a_i$排序，就变成了对$a_i-i$求最长不上升子序列。

好神奇啊。

代码

#include <cstdio>
#include <algorithm>

const int maxn=100000;
const int inf=0x3f3f3f3f;

int read()
{
  int x=0,f=1;
  char ch=getchar();
  while((ch<'0')||(ch>'9'))
    {
      if(ch=='-')
        {
          f=-f;
        }
      ch=getchar();
    }
  while((ch>='0')&&(ch<='9'))
    {
      x=x*10+ch-'0';
      ch=getchar();
    }
  return x*f;
}

struct pps
{
  int x,y;

  bool operator <(const pps &other) const
  {
    if(x==other.x)
      {
        return y>other.y;
      }
    return x<other.x;
  }
};

pps p[maxn+10];
int n,len,f[maxn+10],tot;

int main()
{
  n=read();
  for(register int i=1; i<=n; ++i)
    {
      int a=read();
      if(a<=i)
        {
          p[++tot].x=a;
          p[tot].y=i-a;
        }
    }
  if(!tot)
    {
      puts("0");
      return 0;
    }
  std::sort(p+1,p+tot+1);
  f[1]=p[1].y;
  len=1;
  for(register int i=2; i<=tot; ++i)
    {
      int pos=std::upper_bound(f+1,f+len+1,p[i].y)-f;
      f[pos]=p[i].y;
      if(pos>len)
        {
          len=pos;
        }
    }
  printf("%d\n",len);
  return 0;
}