题解

想出了一个神奇的技巧

我们先把串反过来(因为我们需要起始位置的值而不是终止位置的值),每个点维护一下 fail树上子树里的点,作为正数绝对值最大的两个数,作为负数绝对值最大的两个数

我们发现这个两两相乘的最大值肯定是一个后缀最大值,我们对每个节点求一个最大值存到这个节点的len处理出后缀最大值就行

然后关键是怎么求这个多少个位置两两匹配前缀是0,1,2....n - 1

似乎会想到把后缀自动机上节点的cnt累加进去就好?然而,这样是错的,你连第一个样例都过不了

为什么呢,我们把样例第一个串反过来之后看一看

ponoiiipoi

iopiiionop

你如果要统计两个p的位置作为相同长度为1的两个位置

非常尴尬,没有这个节点

怎么办,你BFS显然会超时……为什么i这个位置就会被统计进去……因为有这个节点

ii这个节点的par就是i

经过一番大胆的猜想,发现我们需要把\([p->par->len + 1,p->len]\)这一段的答案要加上\(p->cnt *(p->cnt - 1) / 2\)

那么为什么这样是正确的呢……显然= =

因为p->par就是p的一段后缀,这段后缀显然会累加上由于最长匹配是个长串带来的贡献,而中间的一段也是这个长串的一部分,没有被加上

代码

#include <bits/stdc++.h>
//#define ivorysi
#define enter putchar('\n')
#define space putchar(' ')
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define eps 1e-8
#define mo 974711
#define MAXN 300005
#define pii pair<int,int>
using namespace std;
typedef long long int64;
typedef double db;
template<class T>
void read(T &res) {
res = 0;char c = getchar();T f = 1;
while(c < '0' || c > '9') {
if(c == '-') f = -1;
c = getchar();
}
while(c >= '0' && c <= '9') {
res = res * 10 + c - '0';
c = getchar();
}
res *= f;
}
template<class T>
void out(T x) {
if(x < 0) {putchar('-');x = -x;}
if(x >= 10) {
out(x / 10);
}
putchar('0' + x % 10);
}
struct node {
int len,cnt;
int val[2][2],vc[2];
node *nxt[26],*par;
}pool[MAXN * 2],*tail = pool,*root,*last,*que[MAXN * 2];
int a[MAXN],N,c[MAXN];
int64 ans[MAXN],sum[MAXN],d[MAXN];
char str[MAXN];
void build_SAM(int e,int len,int v) {
node *nowp = tail++,*p;
nowp->len = len;nowp->cnt = 1;
if(v < 0) nowp->val[0][nowp->vc[0]++] = v;
else nowp->val[1][nowp->vc[1]++] = v;
for(p = last ; p && !p->nxt[e] ; p = p->par) {
p->nxt[e] = nowp;
}
if(!p) nowp->par = root;
else {
node *q = p->nxt[e];
if(q->len == p->len + 1) nowp->par = q;
else {
node *copyq = tail++;
*copyq = *q;
copyq->vc[0] = copyq->vc[1] = 0;
memset(copyq->val,0,sizeof(copyq->val));
copyq->cnt = 0;
copyq->len = p->len + 1;
q->par = nowp->par = copyq;
for(; p && p->nxt[e] == q ; p = p->par) p->nxt[e] = copyq;
}
}
last = nowp;
}
void update(node *p,int on,int v) {
if(p->vc[on] < 2) {p->val[on][1] = v;++p->vc[on];}
else {
if(abs(v) > abs(p->val[on][1])) p->val[on][1] = v;
}
if(abs(p->val[on][1]) > abs(p->val[on][0])) swap(p->val[on][1],p->val[on][0]);
}
void Init() {
read(N);scanf("%s",str + 1);
for(int i = 1 ; i <= N ; ++i) read(a[i]);
for(int i = 1 ; i <= N / 2 ; ++i) {
swap(a[i],a[N - i + 1]);swap(str[i],str[N - i + 1]);
}
root = last = tail++;
for(int i = 1 ; i <= N ; ++i) build_SAM(str[i] - 'a',i,a[i]);
}
void Solve() {
int m = tail - pool;
for(int i = 0 ; i < m ; ++i) c[pool[i].len]++;
for(int i = 1 ; i <= N ; ++i) c[i] += c[i - 1];
for(int i = 0 ; i < m ; ++i) que[c[pool[i].len]--] = &pool[i];
for(int i = 0 ; i <= N ; ++i) ans[i] = -1e18;
for(int i = m ; i >= 1 ; --i) {
if(que[i]->cnt > 1) {
sum[que[i]->len] += 1LL * (que[i]->cnt - 1) * (que[i]->cnt) / 2;
if(que[i]->vc[0] >= 2) ans[que[i]->len] = max(ans[que[i]->len],1LL * que[i]->val[0][1] * que[i]->val[0][0]);
if(que[i]->vc[1] >= 2) ans[que[i]->len] = max(ans[que[i]->len],1LL * que[i]->val[1][0] * que[i]->val[1][1]);
for(int k = 0 ;k < que[i]->vc[0] ; ++k) {
for(int h = 0 ; h < que[i]->vc[1] ; ++h) {
ans[que[i]->len] = max(ans[que[i]->len],1LL * que[i]->val[0][k] * que[i]->val[1][h]);
}
}
}
d[que[i]->len - 1] += 1LL * (que[i]->cnt - 1) * que[i]->cnt / 2;
if(que[i]->par) {
que[i]->par->cnt += que[i]->cnt;
for(int k = 0 ; k < que[i]->vc[0] ; ++k) {
update(que[i]->par,0,que[i]->val[0][k]);
}
for(int k = 0 ; k < que[i]->vc[1] ; ++k) {
update(que[i]->par,1,que[i]->val[1][k]);
}
d[que[i]->par->len] -= 1LL * (que[i]->cnt - 1) * que[i]->cnt / 2;
} }
for(int i = N ; i >= 1 ; --i) d[i] += d[i + 1],sum[i] += d[i];
for(int i = N - 1; i >= 1 ; --i) ans[i] = max(ans[i + 1],ans[i]);
for(int i = 0 ; i <= N - 1 ; ++i) {
if(!sum[i]) ans[i] = 0;
out(sum[i]);space;out(ans[i]);enter;
}
}
int main() {
#ifdef ivorysi
freopen("f1.in","r",stdin);
#endif
Init();
Solve();
return 0;
}

【LOJ】#2133. 「NOI2015」品酒大会的更多相关文章

  1. loj#2129. 「NOI2015」程序自动分析

    题目链接 loj#2129. 「NOI2015」程序自动分析 题解 额... 考你会不会离散化优化常数 代码 #include<queue> #include<cstdio> ...

  2. *LOJ#2134. 「NOI2015」小园丁与老司机

    $n \leq 5e4$个平面上的点,从原点出发,能从当前点向左.右.上.左上或右上到达该方向最近的给定点.问三个问:一.最多经过多少点:二.前一问的方案:三.其所有方案种非左右走的边至少要开几辆挖掘 ...

  3. LOJ#2132. 「NOI2015」荷马史诗

    $n \leq 100000$个数字,放进$k$叉树里,一个点只能放一个数,使所有数字乘以各自深度这个值之和最小的同时,最大深度的数字最小. 哈夫曼.这是我刚学OI那段时间看到的,感觉就是个很无聊的贪 ...

  4. LOJ#2131. 「NOI2015」寿司晚宴

    $n \leq 500$,$2-n$这些数字,两个人挑,可以重复挑,问有几种方案中,一个人选的所有数字与另一个人选的所有数字都互质. 不像前两题那么抠脚.. 如果$n$比较小的话,可以把两个人选的数字 ...

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

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

  6. 【BZOJ4199】【NOI2015】品酒大会(后缀数组)

    [BZOJ4199][NOI2015]品酒大会 题面 BZOJ Uoj 洛谷 题解 考虑最裸的暴力 枚举每次的长度 以及两个开始的位置 检查以下是否满足条件,如果可以直接更新答案 复杂度\(O(n^3 ...

  7. Loj #2192. 「SHOI2014」概率充电器

    Loj #2192. 「SHOI2014」概率充电器 题目描述 著名的电子产品品牌 SHOI 刚刚发布了引领世界潮流的下一代电子产品--概率充电器: 「采用全新纳米级加工技术,实现元件与导线能否通电完 ...

  8. Loj #3096. 「SNOI2019」数论

    Loj #3096. 「SNOI2019」数论 题目描述 给出正整数 \(P, Q, T\),大小为 \(n\) 的整数集 \(A\) 和大小为 \(m\) 的整数集 \(B\),请你求出: \[ \ ...

  9. Loj #3093. 「BJOI2019」光线

    Loj #3093. 「BJOI2019」光线 题目描述 当一束光打到一层玻璃上时,有一定比例的光会穿过这层玻璃,一定比例的光会被反射回去,剩下的光被玻璃吸收. 设对于任意 \(x\),有 \(x\t ...

随机推荐

  1. Codeforces 932.F Escape Through Leaf

    F. Escape Through Leaf time limit per test 3 seconds memory limit per test 256 megabytes input stand ...

  2. springboot集成Guava缓存

    很久没有写博客了,这段时间一直忙于看论文,写论文,简直头大,感觉还是做项目比较舒服,呵呵,闲话不多说,今天学习了下Guava缓存,这跟Redis类似的,但是适用的场景不一样,学习下吧.今天我们主要是s ...

  3. ACE服务端编程2:ACE跨平台之数据类型和宽字符

    ACE网络库的主要优势之一就是跨平台,ACE提供了操作系统API和编译器级别的跨平台解决方法,使开发人员不用再去关心操作系统和编译器的差异,但因此也带来了ACE的复杂性. ACE网络库的组织结构主要分 ...

  4. LeakCanary原理分析

    参考文档 http://blog.csdn.net/wyfei021/article/details/46506521http://vjson.com/wordpress/leakcanary%e6% ...

  5. (转)select、poll、epoll之间的区别

    本文来自:https://www.cnblogs.com/aspirant/p/9166944.html (1)select==>时间复杂度O(n) 它仅仅知道了,有I/O事件发生了,却并不知道 ...

  6. 816B. Karen and Coffee 前缀和思维 或 线段树

    LINK 题意:给出n个[l,r],q个询问a,b,问被包含于[a,b]且这样的区间数大于k个的方案数有多少 思路:预处理所有的区间,对于一个区间我们标记其(左边界)++,(右边界+1)--这样就能通 ...

  7. 重构改善既有代码设计--重构手法04:Replace Temp with Query (以查询取代临时变量)

    所谓的以查询取代临时变量:就是当你的程序以一个临时变量保存某一个表达式的运算效果.将这个表达式提炼到一个独立函数中.将这个临时变量的所有引用点替换为对新函数的调用.此后,新函数就可以被其他函数调用. ...

  8. [转] Linux下程序的加载、运行和终止流程

    TAG: linux, main, _start DATE: 2013-08-08 原文地址: http://blog.csdn.net/tigerscorpio/article/details/62 ...

  9. 【CodeForces】889 B. Restoration of string

    [题目]B. Restoration of string [题意]当一个字符串在字符串S中的出现次数不小于任意子串的出现次数时,定义这个字符串是高频字符串.给定n个字符串,求构造出最短的字符串S满足着 ...

  10. Spring Boot工程结构推荐

    工程结构(最佳实践) Spring Boot框架本身并没有对工程结构有特别的要求,但是按照最佳实践的工程结构可以帮助我们减少可能会遇见的坑,尤其是Spring包扫描机制的存在,如果您使用最佳实践的工程 ...