来自F allDream的博客,未经允许,请勿转载,谢谢。


如果一个字符串可以被拆分为 AABB 的形式,其中 A和 B是任意非空字符串,则我们称该字符串的这种拆分是优秀的。

例如,对于字符串 aabaabaa,如果令 A=aab,B=a,我们就找到了这个字符串拆分成 AABB的一种方式。

一个字符串可能没有优秀的拆分,也可能存在不止一种优秀的拆分。比如我们令 A=a,B=baa,也可以用 AABB表示出上述字符串;但是,字符串 abaabaa 就没有优秀的拆分。

现在给出一个长度为 n的字符串 S,我们需要求出,在它所有子串的所有拆分方式中,优秀拆分的总个数。这里的子串是指字符串中连续的一段。

以下事项需要注意:

1.出现在不同位置的相同子串,我们认为是不同的子串,它们的优秀拆分均会被记入答案。

2.在一个拆分中,允许出现 A=B。例如 cccc 存在拆分 A=B=c。

3.字符串本身也是它的一个子串

你需要计算一个字符串的所有子串的优秀的拆分的个数之和

T(T<=10)组数据,n<=30000

这道题n^2暴力95分..

正解的话反正我是想不到...

如果知道以每个点为左/右端点构成AA的方案数,就可以计算答案

考虑求出正着反着各求出后缀数组之后,枚举长度L,然后每段长度是L把这个字符串分段

每次枚举两个端点i,j(i和j都是L的倍数,j=i+L),并且计算左端点在[i-L+1,i]内的答案。

所以求一下点i,j正着的lcp,反着的lcp,也就是从点i开始能向左右扩展到哪里,如果能扩展到的长度大等于L,那么就有答案

假设倒着时候的lcp是a,正着的是b,也就是[i-a+1,i]和[j-a+1,j]相同(显然a>0),b同理,那么有答案的左端点区间是[i-a+1,i-a+1 + (a+b-1-L) -1] 差分即可。

复杂度Tnlogn

#include<iostream>
#include<cstdio>
#include<cstring>
#define ll long long
#define MN 30000
#define MD 15
using namespace std;
inline int read()
{
int x = ; char ch = getchar();
while(ch < '' || ch > '') ch = getchar();
while(ch >= '' && ch <= ''){x = x * + ch - '';ch = getchar();}
return x;
} int n,k,lt[MN+],rt[MN+],v[MN+],Log[MN+];ll ans;
char st[MN*+],st2[MN*+]; void CalSa(int*SA,int*RK,int*sa,int*rk)
{
for(int i=;i<=n;++i) v[rk[sa[i]]]=i;
for(int i=n;i;--i)
if(sa[i]>k) SA[v[rk[sa[i]-k]]--]=sa[i]-k;
for(int i=n-k+;i<=n;++i) SA[v[rk[i]]--]=i;
for(int i=;i<=n;++i)
RK[SA[i]]=RK[SA[i-]]+(rk[SA[i]]!=rk[SA[i-]]||rk[SA[i]+k]!=rk[SA[i-]+k]);
} void GetH(int*H,int*sa,int*rk,char*s)
{
for(int i=,k=;i<=n;H[rk[i++]-]=k,k?--k:) if(rk[i]>)
for(int j=sa[rk[i]-];s[j+k]==s[i+k];++k); } struct SuffixArray
{
int p,q,sa[][MN+],rk[][MN+],f[MD+][MN+],F[MD+][MN+],fa[MD+][MN+],Fa[MD+][MN+];
inline void Clear(){p=;q=;
memset(rk,,sizeof(rk));
memset(sa,,sizeof(sa));
}
void Build(char *s)
{
memset(v,,sizeof(int)*);
for(int i=;i<=n;++i)++v[s[i]-'a'];
for(int i=;i<;++i) v[i]+=v[i-];
for(int i=;i<=n;++i) sa[q][v[s[i]-'a']--]=i;
for(int i=;i<=n;++i)
rk[q][sa[q][i]]=rk[q][sa[q][i-]]+(s[sa[q][i]]!=s[sa[q][i-]]);
for(k=;k<n;k<<=)
{
CalSa(sa[p],rk[p],sa[q],rk[q]);
swap(p,q);
}
}
void BuildHeight(char *s)
{
GetH(f[],sa[q],rk[q],s);
for(int i=;i<=n;++i) F[][i]=f[][i-];
for(int i=;i<=MD;++i) F[i][]=f[i][]=f[i][n]=F[i][]=MN;
for(int i=;i<n;++i) fa[][i]=i+,Fa[][i+]=i;
for(int i=;i<=MD;++i)
for(int j=;j<=n;++j)
fa[i][j]=fa[i-][fa[i-][j]],
Fa[i][j]=Fa[i-][Fa[i-][j]],
f[i][j]=min(f[i-][j],f[i-][fa[i-][j]]),
F[i][j]=min(F[i-][j],F[i-][Fa[i-][j]]);
}
int query(int i,int j)
{
i=rk[q][i],j=rk[q][j];
if(i>j) swap(i,j);
int L=Log[j-i];
return min(f[L][i],F[L][j]);
}
}P,S; int main()
{
Log[]=-;
for(int i=;i<=MN;++i) Log[i]=Log[i>>]+;
for(int T=read();T;--T)
{
scanf("%s",st+);n=strlen(st+);
for(int i=;i<=n;++i) st2[i]=st[n-i+];
P.Clear();S.Clear();ans=;
P.Build(st);S.Build(st2);
P.BuildHeight(st);
S.BuildHeight(st2);
memset(lt,,sizeof(lt));
memset(rt,,sizeof(rt));
for(int L=;L<=n;++L)
for(int i=L,j=i+L;j<=n;i=j,j+=L)
{
int a=P.query(i,j),b=S.query(n-j+,n-i+);
if(min(a,L)+min(b,L)->=L)
{
int len=min(a,L)+min(b,L)--L;
++lt[i-min(b,L)+],--lt[i-min(b,L)+len+],
++rt[j-min(b,L)+L],--rt[j-min(b,L)+L+len+];
}
}
for(int i=;i<=n;++i) lt[i]+=lt[i-],rt[i]+=rt[i-];
for(int i=;i<n;++i) ans+=1LL*rt[i]*lt[i+];
printf("%lld\n",ans);
}
return ;
}

[Noi2016]优秀的拆分的更多相关文章

  1. [NOI2016]优秀的拆分&&BZOJ2119股市的预测

    [NOI2016]优秀的拆分 https://www.lydsy.com/JudgeOnline/problem.php?id=4650 题解 如果我们能够统计出一个数组a,一个数组b,a[i]表示以 ...

  2. luogu1117 [NOI2016]优秀的拆分

    luogu1117 [NOI2016]优秀的拆分 https://www.luogu.org/problemnew/show/P1117 后缀数组我忘了. 此题哈希可解决95分(= =) 设\(l_i ...

  3. 【BZOJ4560】[NOI2016]优秀的拆分

    [BZOJ4560][NOI2016]优秀的拆分 题面 bzoj 洛谷 题解 考虑一个形如\(AABB\)的串是由两个形如\(AA\)的串拼起来的 那么我们设 \(f[i]\):以位置\(i\)为结尾 ...

  4. [UOJ#219][BZOJ4650][Noi2016]优秀的拆分

    [UOJ#219][BZOJ4650][Noi2016]优秀的拆分 试题描述 如果一个字符串可以被拆分为 AABBAABB 的形式,其中 A 和 B 是任意非空字符串,则我们称该字符串的这种拆分是优秀 ...

  5. [NOI2016]优秀的拆分(SA数组)

    [NOI2016]优秀的拆分 题目描述 如果一个字符串可以被拆分为 \(AABB\) 的形式,其中 A和 B是任意非空字符串,则我们称该字符串的这种拆分是优秀的. 例如,对于字符串 \(aabaaba ...

  6. 题解-NOI2016 优秀的拆分

    NOI2016 优秀的拆分 \(T\) 组测试数据.求字符串 \(s\) 的所有子串拆成 \(AABB\) 形式的方案总和. 数据范围:\(1\le T\le 10\),\(1\le n\le 3\c ...

  7. [NOI2016]优秀的拆分 后缀数组

    题面:洛谷 题解: 因为对于原串的每个长度不一定等于len的拆分而言,如果合法,它将只会被对应的子串统计贡献. 所以子串这个限制相当于是没有的. 所以我们只需要对于每个位置i求出f[i]表示以i为开头 ...

  8. [BZOJ]4650: [Noi2016]优秀的拆分

    Time Limit: 30 Sec  Memory Limit: 512 MB Description 如果一个字符串可以被拆分为 AABBAABB 的形式,其中 AA 和 BB 是任意非空字符串, ...

  9. 【刷题】BZOJ 4650 [Noi2016]优秀的拆分

    Description 如果一个字符串可以被拆分为 AABBAABB 的形式,其中 AA 和 BB 是任意非空字符串,则我们称该字符串的这种拆分是优秀的.例如,对于字符串 aabaabaa,如果令 A ...

随机推荐

  1. verilog学习笔记(1)_两个小module

    第一个小module-ex_module module ex_module( input wire sclk,//声明模块的时候input变量一定是wire变量 input wire rst_n,// ...

  2. Android 4.4 沉浸式透明状态栏

    原文链接:http://www.bkjia.com/Androidjc/913061.html 第一种方法 这里写代码片第一种方法,在代码设置: if(VERSION.SDK_INT >= VE ...

  3. LeetCode & Q1-Two Sum-Easy

    Array Hash Table Question Given an array of integers, return indices of the two numbers such that th ...

  4. c 语言typedef 和 define的使用和区别

    #define是C的指令,用于为各种数据类型定义别名,与typedef 类似,但是有一下几点不同 1,typedef仅限于为类型定义符号名称,而#define不仅可以为类型定义符号名称,也能为数值定义 ...

  5. Groovy入门(2-2)Groovy的eclipse插件安装

    1.安装eclipse插件 启动eclipse,点击help -> Install New Software... 在弹出的窗口中点击:Add... Groovy插件的地址:http://dis ...

  6. HTTP协议扫盲(三)HTTP协议的请求头列表和分类描述

    一.请求报头和响应报头列表 1.Requests 头列表 Header 解释 示例 Accept 指定客户端能够接收的内容类型 Accept: text/plain, text/html Accept ...

  7. Spring Security入门(3-4)Spring Security 异常处理、异常传递和异常获取

  8. 记一次将公司网站http换成https

    看了博客园将近一年了,一直都只是在被动的看,总觉得应该写点什么,但是又不知道该写点什么.今天正好公司要把网站由http换成https,那我就顺便记录一下吧. 由于之前没有弄过,所以就面向百度编程. 首 ...

  9. express实践(一)

    涉及以下这些内容: 主体. cookie.session 数据 模板引擎 服务器基本结构: const express=require('express'); const static=require ...

  10. [Luogu1342] 请柬 - 最短路模板

    Description 在电视时代,没有多少人观看戏剧表演.Malidinesia古董喜剧演员意识到这一事实,他们想宣传剧院,尤其是古色古香的喜剧片.他们已经打印请帖和所有必要的信息和计划.许多学生被 ...