专栏文章

题解:P8322 『JROI-4』少女幻葬

P8322题解参与者 1已保存评论 0

文章操作

快速查看文章及其快照的属性,并进行相关操作。

当前评论
0 条
当前快照
1 份
快照标识符
@miq7xmbj
此快照首次捕获于
2025/12/04 00:25
3 个月前
此快照最后确认于
2025/12/04 00:25
3 个月前
查看原文
感觉前两篇题解的式子太意识流了。
记值域为 mm
由于 gcd(ai,ai+1,ai+2)=k\gcd(a_i,a_{i+1},a_{i+2})=k,所以有 kaik|a_i。不妨令 aiaik,lilik,ririka_i\leftarrow\dfrac{a_i}{k},l_i\leftarrow\lceil\dfrac{l_i}{k}\rceil,r_i\leftarrow\lfloor\dfrac{r_i}{k}\rfloor。则原限制转化为 liairi,gcd(ai,ai+1)1,gcd(ai,ai+1,ai+2)=1l_i\le a_i\le r_i,\gcd(a_i,a_{i+1})\neq 1,\gcd(a_i,a_{i+1},a_{i+2})=1。接下来的 kk 和题目里的没有关系。
设一个状态转移会比较抽象。设 f(i,j,k)f(i,j,k) 表示考虑到第 ii 个数,ai=ka_i=kgcd(ai,ai1)=j\gcd(a_i,a_{i-1})=j 的方案数;设 g(i,j,k)g(i,j,k) 表示考虑到第 ii 个数,ai=ka_i=kgcd(ai,ai+1)=j\gcd(a_i,a_{i+1})=j 的方案数。
首先是 ggff 的转移为 f(i,j,k)=v=li1ri1g(i1,j,v)[gcd(v,k)=j]f(i,j,k)=\sum\limits_{v=l_{i-1}}^{r_{i-1}}g(i-1,j,v)[\gcd(v,k)=j]
除掉 jj 后反演得到 f(i,j,k)=v=li1jri1jg(i1,j,vj)dv,dkjμ(d)f(i,j,k)=\sum\limits_{v=\lceil\frac{l_{i-1}}{j}\rceil}^{\lfloor\frac{r_{i-1}}{j}\rfloor}g(i-1,j,vj)\sum\limits_{d|v,d|\frac{k}{j}}\mu(d)
即:
f(i,j,k)=dkjμ(d)v=li1djri1djg(i1,j,vdj)f(i,j,k)=\sum\limits_{d|\frac{k}{j}}\mu(d)\sum\limits_{v=\lceil\frac{l_{i-1}}{dj}\rceil}^{\lfloor\frac{r_{i-1}}{dj}\rfloor}g(i-1,j,vdj)
不妨 kkjk\leftarrow \dfrac k j,对第三维做狄利克雷后缀和,对位乘 μ\mu 后再做狄利克雷前缀和即可转移,时间复杂度 O(nmlogmloglogm)O(nm\log m\log\log m),前一个 logm\log m 是枚举因数的调和级数,后一个 loglogm\log\log m 是狄利克雷前后缀和。
然后再是 ffgg 的转移为:
g(i,j,k)=xkf(i,x,k)[gcd(x,j)=1]g(i,j,k)=\sum\limits_{x|k}f(i,x,k)[\gcd(x,j)=1]
枚举 kk,把 f(i,x,k)f(i,x,k) 贡献到 xx 的质因数集合的位置然后求高维后缀和,那么 g(i,j,k)g(i,j,k) 就是 jj 的质因数集合的补集的位置。这一部分复杂度是 kmd(k)ω(k)=kmω(k)ik1=ikmω(ik)ikm(ω(i)+ω(k))=2imkmiω(i)=2imO(miloglogvi)=O(mlogmloglogm)\sum\limits_{k\le m}d(k)\omega(k)=\sum\limits_{k\le m}\omega(k)\sum\limits_{i|k}1=\sum\limits_{ik\le m}\omega(ik)\le\sum\limits_{ik\le m}(\omega(i)+\omega(k))=2\sum\limits_{i\le m}\sum\limits_{k\le\frac m i}\omega(i)=2\sum\limits_{i\le m}O(\dfrac m i\log\log\dfrac v i)=O(m\log m\log\log m)。加上枚举第一维,总的复杂度是 O(nmlogmloglogm)O(nm\log m\log\log m)
然后就做完了,时间复杂度 O(nmlogmloglogm)O(nm\log m\log\log m)

评论

0 条评论,欢迎与作者交流。

正在加载评论...