(以下字符串下标从0开始,并定义$2s=s+s$)

考虑$f(S)$,即令$l=\max_{2i<|S|且S[0,i)=S[|S|-i,|S|)]}i$,则$f(S)=S+S[l,|S|-l)$

由于次数足够多,先做一次$S=f(S)$不影响答案,因此假设原串为$2S$(这个$S$不同于初始的$S$)

考虑$f(2S)$,类似的即令$l=\max_{i<|S|且S[0,i)=S[|S|-i,|S|)}i$(最长非自身的border),则$f(2S)=2(S+S[0,|S|-l))$

构造一个新的函数$g(S)=S+S[0,|S|-l)$,不难归纳得到$f_{k}(2S)=2g_{k}(S)$,且由于$|g_{10^{100}}(S)|\ge 10^{18}$,因此不需要考虑2倍,只需要求$g_{10^{100}}(S)$

令$T=S[0,|S|-l)$,那么$g_{2}(S)=g(S+T)$,记$S'=S+T$,即求$S'$时的$l$

根据$l$的定义可得$S[0,|S|-|T|)=S[|T|,|S|)$,代入到每一个字符即$\forall 0\le i<|S|-|T|,S_{i}=S_{i+|T|}$

若$|T|$为$|S|$的约数,由此不难得到$S=TT...T$,否则$S=TT...TP$(其中$P$为$T$的前缀),之后$l$所对应的$S$的border即为去掉第一个$T$后的串

对于第一种情况有$l=|S|$,可得$S'[0,l)=S[|S'|-l,|S'|)=S$,同时若有更大的$l'$,可得$TT...TP=PTT...T$(其中$P=T[0,l'-|S|)$),那么将两边同时去掉串首的$T$即得到$S$一个更长的border

(关于第二个为什么以$P$开头:首先由于其与$T$相同,必然是$T$的前缀,长度又等于$P$,根据前缀的唯一性即为$P$)

在此基础上继续归纳下去,即得$g_{k}(S)=g_{k-1}(S)+T=STT...T$(共$k$个$T$),简单计算一下即可

对于第二种情况有$l=|T|$,可得$S'[0,l)=S'[|S'|-l,|S'|)=T$,同时若存在更大的$l'$,则考虑两者相同的最后$|T|$个字符,前者是$T$的一个循环,后者即为$T$

假设以$T_{i}$开头,即得到$T_{j}=T_{(j+i)mod\ |T|}$,不难得到$T$必然是一个串循环若干次,这与KMP的性质矛盾,因此不存在更大的$l'$

接下来同样归纳,可得$g_{k+1}(S)=g_{k}(S)+g_{k-1}(S)$,考虑如何处理:

先将其差分,求长为$l$的前缀的各个字母数量,而由于$|g_{0}(S)|,|g_{1}(S)|\ge 1$,不难发现其在100项时长度一定超过$10^{18}$,找到最后一个$g_{k}(S)\le l$,之后统计$g_{k}(S)$的字母数量并递归$l-g_{k}(S)$即可

这样的时间复杂度最优可以达到$o(|S|+|\alpha|m)$(其中$m=100$)

 1 #include<bits/stdc++.h>
2 using namespace std;
3 #define N 200005
4 #define M 105
5 #define ll long long
6 int n,m,nex[N];
7 ll l,r,len[M],sum[M][31],ans[31];
8 char s[N];
9 void calc1(ll k,int p){
10 if (k<n)
11 for(int i=0;i<k;i++)ans[s[i]-'a']+=p;
12 else{
13 for(int i=0;i<n;i++)ans[s[i]-'a']+=p;
14 k-=n;
15 for(int i=0;i<nex[n];i++)ans[s[i]-'a']+=p*k/nex[n];
16 for(int i=0;i<k%nex[n];i++)ans[s[i]-'a']+=p;
17 }
18 }
19 void calc2(ll k,int m,int p){
20 if (k<=n){
21 for(int i=0;i<k;i++)ans[s[i]-'a']+=p;
22 return;
23 }
24 while (len[m]>k)m--;
25 for(int i=0;i<26;i++)ans[i]+=sum[m][i]*p;
26 calc2(k-len[m],m,p);
27 }
28 int main(){
29 scanf("%s%lld%lld",s,&l,&r);
30 n=strlen(s);
31 nex[0]=nex[1]=0;
32 for(int i=1,j=0;i<n;i++){
33 while ((j)&&(s[j]!=s[i]))j=nex[j];
34 if (s[j]==s[i])j++;
35 nex[i+1]=j;
36 }
37 int nn=n;
38 while (2*n>=nn)n=nex[n];
39 n=nn-n;
40 nex[n]=n-nex[n];
41 if (n%nex[n]==0){
42 calc1(r,1);
43 calc1(l-1,-1);
44 for(int i=0;i<26;i++)printf("%lld ",ans[i]);
45 return 0;
46 }
47 len[0]=n;
48 for(int i=0;i<n;i++)sum[0][s[i]-'a']++;
49 len[1]=n+nex[n];
50 for(int i=0;i<n;i++)sum[1][s[i]-'a']++;
51 for(int i=0;i<nex[n];i++)sum[1][s[i]-'a']++;
52 m=1;
53 while (len[m]<r){
54 m++;
55 len[m]=len[m-1]+len[m-2];
56 for(int i=0;i<26;i++)sum[m][i]=sum[m-1][i]+sum[m-2][i];
57 }
58 calc2(r,m,1);
59 calc2(l-1,m,-1);
60 for(int i=0;i<26;i++)printf("%lld ",ans[i]);
61 }

[atARC077F]SS的更多相关文章

  1. 极路由2(极贰)在OpenWrt下定制自己的ss服务

    默认刷入的OpenWrt带的ss, 只有ss-redir服务, 但是在实际使用中, 很多时候还是希望访问直接通过正常网关, 只有少部分访问需要通过ss, 所以希望能配置成为ss-local服务. 在保 ...

  2. 记一次ss故障

    本文主要参考: https://github.com/shadowsocks/shadowsocks shadowssocks 分为客户端和服务器端. 我们平时买的服务,使用是要用的是客户端. 如果你 ...

  3. ss命令和Recv-Q和Send-Q状态

    ss 用来显示处于活动状态的套接字信息.ss命令可以用来获取socket统计信息,它可以显示和netstat类似的内容.但ss的优势在于它能够显示更多更详细的有关TCP和连接状态的信息,而且比nets ...

  4. SimpleDateFormat df = new SimpleDateFormat("yyyy-MM-dd HH:mm:ss");//设置日期格式

    java日期格式大全 format SimpleDateFormat(转) SimpleDateFormat df = new SimpleDateFormat("yyyy-MM-dd HH ...

  5. SS - DIY一个前端模板引擎.(一)

    前端MVVM 模式有点很多,完全摆脱了意大利面条式的代码. 个人认为,所有MVVM 的框架基础就是一个高性能的JS模板引擎,它极大简化了 DOM 操作, 使页面渲染和业务逻辑彻底分离. 为了理解模板引 ...

  6. YYYY-mm-dd HH:MM:SS

    备忘:YYYY-mm-dd HH:MM:SS部分解释 d               月中的某一天.一位数的日期没有前导零.    dd             月中的某一天.一位数的日期有一个前导零 ...

  7. js 获取当前日期时间3种格式化方法 yyyy-mm-dd hh:MM:ss

    方法一: Date.prototype.format = function (format) { var args = { "M+": this.getMonth() + 1, & ...

  8. 在VPS上搭建SS访问火星

    前段时间发布了Visual Studio 2017 RC,由于现在VS没有离线的ISO了,只有一个在线安装文件.虽然可以通过这个在线安装文件生成完整的离线安装包(之前的ISO版本在安装过程中仍然需要联 ...

  9. 每天一个linux命令(57):ss命令

    ss是Socket Statistics​的缩写.顾名思义,ss命令可以用来获取socket统计信息,它可以显示和netstat类似的内容.但ss的优势在于它能够显示更多更详细的有关TCP和连接状态的 ...

随机推荐

  1. Space Time Pattern Mining Tools(时空模式挖掘工具)

    时空模式挖掘工具 # Process: 局部异常值分析 arcpy.LocalOutlierAnalysis_stpm("", "", 输出要素, " ...

  2. 洛谷2375 NOI2014动物园(KMP)

    题目链接: 题目. 简单一点来说,这个题就是求一个字符串的\(num\)数组的和,其中有\(num[i]\)表示1~i中有多少个不交叉的相等的前缀和后缀 的数目,要求一个\(O(n)\)的做法 QwQ ...

  3. Python 3.10 正式发布,新增模式匹配,同事用了直呼真香!

    关注微信公众号:K哥爬虫,QQ交流群:808574309,持续分享爬虫进阶.JS/安卓逆向等技术干货! 前几天,也就是 10 月 4 日,Python 发布了 3.10.0 版本,什么?3.9 之后居 ...

  4. C 库函数源码

    github URL git://sourceware.org/git/glibc.git 码云 URL https://gitee.com/jason.R.xie/glibc.git

  5. AES解密尾部出现乱码问题

    说明 在使用AES解密的时候我发现解密出来的字符串尾部一直都有乱码 解决方案 尾部字符串的ascii码就是删除位索引 具体代码: cryptor = AES.new('AES_KEY'.encode( ...

  6. 什么是js事件,冒泡机制,事件捕获,默认行为

    js事件: javascript使我们能够有能力创建动态页面,事件就是可以被js侦测到的行为,网页中每个元素都可以产生某些触发js函数的事件. 例如我们可以在用户点击某个按钮时产生一个click事件来 ...

  7. Java基础之原生JDBC操作数据库

    前言 日常开发中,我们都习惯了使用ORM框架来帮我们操作数据库,本文复习.记录Java如何使用原生JDBC操作数据库 代码编写 封装几个简单方法 find查询方法 findOne查询方法 update ...

  8. 21.10.18 test

    可可大神出题,四款有趣的游戏推荐,第四个好玩/se T1 loopers \(\color{green}{100}\) 考虑钦定 \(a_1,a_i\) 的位置,固定左边一坨,那么剩下的一坨的 \(\ ...

  9. VirtualBox问题解决合集 - [drm:vmw_host_log [vmwgfx]] *ERROR* Failed to send host log message

    转载:https://blog.csdn.net/mychangee/article/details/104954262 问题描述:[drm:vmw_host_log [vmwgfx]] ERROR  ...

  10. 数组中重复的数字 牛客网 剑指Offer

    数组中重复的数字 牛客网 剑指Offer 题目描述 在一个长度为n的数组里的所有数字都在0到n-1的范围内. 数组中某些数字是重复的,但不知道有几个数字是重复的.也不知道每个数字重复几次.请找出数组中 ...