BZOJ 4698: Sdoi2008 Sandy的卡片 后缀数组_RMQ_查分

BZOJ 4698: Sdoi2008 Sandy的卡片 后缀数组_RMQ_查分

Code:

#include <bits/stdc++.h>
#define setIO(s) freopen(s".in", "r" , stdin) 
#define maxn 500000
#define N 2000 
#define inf 10000  
using namespace std;
int n,m; 
int arr[maxn], num[maxn], A[maxn], tag[maxn], cnt[maxn];  
namespace SA
{
    int rk[maxn], tp[maxn], tax[maxn], sa[maxn], height[maxn], Log[maxn], dp[maxn][23];  
    void qsort()
    {
        for(int i=0;i<=m;++i) tax[i]=0;
        for(int i=1;i<=n;++i) ++tax[rk[i]]; 
        for(int i=1;i<=m;++i) tax[i]+=tax[i-1];
        for(int i=n;i>=1;--i) sa[tax[rk[tp[i]]]--]=tp[i]; 
    }
    void construct()
    {
        for(int i=1;i<=n;++i) rk[i]=arr[i],tp[i]=i;
        qsort();
        for(int k=1;k<=n;k<<=1)
        {
            int p=0;
            for(int i=n-k+1;i<=n;++i) tp[++p]=i;
            for(int i=1;i<=n;++i) if(sa[i]>k) tp[++p]=sa[i]-k;
            qsort(), swap(tp,rk), rk[sa[1]]=p=1;
            for(int i=2;i<=n;++i)
                rk[sa[i]]=(tp[sa[i-1]]==tp[sa[i]]&&tp[sa[i-1]+k]==tp[sa[i]+k])?p:++p;
            if(p==n) break;
            m=p; 
        }
        int k=0;
        for(int i=1;i<=n;++i)rk[sa[i]]=i;
        for(int i=1;i<=n;++i)
        {
            if(k)--k;
            int j=sa[rk[i]-1]; 
            while(arr[i+k]==arr[j+k])++k;
            height[rk[i]]=k; 
        }
    }
    void RMQ()
    {    
        for(int i=2;i<=n;++i) Log[i]=Log[i>>1]+1;     
        for(int i=1;i<=n;++i) dp[i][0]=height[i];
        for(int j=1;(1<<j)<=n;++j)
            for(int i=1;i+(1<<j)-1<=n;++i) 
                dp[i][j]=min(dp[i][j-1],dp[i+(1<<j-1)][j-1]); 
    }
    int getmin(int a,int b)
    { 
        int k=Log[b-a+1];     
        return min(dp[a][k],dp[b-(1<<k)+1][k]);
    } 
};
int main()
{
    // setIO("input");
    int tot;
    scanf("%d",&tot);
    for(int i=1,o,h;i<=tot;++i)
    {
        scanf("%d",&o);
        for(int j=1;j<=o;++j) scanf("%d",&A[j]);
        for(int j=2;j<=o;++j) arr[++n]=(A[j]-A[j-1]+inf),tag[n]=i; 
        arr[++n]=(inf<<1)+i; 
    }  
    m=(inf<<1)+(tot<<1), SA::construct(), SA::RMQ();           
    int ty=0,l=1,ans=0;     
    for(int i=1;i<=n;++i)
    {
        int o=SA::sa[i];              
        if(!tag[o]) continue; 
        if(!cnt[tag[o]]) ++ty; 
        ++cnt[tag[o]]; 
        for(;l<i;++l)
        {
            if(!tag[SA::sa[l]]) continue; 
            if(cnt[tag[SA::sa[l]]] > 1) 
                --cnt[tag[SA::sa[l]]];        
            else 
                break; 
        } 
        if(ty==tot) 
        { 
            ans=max(ans,SA::getmin(l+1,i) + 1);   
        }
    }
    printf("%d\n",ans);
    return 0; 
}

　　

posted @ 2019-06-10 13:27 EM-LGH 阅读( ...) 评论( ...) 编辑 收藏