LOJ3123 CTS2019 重复 KMP自动机、DP、多项式求逆
CTS的计数题更完辣(撒花
Orz zx2003,下面的内容在上面的博客基础上进行一定的补充。
考虑计算无限循环之后不存在子串比\(s\)字典序小的串的个数。先对串\(s\)建立KMP自动机,那么对于点\(i\)连出的所有边,只有不是回到起点的字符最大的那条边以及字符比它大的走向起点的边可以走,因为如果走其他边,那么当前匹配的串就会存在一个后缀比\(s\)字典序小。
经过这样的删边后,对于每一个点只有两种边,一种是去往另一个点的边,一种是去往起点的边。如果不考虑去往起点的边,与起点连通的连通块刚好构成一个rho的形状(也就是\(\rho\),即一条从起点开始的链,链底是一个环)。我们考虑在这个rho上计数。
不妨设当前构造的满足条件的串是串\(t\)。有一个引理:设\(a_0 = 0\),\(a_i =\)从KMP自动机上的\(a_{i-1}\)号点开始匹配串\(t\)走到的位置,可以得到一个长度无限的序列\(a\),在序列\(a\)上必定存在一个位置\(x\)满足\(\forall y > x , a_y = a_x\)。
证明:\(a_x\)实质上就是把\(x\)个\(t\)串拼在一起放在KMP自动机上匹配最后到达的节点。因为如果\(t\)合法,那么\(t\)只会走每个点字符最大的一条边和比这个字符更大的回到起点的边,所以将\(x\)个\(t\)串放在现在的KMP自动机上匹配等价于将\(x\)个\(t\)串放在没有删边之前的KMP自动机上匹配,而在后者意义下匹配完到达的节点表示\(x\)个\(t\)串的最长后缀满足它是\(s\)的前缀。当\(x|t|>|s|\),因为在后面加上一个\(t\)相当于在前面加上一个\(t\),而此时加上一个\(t\)不可能对最长后缀的长度产生影响,所以再加上一个\(t\)进行匹配,仍会匹配到之前的节点。
显然这个匹配到的节点是唯一的,所以我们可以枚举这个节点,那么我们需要求的就是长度为\(m\)的路径条数满足起点和终点都是这一个点。
如果从一个节点开始匹配串\(t\),没有经过起点回到自己,那么一定是在rho上转,方案数唯一。
接下来考虑回到起点的情况。考虑一个dp:设\(f_{i,j}\)表示从起点开始走\(i\)步走到节点\(j\)的方案数,这个可以直接做。接下来我们枚举从这个节点沿着\(\rho\)走了多少步回到根,用回到根的方案数乘上从起点经过若干步到达当前节点的方案数贡献答案。
总复杂度\(O(n^2)\)。
接下来使用多项式运算优化上面后半部分计算的复杂度。
上面枚举的路径相当于三个过程:从当前节点开始走到起点->走若干次起点到起点的路径->从起点走到当前节点停下来,而第一个过程和第三个过程合起来等价于先枚举一个从起点开始走回起点的路径,然后再在这条路径上选择一个点作为匹配终点。那么先枚举这条路径,对于一条长度为\(x\)的路径,就有\(x\)种选择匹配终点的方式。
不妨设\(F = \sum\limits_{i=1}^\infty f_ix^i\),其中\(f_i\)表示从起点开始走\(i\)步走回到起点的方案数,这个可以在KMP自动机建立之后计算。那么枚举第一个过程和第三个过程的路径的并,方案数的多项式就是\(x F'\)。接下来再枚举经过了多少次起点,需要乘上多项式\(\sum\limits_{i=0}^\infty F^i = \frac{1}{1 - F}\),那么最终的答案就是\([x^m] (xF' \times \frac{1}{1 - F} + 1)\),最后加上的\(1\)是一个余项。
LOJ3123 CTS2019 重复 KMP自动机、DP、多项式求逆的更多相关文章
- CF848E Days of Floral Colours——DP+多项式求逆/分治NTT
官方题解:http://codeforces.com/blog/entry/54233 就是由简入繁 1.序列处理,只考虑一个半圆 2.环形处理(其实这个就是多了旋转同构) 然后基于分割线邻居的跨越与 ...
- BZOJ 3456: 城市规划(dp+多项式求逆)
传送门 解题思路 这道题就是求带标号的无向连通图个数,首先考虑\(O(n^2)\)的做法,设\(f_i\)表示有\(i\)个节点的无向连通图个数,那么考虑容斥,先把所有的无向图求出,即为\(2^{C( ...
- 【bzoj3456】城市规划 dp+多项式求逆
Description 刚刚解决完电力网络的问题, 阿狸又被领导的任务给难住了. 刚才说过, 阿狸的国家有n个城市, 现在国家需要在某些城市对之间建立一些贸易路线, 使得整个国家的任意两个城市都直接或 ...
- BZOJ 3456: 城市规划 与 多项式求逆算法介绍(多项式求逆, dp)
题面 求有 \(n\) 个点的无向有标号连通图个数 . \((1 \le n \le 1.3 * 10^5)\) 题解 首先考虑 dp ... 直接算可行的方案数 , 容易算重复 . 我们用总方案数减 ...
- 【bzoj3456】城市规划(多项式求逆+dp)
Description 求\(~n~\)个点组成的有标号无向连通图的个数.\(~1 \leq n \leq 13 \times 10 ^ 4~\). Solution 这道题的弱化版是poj1737, ...
- [BZOJ 3625] [Codeforces 438E] 小朋友的二叉树 (DP+生成函数+多项式开根+多项式求逆)
[BZOJ 3625] [Codeforces 438E] 小朋友的二叉树 (DP+生成函数+多项式开根+多项式求逆) 题面 一棵二叉树的所有点的点权都是给定的集合中的一个数. 让你求出1到m中所有权 ...
- Luogu4233 射命丸文的笔记 DP、多项式求逆
传送门 注意到总共有\(\frac{n!}{n}\)条本质不同的哈密顿回路,每一条哈密顿回路恰好会出现在\(2^{\binom{n}{2} - n}\)个图中,所以我们实际上要算的是强连通有向竞赛图的 ...
- [BZOJ3456]城市规划:DP+NTT+多项式求逆
写在前面的话 昨天听吕老板讲课,数数题感觉十分的神仙. 于是,ErkkiErkko这个小蒟蒻也要去学数数题了. 分析 Miskcoo orz 带标号无向连通图计数. \(f(x)\)表示\(x\)个点 ...
- bzoj 3456 城市规划——分治FFT / 多项式求逆 / 多项式求ln
题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3456 分治FFT: 设 dp[ i ] 表示 i 个点时连通的方案数. 考虑算补集:连通的方 ...
随机推荐
- 洛谷P1854 花店橱窗布置
题目 DP,直接递推比记忆化搜索简单. 定义状态\(dp[i][j]\)为前i行最后一个选择第i行第j个数所得到最大值. 易得状态转移方程 \(dp[i][j]=max(dp[i-1][k]+a[i] ...
- 意图Intent
意图点击官方链接 前言 对意图Intent,学习安卓需掌握.以官方链接:http://www.android-doc.com/reference/android/content/Intent.html ...
- javascript 中的方法注入
js 中的方法注入 java中很多框架支持 apo 的注入, js中也可以类似的进行实现 主要是通过扩展js中方法的老祖 Function 对象来进行实现. Function.prototype.af ...
- MATLAB 线性规划实例应用
线性规划 线性规划函数 功能:求解线性规划问题 语法 x = linprog(f,A,b):求解问题 min fx,约束条件为 Ax <= b x = linprog(f,A,b,Aeq,beq ...
- linux(deepin) 下隐藏firefox标题栏
1. 右上角菜单 -> 定制 -> 左下角 "标题栏" 取消打钩 2. 如果上面无法解决,在firefox的启动前插入一个环境变量,具体修改 /usr/share/ap ...
- TOMCAT到底能 承受多少并发,并发量计算你方法
TOMCAT 可以稳定支持的最大并发用户数 https://www.jianshu.com/p/d306826aef7a tomcat并发数优化maxThreads.acceptCount(最 ...
- 系统树图 | Dendrogram construction | Phylogenetic Analysis
Molecular Architecture of the Mouse Nervous System 表示亲缘关系的树状图解 先看文章里是怎么做的: Dendrogram construction A ...
- 安装Rancher容器化管理平台
卸载已有的避免冲突 yum remove docker docker-common docker-selinux docker-engine 安装需要的软件包,yum-util 提供yum-co ...
- typescript - 1.环境搭建
安装 1.安装node.js 安装typescript npm install -g typescript 3.编译 tsc helloworld.ts Typescript开发工具Vscode自动编 ...
- Python3基础 tuple(list) 改变list元素的内容时,元组的id值不变
Python : 3.7.3 OS : Ubuntu 18.04.2 LTS IDE : pycharm-community-2019.1.3 ...