标签:

frcode 2 0 6 0 6 unitedstatesofamerica 3 0 6 0 12 0 21 myxophytamyxopodnabnabbednabbingnabit 6 0 9 9 16 16 19 19 25 25 32 32 37
14 12 42 31 43 40
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
#define min(a,b) (a>b?b:a)
#define max(a,b) (a>b?a:b)
using namespace std;
char str[103030];
int sa[103030],Rank[103030],rank2[103030],height[103030],c[103030],*x,*y,len;
int n;
void cmp(int n,int sz)
{
int i;
memset(c,0,sizeof(c));
for(i=0;i<n;i++)
c[x[y[i]]]++;
for(i=1;i<sz;i++)
c[i]+=c[i-1];
for(i=n-1;i>=0;i--)
sa[--c[x[y[i]]]]=y[i];
}
void build_sa(char *s,int n,int sz)
{
x=Rank,y=rank2;
int i,j;
for(i=0;i<n;i++)
x[i]=s[i],y[i]=i;
cmp(n,sz);
int len;
for(len=1;len<n;len<<=1)
{
int yid=0;
for(i=n-len;i<n;i++)
{
y[yid++]=i;
}
for(i=0;i<n;i++)
if(sa[i]>=len)
y[yid++]=sa[i]-len;
cmp(n,sz);
swap(x,y);
x[sa[0]]=yid=0;
for(i=1;i<n;i++)
{
if(y[sa[i-1]]==y[sa[i]]&&sa[i-1]+len<n&&sa[i]+len<n&&y[sa[i-1]+len]==y[sa[i]+len])
x[sa[i]]=yid;
else
x[sa[i]]=++yid;
}
sz=yid+1;
if(sz>=n)
break;
}
for(i=0;i<n;i++)
Rank[i]=x[i];
}
void getHeight(char *s,int n)
{
int k=0;
for(int i=0;i<n;i++)
{
if(Rank[i]==0)
continue;
k=max(0,k-1);
int j=sa[Rank[i]-1];
while(s[i+k]==s[j+k])
k++;
height[Rank[i]]=k;
}
}
int minv[103010][20],lg[103030];
void init_lg()
{
int i;
lg[1]=0;
for(i=2;i<102020;i++)
{
lg[i]=lg[i>>1]+1;
}
}
void init_RMQ(int n)
{
int i,j,k;
for(i=1;i<=n;i++)
{
minv[i][0]=height[i];
}
for(j=1;j<=lg[n];j++)
{
for(k=0;k+(1<<j)-1<=n;k++)
{
minv[k][j]=min(minv[k][j-1],minv[k+(1<<(j-1))][j-1]);
}
}
}
int lcp(int l,int r)
{
if(l==r)
{
return len-l;
}
l=Rank[l];
r=Rank[r];
if(l>r)
swap(l,r);
l++;
int k=lg[r-l+1];
return min(minv[l][k],minv[r-(1<<k)+1][k]);
}
__int64 fun(__int64 x)
{
__int64 ans=0;
if(x==0)
return 1;
while(x)
{
x/=10;
ans++;
}
return ans;
}
int main()
{
while(scanf("%s",str)!=EOF)
{
int t;
len=strlen(str);
build_sa(str,len+1,256);
getHeight(str,len);
scanf("%d",&t);
init_lg();
init_RMQ(len);
__int64 ans1,ans2;
ans1=ans2=0;
int l,r;
scanf("%d%d",&l,&r);
ans1+=(r-l)+1;
t--;
ans2+=(r-l)+3;
while(t--)
{
int a,b;
scanf("%d%d",&a,&b);
ans1+=(b-a)+1;
__int64 p;
p=min(lcp(l,a),min((b-a),(r-l)));
ans2+=(b-a)-p+2+fun(p);
l=a;
r=b;
}
printf("%I64d %I64d\n",ans1,ans2);
}
return 0;
}版权声明:本文为博主原创文章,未经博主允许不得转载。
HDOJ 题目4691 Front compression(后缀数组+RMQ最长前缀)
标签:
原文地址:http://blog.csdn.net/yu_ch_sh/article/details/48086521