题目链接

很容易想到p和q"r相似"就等价于在后缀数组中q与p之间的height值\(\ge r\),也就是说\(<r\)的那些height值会把排好序后的后缀分割成若干段,可以想到倒序枚举r,每当r减小1时,中间把他们分开的那些"屏障"就会减少一些,两边的集合会合并到一起,现在的问题就是给定一个集合,支持合并操作,每次询问元素乘积最大值(方案数显然是\(C_{siz}^2\))。由于有负值,只需记正值中最大的两个与负值中最小的两个即可。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<vector>
#include<algorithm>
#include<cmath>
#include<queue>
#include<bitset>
#define P puts("lala")
#define cp cerr<<"lala"<<endl
#define ln putchar('\n')
#define pb push_back
#define fi first
#define se second
#define mkp make_pair
using namespace std;
inline int read()
{
char ch=getchar();int g=1,re=0;
while(ch<'0'||ch>'9') {if(ch=='-')g=-1;ch=getchar();}
while(ch<='9'&&ch>='0') re=(re<<1)+(re<<3)+(ch^48),ch=getchar();
return re*g;
}
typedef long long ll;
typedef pair<int,int> pii; const int N=300050;
char s[N];
int n,Rank[N],tt2[N],sa[N],cnt[N],hei[N];
void buildsa()
{
int i,m=0,*x=Rank,*y=tt2;
for(i=1;i<=n;++i) cnt[x[i]=s[i]]++,m=max(x[i],m);
for(i=2;i<=m;++i) cnt[i]+=cnt[i-1];
for(i=n;i;--i) sa[cnt[x[i]]--]=i;
for(int k=1;k<n;k<<=1)
{
int p=0;
for(i=n-k+1;i<=n;++i) y[++p]=i;
for(i=1;i<=n;++i) if(sa[i]>k) y[++p]=sa[i]-k;
for(i=1;i<=m;++i) cnt[i]=0;
for(i=1;i<=n;++i) cnt[x[y[i]]]++;
for(i=2;i<=m;++i) cnt[i]+=cnt[i-1];
for(i=n;i;--i) sa[cnt[x[y[i]]]--]=y[i];
swap(x,y);
p=1; x[sa[1]]=1;
for(i=2;i<=n;++i) x[sa[i]]=y[sa[i-1]]==y[sa[i]]
&&y[sa[i-1]+k>n?0:sa[i-1]+k]==y[sa[i]+k>n?0:sa[i]+k]?p:++p;
m=p;
if(m>=n) break;
}
for(i=1;i<=n;++i) Rank[sa[i]]=i;
}
void buildheight()
{
int i,k=0;
for(i=1;i<=n;++i)
{
if(k) k--;
int j=sa[Rank[i]-1];
while(s[i+k]==s[j+k]) k++;
hei[Rank[i]]=k;
}
} struct dsynh
{
vector<int>neg,pos;
dsynh() {neg.clear();pos.clear();}
inline ll query()
{
if(neg.size()+pos.size()<=1) return 0;
ll ans=-3e18;
if(neg.size()==2) ans=max(ans,1ll*neg[0]*neg[1]);
if(pos.size()==2) ans=max(ans,1ll*pos[0]*pos[1]);
if(neg.size()&&pos.size()) ans=max(ans,1ll*neg[0]*pos[0]);
return ans;
}
inline void push(ll x)
{
if(x<0) neg.pb(x);
else pos.pb(x);
}
}t[N];
vector<int>ve[N];
int fa[N],siz[N];
inline int find(int x)
{
if(fa[x]!=x) return fa[x]=find(fa[x]);
return fa[x];
}
ll ans1=0;
ll merge(int x,int y)
{
x=find(x); y=find(y);
ans1-=(1ll*siz[x]*(siz[x]-1)/2+1ll*siz[y]*(siz[y]-1)/2);
fa[x]=y; siz[y]+=siz[x];
ans1+=1ll*siz[y]*(siz[y]-1)/2; for(int i=0;i<t[x].neg.size();++i) t[y].push(t[x].neg[i]);
for(int i=0;i<t[x].pos.size();++i) t[y].push(t[x].pos[i]);
t[x].neg.clear(); t[x].pos.clear(); sort(t[y].neg.begin(),t[y].neg.end(),greater<int>() );
while(t[y].neg.size()>2) t[y].neg.erase(t[y].neg.begin()); sort(t[y].pos.begin(),t[y].pos.end());
while(t[y].pos.size()>2) t[y].pos.erase(t[y].pos.begin());
return t[y].query();
} int val[N];
ll out1[N],out2[N];
int main()
{
#ifndef ONLINE_JUDGE
freopen("1.in","r",stdin);freopen("1.out","w",stdout);
#endif
int i,j,opt,T;
n=read();
scanf("%s",s+1);
for(i=1;i<=n;++i) val[i]=read();
buildsa(); buildheight(); for(i=2;i<=n;++i) ve[hei[i]].pb(i);
for(i=1;i<=n;++i) siz[i]=1,fa[i]=i,t[i].push(val[sa[i]]); ll ans2=-3e18;
for(int r=n-1;r>=0;--r)
{
int siz=ve[r].size();
for(i=0;i<siz;++i)
ans2=max(merge(ve[r][i]-1,ve[r][i]),ans2);
out1[r]=ans1; out2[r]=ans2==-3e18?0:ans2;
}
for(i=0;i<n;++i) printf("%lld %lld\n",out1[i],out2[i]);
return 0;
}

uoj131【NOI2015】品酒大会的更多相关文章

  1. BZOJ4199/UOJ131 [Noi2015]品酒大会

    本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000 作者博客:http://www.cnblogs.com/ljh2000-jump/ ...

  2. UOJ131 [NOI2015] 品酒大会

    考前挣扎(bu shi 之前留下来的坑 首先注意到 SAM的parent树 是反串的后缀树 也就是原串的前缀树 这个性质很重要 所以说我们在树上统计的时候两个点的lca就是两个后缀串的lcp 于是可以 ...

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

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

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

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

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

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

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

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

  7. [UOJ#131][BZOJ4199][NOI2015]品酒大会

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

  8. BZOJ_4199_[Noi2015]品酒大会_后缀自动机

    BZOJ_4199_[Noi2015]品酒大会_后缀自动机 Description 一年一度的“幻影阁夏日品酒大会”隆重开幕了.大会包含品尝和趣味挑战两个环节,分别向优胜者颁发“首席品 酒家”和“首席 ...

  9. [NOI2015]品酒大会(SA数组)

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

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

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

随机推荐

  1. npm run dev 启动项目报错我的解决办法

    我的报错截屏 解决方案   1.    config文件中 index 文件中的 host 值如果是数字串就将其改为 localhost 2.    再次尝试 如果有遇到其他问题阔以将 node-mo ...

  2. Luogu1083 借教室 (线段树)

    当根都<0,就不行了 #include <iostream> #include <cstdio> #include <cstring> #include &l ...

  3. Ceph创建一个新集群 报错: File "/usr/bin/ceph-deploy", line 18, in..........

    [root@ceph-node1 ceph]# ceph-deploy new node1 Traceback (most recent call last): File "/usr/bin ...

  4. Mac设置外网访问本地项目

    > 官网地址:https://ngrok.com/download 步骤(官网基本已经说明了步骤,但还不完善,以下为亲测步骤): 下载并注册账号 打开终端,进入ngrok.zip所在文件夹(方法 ...

  5. Linux-shell笔记1

    一次执行很多命令,可以用:分割每个命令,依次运行所有命令.但是不是进程列表,要用()包围命令才是进程列表.它们有什么差别呢?进程列表是启动了一个子SHELL来执行的.用echo $BASH_SUBSH ...

  6. 「雅礼集训 2017 Day2」线段游戏(线段树懒标记“启发式下传”,李超树)

    题面 题解 加入一条线段,可以把它转化为在[L,R]区间内加一条线 y=ax+b (如果原线段与y轴平行,就相当于在{x1}处加一条线 y=max(y1,y2)) 我们可以把它加到线段树上,线段树上每 ...

  7. 【MySQL】从入门到精通8-SQL数据库编程

    上期:[MySQL]从入门到精通7-设计多对多数据库 第零章:Mac用户看这里: mac终端写MySQL和windows基本相同,除了配置环境变量和启动有些许不同以外. 先配置环境变量,在终端输入vi ...

  8. 快Key:按一下鼠标【滚轮】,帮你自动填写用户名密码,快速登录,可制作U盘随身(开源免费-附安装文件和源代码)

    * 代码以本文所附下载文件包为准,安装文件和源文件包均在本文尾部可下载. * 快Key及本文所有内容仅供交流使用,使用者责任自负,由快Key对使用者及其相关人员或组织造成的任何损失均由使用者自负,与本 ...

  9. Win32简单图形界面程序逆向

    Win32简单图形界面程序逆向 前言 为了了解与学习底层知识,从 汇编开始 -> C语言 -> C++ -> PE文件 ,直至今天的Win32 API,着实学的令我头皮发麻(笑哭). ...

  10. WinUI 3 踩坑记:从创建项目到发布

    本文是 WinUI 3 踩坑记 的一部分,该系列发布于 GitHub@Scighost/WinUI3Keng,若内容出现冲突以 GitHub 上的为准. 创建项目 现在 WinUI 3 的入门体验比刚 ...