【NOI2015】品酒大会

 

分析：首先可以想到 r 相似的两个子串肯定是满足 1...r-1 相似的，我们先求一遍后缀数组，根据 sa 数组的性质，0相似 即sa[1..n]整个区间，1 相似即将 1..n 划分成若干个区间，每个区间内的最长公共子串为 1，依次类推下去，则 一个 i-1 相似 的区间是由两个或多个 i 相似的区间合并的，类似并查集的合并过程。由此，我们可以以ht 数组的值从大到小遍历，维护并查集的合并，并且维护每个连通块的最大值最小值来判断最大美味度。

 

代码：

#include<vector>
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;

typedef long long LL;
const int N = 3E5+10;

char str[N];
int n,s[N],a[N],sa[N],rk[N],oldrk[N<<1];
int ht[N],px[N],id[N],cnt[N];

bool cmp(int x,int y,int w){
	return oldrk[x]==oldrk[y] && oldrk[x+w]==oldrk[y+w];
}

void da(int s[],int n,int m){
	int i,p=0,w,k;
	
	for(i=1;i<=n;i++) ++cnt[rk[i] = s[i]];
	for(i=1;i<=m;i++) cnt[i] += cnt[i-1];
	for(i=n;i>=1;i--) sa[cnt[rk[i]]--] = i;
	
	for(w=1;w<n;w<<=1,m=p){
		for(p=0,i=n;i>n-w;i--) id[++p]=i;
		for(i=1;i<=n;i++)
		    if(sa[i]>w) id[++p]=sa[i]-w;
		memset(cnt,0,sizeof(cnt));
		for(i=1;i<=n;i++) ++cnt[px[i] = rk[id[i]]];
		for(i=1;i<=m;i++) cnt[i] += cnt[i-1];
		for(i=n;i>=1;i--) sa[cnt[px[i]]--] = id[i];
		memcpy(oldrk,rk,sizeof(rk));
		
		for(p=0,i=1;i<=n;i++)
		    rk[sa[i]]=cmp(sa[i],sa[i-1],w)?p:++p;
	}
	
	for(i=1,k=0;i<=n;i++){
		if(k) --k;
		while(s[i+k]==s[sa[rk[i]-1]+k]) ++k;
		ht[rk[i]]=k;
	} 
}

struct P{
	int Min,Max,fa,size;
	P(){}
	P(int _fa,int _size,int _Min,int _Max){
		fa=_fa,size=_size,Min=_Min,Max=_Max;
	}
}f[N];

vector<int>H[N];
LL ANS[N][2],num,res;

int find(int x){
	if(f[x].fa==x) return x;
	return f[x].fa=find(f[x].fa);
}

void merge(int a,int b){
	int x=find(a),y=find(b);
	num+=(LL)f[x].size*f[y].size;
	res=max(res,max((LL)f[x].Min*f[y].Min,(LL)f[x].Max*f[y].Max));
	f[y].fa=x;
	f[x].size+=f[y].size;
	f[x].Min=min(f[x].Min,f[y].Min);
	f[x].Max=max(f[x].Max,f[y].Max);
}

int main()
{
	scanf("%d",&n);
	scanf("%s",str+1);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
	    s[i]=str[i];	
	}
	da(s,n,128);
	for(int i=1;i<=n;i++){
		f[i]=P(i,1,a[sa[i]],a[sa[i]]); 
	}
	for(int i=2;i<=n;i++) H[ht[i]].push_back(i);
	
	num=0,res=-2e18;
	for(int i=n-1;i>=0;i--){
	    for(int j=0;j<H[i].size();j++){
	    	merge(H[i][j],H[i][j]-1);
		}
		if(num) ANS[i][0]=num,ANS[i][1]=res;
	}
	
	for(int i=0;i<n;i++){
		printf("%lld %lld\n",ANS[i][0],ANS[i][1]);
	}
}