【kmp-循环节】ARC060D 最良表現/Best Representation
题意简述
某谷的题意太绕了qwq
若一个字符串存在长度小于循环节,我们称它是一个好串。
现在要把分割成若干段,要求每一段都是一个好串。
所有满足条件的分割方案中,分割段数最小的那个方案被称为“全场最佳”。
求“全场最佳”的分割段数以及“全场最佳”的方案数。
题目解析
首先考虑怎么把循环节给搞出来。
这个可以用的失配数组来做。
如果一个长度为的串由个相同的长度为的子串拼接而成,并且是可以划分的最小循环串,那么一定有(前个子串拼在一起等于后个子串拼在一起)
并且根据失配数组()的定义可以知道,就是去掉最后一个的前一个位置。
所以
那么只需要判断就知道是否有循环节了
然后考虑这道题。
如果没有循环节,那么很显然答案就是不分割,
如果循环节为,那么就要全部分开 答案就是
如果有不为的循环节,那么可以从某个循环节的中间断开,分成段。
刚开始还以为可以直接算答案,所以我是这么写的:
int k=len/m;
if(k==2) k--;
int ans=len-1-(k-1);
printf("2\n%d\n",ans);
然后掉了才发现有这样的情况:
右边是循环串
发现总方案数不超过(所以是模了个寂寞)可以枚举一个断点来看是否满足左右两边都不循环
所以还需要倒着做一遍,用之前的方法判断左右两个子串是否为循环串。
由于是最小循环节,所以可以证明至少存在一种方案满足条件(分类讨论循环次数为和大于),所以是最小答案。
还是简单说一下:循环次数为时,从两个循环节之间劈开就可以,因为是最小循环节,所以两个子串一定不循环;循环次数大于时,可以把它分成部分,一段字符串个字符的形式,那么显然结尾的个字符不循环,接下来证明开头的那一个字符串一定不循环。假设开头的字符串循环,可以画一个图大概看一下,利用两个循环关系对应位置上的数相等,可以导出存在一个更小的循环节,而这与是最小的循环节相矛盾,所以开头的字符一定不循环。
方框是一个循环节,蓝色线段是另外一个循环节(一个是刚开始的,另一个是去掉最后一个字符的循环节)根据对应关系可以得出这些红色的线段都是相等的:
然后你发现两条蓝色的线段是相等的,而其中一条蓝色线段有两条红色线段,另一条蓝色线段只有一条红色线段,于是可以对应得到剩下那一部分是红色线段,然后在根据剩下的黑框框相等得到左边那个黑框框的下一部分也是红色线段...然后套娃下去。
并且红色线段的长度是黑框框的长度和蓝色线段长度的因数,所以红色线段一定是一个循环节,所以存在一个更短的循环节,与是最小的循环节相矛盾。
updated:
然后路过的大佬给了我一个数学上的思路:(果然数学好和数学不好的人的差距就出来了)
假设两个循环节的长度分别为,那么有
组合一下,,可以取一切整数(包括负数,相当于线性组合)
那么
而,则
所以存在一个更小的循环节。
最后我写的时候倒过来的时候写掉了一行,(注释里有),改了好久,还是要多谢大佬orz
►Code View
#include<cstdio>
#include<algorithm>
#include<queue>
#include<cstring>
#include<string>
#include<map>
#include<cmath>
using namespace std;
#define N 500005
#define MOD 1000000007
#define LL long long
int rd()
{
int x=0,f=1;char c=getchar();
while(c<'0'||c>'9'){if(c=='-')f=-1; c=getchar();}
while(c>='0'&&c<='9'){x=(x<<3)+(x<<1)+(c^48); c=getchar();}
return f*x;
}
char s[N];
int p[N],q[N];
int j,len;
bool check1(int i)
{
int t=i-p[i];
if(i==t||i%t) return 1;
return 0;
}
bool check2(int i)
{
int t=q[i]-i;
if((len-i+1)==t||(len-i+1)%t) return 1;
return 0;
}
int main()
{
scanf("%s",s+1);
p[1]=0;
j=0,len=strlen(s+1);
//表示1...j 和 len-j+1...len 可以匹配
for(int i=2;i<=len;i++)
{//求失配数组
while(j>0&&s[j+1]!=s[i]) j=p[j];
if(s[j+1]==s[i]) p[i]=j+1;
else p[i]=0;
j=p[i];
}
j=len+1,q[len]=len+1;
for(int i=len-1;i>=1;i--)
{
//printf("%d\n",i);
while(j<len+1&&s[j-1]!=s[i]) j=q[j];
if(s[j-1]==s[i]) q[i]=j-1;
else q[i]=len+1;//啊 这里打掉了
j=q[i];
}
int m=len-p[len];
if(m==1)
{
printf("%d\n1\n",len);
return 0;
}
if(len%m||len==m)
{
puts("1"); puts("1");
return 0;
}
int ans=0;
for(int i=1;i<=len-1;i++)
if(check1(i)&&check2(i+1)) ans++;
printf("2\n%d\n",ans);
return 0;
}
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· go语言实现终端里的倒计时
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 使用C#创建一个MCP客户端
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列1:轻松3步本地部署deepseek,普通电脑可用
· 按钮权限的设计及实现
2019-10-31 C++ 结构体重载运算符
2019-10-31 C++ 优先队列 priority_queue
2019-10-31 USACO3.3 A Game【区间dp】
2019-10-31 USACO3.3 Home on the Range【思维】
2019-10-31 【生成树趣题】CF723F st-Spanning Tree