利用SAM建出后缀树,树上每个节点计算一下|right|、right集合中ai的最大、次大、最小、次小值即可。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define N 600010
#define inf 1000000010
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{
int x=0,f=1;char c=getchar();
while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
return x*f;
}
int n,a[N],son[N][26],fail[N],len[N],size[N],mx[N],mx2[N],mn[N],mn2[N],cnt=1,last=1;
ll ans[N],tot[N];
char s[N];
int p[N],t;
struct data{int to,nxt;
}edge[N<<1];
void addedge(int x,int y){t++;edge[t].to=y,edge[t].nxt=p[x],p[x]=t;}
void ins(int c)
{
int x=++cnt,p=last;last=x;len[x]=len[p]+1;size[x]=1;mx[x]=mn[x]=a[len[x]];mx2[x]=-inf,mn2[x]=inf;
while (!son[p][c]) son[p][c]=x,p=fail[p];
if (!p) fail[x]=1;
else
{
int q=son[p][c];
if (len[p]+1==len[q]) fail[x]=q;
else
{
int y=++cnt;mx[y]=mx2[y]=-inf,mn[y]=mn2[y]=inf;
len[y]=len[p]+1;
memcpy(son[y],son[q],sizeof(son[q]));
fail[y]=fail[q],fail[q]=fail[x]=y;
while (son[p][c]==q) son[p][c]=y,p=fail[p];
}
}
}
void dfs(int k)
{
for (int i=p[k];i;i=edge[i].nxt)
{
dfs(edge[i].to);
size[k]+=size[edge[i].to];
tot[len[k]]-=1ll*size[edge[i].to]*(size[edge[i].to]-1)>>1;
if (mx2[edge[i].to]>mx[k]) mx[k]=mx[edge[i].to],mx2[k]=mx2[edge[i].to];
else if (mx[edge[i].to]>mx[k]) mx2[k]=max(mx2[edge[i].to],mx[k]),mx[k]=mx[edge[i].to];
else mx2[k]=max(mx2[k],mx[edge[i].to]);
if (mn2[edge[i].to]<mn[k]) mn[k]=mn[edge[i].to],mn2[k]=mn2[edge[i].to];
else if (mn[edge[i].to]<mn[k]) mn2[k]=min(mn2[edge[i].to],mn[k]),mn[k]=mn[edge[i].to];
else mn2[k]=min(mn2[k],mn[edge[i].to]);
}
tot[len[k]]+=1ll*size[k]*(size[k]-1)>>1;
if (mx2[k]>-inf) ans[len[k]]=max(ans[len[k]],1ll*mx[k]*mx2[k]);
if (mn2[k]<inf) ans[len[k]]=max(ans[len[k]],1ll*mn[k]*mn2[k]);
}
int main()
{
#ifndef ONLINE_JUDGE
freopen("a.in","r",stdin);
freopen("a.out","w",stdout);
const char LL[]="%I64d\n";
#else
const char LL[]="%lld\n";
#endif
n=read();mx[1]=mx2[1]=-inf,mn[1]=mn2[1]=inf;
scanf("%s",s+1);reverse(s+1,s+n+1);
for (int i=1;i<=n;i++) a[i]=read();reverse(a+1,a+n+1);
for (int i=1;i<=n;i++) ins(s[i]-'a');
for (int i=2;i<=cnt;i++) addedge(fail[i],i);
memset(ans,200,sizeof(ans));
dfs(1);
for (int i=n-1;i>=0;i--) ans[i]=max(ans[i],ans[i+1]),tot[i]+=tot[i+1];
for (int i=0;i<n;i++) if (tot[i]==0) ans[i]=0;
for (int i=0;i<n;i++) printf("%lld %lld\n",tot[i],ans[i]);
return 0;
}

  

最新文章

  1. java SSH框架详解(面试和学习都是最好的收藏资料)
  2. spi_flash
  3. 4D数据介绍
  4. JavaScript高级程序设计(第三版)学习笔记6、7章
  5. 用zepto找到当前被选中的option
  6. TCPL 行计数
  7. 【xcode插件介绍】Alcatraz ----The package manager for Xcode
  8. [SOJ] shortest path in unweighted graph
  9. JVM类加载以及执行的实战
  10. 恶补web之八:jQuery(2)
  11. day09 函数学习
  12. 《DOM Scripting》学习笔记-——第八章 充实文档的内容
  13. 2018-12-09 疑似bug_中文代码示例之Programming in Scala笔记第九十章
  14. nginx1.14的安装
  15. Python操作redis系列之 列表(list) (五)
  16. ie6,ie7下a标签无法点击(转载)
  17. 第三百七十四节,Django+Xadmin打造上线标准的在线教育平台—创建课程app,在models.py文件生成4张表,课程表、课程章节表、课程视频表、课程资源表
  18. 设计模式初探-桥接(Bridge)模式
  19. PhoneGap+Cordova+SenchaTouch-01-环境搭建
  20. EventBus的基本使用步骤

热门文章

  1. selenium反爬机制
  2. 【spring源码分析】IOC容器初始化——查漏补缺(五)
  3. CentOS 修改固定IP地址
  4. jQuery插件fontIconPicker配合FontAwesome字体图标库的使用
  5. QDateTime QString
  6. Linux -- 管理锁争用(翻译)
  7. OpenStack社区中的GO语言之争
  8. 123457123457#0#-----com.cym.shuXueWangGuo1--前拼后广--儿童数学
  9. Swift学习 (二)
  10. PostgreSQL学习笔记——事务