HDU4691 Front compression【后缀数组】

题意：先给出一个串，然后下面给一个n表示输入串的个数，然后每行两个数表示输入的串的起点和重点（左闭右开）。。问的是 不压缩和压缩之后的大小。。

每个字符占1大小，空格和回车也占1.。

很基本的后缀数组，求lcp，因为不会写RMQ。。。所以直接写线段树了。。每次求一下这个串与上一个串的lcp。。然后处理就好了

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ls (k<<1)
#define rs (k<<1|1)
using namespace std;
const int MAXN=100010;
typedef long long ll;
int wa[MAXN],wb[MAXN],ws[MAXN],wv[MAXN];
bool 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=1,*x=wa,*y=wb;
	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<n;j<<=1)
	{
		p=0;
		for(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<m;i++)
			ws[i]=0;
		for(i=0;i<n;i++)
			wv[i]=x[y[i]];
		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];
		swap(x,y);
		p=1;
		x[sa[0]]=0;
		for(i=1;i<n;i++)
			x[sa[i]]=cmp(y,sa[i],sa[i-1],j)?p-1:p++;
		m=p;
	}
}
int height[MAXN],RANK[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;i++)
	{
		if(k)
			k--;
		j=sa[RANK[i]-1];
		while(r[i+k]==r[j+k])
			k++;
		height[RANK[i]]=k;
	}
}
struct TREE
{
	int l,r;
	int val;
	int mid()
	{
		return (l+r)>>1;
	}
}tree[MAXN<<2];
void pushup(int k)
{
	tree[k].val=min(tree[ls].val,tree[rs].val);
}
void build(int l,int r,int k)
{
	tree[k].l=l;
	tree[k].r=r;
	if(l==r)
	{
		tree[k].val=height[l];
		return;
	}
	int mid=tree[k].mid();
	build(l,mid,ls);
	build(mid+1,r,rs);
	pushup(k);
}
int query(int l,int r,int k)
{
	if(tree[k].l>=l&&tree[k].r<=r)
		return tree[k].val;
	int mid=tree[k].mid();
	if(r<=mid)
		return query(l,r,ls);
	else if(l>mid)
		return query(l,r,rs);
	else
	{
		return min(query(l,mid,ls),query(mid+1,r,rs));
	}
}
int Query(int l,int r)
{
	if(l==r)
	{
		return 99999999;
	}
	int u=RANK[l];
	int v=RANK[r];
	if(u>v)
		swap(u,v);
	return query(u+1,v,1);
}
char str[MAXN];
int r[MAXN],sa[MAXN];
int main()
{
	int i;
	while(scanf("%s",str)==1)
	{
		int q;
		int n=strlen(str);
		for(i=0;i<n;i++)
			r[i]=str[i];
		r[n]=0;
		da(r,sa,n+1,128);
		calheight(r,sa,n);
		build(1,n,1);
		scanf("%d",&q);
		ll ans1=0,ans2=0;
		ans1+=q;
		ans2+=q+2;
		int a,b,aa,bb;
		q--;
		scanf("%d%d",&aa,&bb);
		ans1+=(ll)bb-aa;
		ans2+=(ll)bb-aa;
		while(q--)
		{
			scanf("%d%d",&a,&b);
			int temp=Query(a,aa);
			temp=min(temp,bb-aa);
			temp=min(temp,b-a);
			ans1+=(ll)b-a;
			ans2+=(ll)b-a-temp;
			int cnt=0;
			int num=temp;
			while(num)
			{
				cnt++;
				num/=10;
			}
			if(temp==0)
				cnt++;
			ans2+=(ll)cnt+1;
			aa=a;
			bb=b;
		}
		printf("%I64d %I64d\n",ans1,ans2);
	}
	return 0;
}