这道题一看就是一种非常恶心的感觉,事实的确如此恶心到我饭都吃不下去。。
我恐怕真的不适合这个竞赛。。
这道题我想的是sa+单调栈然后全世界都比我快,全世界都用的是并查集。。
好吧说下zz的做法:首先我们想想求lcp的过程两个后缀的lcp=rmq min height(i,j),那么我们考虑height[i]并维护一个单调递增栈,如果height[i]<=height[stack[top]],那么height[stack[top]]对后面的height[j](j>i)没有贡献那么就弹出,只要知道它待在栈内的时间新加入了多少个后缀,就能知道它会被累加多少遍。很显然打一个时间戳就能记录height[stack[top]]对多少元素累加的次数。设L[i]表示i左边有多少元素height比它大对于比他的lcm肯定等于自身的height,所以两者相乘就是对于第i个height相似度为height[i]的个数,把所有累加起来就是相似度为height[i]的个数,至于最大值呢呢就用st表维护就好了。说的简单(思考的真是无比困难。。程序也是无比难打。。菜飞了我。。
PS:介于楼主表达能力有限所以有不懂的地方可以评论哟!!
注意:1.long long 2.s[i]可能为负,对于负的情况要维护最小值,所以不一定维护最大值即最大,两个负数乘起来有可能比两个正数乘起来大,神坑。。
代码:
#include<iostream>
#include<stdio.h>
#define ll long long
using namespace std;
const int Maxn=300005;
const long long oo=2e18;
char s[Maxn];
int tax[Maxn],sa[Maxn],ranks[Maxn],height[Maxn],tp[Maxn];
int log2[Maxn],st[Maxn];
int Min[Maxn][20],Max[Maxn][20],a[Maxn],L[Maxn],times[Maxn];
long long sum[Maxn],Maxa[Maxn],Mina[Maxn],ans[Maxn];
int m=300,n,now=1,top;
inline long long max(long long x,long long y)
{
return x>=y?x:y;
}
inline long long min(long long x,long long y)
{
return x>=y?y:x;
}
inline void Rsort()
{
for(int i=1;i<=m;i++)
tax[i]=0;
for(int i=1;i<=n;i++)
tax[ranks[tp[i]]]++;
for(int i=1;i<=m;i++)
tax[i]+=tax[i-1];
for(int i=n;i>=1;i--)
sa[tax[ranks[tp[i]]]--]=tp[i];
}
inline int cmp(int *f,int x,int y,int w)
{
return ((f[x]==f[y])&&(f[x+w]==f[y+w]));
}
inline void suffix()
{
for(int i=1;i<=n;i++)
ranks[i]=s[i],tp[i]=i;
Rsort();
for(int p=1,w=1,i;p<n;w+=w,m=p)
{
for(p=0,i=n-w+1;i<=n;i++)
tp[++p]=i;
for(i=1;i<=n;i++)
if(sa[i]>w)
tp[++p]=sa[i]-w;
Rsort();
swap(ranks,tp);
ranks[sa[1]]=p=1;
for(i=2;i<=n;i++)
ranks[sa[i]]=cmp(tp,sa[i],sa[i-1],w)?p:++p;
}
for(int i=1,p,j;i<=n;height[ranks[i]]=j,i++)
for(j=max(height[ranks[i-1]]-1,0),p=sa[ranks[i]-1];s[p+j]==s[i+j];j++);
}
inline long long rmqmin(int q1,int q2)
{
return q1>q2?0x3f3f3f3f:min(Min[q1][log2[q2-q1+1]],Min[q2-(1<<log2[q2-q1+1])+1][log2[q2-q1+1]]);
}
inline long long rmqmax(int q1,int q2)
{
return q1>q2?-0x3f3f3f3f:max(Max[q1][log2[q2-q1+1]],Max[q2-(1<<log2[q2-q1+1])+1][log2[q2-q1+1]]);
}
int main()
{
scanf("%d",&n);
scanf("%s",s+1);
suffix();
for(int i=0;i<=n;i++)
{
ans[i]=-oo;
Maxa[i]=-oo;
Mina[i]=oo;
}
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
log2[0]=-1;
for(int i=1;i<=n;i++)
log2[i]=log2[i>>1]+1;
for(int i=1;i<=n;i++)
{
Min[i][0]=a[sa[i]];
Max[i][0]=a[sa[i]];
}
for(int j=1;j<=log2[n];j++)
for(int i=1;i+(1<<j)-1<=n;i++)
{
Min[i][j]=min(Min[i][j-1],Min[i+(1<<j-1)][j-1]);
Max[i][j]=max(Max[i][j-1],Max[i+(1<<j-1)][j-1]);
}
for(int i=2;i<=n;i++)
{
L[i]=i-1;
Mina[i]=min(Mina[i],a[sa[i-1]]);
Maxa[i]=max(Maxa[i],a[sa[i-1]]);
while(top!=0&&height[st[top]]>=height[i])
{
int len=st[top]-L[st[top]],time_now=now-times[st[top]];
L[i]=L[st[top]];
sum[height[st[top]]]+=(long long)len*time_now;
long long Ma=rmqmax(times[st[top]]+1,now),Mi=rmqmin(times[st[top]]+1,now);
ans[height[st[top]]]=max(ans[height[st[top]]],Ma*Maxa[st[top]]);
ans[height[st[top]]]=max(ans[height[st[top]]],Mi*Mina[st[top]]);
Mina[i]=min(Mina[st[top]],Mina[i]);
Maxa[i]=max(Mina[st[top]],Maxa[i]);
Mina[i]=min(Maxa[st[top]],Mina[i]);
Maxa[i]=max(Maxa[st[top]],Maxa[i]);
top--;
}
st[++top]=i;
times[i]=now;
now++;
}
while(top!=0)
{
int len=st[top]-L[st[top]],time_now=now-times[st[top]];
sum[height[st[top]]]+=(long long)len*time_now;
long long Ma=rmqmax(times[st[top]]+1,now),Mi=rmqmin(times[st[top]]+1,now);
ans[height[st[top]]]=max(ans[height[st[top]]],Ma*Maxa[st[top]]);
ans[height[st[top]]]=max(ans[height[st[top]]],Mi*Mina[st[top]]);
top--;
}
for(int i=n-1;i>=0;i--)
{
ans[i]=max(ans[i+1],ans[i]);
sum[i]+=sum[i+1];
}
for(int i=0;i<n;i++)
if(ans[i]==-oo)
{
ans[i]=0;
}
for(int i=0;i<n;i++)
printf("%lld %lld\n",sum[i],ans[i]);
return 0;
}