题目大意:略

40分暴力还是很好写的,差分再跑个后缀和 和 后缀最大值就行了

一种正解是后缀数组+并查集

但据说还有后缀数组+单调栈的高端操作蒟蒻的我当然不会

后缀数组求出height,然后从大到小排个序,那么每次把最大height的两个后缀放到一个并查集里

这样,每次根据height合并两个并查集的时候,一个并查集任何一个后缀和另一个并查集任何一个后缀的LCP长度都是height[i]

并查集维护集合内a[i]最大值,最小值(负负得正),然后每次合并都打差分,统计答案的时候跑一遍后缀和和后缀最大值即可

 #include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define ll long long
#define N 301000
#define rint register int
#define inf 1000000000000000001ll
using namespace std;
//re
int len;
char str[N];
int tr[N],rk[N],hs[N],sa[N],fa[N];
ll a[N],ma[N],mi[N],ans[N],sum[N],p[N],sz[N];
struct Height{int id,w;}h[N];
int cmp(Height s1,Height s2){return s1.w>s2.w;}
bool check(int k,int x,int y){
if(x+k>len||y+k>len) return ;
else return (rk[x]==rk[y]&&rk[x+k]==rk[y+k])?:;
}
void get_sa()
{
int cnt=;rint i;
for(i=;i<=len;++i) hs[str[i]]++;
for(i=;i<=;++i) if(hs[i]) tr[i]=++cnt;
for(i=;i<=;++i) hs[i]+=hs[i-];
for(i=;i<=len;++i) rk[i]=tr[str[i]],sa[hs[str[i]]--]=i;
for(int k=;cnt<len;k<<=)
{
for(i=;i<=cnt;++i) hs[i]=;
for(i=;i<=len;++i) hs[rk[i]]++;
for(i=;i<=cnt;++i) hs[i]+=hs[i-];
for(i=len;i>=;i--) if(sa[i]>k) tr[sa[i]-k]=hs[rk[sa[i]-k]]--;
for(i=;i<=k;++i) tr[len-i+]=hs[rk[len-i+]]--;
for(i=;i<=len;++i) sa[tr[i]]=i;
for(i=,cnt=;i<=len;++i) tr[sa[i]]=check(k,sa[i],sa[i-])?cnt:++cnt;
for(i=;i<=len;++i) rk[i]=tr[i];
}
for(i=;i<=len;++i){
h[i].id=i;
if(rk[i]==) continue;
for(int j=max(,h[rk[i-]].w-);;j++)
if(str[i+j-]==str[sa[rk[i]-]+j-]) h[rk[i]].w=j;
else break;
}
}
int find_fa(int x){
int y=x;
while(y!=fa[y]) y=fa[y];
while(fa[x]!=y){
int pre=fa[x];
fa[x]=y;
x=pre;
}return y;
}
void solve()
{
rint i;int x,y,fx,fy;
for(i=;i<=len;++i) sz[i]=,fa[i]=i,mi[i]=ma[i]=a[sa[i]];
for(i=;i<=len;)
{
int j=h[i].w;
for(;h[i].w>=j&&i<=len;i++){
x=h[i].id-;fx=find_fa(x);
y=h[i].id;fy=find_fa(y);
p[j]+=sz[fx]*sz[fy];
ans[j]=max(ans[j],max(mi[fx]*mi[fy],ma[fx]*ma[fy]));
fa[fy]=fx,sz[fx]+=sz[fy];
ma[fx]=max(ma[fx],ma[fy]),mi[fx]=min(mi[fx],mi[fy]);
}
}
} int main()
{
scanf("%d",&len);
scanf("%s",str+);
for(int i=;i<=len;++i)
scanf("%lld",&a[i]);
memset(ans,-0x3f,sizeof(ans));
get_sa();
sort(h+,h+len+,cmp);
solve();
for(int i=len;i>=;i--)
sum[i]=sum[i+]+p[i],ans[i]=max(ans[i],ans[i+]);
for(int i=;i<len;i++)
printf("%lld %lld\n",sum[i],(ans[i]<=-inf)?:ans[i]);
return ;
}

NOI 2015 品酒大会 (后缀数组+并查集)的更多相关文章

  1. [UOJ#131][BZOJ4199][NOI2015]品酒大会 后缀数组 + 并查集

    [UOJ#131][BZOJ4199][NOI2015]品酒大会 试题描述 一年一度的“幻影阁夏日品酒大会”隆重开幕了.大会包含品尝和趣味挑战两个环节,分别向优胜者颁发“首席品酒家”和“首席猎手”两个 ...

  2. 【BZOJ4199】[Noi2015]品酒大会 后缀数组+并查集

    [BZOJ4199][Noi2015]品酒大会 题面:http://www.lydsy.com/JudgeOnline/wttl/thread.php?tid=2144 题解:听说能用SAM?SA默默 ...

  3. [NOI2015] 品酒大会 - 后缀数组,并查集,STL,启发式合并

    [NOI2015] 品酒大会 Description 对于每一个 \(i \in [0,n)\) 求有多少对后缀满足 LCP 长度 \(\le i\) ,并求满足条件的两个后缀权值乘积的最大值. So ...

  4. BZOJ 4199: [Noi2015]品酒大会( 后缀数组 + 并查集 )

    求出后缀数组后, 对height排序, 从大到小来处理(r相似必定是0~r-1相似), 并查集维护. 复杂度O(NlogN + Nalpha(N)) ------------------------- ...

  5. 【学术篇】NOI2015 品酒大会 后缀数组+并查集

    省选前大致是刷不了几道题了... 所以就找一些裸一点的题目练练板子算了= = 然而这题一点都不裸, 也并不怎么好写... 于是就浪费了将近一下午的时间... 然而还不是因为后缀数组板子不熟= = 首先 ...

  6. Uoj #131. 【NOI2015】品酒大会 后缀数组,并查集

    #131. [NOI2015]品酒大会 统计 描述 提交 自定义测试 一年一度的“幻影阁夏日品酒大会”隆重开幕了.大会包含品尝和趣味挑战两个环节,分别向优胜者颁发“首席品酒家”和“首席猎手”两个奖项, ...

  7. 【BZOJ-4199】品酒大会 后缀数组 + 并查集合并集合

    4199: [Noi2015]品酒大会 Time Limit: 10 Sec  Memory Limit: 512 MBSubmit: 436  Solved: 243[Submit][Status] ...

  8. [LOJ 2133][UOJ 131][BZOJ 4199][NOI 2015]品酒大会

    [LOJ 2133][UOJ 131][BZOJ 4199][NOI 2015]品酒大会 题意 给定一个长度为 \(n\) 的字符串 \(s\), 对于所有 \(r\in[1,n]\) 求出 \(s\ ...

  9. BZOJ 4566 JZYZOJ 1547 [haoi2016T5]找相同子串 后缀数组 并查集

    http://172.20.6.3/Problem_Show.asp?id=1547 http://www.lydsy.com/JudgeOnline/problem.php?id=4566 单纯后缀 ...

随机推荐

  1. 【LibreOJ 6277】数列分块入门 1 (分块)

    emmm-学下分块~ 区间:数列中连续一段的元素 区间操作:将某个区间[a,b]的所有元素进行某种改动的操作 块:我们将数列划分成若干个不相交的区间,每个区间称为一个块 整块:在一个区间操作时,完整包 ...

  2. iptables防火墙和selinux

    iptables 存在以下两种方式: 一.service方式 查看防火墙状态: [root@centos6 ~]# service iptables status iptables:未运行防火墙 开启 ...

  3. Problem 14

    Problem 14 # Problem_14.py """ The following iterative sequence is defined for the se ...

  4. Centos6.5安装Seafile,遇到的问题处理记录。

    问题1:启动Seafile安装脚本时,提示找不到MySQL-python模块,使用yum安装成功也提示未安装该软件包 问题1解决方法:需要通过 python 的工具pip来安装MySQL-python ...

  5. J2EE异常问题总结

    http://blog.csdn.net/rchm8519/article/details/41624381 http://www.51testing.com/html/49/240349-14393 ...

  6. socket 客户端的认证

    一:使用 hashlib 进行加密验证: # server.py 服务端 import os import socket import hashlib ​ def check_conn(conn): ...

  7. MAVEN项目的搭建

    MAVEN能为我们做什么? 1.Jar的声明式依赖性管理 2.项目的自动构建 搭建流程 环境配置 http://maven.apache.org/download.html 下载最新版本Maven 3 ...

  8. 解析xml文件,遍历输出xml文件中的所有节点, 最终模仿实现struts2框架自动封装参数的功能

    程序结构:src文件夹下存放xml文件 该文件内容: <?xml version="1.0" encoding="UTF-8"?> <myst ...

  9. JavaScript之Ajax技术

  10. ASP.NET--identity笔记及截图

    aspnetUsers aspnetRoles aspnetUserroles aspnetUserClaims 存储用户user额外信息的键值对 aspnetUserLogins 第三方登陆的控件 ...