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

就像这样,但如果是上述情况的话,是无法拼成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 ;
}
股市的预测
#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股市的预测的更多相关文章
- bzoj 4650(洛谷 1117) [Noi2016]优秀的拆分——枚举长度的关键点+后缀数组
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=4650 https://www.luogu.org/problemnew/show/P1117 ...
- luogu1117 [NOI2016]优秀的拆分
luogu1117 [NOI2016]优秀的拆分 https://www.luogu.org/problemnew/show/P1117 后缀数组我忘了. 此题哈希可解决95分(= =) 设\(l_i ...
- 【BZOJ4560】[NOI2016]优秀的拆分
[BZOJ4560][NOI2016]优秀的拆分 题面 bzoj 洛谷 题解 考虑一个形如\(AABB\)的串是由两个形如\(AA\)的串拼起来的 那么我们设 \(f[i]\):以位置\(i\)为结尾 ...
- [UOJ#219][BZOJ4650][Noi2016]优秀的拆分
[UOJ#219][BZOJ4650][Noi2016]优秀的拆分 试题描述 如果一个字符串可以被拆分为 AABBAABB 的形式,其中 A 和 B 是任意非空字符串,则我们称该字符串的这种拆分是优秀 ...
- [NOI2016]优秀的拆分(SA数组)
[NOI2016]优秀的拆分 题目描述 如果一个字符串可以被拆分为 \(AABB\) 的形式,其中 A和 B是任意非空字符串,则我们称该字符串的这种拆分是优秀的. 例如,对于字符串 \(aabaaba ...
- 题解-NOI2016 优秀的拆分
NOI2016 优秀的拆分 \(T\) 组测试数据.求字符串 \(s\) 的所有子串拆成 \(AABB\) 形式的方案总和. 数据范围:\(1\le T\le 10\),\(1\le n\le 3\c ...
- [NOI2016]优秀的拆分 后缀数组
题面:洛谷 题解: 因为对于原串的每个长度不一定等于len的拆分而言,如果合法,它将只会被对应的子串统计贡献. 所以子串这个限制相当于是没有的. 所以我们只需要对于每个位置i求出f[i]表示以i为开头 ...
- [BZOJ]4650: [Noi2016]优秀的拆分
Time Limit: 30 Sec Memory Limit: 512 MB Description 如果一个字符串可以被拆分为 AABBAABB 的形式,其中 AA 和 BB 是任意非空字符串, ...
- [Noi2016]优秀的拆分
来自F allDream的博客,未经允许,请勿转载,谢谢. 如果一个字符串可以被拆分为 AABB 的形式,其中 A和 B是任意非空字符串,则我们称该字符串的这种拆分是优秀的. 例如,对于字符串 aab ...
随机推荐
- [转帖]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 ...
- [转帖]firewall-cmd
firewall-cmd https://wangchujiang.com/linux-command/c/firewall-cmd.html 高手大作 等哪天需要防火墙了 再练习一下. Linux上 ...
- 极验3.0滑动拼图验证的使用--java
[ 前言: 在登录其他网站的时候,看到有个滑动拼图的验证觉得挺好玩的,以前做一个图片验证的小demo,现在发现很多网站都开始流行滑动拼图的验证了,今天也想自己动手来弄一个. 废话不多说,开始撸起来! ...
- fatal: HttpRequestException encountered解决方法
最近在windows下git push提交就会弹出如下错误: 网上查了一下发现是Github 禁用了TLS v1.0 and v1.1,必须更新Windows的git凭证管理器,才行. https:/ ...
- JSON Support in PostgreSQL and Entity Framework
JSON 和JSONB的区别(What's difference between JSON and JSONB data type in PosgresSQL?) When should be use ...
- 手机端图像编辑上传-cropper
编辑头像,实现相册,照像功能,并能缩放裁剪功能,可自定义UI,引用'cropper.js', 'exif.js' /*初始化裁剪插件*/ var screenWidth = $(window).wid ...
- qtp 自动化测试--点滴 菜单没有了,有些控件运行时找不到
test项目页签下-没有了 菜单栏:file edit view insert 看不到了 1 解决:在startpage标签下-tool-option-点击 restore layout-确定 2 菜 ...
- Stream、FileStream、MemoryStream的区别
1.Stream:流,在msdn的定义:提供字节序列的一般性视图,Stream提供了读写流的方法是以字节的形式从流中读取内容.而我们经常会用到从字节流中读取文本或者写入文本,微软提供了StreamRe ...
- java中的几个概念
1.JDK(Java Development Kit ):编写Java程序的程序员使用的软件(它是编写java程序,使用到的工具包,为程序员提供一些已经封装好的 java 类库) 2.JRE(Java ...
- Mac下搭建PHP服务器
打开终端 1. 输入 sudo vi /etc/apache2/httpd.conf 2.把167-170的前面#去掉即加载下面几个模块 1.LoadModule alias_module libe ...