[BZOJ2253][2010BJWC]纸箱堆叠（dp+cdq分治）

题目：

我是超链接

题解：

这个东西$O(n^2)$的最长上升子序列嘛
不过过不了啊，回忆一下求最长上升子序列的过程，我们会发现其实就是枚举n，然后寻找三项都小于ta的延续个数最长的
而这是三个元素，三维偏序cdq分治没跑了
其实我们只需要修改一下树状数组，用于求i < q[i].c的连续个数最长是几就行啦
并且因为前面的对后面的一定是有影响的，所以我们并不能变换顺序，必须是：处理左区间，处理左区间对右区间的影响，处理右区间，这样的步骤，不然对右区间的修改会有影响
由于此时是严格上升，我们还需要对a这一维进行特殊的处理才行，比如说找到ta可以延伸的最左点什么的

代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define LL long long
using namespace std;
const int N=50005;
struct hh{LL a,b,c;int id;}q[N],jl[N];
LL a,mod;int c[N*3],f[N],n,has[N*3],s,L[N],R[N];
int cmp(hh a,hh b){return ((a.a<b.a) || (a.a==b.a&&a.b<b.b) || (a.a==b.a && a.b==b.b && a.c<b.c));}
void add(int loc,int v)
{
    for (int i=loc;i<=s;i+=i&(-i)) c[i]=max(c[i],v);
}
int qurry(int loc)
{
    int ans=0;
    for (int i=loc;i>=1;i-=i&(-i)) ans=max(ans,c[i]);
    return ans;
}
void cover(int loc){for (int i=loc;i<=s;i+=i&(-i)) c[i]=0;}
int cmpb(hh a,hh b){return a.b<b.b;}
void cdq(int l,int r)
{
    if (q[l].a==q[r].a) return;
    int mid=(l+r)>>1;
    if (L[mid]!=l) mid=L[mid]-1;
    else mid=R[mid];
    cdq(l,mid); 
    sort(q+l,q+mid+1,cmpb);
    for (int i=mid+1;i<=r;i++) jl[i]=q[i];
    sort(jl+mid+1,jl+r+1,cmpb); 
    int t1=l,t2=mid+1;
    for (int i=l;i<=r;i++)
      if ((t1<=mid && q[t1].b<jl[t2].b) || t2>r) add(q[t1].c,f[q[t1].id]),t1++;
      else f[jl[t2].id]=max(f[jl[t2].id],qurry(jl[t2].c-1)+1),t2++;
    for (int i=l;i<=mid;i++) cover(q[i].c);
    cdq(mid+1,r);
}
int main()
{
    scanf("%lld%lld%d",&a,&mod,&n);
    LL b=1;
    for (int i=1;i<=n;i++)
    {
        b=b*a%mod; q[i].a=b;
        b=b*a%mod; q[i].b=b;
        b=b*a%mod; q[i].c=b;
        if (q[i].b>q[i].a) swap(q[i].b,q[i].a);
        if (q[i].c>q[i].a) swap(q[i].c,q[i].a);
        if (q[i].c>q[i].b) swap(q[i].b,q[i].c);
        has[++s]=q[i].a; has[++s]=q[i].b; has[++s]=q[i].c;
    }
    sort(has+1,has+s+1);
    s=unique(has+1,has+s+1)-has-1;
    for (int i=1;i<=n;i++) 
      q[i].a=lower_bound(has+1,has+s+1,q[i].a)-has,q[i].b=lower_bound(has+1,has+s+1,q[i].b)-has,q[i].c=lower_bound(has+1,has+s+1,q[i].c)-has;
    sort(q+1,q+n+1,cmp);
    for (int i=1;i<=n;i++) q[i].id=i,f[i]=1;

    for (int i=1;i<=n;i++)
      if (i==1 || (q[i].a!=q[i-1].a)) L[i]=i;
      else L[i]=L[i-1];

    for (int i=n;i>=1;i--)
      if (i==n || (q[i].a!=q[i+1].a)) R[i]=i;
      else R[i]=R[i+1];
    cdq(1,n);
    int ans=0;
    for (int i=1;i<=n;i++) ans=max(ans,f[i]);
    printf("%d",ans);
}