[NOI2016]优秀的拆分
https://www.lydsy.com/JudgeOnline/problem.php?id=4650
题解
如果我们能够统计出一个数组a,一个数组b,a[i]表示以i为结尾的AA串个数,b[i]表示以i开头的AA串个数,我们就可以O(n)的统计答案了。
我们开看这么一个算法。
先枚举A的长度len,然后按照A的长度把序列分成很多段,每段有一个代表点即这一段的第一个点。
然后枚举所有代表点l,再枚举r=l+len,求一下(1,l)和(1,r)的最长公共后缀,(l,n)和(r,n)的最长公共前缀。

就像这样,但如果是上述情况的话,是无法拼成AA串的,因为蓝色和橙色并没有完全覆盖中间的红色段。

如果是这样的话就可以了,有多少个呢?蓝色和橙色组成的大的线段的左右端点就是第一条线段的左端点和第二条线段的右端点,这条线段可以在中间任意滑动。

然后就可以差分统计答案了。

这种算法相当于利用两个A的开头不会在同一个块内出现,只会在相邻两个块出现的性质做的。

注意,在求公共前后缀时要和和len取min,否则就会算多。

代码

#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#include<cstring>
#define N 30002
using namespace std;
int n,m,y[N],tong[N],T;
long long tag1[N],tag2[N],ans;
inline int rd(){
int x=;char c=getchar();bool f=;
while(!isdigit(c)){if(c=='-')f=;c=getchar();}
while(isdigit(c)){x=(x<<)+(x<<)+(c^);c=getchar();}
return f?-x:x;
}
struct suffixarray{
int rnk[N],sa[N],height[N],p[][N];
char s[N];
inline void qsort(){
for(int i=;i<=n;++i)tong[rnk[i]]++;
for(int i=;i<=m;++i)tong[i]+=tong[i-];
for(int i=n;i>=;--i)sa[tong[rnk[y[i]]]--]=y[i];
for(int i=;i<=m;++i)tong[i]=;
}
inline void SA(){
m=;
for(int i=;i<=n;++i)rnk[i]=s[i],y[i]=i;
qsort();
for(int w=,p=;p<n;m=p,w<<=){
p=;
for(int i=n-w+;i<=n;++i)y[++p]=i;
for(int i=;i<=n;++i)if(sa[i]>w)y[++p]=sa[i]-w;
qsort();swap(rnk,y);
rnk[sa[]]=p=;
for(int i=;i<=n;++i)
rnk[sa[i]]=((y[sa[i]]==y[sa[i-]])&&(y[sa[i]+w]==y[sa[i-]+w]))?p:++p;
}
for(int i=;i<=n;++i){
if(rnk[i]==)continue;
int j=max(,height[rnk[i-]]-);
while(s[i+j]==s[sa[rnk[i]-]+j])++j;
height[rnk[i]]=j;
p[][rnk[i]]=j;
}
for(int i=;(<<i)<=n;++i)
for(int j=;j+(<<i)-<=n;++j)
p[i][j]=min(p[i-][j],p[i-][j+(<<i-)]);
}
inline int query(int l,int r){
if(l>r)swap(l,r);l++;
if(l>r)return 2e9;
int lo=log2(r-l+);
return min(p[lo][l],p[lo][r-(<<lo)+]);
}
inline void clear(){
// memset(p,0,sizeof(p));
memset(height,,sizeof(height));
memset(s,,sizeof(s));
memset(rnk,,sizeof(rnk));
memset(sa,,sizeof(sa));
}
}z,f;
inline void init(){
memset(y,,sizeof(y));
ans=;
memset(tag1,,sizeof(tag1));memset(tag2,,sizeof(tag2));
}
int main(){
T=rd();
while(T--){
z.clear();f.clear();init();
scanf("%s",z.s+);n=strlen(z.s+);
for(int i=;i<=n;++i)f.s[i]=z.s[n-i+];f.s[n+]=z.s[n+]=;
f.SA();z.SA();
for(int i=;i<=n/;++i)
for(int j=;j+i<=n;j+=i){
int l=j,r=j+i;
int lcp=min(i,z.query(z.rnk[l],z.rnk[r])),lcs=min(i,f.query(f.rnk[n-l+],f.rnk[n-r+]));
if(lcp+lcs>=i+){
tag1[l-lcs+(i<<)]++;tag1[r+lcp]--;
tag2[l-lcs+]++;tag2[r+lcp-(i<<)+]--;
}
//cout<<l<<" "<<r<<" "<<lcp<<" "<<lcs<<endl;
}
for(int i=;i<=n;++i)tag1[i]+=tag1[i-],tag2[i]+=tag2[i-];
for(int i=;i<=n;++i)ans+=tag1[i-]*tag2[i];
cout<<ans<<endl;
} return ;
}

股市的预测

墨墨的妈妈热爱炒股,她要求墨墨为她编写一个软件,预测某只股票未来的走势。股票折线图是研究股票的必备工
具,它通过一张时间与股票的价位的函数图像清晰地展示了股票的走势情况。经过长时间的观测,墨墨发现很多股
票都有如下的规律:之前的走势很可能在短时间内重现!如图可以看到这只股票A部分的股价和C部分的股价的走势
如出一辙。通过这个观测,墨墨认为他可能找到了一个预测股票未来走势的方法。进一步的研究可是难住了墨墨,
他本想试图统计B部分的长度与发生这种情况的概率关系,不过由于数据量过于庞大,依赖人脑的力量难以完成,
于是墨墨找到了善于编程的你,请你帮他找一找给定重现的间隔(B部分的长度),有多少个时间段满足首尾部分
的走势完全相同呢?当然,首尾部分的长度不能为零。
题解
我们把序列差分之后,相当于是求形如ABA的串,这和上一题就差不多了。
每次取得l=i,r=i+m+len。然后就和上一题一样了。
每次新增加的答案为max(0,公共前缀+公共后缀-len)
代码
#include<iostream>
#include<cstdio>
#include<cmath>
#include<algorithm>
#define N 50002
using namespace std;
int n,m,y[N],ans,tong[N],M,b[N],c[N];
inline int rd(){
int x=;char c=getchar();bool f=;
while(!isdigit(c)){if(c=='-')f=;c=getchar();}
while(isdigit(c)){x=(x<<)+(x<<)+(c^);c=getchar();}
return f?-x:x;
}
struct suffixarray{
int rnk[N],sa[N],a[N],height[N],p[][N];
inline void qsort(){
for(int i=;i<=m;++i)tong[i]=;
for(int i=;i<=n;++i)tong[rnk[i]]++;
for(int i=;i<=m;++i)tong[i]+=tong[i-];
for(int i=n;i>=;--i)sa[tong[rnk[y[i]]]--]=y[i];
}
inline void SA(){
m=;
for(int i=;i<=n;++i)rnk[i]=a[i],y[i]=i;
qsort();
for(int w=,p=;p<n;m=p,w<<=){
p=;
for(int i=n-w+;i<=n;++i)y[++p]=i;
for(int i=;i<=n;++i)if(sa[i]>w)y[++p]=sa[i]-w;
qsort();swap(rnk,y);
rnk[sa[]]=p=;
for(int i=;i<=n;++i)
rnk[sa[i]]=((y[sa[i]]==y[sa[i-]])&&(y[sa[i]+w]==y[sa[i-]+w]))?p:++p;
}
for(int i=;i<=n;++i){
if(rnk[i]==)continue;
int j=max(,height[rnk[i-]]-);
while(a[i+j]==a[sa[rnk[i]-]+j])++j;
height[rnk[i]]=j;
p[][rnk[i]]=j;
}
for(int i=;(<<i)<=n;++i)
for(int j=;j+(<<i)-<=n;++j)
p[i][j]=min(p[i-][j],p[i-][j+(<<i-)]);
}
inline int query(int l,int r){
if(l>r)swap(l,r);l++;
if(l>r)return 2e9;
int lo=log2(r-l+);
return min(p[lo][l],p[lo][r-(<<lo)+]);
}
}z,f;
int main(){
n=rd();M=rd();
for(int i=;i<=n;++i)b[i]=rd();
for(int i=n-;i>=;--i)z.a[i]=b[i+]-b[i],c[i]=z.a[i];n--;
sort(c+,c+n+);int tt=unique(c+,c+n+)-c-;
for(int i=;i<=n;++i)z.a[i]=lower_bound(c+,c+tt+,z.a[i])-c;
for(int i=;i<=n;++i)f.a[n-i+]=z.a[i];
z.SA();f.SA();
for(int i=;i<=(n-M)/;++i){
for(int j=;j<=n&&j+i+M<=n;j+=i){
int l=j,r=i+j+M;
int lcp=min(i,z.query(z.rnk[l],z.rnk[r])),lcs=min(i,f.query(f.rnk[n-l+],f.rnk[n-r+]));
int x=lcs+lcp-i;
if(x>)ans+=x;
}
}
cout<<ans;
return ;
}

[NOI2016]优秀的拆分&&BZOJ2119股市的预测的更多相关文章

  1. bzoj 4650(洛谷 1117) [Noi2016]优秀的拆分——枚举长度的关键点+后缀数组

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4650 https://www.luogu.org/problemnew/show/P1117 ...

  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. [Noi2016]优秀的拆分

    来自F allDream的博客,未经允许,请勿转载,谢谢. 如果一个字符串可以被拆分为 AABB 的形式,其中 A和 B是任意非空字符串,则我们称该字符串的这种拆分是优秀的. 例如,对于字符串 aab ...

随机推荐

  1. [转帖]SAP S/4 HANA与SAP Business Suite/R3(ECC)的区别

    SAP S/4 HANA与SAP Business Suite/R3(ECC)的区别 https://blog.csdn.net/zhongguomao/article/details/5351520 ...

  2. [转帖]firewall-cmd

    firewall-cmd https://wangchujiang.com/linux-command/c/firewall-cmd.html 高手大作 等哪天需要防火墙了 再练习一下. Linux上 ...

  3. 极验3.0滑动拼图验证的使用--java

    [ 前言: 在登录其他网站的时候,看到有个滑动拼图的验证觉得挺好玩的,以前做一个图片验证的小demo,现在发现很多网站都开始流行滑动拼图的验证了,今天也想自己动手来弄一个. 废话不多说,开始撸起来! ...

  4. fatal: HttpRequestException encountered解决方法

    最近在windows下git push提交就会弹出如下错误: 网上查了一下发现是Github 禁用了TLS v1.0 and v1.1,必须更新Windows的git凭证管理器,才行. https:/ ...

  5. JSON Support in PostgreSQL and Entity Framework

    JSON 和JSONB的区别(What's difference between JSON and JSONB data type in PosgresSQL?) When should be use ...

  6. 手机端图像编辑上传-cropper

    编辑头像,实现相册,照像功能,并能缩放裁剪功能,可自定义UI,引用'cropper.js', 'exif.js' /*初始化裁剪插件*/ var screenWidth = $(window).wid ...

  7. qtp 自动化测试--点滴 菜单没有了,有些控件运行时找不到

    test项目页签下-没有了 菜单栏:file edit view insert 看不到了 1 解决:在startpage标签下-tool-option-点击 restore layout-确定 2 菜 ...

  8. Stream、FileStream、MemoryStream的区别

    1.Stream:流,在msdn的定义:提供字节序列的一般性视图,Stream提供了读写流的方法是以字节的形式从流中读取内容.而我们经常会用到从字节流中读取文本或者写入文本,微软提供了StreamRe ...

  9. java中的几个概念

    1.JDK(Java Development Kit ):编写Java程序的程序员使用的软件(它是编写java程序,使用到的工具包,为程序员提供一些已经封装好的 java 类库) 2.JRE(Java ...

  10. Mac下搭建PHP服务器

    打开终端 1. 输入  sudo vi /etc/apache2/httpd.conf 2.把167-170的前面#去掉即加载下面几个模块 1.LoadModule alias_module libe ...