hdu 4691 Front compression (后缀数组+rmq)

最新推荐文章于 2019-02-19 18:03:23 发布

My_ACM_Dream

最新推荐文章于 2019-02-19 18:03:23 发布

阅读量597

点赞数

CC 4.0 BY-SA版权

分类专栏：字符串字符串—后缀数组

本文链接：https://blog.youkuaiyun.com/My_ACM_Dream/article/details/44157115

字符串同时被 2 个专栏收录

72 篇文章

订阅专栏

字符串—后缀数组

16 篇文章

订阅专栏

本文深入探讨了如何面对看似简单的题目，却因表述晦涩导致理解困难的问题。通过实例分析，揭示了解题的关键在于细致阅读、准确把握题意，并提供了有效的策略帮助读者克服这一障碍。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

明知道水题，就是读不懂题目。

#include<iostream>
#include<math.h>
#include<stdio.h>
#include<algorithm>
#include<string.h>
#include<vector>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define B(x) (1<<(x))
typedef long long ll;
const int oo=0x3f3f3f3f;
const ll OO=1LL<<61;
const int MOD=10007;
const int maxn=100000+100;
int Rank[maxn],SA[maxn],height[maxn];
int t1[maxn],t2[maxn],t3[maxn],t4[maxn];
char str[maxn];
int len;

void Swap(int*& x,int*& y){

    int *temp=x;
    x=y;
    y=temp;
}

int cmp(int t[],int a,int b,int l){

    return t[a]==t[b]&&t[a+l]==t[b+l];
}

void build_SA(char s[],int len,int up){

    int *k1=t1,*k2=t2,*r=t3,*cnt=t4;
    for(int i=0;i<up;i++)cnt[i]=0;
    for(int i=0;i<len;i++)cnt[k1[i]=s[i]]++;
    for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
    for(int i=len-1;i>=0;i--)SA[--cnt[k1[i]]]=i;

    for(int d=1,p=1;p<len;d<<=1,up=p){

        p=0;
        for(int i=len-d;i<len;i++)k2[p++]=i;
        for(int i=0;i<len;i++)if(SA[i]>=d)k2[p++]=SA[i]-d;
        for(int i=0;i<len;i++)r[i]=k1[k2[i]];

        for(int i=0;i<up;i++)cnt[i]=0;
        for(int i=0;i<len;i++)cnt[r[i]]++;
        for(int i=1;i<up;i++)cnt[i]+=cnt[i-1];
        for(int i=len-1;i>=0;i--)SA[--cnt[r[i]]]=k2[i];

        Swap(k1,k2);
        k1[SA[0]]=0;
        p=1;
        for(int i=1;i<len;i++){
            k1[SA[i]]= cmp(k2,SA[i-1],SA[i],d) ? p-1 : p++;
        }
    }
}

void get_height(char s[],int len){

    for(int i=1;i<=len;i++)Rank[SA[i]]=i;
    for(int i=0,p=0;i<len;i++){

        int j=SA[Rank[i]-1];
        while(s[i+p]==s[j+p])p++;
        height[Rank[i]]=p;
        if(p)p--;
    }
}

int dp[maxn][20];
void rmq_init(int n){

    int m=floor(log(n+0.0)/log(2.0));
    for(int i=1;i<=n;i++)dp[i][0]=height[i];
    for(int j=1;j<=m;j++){
        for(int i=n;i;i--){
            dp[i][j]=dp[i][j-1];
            if(i+B(j-1)<=n){
                dp[i][j]=min(dp[i][j],dp[i+B(j-1)][j-1]);
            }
        }
    }
}

int rmq_query(int l,int r){

    if(l==r) return len-l;
    int a=Rank[l],b=Rank[r];
    if(a>b)swap(a,b);
    a++;
    int k=floor(log(b-a+1.0)/log(2.0));
    return min(dp[a][k],dp[b-B(k)+1][k]);
}

int bit(int n){
    int i=0;
    if(n==0)i=1;
    while(n){
        i++;
        n/=10;
    }
    return i;
}

int main(){

    int n;
    int a,b,l,r;
    ll ans1,ans2;
    while(scanf("%s",str)!=EOF){

        len=strlen(str);
        build_SA(str,len+1,130);
        get_height(str,len);
        rmq_init(len);

        ans1=ans2=0;
        scanf("%d",&n);
        if(n){
            scanf("%d %d",&a,&b);
            ans1+=(b-a+1);
            ans2+=(b-a+3);
        }
        while(--n){
            scanf("%d %d",&l,&r);
            int lcp=rmq_query(a,l);
            lcp=min(lcp,min(r-l,b-a));
            ans1+=(r-l+1);
            ans2+=(r-l)-lcp+2+bit(lcp);
            a=l;
            b=r;
        }
        printf("%I64d %I64d\n",ans1,ans2);
    }
    return 0;
}
/**
3
abacabab
abababacb
aaaabacababc
*/