NOI 2015 品酒大会 (后缀数组+并查集)
题目大意:略
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 品酒大会 (后缀数组+并查集)的更多相关文章
- [UOJ#131][BZOJ4199][NOI2015]品酒大会 后缀数组 + 并查集
[UOJ#131][BZOJ4199][NOI2015]品酒大会 试题描述 一年一度的“幻影阁夏日品酒大会”隆重开幕了.大会包含品尝和趣味挑战两个环节,分别向优胜者颁发“首席品酒家”和“首席猎手”两个 ...
- 【BZOJ4199】[Noi2015]品酒大会 后缀数组+并查集
[BZOJ4199][Noi2015]品酒大会 题面:http://www.lydsy.com/JudgeOnline/wttl/thread.php?tid=2144 题解:听说能用SAM?SA默默 ...
- [NOI2015] 品酒大会 - 后缀数组,并查集,STL,启发式合并
[NOI2015] 品酒大会 Description 对于每一个 \(i \in [0,n)\) 求有多少对后缀满足 LCP 长度 \(\le i\) ,并求满足条件的两个后缀权值乘积的最大值. So ...
- BZOJ 4199: [Noi2015]品酒大会( 后缀数组 + 并查集 )
求出后缀数组后, 对height排序, 从大到小来处理(r相似必定是0~r-1相似), 并查集维护. 复杂度O(NlogN + Nalpha(N)) ------------------------- ...
- 【学术篇】NOI2015 品酒大会 后缀数组+并查集
省选前大致是刷不了几道题了... 所以就找一些裸一点的题目练练板子算了= = 然而这题一点都不裸, 也并不怎么好写... 于是就浪费了将近一下午的时间... 然而还不是因为后缀数组板子不熟= = 首先 ...
- Uoj #131. 【NOI2015】品酒大会 后缀数组,并查集
#131. [NOI2015]品酒大会 统计 描述 提交 自定义测试 一年一度的“幻影阁夏日品酒大会”隆重开幕了.大会包含品尝和趣味挑战两个环节,分别向优胜者颁发“首席品酒家”和“首席猎手”两个奖项, ...
- 【BZOJ-4199】品酒大会 后缀数组 + 并查集合并集合
4199: [Noi2015]品酒大会 Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 436 Solved: 243[Submit][Status] ...
- [LOJ 2133][UOJ 131][BZOJ 4199][NOI 2015]品酒大会
[LOJ 2133][UOJ 131][BZOJ 4199][NOI 2015]品酒大会 题意 给定一个长度为 \(n\) 的字符串 \(s\), 对于所有 \(r\in[1,n]\) 求出 \(s\ ...
- BZOJ 4566 JZYZOJ 1547 [haoi2016T5]找相同子串 后缀数组 并查集
http://172.20.6.3/Problem_Show.asp?id=1547 http://www.lydsy.com/JudgeOnline/problem.php?id=4566 单纯后缀 ...
随机推荐
- 【LibreOJ 6277】数列分块入门 1 (分块)
emmm-学下分块~ 区间:数列中连续一段的元素 区间操作:将某个区间[a,b]的所有元素进行某种改动的操作 块:我们将数列划分成若干个不相交的区间,每个区间称为一个块 整块:在一个区间操作时,完整包 ...
- iptables防火墙和selinux
iptables 存在以下两种方式: 一.service方式 查看防火墙状态: [root@centos6 ~]# service iptables status iptables:未运行防火墙 开启 ...
- Problem 14
Problem 14 # Problem_14.py """ The following iterative sequence is defined for the se ...
- Centos6.5安装Seafile,遇到的问题处理记录。
问题1:启动Seafile安装脚本时,提示找不到MySQL-python模块,使用yum安装成功也提示未安装该软件包 问题1解决方法:需要通过 python 的工具pip来安装MySQL-python ...
- J2EE异常问题总结
http://blog.csdn.net/rchm8519/article/details/41624381 http://www.51testing.com/html/49/240349-14393 ...
- socket 客户端的认证
一:使用 hashlib 进行加密验证: # server.py 服务端 import os import socket import hashlib def check_conn(conn): ...
- MAVEN项目的搭建
MAVEN能为我们做什么? 1.Jar的声明式依赖性管理 2.项目的自动构建 搭建流程 环境配置 http://maven.apache.org/download.html 下载最新版本Maven 3 ...
- 解析xml文件,遍历输出xml文件中的所有节点, 最终模仿实现struts2框架自动封装参数的功能
程序结构:src文件夹下存放xml文件 该文件内容: <?xml version="1.0" encoding="UTF-8"?> <myst ...
- JavaScript之Ajax技术
- ASP.NET--identity笔记及截图
aspnetUsers aspnetRoles aspnetUserroles aspnetUserClaims 存储用户user额外信息的键值对 aspnetUserLogins 第三方登陆的控件 ...