[nowcoder5668H]Sort the Strings Revision

考虑对于$p_{i}=0$，那么可以快速比较出$s_{0},s_{1},...,s_{i-1}$与$s_{i},s_{i+1},...,s_{n}$之间的大小关系，然后对两边分别找到最小的$p_{i}$即可，用线段树维护复杂度多了一个log无法通过，因此需要用笛卡尔树来维护

笛卡尔树：https://oi-wiki.org/ds/cartesian-tree/

搜索时维护在其前面的数量，即要求出笛卡尔树的子树大小，另外需要特殊处理$d[i]=s_{0}[p[i]]$的情况，可以令$p[i]=n$，还有当两个字符串相同时编号小的优先

#include<bits/stdc++.h>
using namespace std;
#define N 2000005
#define mod 1000000007
int t,n,p[N],d[N],st[N],ls[N],rs[N],rk[N];
void dfs(int k,int l,int r,int s){
    if ((p[k]==n)||(l>=r)){
        for(int i=l;i<=r;i++)rk[i]=s+(i-l);
        return;
    }
    dfs(ls[k],l,k,s+(p[k]%10>d[k])*(r-k));
    dfs(rs[k],k+1,r,s+(p[k]%10<d[k])*(k-l+1));
}
int main(){
    scanf("%d",&t);
    while (t--){
        scanf("%d",&n);
        int p1,p2,p3,p4;
        scanf("%d%d%d%d",&p1,&p2,&p3,&p4);
        for(int i=0;i<n;i++)p[i]=i;
        for(int i=1;i<n;i++){
            swap(p[p1%(i+1)],p[i]);
            p1=(1LL*p1*p2+p3)%p4;
        }
        int d1,d2,d3,d4;
        scanf("%d%d%d%d",&d1,&d2,&d3,&d4);
        for(int i=0;i<n;i++){
            d[i]=d1%10;
            d1=(1LL*d1*d2+d3)%d4;
            if (p[i]%10==d[i])p[i]=n;
        }
        st[0]=0;
        for(int i=0;i<n;i++){
            int k=st[0];
            while ((k>0)&&(p[st[k]]>p[i]))k--;
            if (k)rs[st[k]]=i;
            if (k<st[0])ls[i]=st[k+1];
            st[++k]=i;
            st[0]=k;
        }
        dfs(st[1],0,n,0);
        int s=1,ans=0;
        for(int i=0;i<=n;i++){
            ans=(ans+1LL*s*rk[i])%mod;
            s=s*10000019LL%mod;
        }
        printf("%d
",ans);
    }
}