正题

题目链接:https://www.luogu.com.cn/problem/P2178


题目大意

长度为nnn的字符串,每个位置有一个权值。
如果对于q,pq,pq,p有sq,q+r−1=sp,p+r−1s_{q,q+r-1}=s_{p,p+r-1}sq,q+r−1​=sp,p+r−1​那么就称q,pq,pq,p为rrr相似。
对于每个rrr求有多少对q,pq,pq,p是rrr相似并且求一对rrr相似的p,qp,qp,q似的权值的乘积最大。


解题思路

如果两个后缀的LCP>=rLCP>=rLCP>=r那么这两个后缀就是rrr相似的。

然后考虑排序后的heightheightheight,如果我们按照heightheightheight从大到小合并那么就说明这次合并时可以使得左右两个集合heightiheight_iheighti​相似。

开一个并查集然后按照heightheightheight从大到小排序然后一个一个合并即可


codecodecode

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=3e5+10;
ll n,m,p[N],siz[N],val[N],wal[N];
ll c[N],x[N],sa[N],y[N],rk[N],h[N];
ll ans[N],ans1[N],ans2[N],fa[N];
char s[N];
void Qsort(){for(ll i=1;i<=m;i++)c[i]=0;for(ll i=1;i<=n;i++)c[x[i]]++;for(ll i=1;i<=m;i++)c[i]+=c[i-1];for(ll i=n;i>=1;i--)sa[c[x[y[i]]]--]=y[i],y[i]=0;return;
}
void Get_SA(){m=26;for(ll i=1;i<=n;i++)x[i]=s[i]-'a'+1,y[i]=i;Qsort();for(ll w=1;w<=n;w<<=1){ll p=0;for(ll i=n-w+1;i<=n;i++)y[++p]=i;for(ll i=1;i<=n;i++)if(sa[i]>w)y[++p]=sa[i]-w;Qsort();swap(x,y);x[sa[1]]=p=1;for(ll i=2;i<=n;i++)x[sa[i]]=(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+w]==y[sa[i-1]+w])?p:++p;if(p==n)break;m=p;}return;
}
void Get_Height(){ll k=0;for(ll i=1;i<=n;i++)rk[sa[i]]=i;for(ll i=1;i<=n;i++){if(rk[i]==1)continue;if(k)k--;ll j=sa[rk[i]-1];while(i+k<=n&&j+k<=n&&s[i+k]==s[j+k])k++;h[rk[i]]=k;}return;
}
bool cmp(ll x,ll y)
{return h[x]>h[y];}
ll find(ll x)
{return fa[x]==x?(x):(fa[x]=find(fa[x]));}
void unionn(ll x,ll y,ll h){ans1[h]+=siz[x]*siz[y];ans[x]=max(ans[x],ans[y]);ans[x]=max(ans[x],max(val[x]*val[y],wal[x]*wal[y]));ans[x]=max(ans[x],max(val[x]*wal[y],wal[x]*val[y]));val[x]=min(val[x],val[y]);wal[x]=max(wal[x],wal[y]);fa[y]=x;siz[x]+=siz[y];ans2[h]=max(ans2[h],ans[x]);return;
}
int main()
{scanf("%lld",&n);scanf("%s",s+1);for(ll i=1;i<=n;i++)scanf("%lld",&val[i]),p[i]=i,siz[i]=1,wal[i]=val[i],fa[i]=i;Get_SA();Get_Height();sort(p+2,p+1+n,cmp);memset(ans2,0xcf,sizeof(ans2));memset(ans,0xcf,sizeof(ans));for(ll i=2;i<=n;i++){ll x=find(sa[p[i]-1]),y=find(sa[p[i]]);if(x==y)continue;unionn(x,y,h[p[i]]);}for(ll i=n-1;i>=0;i--)ans1[i]+=ans1[i+1],ans2[i]=max(ans2[i],ans2[i+1]);for(ll i=0;i<n;i++)printf("%lld %lld\n",ans1[i],(!ans1[i])?0:ans2[i]);
}

P2178-[NOI2015]品酒大会【SA,并查集】相关推荐

  1. 洛谷 P2178 [NOI2015]品酒大会 解题报告

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

  2. 洛谷P2178 [NOI2015]品酒大会 后缀数组+单调栈

    P2178 [NOI2015]品酒大会 题目链接 https://www.luogu.org/problemnew/show/P2178 题目描述 一年一度的"幻影阁夏日品酒大会" ...

  3. P2178 [NOI2015] 品酒大会(并查集+后缀数组)

    LINK T1\color{Red}{T1}T1 第ppp杯酒和第qqq杯酒的最大相似度就是LCP(suffix(p),suffix(q))LCP(suffix(p),suffix(q))LCP(su ...

  4. BZOJ 4199 [Noi2015]品酒大会(后缀自动机 + parent树上统计)

    4199: [Noi2015]品酒大会 Time Limit: 10 Sec  Memory Limit: 512 MB Submit: 1598  Solved: 905 [Submit][Stat ...

  5. [NOI2015]品酒大会

    题目描述 一年一度的"幻影阁夏日品酒大会"隆重开幕了.大会包含品尝和趣味挑战 两个环节,分别向优胜者颁发"首席品酒家"和"首席猎手"两个奖项 ...

  6. [BZOJ4199][NOI2015]品酒大会-后缀数组

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

  7. 【bzoj4195】[Noi2015]程序自动分析 离散化+并查集

    题目描述 在实现程序自动分析的过程中,常常需要判定一些约束条件是否能被同时满足. 考虑一个约束满足问题的简化版本:假设x1,x2,x3,-代表程序中出现的变量,给定n个形如xi=xj或xi≠xj的变量 ...

  8. BZOJ.4199.[NOI2015]品酒大会(后缀自动机 树形DP)

    BZOJ 洛谷 后缀数组做法. 洛谷上SAM比SA慢...BZOJ SAM却能快近一倍... 只考虑求极长相同子串,即所有后缀之间的LCP. 而后缀的LCP在后缀树的LCA处.同差异这道题,在每个点处 ...

  9. Bzoj4199:[NOI2015]品酒大会

    题面 Bzoj4199 Sol 后缀数组 显然的暴力就是求\(LCP\)+差分 \(40\)分 # include <bits/stdc++.h> # define RG register ...

最新文章

  1. Java的File类
  2. js进阶 12-1 jquery的鼠标事件有哪些
  3. python网络爬虫入门小程序_Python 实现网络爬虫小程序
  4. C排序算法:(二)冒泡排序
  5. linux命令行sip电话,基于Linux和MiniGUI的SIP电话终端设计
  6. BZOJ 4143 The Lawyer
  7. Linux学习笔记003----linux yum命令详解
  8. 索引深入浅出(6/10):选择正确并合适的聚集索引键
  9. 【ArcGIS教程01】前言
  10. 新版本glib使用epoll代替poll
  11. 联想小新 Pad Pro 12.6 怎么样
  12. 干货-做奶白的鲫鱼汤
  13. 小批量、多品种生产模式如何快速响应客户交期
  14. HDU 6217 BBP Formula (数学)
  15. LCD LVDS BT656 SDI
  16. c语言中时间字符化,C语言字符画,字符闪画.doc
  17. BeautifulSoup爬取贝壳网成都二手房源信息(附:完整源代码)
  18. Hive 优化--SQL执行顺序、Hive参数、数据倾斜 、小文件优化
  19. 迅雷服务器未响应,给你传授win10系统打开迅雷就卡死无响应的图文方法
  20. SMS106 短信验证码接口测试

热门文章

  1. javq接口_Java为什么要使用接口_java接口怎么使用
  2. myclipes 配置php,myclipse使用技巧
  3. flask mysql项目模板渲染_21. Flask 模板 - 宏、继承、包含
  4. php xlsx格式转换csv,我想使用C将.xls或.xlsx文件转换为.csv格式#
  5. leedcode04:转换字符串的最少操作次数
  6. vue2的响应式原理学“废”了吗?继续观摩vue3响应式原理Proxy
  7. C++ vector容器中常见的三种遍历方式
  8. 4位无符号比较器设计
  9. C++实现链式存储线索二叉树
  10. 圆桌会议 HDU - 1214(规律+模拟队列)