省选前大致是刷不了几道题了... 所以就找一些裸一点的题目练练板子算了= =

然而这题一点都不裸, 也并不怎么好写... 于是就浪费了将近一下午的时间... 然而还不是因为后缀数组板子不熟= =

首先这个"r相似"很显然就是lcp的值, 也就能想到后缀数组上的height... 不会后缀数组的先左转百度~

那么我们考虑如果有一个连续的区间, 它们的height值都是大于等于r的, 那么这段区间中的后缀两两"r相似".

而"r相似"的话, 也肯定有"r-1相似", "r-2相似", ... "0相似". 这样我们就会重复统计, 就会浪费时间. 所以我们不妨将这个连续的区间表示成一个点, 并查集!!

这样我们把id按照对应位置的height降序排序, 然后对于每个id, 我们把id-1这个点所在的区间和id所在的区间合并(根据height的含义, 就是表示sa[id]和sa[id-1]所对应的后缀的lcp长度..)

合并的同时维护信息即可. 说起来挺轻巧的, 其实不是很好懂.. (当然也可能是我太蒻了 理解能力差)

最后不要忘了统计的时候做一个后缀和, 比r大的答案都要统计一下.

说的很不清楚(然而其实只是用来练板子谁曾想到这破题我写了一下午呢...)

有不懂的可以去看代码看了就更不懂了Emmmm

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const int N=300050;
const long long INF=1ll<<62;
char s[N];
int fa[N],sz[N],mn[N],mx[N],w[N],id[N];
long long cnt[N],ans[N];
int x[N],y[N],sa[N],rnk[N],cc[N],height[N],len;
inline int gn(int a=0,char c=0,int f=1){
for(;(c<'0'||c>'9')&&c!='-';c=getchar());
if(c=='-') c=getchar(),f=-1;
for(;c>47&&c<58;c=getchar()) a=a*10+c-48;
return a*f;}
bool hcmp(int x,int y){
if(height[x]==height[y]) return x<y;
return height[x]>height[y];
}
bool cmp(int *y,int a,int b,int k){
int ra=a+k>=len?-1:y[a+k],rb=b+k>=len?-1:y[b+k];
return y[a]==y[b]&&ra==rb;
}
void make_sa(){ int m=26;
for(int i=0;i<m;++i) cc[i]=0;
for(int i=0;i<len;++i) ++cc[x[i]=s[i]-'a'];
for(int i=1;i<m;++i) cc[i]+=cc[i-1];
for(int i=len-1;i>=0;--i) sa[--cc[x[i]]]=i; for(int k=1;k<=len;k<<=1){ int p=0;
for(int i=len-k;i<len;++i) y[p++]=i;
for(int i=0;i<len;++i) if(sa[i]>=k) y[p++]=sa[i]-k; for(int i=0;i<m;++i) cc[i]=0;
for(int i=0;i<len;++i) ++cc[x[y[i]]];
for(int i=1;i<m;++i) cc[i]+=cc[i-1];
for(int i=len-1;i>=0;--i) sa[--cc[x[y[i]]]]=y[i]; std::swap(x,y); m=1; x[sa[0]]=0;
for(int i=1;i<len;++i)
x[sa[i]]=cmp(y,sa[i],sa[i-1],k)?m-1:m++;
if(m>=len) break;
}
for(int i=0;i<len;++i) rnk[sa[i]]=i;
}
void make_height(){ int k=0;
for(int i=0;i<len;++i){
if(!rnk[i]) continue;
int j=sa[rnk[i]-1];
if(k) --k;
while(s[i+k]==s[j+k]) ++k;
height[rnk[i]]=k;
}
}
int find(int x){ if(fa[x]!=x) fa[x]=find(fa[x]); return fa[x];}
void merge(int x,int y){
sz[y]+=sz[x]; fa[x]=y;
mn[y]=min(mn[x],mn[y]);
mx[y]=max(mx[x],mx[y]);
}
int main(){
len=gn();
scanf("%s",s); len=strlen(s);
make_sa(); make_height();
for(int i=0;i<len;++i)
w[i]=gn();
for(int i=0;i<len;++i){
fa[i]=i; id[i]=i;
mx[i]=w[sa[i]]; mn[i]=w[sa[i]];
sz[i]=1; ans[i]=-INF;
} ::sort(id+1, id+len, hcmp);
for(int i=1;i<len;++i){
int x=find(id[i]-1),y=find(id[i]);
cnt[height[id[i]]]+=1ll*sz[x]*sz[y];
ans[height[id[i]]]=max(ans[height[id[i]]],1ll*mn[x]*mn[y]);
ans[height[id[i]]]=max(ans[height[id[i]]],1ll*mn[x]*mx[y]);
ans[height[id[i]]]=max(ans[height[id[i]]],1ll*mx[x]*mn[y]);
ans[height[id[i]]]=max(ans[height[id[i]]],1ll*mx[x]*mx[y]);
merge(x,y);
}
for(int i=len-2;i>=0;--i)
cnt[i]+=cnt[i+1],ans[i]=max(ans[i],ans[i+1]);
for(int i=0;i<len;++i)
printf("%lld %lld\n",cnt[i],cnt[i]?ans[i]:0);
}

【学术篇】NOI2015 品酒大会 后缀数组+并查集的更多相关文章

  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. NOI 2015 品酒大会 (后缀数组+并查集)

    题目大意:略 40分暴力还是很好写的,差分再跑个后缀和 和 后缀最大值就行了 一种正解是后缀数组+并查集 但据说还有后缀数组+单调栈的高端操作蒟蒻的我当然不会 后缀数组求出height,然后从大到小排 ...

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

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

  7. BZOJ 4199: [Noi2015]品酒大会 [后缀数组 带权并查集]

    4199: [Noi2015]品酒大会 UOJ:http://uoj.ac/problem/131 一年一度的“幻影阁夏日品酒大会”隆重开幕了.大会包含品尝和趣味挑战两个环节,分别向优胜者颁发“首席品 ...

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

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

  9. BZOJ.4199.[NOI2015]品酒大会(后缀数组 单调栈)

    BZOJ 洛谷 后缀自动机做法. 洛谷上SAM比SA慢...BZOJ SAM却能快近一倍... 显然只需要考虑极长的相同子串的贡献,然后求后缀和/后缀\(\max\)就可以了. 对于相同子串,我们能想 ...

随机推荐

  1. C# 编程—字符串(方法应用)、数学Math函数、DateTime、异常捕获、其他

    其他:        #--任意位数字,有几位显示几位        0--至少以为数字,不足则补0        例如:#.00--必须保留两位小数 字符串(string): Length  长度 ...

  2. mac OS 安装 Homebrew及常用命令

    Homebrew  是由国外大神 Max Howell 开发的一款包管理工具,类似Debian的apt,他可以安装任何你想安装的东西. 安装方法 命令行输入 /usr/bin/ruby -e &quo ...

  3. mysql分表详解

    经常听到有人说“数据表太大了,需要分表”,“xxxx了,要分表”的言论,那么,到底为什么要分表? 难道数据量大就要分表? mysql数据量对索引的影响 本人mysql版本为5.7 新增数据测试 为了测 ...

  4. 基于Windows 7(本地)和CentOS7.6(云端)的Minecraft服务器(无Forge/有Forge)搭建方法

    炎炎夏日中想和小伙伴们开黑的同学可以进来看一下了,本教程教你搭建基于两个平台的Minecraft服务器,这里我以Minecraft 1.11.2版本为例给大家讲解搭建流程.其中有Forge版本可以加入 ...

  5. 使用 sar 查看网卡的流量

    1.常用命令 sar -n DEV #查看当天从零点到当前时间的网卡流量信息 sar -n DEV 1 10 #每秒显示一次,共显示10次 sar -n DEV -f /var/log/sa/saxx ...

  6. mybatis generator 使用方法

    环境: ubuntu   eclipse maven 一. 简介 mybatis-geneator是一款mybatis自动代码生成工具,可以通过配置,快速生成mapper和xml文件以及pojo 二. ...

  7. git小乌龟图标不显示 windows系统

    1.修改注册表[windows + R 输入 regedit]进入注册信息,按照步骤找到 计算机\HKEY_LOCAL_MACHINE\SOFTWARE\Microsoft\Windows\Curre ...

  8. Linux的软件包管理

    此博客的环境任意. 主题Linux的软件包管理 一软件管理工具 1编译安装      2rpm包管理          3yum管理 二软件运行和编译 1ABI 应用程序的二进制接口 ABI:Appl ...

  9. LCD驱动程序架构和分析

    一.LCD驱动程序架构 1.裸机驱动代码分析 ①LCD初始化:控制器初始化,端口初始化,指明帧缓冲 ②LCD图形显示:将图形数据写入帧缓冲 void lcd_init() { lcd_port_ini ...

  10. 极简Vue的异步组件函数

    export default new Router({ routes: [ { path: '/live', name: 'live', component: () => import('@/v ...