Description

  

  如果某个串可以由两个一样的串前后连接得到,我们就称之为“偶串”。比如说“xyzxyz”和“aaaaaa”是偶串,而“ababab”和“xyzxy”则不是偶串。

  

​   对于一个非空串\(S\),我们定义\(f(S)\)是在\(S\)后面添加一些字符得到的最短偶串。比如\(f(\)'abaaba'\()=\)'abaababaab'。容易证明,对于一个非空串\(S\),\(f(S)\)是唯一的 。

  

​   现在给定一个由小写英文字母构成的偶串\(S\),你需要求出\(f^{10^{100}}(S)\),并统计计算结果的第\(l\)个字符到第\(r\)个字符中,每个字母出现了多少次 。

  

​   其中,\(f^{10^{100}}(S)\)是指\(f(f(f(...f(S)...)))\),式子中共有\(10^{100}\)个\(f\) 。

  

  

  

Solution

  

​   打表找规律。从整个串来寻找规律看不出什么,反而是分别考虑两半的变化比较有效。

  

​   我们称\(T\)为\(S\)的"Border",当且仅当\(S\)作为\(T\)的循环串的前缀出现,且\(T\)长度最小。

  

​   现有输入串\(SS\),令\(S\)的"Border"为\(T\),则会发现:

\[SS\rightarrow (ST)(ST)\rightarrow (STS)(STS)\rightarrow(STSST)(STSST)...
\]

​   由于操作次数过多,只考虑前半部分的变化,就能回答所有询问:

\[S\rightarrow ST\rightarrow STS\rightarrow STSST\rightarrow STSSTSTS
\]

​   我们发现,从第三个串开始,每一个串\(i\),都是由串\(i-1\)+串\(i-2\)得来。姑且称它们为\(\text{fib}\)串。

  

​   特别地:若\(|S|-|T|\)整除\(|S|\),则变化为

\[S\rightarrow ST\rightarrow STT\rightarrow STTT...
\]

​   特判掉这种情况。

  

​   记\(len_i\)为第\(i\)个\(\text{fib}\)串的长度,\(g_{i,j}\)为第\(i\)个\(\text{fib}\)串中\(j\)字符出现了多少次。由于\(\text{fib}\)在第85项左右就已经超出了\(1e18\),因此我们可以暴力计算出这两个数组。

  

​   对答案差分,假设要求\(1...r\)中,每个字符的出现次数。

  

​   有结论是,只要有\(len_{x_1}+len_{x_2}+...+len_{x_k}=r\),且\(len_{x_1}>len_{x_2}>...>len_{x_k}\),则\(1...r\)这个字符串就可以用\(x_1,x_2,...,x_k\)这\(k\)个\(\text{fib}\)串首尾相接得到。

  

​   所以对于询问我们从大到小枚举\(\text{fib}\)串累加每个字符的答案。最后可能会剩余一定长度无法消去,但剩余的长度一定不超过\(|S|\)(在\(\text{fib}\)串中甚至没有组成一个元素\(S\)),且是\(S\)的一个前缀,额外计算上这些部分即可。

  

​  

  

Code

  

#include <cstdio>
#include <cstring>
using namespace std;
typedef long long ll;
const int N=200005;
int n,up,nex[N],type;
char str[N];
int pre[N][26];
ll l,r,cut,len[100],ans[26];
ll f[100][26];
void readData(){
scanf("%s%lld%lld",str+1,&l,&r);
n=strlen(str+1);
n>>=1;
for(int i=1;i<=n;i++){
for(int j=0;j<26;j++) pre[i][j]=pre[i-1][j];
pre[i][str[i]-'a']++;
}
}
void getNex(){
nex[1]=0;
for(int i=2,j;i<=n;i++){
j=nex[i-1];
while(j&&str[j+1]!=str[i]) j=nex[j];
nex[i]=(str[j+1]==str[i])?j+1:0;
}
}
void preWork(){
len[1]=n;
for(int i=1;i<=n;i++){
f[1][str[i]-'a']++;
f[2][str[i]-'a']++;
}
len[2]=n+cut;
for(int i=0;i<26;i++)
f[2][i]+=pre[cut][i];
for(int i=3;;i++){
for(int j=0;j<26;j++) f[i][j]=f[i-2][j]+f[i-1][j];
len[i]=len[i-2]+len[i-1];
if(len[i]>1e18){
up=i;
break;
}
}
}
void calc1(ll m,ll a){
if(m<=0) return;
ll t=m/cut;
for(int i=0;i<26;i++) ans[i]+=a*t*pre[cut][i];
t=m%cut;
for(int i=0;i<26;i++) ans[i]+=a*pre[t][i];
}
void calc2(ll m,ll a){
for(int i=up;i>=1;i--)
if(len[i]<=m){
m-=len[i];
for(int j=0;j<26;j++) ans[j]+=a*f[i][j];
}
for(int j=0;j<26;j++) ans[j]+=a*pre[m][j];
}
int main(){
readData();
getNex();
cut=n-nex[n];
if(n%cut==0){
calc1(r,1);
calc1(l-1,-1);
}
else{
preWork();
calc2(r,1);
calc2(l-1,-1);
}
for(int i=0;i<26;i++) printf("%lld ",ans[i]);
return 0;
}

【ARC077F】SS的更多相关文章

  1. 【ARC077F】SS kmp+打表找规律

    Description ​ 如果某个串可以由两个一样的串前后连接得到,我们就称之为"偶串".比如说"xyzxyz"和"aaaaaa"是偶串, ...

  2. 【arc077f】AtCoder Regular Contest 077 F - SS

    题意 给你一个形如"SS"的串S,以及一个函数\(f(x)\),\(x\)是一个形如"SS"的字符串,\(f(x)\)也是一个形如"SS"的 ...

  3. 【arc077f】AtCoder Regular Contest 074 F - Lotus Leaves

    题意 给定一个n*m的池塘,每个格子上可能有叶子. 从一个叶子出发,可以跳到相同行或相同列的叶子. 问至少去掉多少叶子,使得起点不能到达终点. \(n,m<=100\) 解法 很显然的最小割模型 ...

  4. 【转载】运维小技巧:使用ss命令代替 netstat

    转自:https://www.91ri.org/12470.html ss是Socket Statistics的缩写. 顾名思义,ss命令可以用来获取socket统计信息,它可以显示和netstat类 ...

  5. 一步步实现ABAP后台导入EXCEL到数据库【1】

    在SAP的应用当中,导入.导出EXCEL文件的情况是一个常见的需求,有时候用户需要将大量数据定期导入到SAP的数据库中.这种情况下,使用导入程序在前台导入可能要花费不少的时间,如果能安排导入程序为后台 ...

  6. 【Django】--Model字段

    参考地址:http://www.cnblogs.com/wupeiqi/articles/6216618.html 所有字段 AutoField(Field) --int自增列,必须填入参数prima ...

  7. 通用js函数集锦<来源于网络> 【二】

    通用js函数集锦<来源于网络> [二] 1.数组方法集2.cookie方法集3.url方法集4.正则表达式方法集5.字符串方法集6.加密方法集7.日期方法集8.浏览器检测方法集9.json ...

  8. 通用js函数集锦<来源于网络/自己> 【一】

    通用js函数集锦<来源于网络/自己>[一] 1.返回一个全地址2.cookie3.验证用户浏览器是否是微信浏览器4.验证用户浏览器是否是微博内置浏览器5.query string6.验证用 ...

  9. 【五】将博客从jekyll迁移到了hexo

    本系列有五篇:分别是  [一]Ubuntu14.04+Jekyll+Github Pages搭建静态博客:主要是安装方面  [二]jekyll 的使用 :主要是jekyll的配置  [三]Markdo ...

随机推荐

  1. Unity3D — — UGUI之RectTransform

    Mask.GetComponent<RectTransform>().anchoredPosition(子物体) = hotKey_image.rectTransform.anchored ...

  2. HackRF 升级固件到新版本

    本文内容.开发板及配件仅限用于学校或科研院所开展科研实验! 淘宝店铺名称:开源SDR实验室 HackRF链接:https://item.taobao.com/item.htm?spm=a1z10.1- ...

  3. GlusterFS分布式存储集群-1. 部署

    参考文档: Quick Start Guide:http://gluster.readthedocs.io/en/latest/Quick-Start-Guide/Quickstart/ Instal ...

  4. SpringMVC 常用注释

    @Controller  Controller控制器是通过服务接口定义的提供访问应用程序的一种行为 @Repository  Dao层的标志 @RequestMapping  标注控制层函数的访问路径 ...

  5. DataTime日期格式化

    C# DateTime日期格式化 在C#中DateTime是一个包含日期.时间的类型,此类型通过ToString()转换为字符串时,可根据传入给Tostring()的参数转换为多种字符串格式. 目录 ...

  6. BATA冲刺准备

    目录 第一部分 调研,评测 福大助手的bug IOS端 Android端 福大助手结构体系的思维导图 为什么开发人员没有发现这个bug 假设团队开发这款app,应注意哪些方面(架构.部署运维.微服务等 ...

  7. java小学生四则运算带面板版 但我不知道为什么同类变量却进不了动作监听中去

    ---恢复内容开始--- package yun; import java.util.*; import java.awt.*; import java.awt.event.ActionEvent; ...

  8. APP案例分析之华为浏览器

    第一部分  对华为浏览器的调研,评测 1.对华为浏览器的第一次上手体验 我会使用华为浏览器呢,是因为我买的手机是华为nova,该浏览器也是手机里面预装的.刚开始用的时候也没太注意,感觉跟以前用的其他浏 ...

  9. vue-cli 安装步骤(转载)

    参考资料:Vue2.0 新手完全填坑攻略—从环境搭建到发布 1.Node.js安装 https://nodejs.org/en/download/ 2.安装vue-cli npm install -g ...

  10. My复利计算程序测试报告

    My复利计算程序测试报告 4.0 单元测试----------------------------- 要求: 对我们和复利计算程序,写单元测试. 有哪些场景? 期待的返回值 写测试程序. 运行测试. ...