洛谷 P3732 [HAOI2017]供给侧改革【trie树】

参考：http://blog.csdn.net/di4covery/article/details/73065684

我以为是后缀数组+某某数据结构，结果居然是01trie！！题解说“因为是随机的所以大概不会超过40个~

好了正解是这样的，先把询问离线按r排序，然后按照每个r的位置把这个后缀的前40个（我比较非就写了50个）放进trie树里，每扫到trie深度为dep的一个节点时，对这个dep记录一个tg代表最后一次更新到这个点的后缀的开始位置，于是把这个后缀塞进去就能得到一个tg数组，tg[i]记录了最晚的后缀开始位置与当前后缀的LCP是i这个深度，于是就可以更新了。

注意这里的LCP并不是最长，所以要贪心的更新，即i*(tg[i]-max(q[ti].l-1,tg[i+1]))，如果i不是最长后缀就留到i+1更新

#include<iostream>

#include<cstdio>

#include<cstring>

#include<string>

#include<algorithm>

using namespace std;

const int N=5000005;

int n,m,tg[N],c[N][2],f[N],tot;

long long ans[N];

char s[N];

struct qwe

{

	int l,r,id;

}q[N];

int read()

{

	int r=0,f=1;

	char p=getchar();

	while(p>'9'||p<'0')

	{

		if(p=='-')

			f=-1;

		p=getchar();

	}

	while(p>='0'&&p<='9')

	{

		r=r*10+p-48;

		p=getchar();

	}

	return r*f;

}

bool cmp(const qwe &a,const qwe &b)

{

	return a.r<b.r;

}

void add(int x)

{

	int now=0;

	for(int i=0;i+x<=n&&i+1<=50;i++)

	{

		int a=s[i+x]-'0';

		if(!c[now][a])

			c[now][a]=++tot;

		tg[i+1]=max(f[c[now][a]],tg[i+1]);

		f[c[now][a]]=x;

		now=c[now][a];

	}

}

int main()

{

	n=read(),m=read();

	scanf("%s",s+1);

	for(int i=1;i<=m;i++)

		q[i].l=read(),q[i].r=read(),q[i].id=i;

	sort(q+1,q+1+m,cmp);

	int ti=1;

	for(int i=1;i<=n;i++)

	{

		add(i);

		while(ti<=m&&q[ti].r==i)

		{

			for(int i=1;i<=50;i++)

			{

				if(tg[i]>=q[ti].l)

					ans[q[ti].id]+=1ll*i*(tg[i]-max(q[ti].l-1,tg[i+1]));

				else

					break;

			}

			ti++;

		}

	}

	for(int i=1;i<=m;i++)

		printf("%lld\n",ans[i]);

	return 0;

}