当前位置: 首页 > news >正文

成为 aoliao 粉丝的第一天

这个博客园渲染做了些什么玩意。

$$
f(x)=f(x_0)+f'(x_0)(x-x_0)+\frac{f''(x_0)}{2!}(x-x_0)^2+\dots\\
f'(x)=\sum a_nnx^{n-1}\\
\int f(x)=\sum\frac{a_n}{n+1}x^{n+1}+C
$$$$
\exp'(x)=\exp(x),\exp(0)=1\\
\text{取麦克劳林展开:}\exp(x)=1+x+\frac{x^2}{2!}+\frac{x^3}{3!}+\dots\\
\ln'(x)=\frac1x,\ln(1)=0\\
\text{取麦克劳林展开:}\ln(1+x)=x-\frac{x^2}2+\frac{x^3}3-\frac{x^4}4+\dots
$$对于可导函数 $f(x)$,求零点 $x_0$。任取一点 $x_1$,泰勒展开之:
$$
f(x)=f(x_1)+f'(x_1)(x-x_1)+R\\
0=f(x_1)+f'(x_1)(x_0-x_1)+R\\
x_0=x_1-\frac{f(x_1)}{f'(x_1)}-R
$$
考虑做若干次后,$x_1\equiv x_0\pmod {x^A}$,根据 $R=\frac{f''(\xi)}{2}(x_0-x_1)^2$,易得精度翻倍。求 $f^{-1}$,即要求关于 $G$ 的多项式 $G^{-1}-f$ 的零点,$G$ 为多项式。
$$
G_{new	}=G+\frac{G^{-1}-f}{G^{-2}}=2G-G^2f
$$
求 $f^{\frac1a}$,$a$ 为正整数,要求关于 $G$ 的多项式 $G^a-f$ 的零点。
$$
G_{new}=G-\frac{G^a-f}{aG^{a-1}}
$$
注意上面都是关于 $G$ 的多项式,因此求导的时候 $f$ 是常数,即微分的应该是 $dG$。对于 $a=2$,$T(n)=T(\frac n2)+O(n\log n)=O(n\log n)$。所以尽管看上去很荒谬,多项式乘法、需要做 $O(\log n)$ 次多项式乘法的多项式求逆、需要做 $O(\log n)$ 次多项式求逆的多项式开方都是 $O(n\log n)$ 的。初始值需要考虑二次剩余。求 $\ln f$,保证常数项为 $1$,直接用麦克劳林展开计算比较低效,考虑 $\ln f=[\int(\ln f)']+C$,这里微分的是 $dx$,根据 $\frac{d}{dx}f[g(x)]=f'[g(x)]\times g'(x)$,得到 $\ln f=(\int\frac{f'}f)+C$。常数项不一定对,代入 $x=0$ 得到 $C=0$。求 $\exp f$,保证常数项为 $0$,使用牛顿迭代,求 $\ln G-f$ 的零点。
$$
G_{new}=G-\frac{\ln G-f}{\frac1G}=G(1-\ln G+f)
$$
这样看来,若要求一个多项式函数 $h(f)$,如果 $h$ 的复合逆比较容易求且可导,那么就能牛顿迭代。如果 $h'$ 比较好求,可以先求导再积分。还可以构造复合 $z^{-1}[z(h)]$,这里 $z^{-1}$ 指 $z$ 的复合逆,但多项式复合逆比较困难,一般使用的也就是 $\exp$ 和 $\ln$ 这一对。求 $f^a$,$a$ 为正整数,有 $f^a=\exp(a\ln f)$,当常数项不为 $1$ 时比较难办。找到第一个非零项 $f_kx^k$,让 $f$ 除掉 $f_kx^k$ 后即可。时间复杂度 $O(n\log n)$。

诶,其实今天还有一个发现。所有容斥题殊途同归都是 \([n=0]=\sum\limits_{i=0}^n(-1)^i\binom ni\) 这个式子,原因在于这个式子本质上是先枚举所有方案,再枚举该方案下所有不合法的子集;我们把它改成先枚举所有不合法的子集,再在此基础上枚举方案,就是一般的容斥做法。本质上仍然是恰好转钦定,不过之前我只会对着二项式反演的式子硬做,今天终于领悟到了容斥是二项式反演特殊情形的道理,此事在去年五月时被 dcytrl 跟我反复提及,顿感自身之不足,当时不应轻视这句话的分量。不过这个东西好像是 classic,只有我现在才领悟罢了。

这样说来二项式反演也是在交换两个 \(\sum\)。本质上在利用 \([n=x]=\sum\limits_{i=x}^n(-1)^{x-i}\binom ix\),去年八月份的时候有人问我二项式反演的问题,当时我只会生硬地求出 \(g(x)\) 去推 \(f(x)\),如今确有一些更好的视角。感觉交换两个 \(\sum\) 这个理解真是太牛了,之前怎么没人跟我说过。

但好像说所有容斥题殊途同归都是这个也不太对,有一些题并不是这样简单地对限制容斥。你说对吧?zr hlt 场 t2。

http://www.jsqmd.com/news/343732/

相关文章:

  • translategemma-4b-it新手教程:不写代码也能用浏览器完成图文翻译任务
  • 深度学习环境配置:Ubuntu系统优化与常用命令大全
  • 世毫九实验室三大核心理论综合研究报告
  • Qwen3-ASR-0.6B语音识别:多语言支持实测分享
  • 如何将电池材料研磨至纳米级,具体方法、注意事项?
  • 灵感画廊开发者案例:基于SDXL 1.0构建跨平台(Web/iOS/Android)艺术APP
  • SiameseUIE中文-base部署避坑指南:常见连接失败与空结果解决方案
  • 使用Cursor开发EasyAnimateV5-7b-zh-InP模型插件
  • java二叉树的简单创建和遍历
  • 24G显存也能玩!BEYOND REALITY Z-Image高效部署方案
  • 2026宿州职业装定制企业评测:谁在引领娱乐行业新风尚? - 2026年企业推荐榜
  • Phi-4-mini-reasoning实测:128K长文本推理能力惊艳展示
  • 2026年音乐留学中介选择:如何识别并选择真正的4E教学体系机构? - 2026年企业推荐榜
  • 本杰明·格雷厄姆的资产负债表分析方法
  • 看不懂 ComfyUI?把它当“生产线”来设计工作流(新手也能上手的版本)
  • 如何识别企业的边缘AI视频分析芯片优势
  • 大数据领域中 RabbitMQ 的安全配置要点
  • 世盟股份深交所上市:市值48亿 预计年营收9亿同比降10%
  • 德业股份冲刺港股:9个月营收88亿 利润23亿 张和君控制60%表决权
  • 基于阿里云v2018数据集的机器负载概率预测:外生特征工程与Python实现
  • eVTOL主机厂沃飞长空完成近10亿融资:光合创投联合投资
  • 摆脱局域网!GoLand+cpolar 解锁 Go 开发远程协作新玩法
  • 多语言 RAG、算法思维、异常值检测与其他问题解决亮点
  • 多层感知器,解释:带有迷你二维数据集的视觉指南
  • 【资讯】CMSIS-DSP库发布V1.7.0版本,ARM发消息Mbed平台将于2026年7月终止服务
  • 没有检索模型构建 RAGs 是一个严重的错误
  • 选品更稳的新打法:用“用户决策阻力”挑品——阻力越大,越有机会做出溢价
  • 计算机毕业设计springboot“聚力”考研互助系统 基于 SpringBoot 的 “研途帮” 考研资源共享平台 基于 SpringBoot 的 “研聚力” 考研信息互助系统
  • 循环网络RNN--评论内容情感分析
  • STM32F103C8T6最小系统板控制Nano-Banana生成模型展示