NOIP2020 T2 字符串匹配题解
首先考虑O(n^3)的暴力怎么写。
显然,可以枚举字符串\(A\)+\(B\)的右端点,左端点显然是1,暴力判断是否能与后面的字符构成循环节,对于满足
\(k*(A+B)+C=\) 整个字符串\((k \in Z)\)
的情况暴力枚举\(A\),\(B\)分界点,对于\(L(x)\le R(x)\)的情况统计答案即可,\(L(x)\)为\(1\)到\(n\)出现奇数次的字符数量,\(R(x)\)为\(x\)到\(n\)出现奇数次的字符数量。
其实这样由于\(A+B\)很难找到循环节,均摊下来时间复杂度在随机数据下远小于\(O(n^{3})\),可以拿到\(48pts\)了,美滋滋。
把暴力贴上来吧:
#include<bits/stdc++.h>
#define LL long long
#define N (1<<20)+50
using namespace std;
LL ans;
int t,len;
char a[N];
int toa[30],toc[30],jia[N],jic[N];
inline int qr()
{
char a=0;int x=0,w=1;
while(a<'0'||a>'9'){if(a=='-')w=-1;a=getchar();}
while(a<='9'&&a>='0'){x=(x<<3)+(x<<1)+(a^48);a=getchar();}
return x*w;
}
int main()
{
t=qr();
while(t--)
{
scanf("%s",(a+1));
len=strlen(a+1);
memset(toa,0,sizeof(toa));
memset(toc,0,sizeof(toc));
memset(jic,0,sizeof(jic));
memset(jia,0,sizeof(jia));
ans=0;
for(register int i=1;i<=len;i++)//计算L(i)
{
toa[a[i]-'a']++;
toa[a[i]-'a']&1 ? jia[i]=jia[i-1]+1 :jia[i]=jia[i-1]-1;
}
for(register int i=len;i>=1;i--)//计算R(i)
{
toc[a[i]-'a']++;
toc[a[i]-'a']&1 ? jic[i]=jic[i+1]+1 :jic[i]=jic[i+1]-1;
}
for(register int i=2;i<len;i++)
{
int flag=1;
for(register int k=1;k<i;k++)//特判k=1的情况
if(jia[k]<=jic[i+1])
ans++;
for(register int j=2;j*i<len;j++)//枚举循环次数
{
int le=(j-1)*i+1;
int re=j*i;
for(register int k=le;k<=re;k++)//判断循环节
if(a[k]!=a[k-i])
{
flag=0;
break;
}
if(flag)//循环成立枚举A,B分界点统计次数
{
for(register int k=1;k<i;k++)
if(jia[k]<=jic[re+1])
ans++;
}
else
break;//循环不成立退出
}
}
printf("%lld\n",ans);
continue;
}
return 0;
}
下面考虑优化时间复杂度。
前置知识
- 字符串hash
没了
我们可以想到,时间主要浪费在以下两点:
- 在查找循环节时暴力匹配显然不是一个明智之举
- 查找\(A\),\(B\)分界时单一操作重复多次。
那么考虑优化:
显然,一个字符串出现奇数次的字符的数量\(\le\) \(26\) 废话,那么可以
维护一个类似前缀和的东西\(R\)(\(i\),\(j\))表示从\(1-i\)奇数次的字符的数量满足\(\le\) \(j\)的位置,用\(L\)(\(i\))表示从\(i-n\)奇数次的字符的数量。
于是,对可成立的A+B字符串判断复杂度变成了预处理\(O\)(\(27n\)),查询\(O\)(\(1\))。
其次,可以用字符串哈希或KMP判循环节,将判循环节复杂度优化为调和级数级即\(O\)(\(n\) \(\ln\) \(n\))。
加上这两个优化后复杂度大概在\(O( T(27n+n \ln n) )\)。
喜闻乐见的 Code :
#include<bits/stdc++.h>
#define LL long long
#define ULL unsigned long long
#define N (1<<20)+50
using namespace std;
LL ans;
int t,len;
char a[N];
int sm[N][27],sum[N];//分别对应题解中的L(i,j),R(i)
ULL Hash[N],jc[N],p=13331;//hash采取p进制,不过最好用双模数
//因为我太懒了,用了单模数哈希QAQ
inline int qr()//平平无奇的快读
{
char a=0;int x=0,w=1;
while(a<'0'||a>'9'){if(a=='-')w=-1;a=getchar();}
while(a<='9'&&a>='0'){x=(x<<3)+(x<<1)+(a^48);a=getchar();}
return x*w;
}
inline bool check(int l,int r,int L,int R)//字符串哈希判断相等
{
ULL hash1=(Hash[r]-Hash[l-1]*jc[r-l+1]);
ULL hash2=(Hash[R]-Hash[L-1]*jc[R-L+1]);
if(hash1==hash2)
return 1;
return 0;
}
int main()
{
t=qr();
jc[0]=1;
for(register int i=1;i<=(1<<20)+5;i++)
jc[i]=jc[i-1]*p;//处理p^i
while(t--)
{
scanf("%s",(a+1));
len=strlen(a+1);
ans=0;
int opk[27]={},val=0;
for(register int i=1;i<=len;i++)
Hash[i]=(Hash[i-1]*p+a[i]);//计算hash值
for(register int i=1;i<=len;i++)//预处理L(i,j)
{
for(register int j=0;j<=26;j++)
sm[i][j]=sm[i-1][j];//首先继承i-1时的前缀和
opk[(int)(a[i]-'a')+1]++;
opk[(int)(a[i]-'a')+1]&1?val++:val--;
sm[i][val]++;//加上统计的最新答案
}//这时sm(i,j)表示的仅为 1~i 中出现奇数次的字符数量=j的数量,于是需要再做一下前缀和
for(register int i=1;i<=len;i++)
for(register int j=1;j<=26;j++)
sm[i][j]+=sm[i][j-1];//现在sm(i,j)表示的就是 1~i 中出现奇数次的字符数量<=j的数量了
val=0;
memset(opk,0,sizeof(opk));
for(register int i=len;i>=1;i--)//预处理R(i,j)
{
opk[(int)(a[i]-'a')+1]++;
opk[(int)(a[i]-'a')+1]&1?val++:val--;
sum[i]=val;
}//没啥可说的
for(register int i=2;i<len;i++)
for(register int j=1;j*i<len;j++)//开始判循环节
{
int re=j*i+1;
if(check(1,i,i*(j-1)+1,j*i))//如果循环节成立
ans+=sm[i-1][sum[re]];//统计答案
else
break;//不成立直接退出即可
}
printf("%lld\n",ans);
continue;
}
return 0;
}
Update:
发现\(R\)(\(i\),\(j\))可以边进行计算边求,优化了一些常数(不开O2有96分,C++11就过了),代码如下:
#include<bits/stdc++.h>
#define LL long long
#define ULL unsigned long long
#define N (1<<20)+50
using namespace std;
LL ans;
int t,len;
char a[N];
int sm[27],sum[N];//分别对应题解中的L(i,j),R(i)
ULL Hash[N],jc[N],p=13331;//hash采取p进制,不过最好用双模数
//因为我太懒了,用了单模数哈希QAQ
inline int qr()//平平无奇的快读
{
char a=0;int x=0,w=1;
while(a<'0'||a>'9'){if(a=='-')w=-1;a=getchar();}
while(a<='9'&&a>='0'){x=(x<<3)+(x<<1)+(a^48);a=getchar();}
return x*w;
}
inline bool check(int l,int r,int L,int R)//字符串哈希判断相等
{
ULL hash1=(Hash[r]-Hash[l-1]*jc[r-l+1]);
ULL hash2=(Hash[R]-Hash[L-1]*jc[R-L+1]);
if(hash1==hash2)
return 1;
return 0;
}
int main()
{
t=qr();
jc[0]=1;
for(register int i=1;i<=(1<<20)+5;i++)
jc[i]=jc[i-1]*p;//处理p^i
while(t--)
{
scanf("%s",(a+1));
len=strlen(a+1);
ans=0;
int opk[27]={},val=0;
for(register int i=1;i<=len;i++)
Hash[i]=(Hash[i-1]*p+a[i]);//计算hash值
val=0;
for(register int i=len;i>=1;i--)//预处理R(i,j)
{
opk[(int)(a[i]-'a')+1]++;
opk[(int)(a[i]-'a')+1]&1?val++:val--;
sum[i]=val;
}//没啥可说的
val=0;
memset(sm,0,sizeof(sm));
memset(opk,0,sizeof(opk));
for(register int i=2;i<len;i++)
{
opk[(int)(a[i-1]-'a')+1]++;
opk[(int)(a[i-1]-'a')+1]&1?val++:val--;
for(register int j=val;j<=26;j++)
sm[j]++;//求R(i,j)
for(register int j=1;j*i<len;j++)//开始判循环节
{
int re=j*i+1;
if(check(1,i,i*(j-1)+1,j*i))//如果循环节成立
ans+=sm[sum[re]];//统计答案
else
break;//不成立直接退出即可
}
}
printf("%lld\n",ans);
continue;
}
return 0;
}
NOIP2020 T2 字符串匹配题解的更多相关文章
- KMP算法详解&&P3375 【模板】KMP字符串匹配题解
KMP算法详解: KMP算法是一种改进的字符串匹配算法,由D.E.Knuth,J.H.Morris和V.R.Pratt(雾)提出的. 对于字符串匹配问题(such as 问你在abababb中有多少个 ...
- 洛谷 P3375 【模板】KMP字符串匹配 题解
KMP模板,就不解释了 #include<iostream> #include<cstdio> #include<cstring> #include<algo ...
- 题解-洛谷P7114 字符串匹配
题面 洛谷P7114 字符串匹配 \(T\) 组测试数据.给定字符串 \(S\),问有多少不同的非空字符串 \(A\),\(B\),\(C\) 满足 \(S=ABABAB...ABC\) 且 \(A\ ...
- 字符串匹配之KMP---全力解析
近日,一同学面试被问到字符串匹配算法,结果因为他使用了暴力法,直接就跪了(如今想想这种面试官真的是不合格的,陈皓的一篇文章说的非常好,点击阅读).字符串匹配方法大概有:BF(暴力破解法), 简化版的B ...
- COJN 0558 800600带通配符的字符串匹配
800600带通配符的字符串匹配 难度级别:B: 运行时间限制:1000ms: 运行空间限制:51200KB: 代码长度限制:2000000B 试题描述 通配符是一类键盘字符,当我们不知道真正字符或者 ...
- HDU 5716 带可选字符的多字符串匹配(ShiftAnd)
[题目链接] http://acm.hdu.edu.cn/showproblem.php?pid=5716 [题目大意] 给出一个字符串,找出其中所有的符合特定模式的子串位置,符合特定模式是指,该子串 ...
- 2018 ACM-ICPC 中国大学生程序设计竞赛线上赛 H题 Rock Paper Scissors Lizard Spock.(FFT字符串匹配)
2018 ACM-ICPC 中国大学生程序设计竞赛线上赛:https://www.jisuanke.com/contest/1227 题目链接:https://nanti.jisuanke.com/t ...
- CCF CSP 201409-3 字符串匹配
CCF计算机职业资格认证考试题解系列文章为meelo原创,请务必以链接形式注明本文地址 CCF CSP 201409-3 字符串匹配 问题描述 给出一个字符串和多行文字,在这些文字中找到字符串出现的那 ...
- 字符串匹配&Rabin-Karp算法讲解
问题描述: Rabin-Karp的预处理时间是O(m),匹配时间O( ( n - m + 1 ) m )既然与朴素算法的匹配时间一样,而且还多了一些预处理时间,那为什么我们还要学习这个算法呢?虽然Ra ...
随机推荐
- APEX-数据导出/打印
前言: 由于公司使用了Oracle APEX构建应用,且在APEX新版本v20.2版本中增强了相关报表导出数据相关功能:正好现在做的事情也需要类似的功能,就先来学习一下Oracle的APEX相关功能及 ...
- NET 5 Cron表达式
cron表达式通过特定的规则指定时间,用于定时任务 1. 整体结构 cron表达式是一个字符串,分为6或7个域,每两个域之间用空格分隔,其语法格式为: "秒域 分域 时域 日域 月域 周域 ...
- HTTP Error 405.0 - Method Not Allowed 无法显示您正在查找的页面,因为使用了无效方法(HTTP 谓词)。
将submit改成button即可 因为触发了form表单
- (一)必须掌握的linux命令行快捷键
1.序 使用linux时,接触最多的莫过于命令行,参差不齐,形形色色,对于短的命令行使用脑残的上下左右,back,del就够用了:而对于带有很多参数的长的命令行,再使用上下左右,del,back显得那 ...
- BF,BM,KMP,就这?
为保证代码严谨性,文中所有代码均在 leetcode 刷题网站 AC ,大家可以放心食用. 皇上生辰之际,举国同庆,袁记菜馆作为天下第一饭店,所以被选为这次庆典的菜品供应方,这次庆典对于袁记菜馆是一项 ...
- [Machine Learning] 逻辑回归 (Logistic Regression) -分类问题-逻辑回归-正则化
在之前的问题讨论中,研究的都是连续值,即y的输出是一个连续的值.但是在分类问题中,要预测的值是离散的值,就是预测的结果是否属于某一个类.例如:判断一封电子邮件是否是垃圾邮件:判断一次金融交易是否是欺诈 ...
- JAVA数据结构(十一)—— 堆及堆排序
堆 堆基本介绍 堆排序是利用堆这种数据结构而设计的一种排序算法,堆排序是一种选择排序,最坏,最好,平均时间复杂度都是O(nlogn),不稳定的排序 堆是具有以下性质的完全二叉树:每个节点的值都大于或等 ...
- Java学习日报7.7
今天进一步学习了eclipse软件,遇到了几次程序运行不成功的问题,检查之后运行成功!明天继续学习程序逻辑控制!
- Java学习日报7.27
笔记
- Java学习日报7.18
/** * *//** * @author 86152 * */ package trangle;import java.util.Scanner; public class Trangle{ pub ...