HDU 6987 - Cycle Binary(找性质+杜教筛)
首先 mol 一发现场 AC 的 csy 神仙
为什么这题现场这么多人过啊啊啊啊啊啊
继续搬运官方题解(
首先对于题目中的 \(k,P\),我们有若存在字符串 \(k,P,P’\) 满足 \(S=kP+P’\),且 \(P’\) 为 \(P\) 前缀,那么称 \(P\) 为 \(S\) 的一个周期,显然 \(k\) 最大等价于 \(|P|\) 最小,也就是说我们要找到对于每个 \(|P|\),最短周期长度为 \(|P|\) 的字符串个数,那么显然对于长度为 \(x\) 的周期,如果一个字符串最短周期长度为 \(x\),那么其对应的 \(k=\lfloor\dfrac{|S|}{n}\rfloor\),也就是说,假设 \(f(i)\) 表示最短周期长度为 \(i\) 的字符串个数,那么
\]
那么怎么求 \(f(i)\) 呢?我们不妨做以下猜想:
Observation \(-1\). 一个串 \(S\) 最小周期长度为 \(x\),当且仅当长度为 \(x\) 的前缀是 \(S\) 的周期且不存在某个 \(y\mid x\),满足长度为 \(y\) 的前缀是 \(S\) 的周期
为什么是 \(-1\) 呢?因为显然这个结论是错的,反例随便举,譬如 \(S='\text{abaa}',x=3\),u1s1 我现场开场 30min 就在肝这道题并猜了这个结论,然鹅意识到这个结论是错误的,并且注意到现场 30min 没人通过这道题,感觉不对劲就放弃了(
不过注意到违背这个结论的都是 \(x\) 比较大的情况,对于 \(x\) 比较小的情况这个结论是不是就对了呢?
事实果真如此。
Observation. 一个串 \(S\) 最小周期长度为 \(x(x\le\lfloor\dfrac{n}{2}\rfloor)\),当且仅当长度为 \(x\) 的前缀是 \(S\) 的周期且不存在某个 \(y\mid x\),满足长度为 \(y\) 的前缀是 \(S\) 的周期
证明:首先关于周期我们有一个周期引理(Periodicity Lemma),大概说的是如果一个字符串存在两个周期长度分别为 \(x,y\),且 \(x+y\le n+\gcd(x,y)\),那么长度为 \(\gcd(x,y)\) 的前缀也是 \(S\) 的周期。
这样就可以反证法了,假设存在一个 \(y<x\) 满足 \(y\nmid x\) 且长度为 \(y\) 的前缀也是 \(S\) 的周期,那么由于 \(x\le\lfloor\dfrac{n}{2}\rfloor\),必然有 \(x+y\le n<n+\gcd(x,y)\),套用 Periodicity Lemma 可得 \(\gcd(x,y)\) 也是 \(S\) 的周期,而 \(\gcd(x,y)\mid x\),与长度为 \(x\) 的前缀为 \(S\) 的最小周期矛盾,因此假设不成立。
注意到对于 \(x>\lfloor\dfrac{n}{2}\rfloor\),必然有 \(\lfloor\dfrac{n}{x}\rfloor=1\),因此
res&=\sum\limits_{i=1}^n\lfloor\dfrac{n}{i}\rfloor f(i)\\
&=\sum\limits_{i=1}^{n/2}\lfloor\dfrac{n}{i}\rfloor f(i)+2^n-\sum\limits_{i=1}^{n/2}f(i)\\
&=2^n-\sum\limits_{i=1}^{n/2}(\lfloor\dfrac{n}{i}\rfloor-1)f(i)
\end{aligned}
\]
而对于 \(x\le\lfloor\dfrac{n}{2}\rfloor\),根据上面的性质有
\]
即枚举 \(S\) 的前 \(x\) 位是什么,这样即可确定整个 \(S\),同时需扣掉最小周期 \(<x\) 且是 \(x\) 的约数的串的个数。
这看起来和数论有点关系了,移个项:
\]
也即
\]
即
\]
于是上面那个求答案的式子就可以整除分块+杜教筛了。时间复杂度 \(Tn^{2/3}\)
const int MAXV=5e6;
const int MOD=998244353;
int qpow(int x,int e){
int ret=1;
for(;e;e>>=1,x=1ll*x*x%MOD) if(e&1) ret=1ll*ret*x%MOD;
return ret;
}
int f[MAXV+5];
void init(){
for(int i=1,pw=2;i<=MAXV;i++,(pw<<=1)%=MOD) f[i]=pw;
for(int i=1;i<=MAXV;i++) for(int j=i<<1;j<=MAXV;j+=i) f[j]=(f[j]-f[i]+MOD)%MOD;
for(int i=1;i<=MAXV;i++) f[i]=(f[i-1]+f[i])%MOD;
}
unordered_map<int,int> sm;
int calc(int x){
if(x<=MAXV) return f[x];
if(sm.count(x)) return sm[x];int res=(qpow(2,x+1)-2+MOD)%MOD;
for(int l=2,r;l<=x;l=r+1){
r=x/(x/l);res=(res-1ll*(r-l+1)*calc(x/l)%MOD+MOD)%MOD;
} return res;
}
int main(){
init();int qu;scanf("%d",&qu);
while(qu--){
int n;scanf("%d",&n);int res=qpow(2,n);
for(int l=1,r;l<=(n>>1);l=r+1){
r=min(n/(n/l),n>>1);
res=(res+1ll*(calc(r)-calc(l-1)+MOD)*(n/l-1))%MOD;
} printf("%d\n",res);
}
return 0;
}
HDU 6987 - Cycle Binary(找性质+杜教筛)的更多相关文章
- 洛谷 P6860 - 象棋与马(找性质+杜教筛)
题面传送门 首先我们来探究一下什么样的 \((a,b)\) 满足 \(p(a,b)=1\).不难发现只要点 \((1,0)\) 能够到达,那么网格上所有点都能到达,因为由于 \((1,0)\) 能够到 ...
- HDU 5608 function(莫比乌斯反演 + 杜教筛)题解
题意: 已知\(N^2-3N+2=\sum_{d|N}f(d)\),求\(\sum_{i=1}^nf(i) \mod 1e9+7\),\(n\leq1e9\) 思路: 杜教筛基础题? 很显然这里已经设 ...
- CCPC 2019 网络赛 HDU huntian oy (杜教筛)
1005 huntian oy (HDU 6706) 题意: 令,有T次询问,求 f(n, a, b). 其中 T = 10^4,1 <= n,a,b <= 1e9,保证每次 a,b互质. ...
- [HDU 5608]Function(莫比乌斯反演 + 杜教筛)
题目描述 有N2−3N+2=∑d∣Nf(d)N^2-3N+2=\sum_{d|N} f(d)N2−3N+2=∑d∣Nf(d) 求∑i=1Nf(i)\sum_{i=1}^{N} f(i)∑i=1Nf ...
- HDU 6706 huntian oy(杜教筛 + 一些定理)题解
题意: 已知\(f(n,a,b)=\sum_{i=1}^n\sum_{j=1}^igcd(i^a-j^a,i^b-j^b)[gcd(i,j)=1]\mod 1e9+7\),\(n\leq1e9\),且 ...
- BZOJ4916 神犇和蒟蒻 【欧拉函数 + 杜教筛】
题目 很久很久以前,有一只神犇叫yzy; 很久很久之后,有一只蒟蒻叫lty; 输入格式 请你读入一个整数N;1<=N<=1E9,A.B模1E9+7; 输出格式 请你输出一个整数A=\sum ...
- Mobius 反演与杜教筛
积性函数 积性函数 指对于所有互质的整数 aaa 和 bbb 有性质 f(ab)=f(a)f(b)f(ab)=f(a)f(b)f(ab)=f(a)f(b) 的数论函数. 特别地,若所有的整数 aaa ...
- LG4213 【模板】杜教筛(Sum)和 BZOJ4916 神犇和蒟蒻
P4213 [模板]杜教筛(Sum) 题目描述 给定一个正整数$N(N\le2^{31}-1)$ 求 $$ans_1=\sum_{i=1}^n\varphi(i)$$ $$ans_2=\sum_{i= ...
- 我也不知道什么是"莫比乌斯反演"和"杜教筛"
我也不知道什么是"莫比乌斯反演"和"杜教筛" Part0 最近一直在搞这些东西 做了将近超过20道题目吧 也算是有感而发 写点东西记录一下自己的感受 如果您真的 ...
随机推荐
- csh
在*unix系统中,常用的shell有sh,bash,csh/tcsh, ksh. sh来自于systemV的Unix,是传统的Unix的shell,直到现在很多的系统管理员仍然喜欢使用sh. ba ...
- 【机器学习基础】逻辑回归——LogisticRegression
LR算法作为一种比较经典的分类算法,在实际应用和面试中经常受到青睐,虽然在理论方面不是特别复杂,但LR所牵涉的知识点还是比较多的,同时与概率生成模型.神经网络都有着一定的联系,本节就针对这一算法及其所 ...
- DOS命令和快捷键
- javascript-jquery-更改jquery对象
在许多情况下,jquery代码所做的事情变成了:生成jquery对象A,操作对jquery象A:更改为jquery对象B,操作jquery对象B:更改为jqueryC,操作jquery对象C..... ...
- 【UE4 C++】学习笔记汇总
UE4 概念知识 基础概念--文件结构.类型.反射.编译.接口.垃圾回收.序列化[导图] GamePlay架构[导图] 类的继承层级关系[导图] 反射机制 垃圾回收机制/算法 序列化 Actor 的生 ...
- 【UE4 C++】获取运行时间、设置时间流速、暂停游戏
基于UGameplayStatics 获取运行时间 /** Returns the frame delta time in seconds, adjusted by time dilation. */ ...
- Python中的括号()、[]、{}
长时间不用容易混淆,仅记! 在Python语言中最常见的括号有三种,分别是:小括号().中括号[].花括号{} . Python中的小括号(): 代表tuple元祖数据类型,元祖是一种不可变序列.大多 ...
- 按之字形顺序打印二叉树 牛客网 剑指Offer
按之字形顺序打印二叉树 牛客网 剑指Offer 题目描述 请实现一个函数按照之字形打印二叉树,即第一行按照从左到右的顺序打印,第二层按照从右至左的顺序打印,第三行按照从左到右的顺序打印,其他行以此类推 ...
- Luogu 520题纪念
一入OI深似海......
- TDengine在浙商银行微服务监控中的实践
作者:楼永红 王轩宇|浙商银行 浙商银行股份有限公司(简称"浙商银行")是 12 家全国性股份制商业银行之一,总部设在浙江杭州,全国第13家"A+H"上市 ...