给出一个字符串,对于每次的k,求所有子串(去除重复后的)中排名为k的L和R,其中如果有多组解,输出L最小的。
后缀数组处理过后,sum[i]表示 sa[1] 到 sa[i] 一共有多少个非重复字串,在sum数组中二分求出k的lower_bound=pos,那么第一次出现这个子串的后缀就是sa[pos],但是有可能有多组解,所以还要往下找所有height[i]>=目标字符串长度 的所有后缀中sa值最小的。其实这步应该用rmq。。。
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define maxn 1000001
typedef long long ll;
int wa[maxn],wb[maxn],wv[maxn],ws[maxn];
int sa[maxn];
char str[maxn];
int r[maxn];
ll sum[maxn];
int cmp(int *r,int a,int b,int l)
{return r[a]==r[b]&&r[a+l]==r[b+l];}
void da(int *r,int *sa,int n,int m)
{
int i,j,p,*x=wa,*y=wb,*t;
for(i=0;i<m;i++) ws[i]=0;
for(i=0;i<n;i++) ws[x[i]=r[i]]++;
for(i=1;i<m;i++) ws[i]+=ws[i-1];
for(i=n-1;i>=0;i--) sa[--ws[x[i]]]=i;
for(j=1,p=1;p<n;j*=2,m=p)
{
for(p=0,i=n-j;i<n;i++) y[p++]=i;
for(i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;
for(i=0;i<n;i++) wv[i]=x[y[i]];
for(i=0;i<m;i++) ws[i]=0;
for(i=0;i<n;i++) ws[wv[i]]++;
for(i=1;i<m;i++) ws[i]+=ws[i-1];
for(i=n-1;i>=0;i--) sa[--ws[wv[i]]]=y[i];
for(t=x,x=y,y=t,p=1,x[sa[0]]=0,i=1;i<n;i++)
x[sa[i]]=cmp(y,sa[i-1],sa[i],j)?p-1:p++;
}
return;
}
int rank[maxn],height[maxn];
void calheight(int *r,int *sa,int n)
{
int i,j,k=0;
for(i=1;i<=n;i++) rank[sa[i]]=i;
for(i=0;i<n;height[rank[i++]]=k)
for(k?k--:0,j=sa[rank[i]-1];r[i+k]==r[j+k];k++);
return;
}
int RMQ[maxn];
int mm[maxn];
int best[20][maxn];
void initRMQ(int n)
{
int i,j,a,b;
for(mm[0]=-1,i=1;i<=n;i++)
mm[i]=((i&(i-1))==0)?mm[i-1]+1:mm[i-1];
for(i=1;i<=n;i++) best[0][i]=i;
for(i=1;i<=mm[n];i++)
for(j=1;j<=n+1-(1<<i);j++)
{
a=best[i-1][j];
b=best[i-1][j+(1<<(i-1))];
if(RMQ[a]<RMQ[b]) best[i][j]=a;
else best[i][j]=b;
}
return;
}
int askRMQ(int a,int b)
{
int t;
t=mm[b-a+1];b-=(1<<t)-1;
a=best[t][a];b=best[t][b];
return RMQ[a]<RMQ[b]?a:b;
}
int lcp(int a,int b)
{
int t;
a=rank[a];b=rank[b];
if(a>b) {t=a;a=b;b=t;}
return(height[askRMQ(a+1,b)]);
}
int main()
{
while(scanf("%s",str)!=EOF)
{
int len,i;
ll q,x;
len=strlen(str);
for(i=0;i<len;i++) r[i]=str[i]-'a'+2;
r[len]=0;
da(r,sa,len+1,30);
calheight(r,sa,len);
scanf("%I64d",&q);
memset(sum,0,sizeof(sum));
for(i=1;i<=len;i++) sum[i]=len-sa[i]-height[i];
for(i=2;i<=len;i++) sum[i]+=sum[i-1];
int l=0,r=0;
while(q--)
{
scanf("%I64d",&x);
ll k=(l^r^x)+1;
int pos=lower_bound(sum+1,sum+len+1,k)-sum;
if(pos>len)
{
l=r=0;
printf("%d %d\n",l,r);
continue;
}
l=sa[pos];
r=len-(sum[pos]-k+1);
int ans=r-l+1;
pos++;
while(pos<=len&&height[pos]>=ans)
{
l=min(l,sa[pos]);
pos++;
}
r=l+ans-1;
l++;r++;
printf("%d %d\n",l,r);
}
}
return 0;
}