UOJ #219 BZOJ 4650 luogu P1117 [NOI2016]优秀的拆分 (后缀数组、ST表)
连NOI Day1T1都不会做。。。看了题解都写不出来还要抄Claris的代码。。
题目链接: (luogu)https://www.luogu.org/problemnew/show/P1117
(bzoj)https://www.lydsy.com/JudgeOnline/problem.php?id=4650
(uoj)http://uoj.ac/problem/219
题解:
\(f[i]\)表示以\(i\)结束的\(AA\)型子串个数,\(g[i]\)表示以\(i\)开始的\(AA\)型子串个数
怎么求\(f,g\)?
打破思维定势,谁说必须要一个一个求呢
分长度来求
枚举长度\(L\), 处理所有长度为\(2L\)的\(AA\)型子串对\(f\)和\(g\)的贡献
如果每隔\(L\)的长度放一个打点计时器,呸,关键点
那么任何\(AA\)型子串都会经过两个相邻关键点
首先肯定要满足这两个关键位置上的字符一样
在这个基础上求出往前往后最多多少个一样的
这个就转化成了LCP和LCS问题,并且两个相邻关键点对\(f,g\)数组的影响是区间+1,使用差分前缀和解决
然后推一推就行了,注意+1-1不要推错
时间复杂度为调和级数,\(O(n\log n)\)
代码
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#define llong long long
using namespace std;
const int N = 1<<15;
const int LGN = 15;
const int S = 26;
int log2[N+3];
struct SparseTable
{
int n;
int str[N+3];
int rk[N+3];
int tmp[N+3];
int height[N+3];
int h[N+3];
int sa[N+3];
int wb[N+3];
int mini[N+3][LGN+3];
void get_sa()
{
int *x = rk,*y = tmp;
for(int i=0; i<=S; i++) wb[i] = 0;
for(int i=1; i<=n; i++) wb[x[i]=str[i]]++;
for(int i=1; i<=S; i++) wb[i] += wb[i-1];
for(int i=n; i>=1; i--) sa[wb[x[i]]--] = i;
int s = S,p = 0;
for(int j=1; p<n; j<<=1)
{
p = 0;
for(int i=n-j+1; i<=n; i++) y[++p] = i;
for(int i=1; i<=n; i++) if(sa[i]>j) y[++p] = sa[i]-j;
for(int i=1; i<=s; i++) wb[i] = 0;
for(int i=1; i<=n; i++) wb[x[y[i]]]++;
for(int i=1; i<=s; i++) wb[i] += wb[i-1];
for(int i=n; i>=1; i--) sa[wb[x[y[i]]]--] = y[i];
swap(x,y);
p = 1; x[sa[1]] = 1;
for(int i=2; i<=n; i++) x[sa[i]] = (y[sa[i]]==y[sa[i-1]] && y[sa[i]+j]==y[sa[i-1]+j]) ? p : ++p;
s = p;
}
for(int i=1; i<=n; i++) rk[sa[i]] = i;
for(int i=1; i<=n; i++)
{
h[i] = h[i-1]==0 ? 0 : h[i-1]-1;
while(i+h[i]<=n && sa[rk[i-1]]+h[i]<=n && str[i+h[i]]==str[sa[rk[i]-1]+h[i]])
{
h[i]++;
}
}
for(int i=1; i<=n; i++) height[i] = h[sa[i]];
for(int i=1; i<=n; i++) mini[i][0] = height[i];
for(int j=1; j<=LGN; j++)
{
for(int i=1; i+(1<<j)-1<=n; i++)
{
mini[i][j] = min(mini[i][j-1],mini[i+(1<<j-1)][j-1]);
}
}
}
int querymin(int lb,int rb)
{
int g = log2[rb-lb+1];
return min(mini[lb][g],mini[rb-(1<<g)+1][g]);
}
int LCP(int x,int y)
{
if(x==y) return n-x+1;
if(rk[x]>rk[y]) swap(x,y);
return querymin(rk[x]+1,rk[y]);
}
void clear()
{
for(int i=1; i<=n; i++) str[i] = rk[i] = tmp[i] = height[i] = h[i] = sa[i] = wb[i] = 0;
for(int i=1; i<=n; i++)
{
for(int j=0; j<=LGN; j++)
{
mini[i][j] = 0;
}
}
}
} s1,s2;
llong f[N+3];
llong g[N+3];
char a[N+3];
int n;
int LCP(int x,int y) {return s1.LCP(x,y);}
int LCS(int x,int y) {return s2.LCP(n+1-x,n+1-y);}
void preprocess()
{
log2[1] = 0; for(int i=2; i<=N; i++) log2[i] = log2[i>>1]+1;
}
void clear()
{
s1.clear(); s2.clear();
for(int i=0; i<=n+1; i++) a[i] = 0,f[i] = g[i] = 0ll;
}
int main()
{
preprocess();
int T; scanf("%d",&T);
while(T--)
{
scanf("%s",a+1); n = strlen(a+1); for(int i=1; i<=n; i++) a[i]-=96;
for(int i=1; i<=n; i++) s1.str[i] = a[i]; s1.n = n;
s1.get_sa();
for(int i=1; i<=n; i++) s2.str[i] = a[n+1-i]; s2.n = n;
s2.get_sa();
for(int i=1; i+i<=n; i++)
{
for(int j=i+i; j<=n; j+=i)
{
if(a[j]==a[j-i])
{
int lb = j-LCS(j,j-i)+1,rb = j+LCP(j,j-i)-1;
lb = max(lb+i-1,j); rb = min(rb,j+i-1);
if(lb<=rb)
{
f[lb]++; f[rb+1]--;
g[lb-i-i+1]++; g[rb+1-i-i+1]--;
}
}
}
}
for(int i=1; i<=n; i++) f[i] += f[i-1],g[i] += g[i-1];
llong ans = 0ll;
for(int i=1; i<n; i++)
{
llong tmp = f[i]*g[i+1];
ans += tmp;
}
printf("%lld\n",ans);
clear();
}
return 0;
}
UOJ #219 BZOJ 4650 luogu P1117 [NOI2016]优秀的拆分 (后缀数组、ST表)的更多相关文章
- UOJ#219. 【NOI2016】优秀的拆分 [后缀数组 ST表]
#219. [NOI2016]优秀的拆分 题意:求有多少AABB样子的子串,拆分不同的同一个子串算多个 一开始一直想直接求,并不方便 然后看了一眼Claris的题解的第一行就有思路了 如果分开,求\( ...
- UOJ#219/BZOJ4650 [NOI2016]优秀的拆分 字符串 SA ST表
原文链接http://www.cnblogs.com/zhouzhendong/p/9025092.html 题目传送门 - UOJ#219 (推荐,题面清晰) 题目传送门 - BZOJ4650 题意 ...
- [NOI2016]优秀的拆分 后缀数组
题面:洛谷 题解: 因为对于原串的每个长度不一定等于len的拆分而言,如果合法,它将只会被对应的子串统计贡献. 所以子串这个限制相当于是没有的. 所以我们只需要对于每个位置i求出f[i]表示以i为开头 ...
- BZOJ.4650.[NOI2016]优秀的拆分(后缀数组 思路)
BZOJ 洛谷 令\(st[i]\)表示以\(i\)为开头有多少个\(AA\)这样的子串,\(ed[i]\)表示以\(i\)结尾有多少个\(AA\)这样的子串.那么\(Ans=\sum_{i=1}^{ ...
- Luogu P1117 [NOI2016]优秀的拆分
题目链接 \(Click\) \(Here\) 这题质量不错,就是暴力分有点足\(hhhhhhhh\),整整有\(95\)分. (搞得我写完暴力都不想写正解直接理解思路之后就直接水过去了\(QwQ\) ...
- BZOJ 4650 [Noi2016]优秀的拆分 ——后缀数组
我们只需要统计在某一个点开始的形如$AA$字符串个数,和结束的个数相乘求和. 首先枚举循环节的长度L.即$\mid (A) \mid=L$ 然后肯定会经过s[i]和[i+L]至少两个点. 然后我们可以 ...
- BZOJ 4453: cys就是要拿英魂![后缀数组 ST表 单调栈类似物]
4453: cys就是要拿英魂! Time Limit: 3 Sec Memory Limit: 128 MBSubmit: 90 Solved: 46[Submit][Status][Discu ...
- bzoj 3796: Mushroom追妹纸【二分+后缀数组+st表】
把三个串加上ASCII大于z的分隔符连起来,然后求SA 显然每个相同子串都是一个后缀的前缀,所以枚举s1的每个后缀的最长和s2相同的前缀串(直接在排序后的数组里挨个找,最近的两个分别属于s1和s2的后 ...
- [LOJ 2083][UOJ 219][BZOJ 4650][NOI 2016]优秀的拆分
[LOJ 2083][UOJ 219][BZOJ 4650][NOI 2016]优秀的拆分 题意 给定一个字符串 \(S\), 求有多少种将 \(S\) 的子串拆分为形如 AABB 的拆分方案 \(| ...
随机推荐
- DBCP数据源使用
DBCP:DataBase Connection Pool 1.须要的jar:commons-dbcp.jar commons-pool.jar 2.把DBCP的配置文件(dbcpconfig.pr ...
- java 集合交并补
通过使用泛型方法和Set来表达数学中的表达式:集合的交并补.在下面三个方法中都将第一个參数Set复制了一份,并未直接改动參数中Set. package Set; import java.util.Ha ...
- c++ 基于Policy 的 模板编程
在没真正接触c++ 模板编程之前.真的没有想到c++ 还能够这么用.最大的感触是:太灵活了,太强大了. 最初接触模板威力还是在Delta3d中,感觉里面的模板使用实在是灵活与方便,特别是dtAI中使 ...
- 截取字符(substr)检索字符位置(instr)
1.SUBSTR(string,start_position,[length]) 求子字符串,返回字符串注释: string 元字符串start_position 开始位置(从0开始)length 可 ...
- Wall(凸包)
http://poj.org/problem?id=1113 题意:给出一些点的坐标,和一个半径r,求出这些点围成的凸包的周长再加上一个半径为r的圆的周长. #include <stdio.h& ...
- 14招搞定JavaScript调试
14招搞定JavaScript调试 译者按: 很多时候,大家可能只是依靠console.log来调试JavaScript代码,这样做的局限性不言而喻,这篇博客将教你几招实用的调试技巧. 原文: The ...
- go之结构体
一.关于结构体 简述 1.go 语言的切片可以存储同一类型的数据,但是结构体可以为不同项定义不同的数据类型 2.结构体是有一系列具有相同类型或不同类型的数据构成的数据集合 3.因为go 没有类似于类的 ...
- 强迫症!一行代码拿到url特定query的值
简单的说一下背景,看到小伙伴给我发的项目中有一段获取当前url特定query值的代码,本着能写1行代码就不写5行代码的原则,我把这个获取方法给改了一下 之前的代码如下: const queryArr ...
- python重定向原理及实例
1. 前言 为了在Python编程中, 利用控制台信息, 我们需要对控制台输出进行接管(重定向).在Python中,控制台输出的接口是sys.stdout,通过分析print与sys.stdout之间 ...
- JS 有趣的eval优化输入验证
//eval就是计算字符串[可以放任何js代码]里的值 . var str1='12+3'; eval(str1); . var str2='[1,2,3]'; eval(str2[]); .eval ...