Codeforces 356E - Xenia and String Problem(哈希)
首先显然一个 gray 串的长度只可能是 \(2^k-1\),其中 \(k\in\mathbb{Z}\)。
考虑将一个字符改成另外一个字符会对答案产生什么样的影响。显然原来包含这个字符的 gray 串的贡献都应被减掉,新出现的 gray 串的贡献又应被加上。我们分别计算两部分的贡献。
对于第一部分,由于 gray 串最多只可能有 \(n\log n\),因此我们逐一枚举它们并差分一下它们贡献。那么如何检验一个长度为 \(2^k-1\) 的字符串是否符合要求呢?考虑倍增。设 \(is_{i,j}\) 表示以 \(s_j\) 开头的长度为 \(2^i-1\) 的字符串是否符合要求,那么 \(is_{i,j}\) 为 \(1\) 当且仅当 \(is_{i-1,j},is_{i-1,j+2^{i-1}}\) 均为 \(1\),并且 \(s[j...j+2^{i-1}-2]=s[j+2^{i-1}...j+2^i-2]\),且 \(s_{j+2^{i-1}-1}\) 在 \(s[j...j+2^i-2]\) 中只出现过一次。这个很好判断。
对于第二部分,我们考虑一个长度为 \(2^{k}-1\) 字符串在修改哪些字符的情况下会变成 gray 串。分三种情况:改左边、改中间、改右边。第一种和第三种本质上是相同的,因此这里只考虑第一种和第二种。如果修改中间字符,那么要求左右两边的子串相同并且都是 gray 串,并且修改过的字符在子串中只出现一次。如果修改左边字符串中的字符,那么修改方式显然是唯一的——就是左右两部分唯一不同的字符串,可以通过二分+hash 求出 LCP 和 LCS 后知道要修改哪个字符,然后再检验下改过之后的字符串是否符合要求即可。
总之是一道不算太难的哈希练手题。
const int MAXN=1e5;
const int LOG_N=17;
int n;char s[MAXN+5];
ll d[MAXN+5],add[MAXN+5][28];
struct hsh{
int BS,MOD;
int hs[MAXN+5],pw[MAXN+5];
void init(){
for(int i=(pw[0]=1);i<=n;i++) pw[i]=1ll*pw[i-1]*BS%MOD;
for(int i=1;i<=n;i++) hs[i]=(1ll*hs[i-1]*BS+s[i])%MOD;
}
int gethash(int l,int r){return (hs[r]-1ll*pw[r-l+1]*hs[l-1]%MOD+MOD)%MOD;}
} h1,h2;
bool check(int l1,int r1,int l2,int r2){
return h1.gethash(l1,r1)==h1.gethash(l2,r2)&&
h2.gethash(l1,r1)==h2.gethash(l2,r2);
}
int occ[MAXN+5][28];
bool is[LOG_N+2][MAXN+5];
int getocc(int c,int l,int r){return occ[r][c]-occ[l-1][c];}
int getlcp(int x,int y){
int l=1,r=n-max(x,y)+1,p=0;
while(l<=r){
int mid=l+r>>1;
if(check(x,x+mid-1,y,y+mid-1)) p=mid,l=mid+1;
else r=mid-1;
} return p;
}
int getlcs(int x,int y){
int l=1,r=min(x,y),p=0;
while(l<=r){
int mid=l+r>>1;
if(check(x-mid+1,x,y-mid+1,y)) p=mid,l=mid+1;
else r=mid-1;
} return p;
}
int main(){
scanf("%s",s+1);n=strlen(s+1);
h1.BS=191;h1.MOD=998244853;h2.BS=193;h2.MOD=1004535809;
h1.init();h2.init();
for(int i=1;i<=n;i++){
for(int j=0;j<26;j++) occ[i][j]=occ[i-1][j];
occ[i][s[i]-'a']++;
}
for(int i=1;i<=n;i++) is[1][i]=1;ll sum=n;
for(int i=2;i<=LOG_N;i++) for(int j=1;j+(1<<i)-2<=n;j++){
is[i][j]=is[i-1][j]&is[i-1][j+(1<<i-1)];
if(is[i][j]){
is[i][j]&=check(j,j+(1<<i-1)-2,j+(1<<i-1),j+(1<<i)-2);
is[i][j]&=(getocc(s[j+(1<<i-1)-1]-'a',j,j+(1<<i)-2)==1);
if(is[i][j]){
int len=(1<<i)-1;sum+=1ll*len*len;
d[j]+=1ll*len*len;d[j+(1<<i)-1]-=1ll*len*len;
}
}
// printf("%d %d %d\n",i,j,is[i][j]);
}
for(int i=1;i<=n;i++) d[i]+=d[i-1];
ll dlt=0;
for(int i=2;i<=LOG_N;i++) for(int j=1;j+(1<<i)-2<=n;j++){
if(!is[i-1][j]&&!is[i-1][j+(1<<i-1)]) continue;
int len=(1<<i)-1;ll bnf=1ll*len*len;
if(is[i-1][j]&&is[i-1][j+(1<<i-1)]){
if(check(j,j+(1<<i-1)-2,j+(1<<i-1),j+(1<<i)-2)){
for(int c=0;c<26;c++) if(!getocc(c,j,j+(1<<i)-2)){
add[j+(1<<i-1)-1][c]+=bnf;
}
}
} if(is[i-1][j+(1<<i-1)]){
int lcp=getlcp(j,j+(1<<i-1)),lcs=getlcs(j+(1<<i-1)-2,j+(1<<i)-2);
if(lcp+lcs==(1<<i-1)-2){
int ori=s[j+lcp]-'a',nd=s[j+(1<<i-1)+lcp]-'a';
if(getocc(s[j+(1<<i-1)-1]-'a',j,j+(1<<i)-2)-(ori==s[j+(1<<i-1)-1]-'a')+(nd==s[j+(1<<i-1)-1]-'a')!=1);
else add[j+lcp][nd]+=bnf;
}
} if(is[i-1][j]){
int lcp=getlcp(j,j+(1<<i-1)),lcs=getlcs(j+(1<<i-1)-2,j+(1<<i)-2);
if(lcp+lcs==(1<<i-1)-2){
int ori=s[j+(1<<i-1)+lcp]-'a',nd=s[j+lcp]-'a';
if(getocc(s[j+(1<<i-1)-1]-'a',j,j+(1<<i)-2)-(ori==s[j+(1<<i-1)-1]-'a')+(nd==s[j+(1<<i-1)-1]-'a')!=1);
else add[j+(1<<i-1)+lcp][nd]+=bnf;
}
}
}
// printf("%lld\n",sum);
// for(int i=1;i<=n;i++) printf("%lld%c",d[i]," \n"[i==n]);
for(int i=1;i<=n;i++) for(int j=0;j<26;j++)
if(s[i]-'a'!=j) chkmax(dlt,-d[i]+add[i][j]);
printf("%lld\n",sum+dlt);
return 0;
}
Codeforces 356E - Xenia and String Problem(哈希)的更多相关文章
- CF356E - Xenia and String Problem
也许更好的阅读体验 \(\mathcal{Description}\) 定义一种字符串\(gray\)串满足: 长度为奇数 正中间的字母只出现一次 左右两端相同,左右两端也是gray串 一个\(gra ...
- Codeforces #541 (Div2) - E. String Multiplication(动态规划)
Problem Codeforces #541 (Div2) - E. String Multiplication Time Limit: 2000 mSec Problem Descriptio ...
- Mediocre String Problem (2018南京M,回文+LCP 3×3=9种做法 %%%千年好题 感谢"Grunt"大佬的细心讲解)
layout: post title: Mediocre String Problem (2018南京M,回文+LCP 3×3=9种做法 %%%千年好题 感谢"Grunt"大佬的细 ...
- hdu String Problem(最小表示法入门题)
hdu 3374 String Problem 最小表示法 view code#include <iostream> #include <cstdio> #include &l ...
- HDU 3374 String Problem(KMP+最大/最小表示)
String Problem Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) T ...
- 【HDU3374】 String Problem (最小最大表示法+KMP)
String Problem Description Give you a string with length N, you can generate N strings by left shift ...
- HDOJ3374 String Problem 【KMP】+【最小表示法】
String Problem Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) T ...
- HDU 3374 String Problem (KMP+最大最小表示)
HDU 3374 String Problem (KMP+最大最小表示) String Problem Time Limit: 2000/1000 MS (Java/Others) Memory ...
- String Problem hdu 3374 最小表示法加KMP的next数组
String Problem Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)To ...
随机推荐
- NXOpen.CAM.CAMSetup.CopyObjects的使用
复制CAM对象 Public Function CopyObjects(ByVal view As NXOpen.CAM.CAMSetup.View, ByVal objectsToBeMoved() ...
- Visual Studio CMake 项目和 WSL
Visual Studio CMake 项目和 WSL https://devblogs.microsoft.com/cppblog/c-with-visual-studio-2019-and-win ...
- 欧姆龙plc通讯协议格式
欧姆龙CPM1A型plc与上位计算机通信的顺序是上位机先发出命令信息给PLC,PLC返回响应信息给上位 机.每次通信发送/接受的一组数据称为一"帧".帧由少于131个字符的数据构成 ...
- AIApe问答机器人Scrum Meeting 5.5&5.6
Scrum Meeting 7 日期:2021年5月5日&2021年5月6日 会议主要内容概述:汇报近日工作,确定下一步计划,放弃"关键词提取"和"改进关键词&q ...
- [Beta]the Agiles Scrum Meeting 11
会议时间:2020.5.26 21:00 1.每个人的工作 今天已完成的工作 成员 已完成的工作 issue yjy 帮助解决技术问题 tq 完成评测机新增评测指标 评测部分增加更多评测指标 wjx ...
- Java版流媒体编解码和图像处理(JavaCPP+FFmpeg)
欢迎访问我的GitHub https://github.com/zq2599/blog_demos 内容:所有原创文章分类汇总及配套源码,涉及Java.Docker.Kubernetes.DevOPS ...
- Qt坐标转换系统的理解
转 https://blog.csdn.net/hgcprg/article/details/53537106 今天又看了一篇对Qt坐标转换系统以及QTransform的博客,作者讲的非常透彻,链接如 ...
- 算法:Z字型(Zigzag)编排
问题:给定 n 行和 m 列的二维数组矩阵.如图所示,以 ZIG-ZAG 方式打印此矩阵. 从对称的角度来看,通过反复施加滑行反射可以从简单的图案如线段产生规则的之字形. 主要思想:算法从(0, 0) ...
- 【Azure 应用服务】App Service For Linux 部署Java Spring Boot应用后,查看日志文件时的疑惑
编写Java Spring Boot应用,通过配置logging.path路径把日志输出在指定的文件夹中. 第一步:通过VS Code创建一个空的Spring Boot项目 第二步:在applicat ...
- 贪心-Saruman‘s Army POJ - 3069
万恶之源 目录 题意 思路 贪心的原则是什么呢? 错解 正解 代码实现 书上的代码 我的代码 比较一下 问题 题意 给定若干个点的坐标,与范围R.每个点可以选择是否标记,标记后这个点的左右范围R内的所 ...