【链接】 链接

【题意】

给你一个字符串s.
定义一个子串的出现值为它出现的次数*字符串的长度。
让你求里面的回文子串的最大出现值
|s|

【题解】

马拉车算法里面.
只有在回文往外扩展超过mx的时候,才出现了不同的回文。
只有这个时候需要计算。
快速计算一个子串在字符串中出现了多少次。
可以用二分。
比如子串为s[x..y]
则先获取temp = Rank[x];
然后看看往左height>=y-x+1的能走多远。
再看看往右height>=y-x+1的有多少。
用RMQ看看height在temp,mid(二分的中点)之间的height最小值是多少
显然有单调性.
找到最左和最右,就是出现的次数了。
马拉车算法虽然改造了原串。
但可以记录每个位置原本是属于字符串的哪一个位置的。

【错的次数】

在这里输入错的次数

【反思】

在这里输入反思

【代码】

#include<bits/stdc++.h>
using namespace std; const int N = 3e5;
const int MAX_CHAR = 255;//每个数字的最大值。
char s[N + 10];//如果是数字,就写成int s[N+10]就好,从0开始存
int Sa[N + 10], T1[N + 10], T2[N + 10], C[N+10];
int Height[N + 10], Rank[N + 10]; void build_Sa(int n, int m) {
int i, *x = T1, *y = T2;
for (i = 0; i<m; i++) C[i] = 0;
for (i = 0; i<n; i++) C[x[i] = s[i]]++;
for (i = 1; i<m; i++) C[i] += C[i - 1];
for (i = n - 1; i >= 0; i--) Sa[--C[x[i]]] = i;
for (int k = 1; k <= n; k <<= 1)
{
int p = 0;
for (i = n - k; i<n; i++) y[p++] = i;
for (i = 0; i<n; i++) if (Sa[i] >= k) y[p++] = Sa[i] - k;
for (i = 0; i<m; i++) C[i] = 0;
for (i = 0; i<n; i++) C[x[y[i]]]++;
for (i = 1; i<m; i++) C[i] += C[i - 1];
for (i = n - 1; i >= 0; i--) Sa[--C[x[y[i]]]] = y[i];
swap(x, y);
p = 1; x[Sa[0]] = 0;
for (i = 1; i<n; i++)
x[Sa[i]] = y[Sa[i - 1]] == y[Sa[i]] && y[Sa[i - 1] + k] == y[Sa[i] + k] ? p - 1 : p++;
if (p >= n) break;
m = p;
}
} void getHeight(int n)
{
int i, j, k = 0;
for (i = 1; i <= n; i++) Rank[Sa[i]] = i;
for (i = 0; i<n; i++) {
if (k) k--;
j = Sa[Rank[i] - 1];
while (s[i + k] == s[j + k]) k++;
Height[Rank[i]] = k;
}
} const int MAXL = 19;//log2数组的最大长度
const int INF = 0x3f3f3f3f;//数值绝对值的最大值 struct abc{
int pre2[MAXL+5],need[N+10];
int fmin[N+10][MAXL+5]; void init(int n)
{
pre2[0] = 1;
for (int i = 1;i <= MAXL;i++)
{
pre2[i] = pre2[i-1]<<1;
}
need[1] = 0; need[2] = 1;
int temp = 2;
for (int i = 3; i <= n; i++)//need[i]表示长度为i是2的多少次方,可以理解为[log2i]
if (pre2[temp] == i)
need[i] = need[i - 1] + 1, temp++;
else
need[i] = need[i - 1];
} void getst(int *a,int n)
{
memset(fmin,INF,sizeof fmin);
for (int i = 1;i <= n;i++)//下标从0开始就改成对应的就好
fmin[i][0] = a[i]; for (int l = 1;pre2[l]<=n;l++)
for (int i = 1;i <= n;i++)
if (i+pre2[l]-1<=n)
fmin[i][l] = min(fmin[i][l-1],fmin[i+pre2[l-1]][l-1]);
} int getmin(int l,int r)
{
int len = need[r-l+1];
return min(fmin[l][len],fmin[r-pre2[len]+1][len]);
} }ST; char S[2*N+100];
int belong[2*N+100],p[2*N+100];
int n;
long long ans = 0; int lcp(int x,int y)
{
return ST.getmin(x+1,y);
} void cal(int x,int y)
{
int len = y - x + 1;//子串s[x..y]的长度
int temp = Rank[x];//这个后缀x的排名 int pr = temp;//最右能到哪里
int l = temp+1,r = n;
while (l <= r)
{
int mid = (l+r)>>1;
if (lcp(temp,mid)>=len)
{
pr = mid;
l = mid + 1;
}else
r = mid - 1;
}
int pl = temp;
l = 1,r = temp - 1;
while (l <= r)
{
int mid = (l+r)>>1;
if (lcp(mid,temp)>=len)
{
pl = mid;
r = mid - 1;
}else
l = mid + 1;
}
ans = max(ans,1LL*(pr-pl+1)*len);
} int main() {
//freopen("F:\\rush.txt","r",stdin);
scanf("%s", s);
n = strlen(s); s[n] = 0;
build_Sa(n + 1, MAX_CHAR);//注意调用n+1
getHeight(n);
ST.init(n);
ST.getst(Height,n); int len = 0;
S[len++] = '$';
S[len++] = '#';
for (int i = 0;i < n;i++)
{
S[len++] = s[i];
belong[len-1] = i;
S[len++] = '#';
} int mx = 0,id = 0;
for (int i = 1;i < len;i++)
{
if (i > mx)
p[i] = 1;
else
p[i] = min(p[2*id-i],mx-i+1);
while (S[i-p[i]]==S[i+p[i]])
{
if (i+p[i] > mx && S[i-p[i]+1]!='#')
{
cal(belong[i-p[i]+1],belong[i+p[i]-1]);
}
p[i]++;
}
if (i + p[i]-1 > mx) mx = i + p[i] - 1,id = i;
}
printf("%lld\n",ans);
return 0;
}

【BZOJ 3676】[Apio2014]回文串的更多相关文章

  1. BZOJ 3676: [Apio2014]回文串

    3676: [Apio2014]回文串 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 2013  Solved: 863[Submit][Status ...

  2. bzoj 3676: [Apio2014]回文串 回文自动机

    3676: [Apio2014]回文串 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 844  Solved: 331[Submit][Status] ...

  3. 字符串(马拉车算法,后缀数组,稀疏表):BZOJ 3676 [Apio2014]回文串

    Description 考虑一个只包含小写拉丁字母的字符串s.我们定义s的一个子串t的“出 现值”为t在s中的出现次数乘以t的长度.请你求出s的所有回文子串中的最 大出现值. Input 输入只有一行 ...

  4. ●BZOJ 3676 [Apio2014]回文串

    题链: http://www.lydsy.com/JudgeOnline/problem.php?id=3676 题解: 后缀数组,Manacher,二分 首先有一个结论:一个串的本质不同的回文串的个 ...

  5. BZOJ 3676 [Apio2014]回文串(回文树)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3676 [题目大意] 考虑一个只包含小写拉丁字母的字符串s. 我们定义s的一个子串t的& ...

  6. bzoj 3676: [Apio2014]回文串【回文自动机】

    回文自动机板子 或者是SAM+manacher+倍增,就是manacher求本质不同回文串(让f++的串),然后在SAM倍增查询对应点出现次数 #include<iostream> #in ...

  7. BZOJ 3676 [Apio2014]回文串 (后缀自动机+manacher/回文自动机)

    题目大意: 给你一个字符串,求其中回文子串的长度*出现次数的最大值 明明是PAM裸题我干嘛要用SAM做 回文子串有一个神奇的性质,一个字符串本质不同的回文子串个数是$O(n)$级别的 用$manach ...

  8. bzoj 3676 [Apio2014]回文串(Manacher+SAM)

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=3676 [题意] 给定一个字符串,定义一个串的权值为长度*出现次数,求最大权的回文子串. ...

  9. BZOJ.3676.[APIO2014]回文串(回文树)

    BZOJ 洛谷 很久之前写(抄)过一个Hash+Manacher的做法,当时十分懵逼=-= 然而是道回文树模板题. 回文树教程可以看这里(真的挺妙的). 顺便再放上MilkyWay的笔记~ //351 ...

  10. bzoj 3676: [Apio2014]回文串【后缀自动机+manacher】

    用manacher找出本质不同的回文子串放在SAM上跑 #include<iostream> #include<cstdio> #include<cstring> ...

随机推荐

  1. 驱动学习3-make

    在向内核中添加驱动的时候要完成3项工作 (1)在Kconfig中添加新代码对应项目的编译条件(下面Makefile文件中需要用到它定义的的宏变量) (2)将驱动源码添加到对应的目录中 (3)在目录Ma ...

  2. Android, IOS 史上最强多语言国际化,不仅第一次会尾随系统,并且会保存用户的语言设置

    劲爆消息,我提供源代码了.你能够先看完再下载.也能够先下载再看完, android源代码地址: https://github.com/hebiao6446/------Bluetooth-Androi ...

  3. Android LoaderManager与CursorLoader用法

    一.基本概念 1.LoaderManager LoaderManager用来负责管理与Activity或者Fragment联系起来的一个或多个Loaders对象. 每个Activity或者Fragme ...

  4. JavaScript系列--JavaScript数组高阶函数reduce()方法详解及奇淫技巧

    一.前言 reduce() 方法接收一个函数作为累加器,数组中的每个值(从左到右)开始缩减,最终计算为一个值. reduce() 可以作为一个高阶函数,用于函数的 compose. reduce()方 ...

  5. JSONP的使用示例(以及jquery版jsonp)超简单

    前言: 平时工作中很少跨域,很少用到jsonp,但是几乎每次面试都会被问到这个问题.很崩溃. 菜鸟教程上的jsonp教程就很好.这里做个笔记,自己捋一遍. Jsonp(JSON with Paddin ...

  6. UVA 11039 - Building designing 水题哇~

    水题一题,按绝对值排序后扫描一片数组(判断是否异号,我是直接相乘注意中间值越界)即可. 感觉是让我练习sort自定义比较函数的. #include<cstdio> #include< ...

  7. mysql :Native table 'performance_schema'.'cond_instances' has the wrong structure

    err: 150418 13:25:06 [ERROR] Native table 'performance_schema'.'cond_instances' has the wrong struct ...

  8. [Angular2 Router] Index router

    Index router as default router. import {RouterModule} from "@angular/router"; import {NotF ...

  9. UVA 11732 - strcmp() Anyone? 字典树

    传送门:http://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&a ...

  10. 【b702】字符串的展开

    Time Limit: 1 second Memory Limit: 50 MB [问题描述] 在初赛普及组的"阅读程序写结果"的问题中,我们曾给出一个字符串展开的例子:如果在输入 ...