首页 > 其他分享 >闲话 24.6.23

闲话 24.6.23

时间:2024-06-23 14:09:23浏览次数:21  
标签:right frac 23 bm 24.6 闲话 partial aligned left

闲话

推歌:Empurple by 春卷饭 feat.初音未来

虽然是商曲吧,但是确实仙品
这种破题方法确实很巧妙,从发色出发,分解成红+蓝再寻找其抽象的指代意
但不停留于此,而是接着探讨更深刻的问题(例如童年、家庭、控制
过程中仍然是大量的隐喻和意象。
我得到的结论就是,紫化不是无法变成红色的妥协,而应当是一种必然的结果——好孩子是难以幸存的。
对埃及文化不太了解,而且只是想推个歌,所以我就不说太多了
请大家听一听歌,展开自己的想象吧(

饭终于用别的歌姬了 很独特(
感觉供给pjsk的那一版 少了 super 春卷饭 time 不是很好听啊(

Please, Forgive me and "Purple"
摇摆不定的我染上紫色(Empurple)

多元拉格朗日反演

因为大家都进化了。而我还没有。所以来进化一下自己。

先来点记号。
由于这里是多元拉格朗日反演,我们不妨先确定占位元的数量 \(m\)。

对一个列 \(a_1, a_2, \dots, a_m\),定义其顺序组成的向量为 \(\bm a\)。这里的列的元素可为占位元、形式幂级数或数字。当 \(\bm n\) 由数字组成时,定义 \(\bm x^{\bm n} = (x_1^{n_1}, x_2^{n_2}, \dots, x_m^{n_m}) = x_1^{n_1} x_2^{n_2}\cdots x_m^{n_m}\)。它既可以作为 \(m\) 元形式幂级数的一组自变量,也可以作为其的某一项。

记 \(n\) 行 \(m\) 列,第 \(i\) 行、第 \(j\) 列值为 \(a_{i, j}\) 的矩阵为\(A = \left\{a_{i, j}\right\}_{n, m}\)。

\(\textbf{定理 1 } \text{(多元拉格朗日反演)}\)

对 \(m\) 个 \(m\) 元形式幂级数 \(f_1, f_2, \dots, f_m\),令其复合逆分别为 \(g_1, g_2, \dots, g_n\),满足对任意 \(k\) 有

\[f_k(\bm x) = x_k g_k(\bm f) \]

则对任意 \(m\) 元形式幂级数 \(h(\bm x)\),我们有

\[[\bm t^{\bm n}] h(\bm f(\bm t)) = [\bm x^{\bm n}] h(\bm x) \bm g(\bm x)^{\bm n} \det \left\{[i = j] - \frac{x_i}{g_j(\bm x)}\frac{\partial g_j(\bm x)}{\partial x_i} \right\}_{m, m} \]

证明:留做习题。\(\square\)

如果展开的话就是

\[\begin{aligned} & [t_1^{n_1} t_2^{n_2}\cdots t_m^{n_m}] h(f_1(t_1, t_2, \dots, t_m), \dots, f_m(t_1, t_2, \dots, t_m)) \\ = \ & [x_1^{n_1} x_2^{n_2}\cdots x_m^{n_m}] h(x_1, x_2, \dots, x_m) \left(\prod_{i = 1}^m g_i(x_1, x_2, \dots, x_m) ^{n_i} \right) \\ &\det \left\{[i = j] - \frac{x_i}{g_j(x_1, x_2, \dots, x_m)}\frac{\partial g_j(x_1, x_2, \dots, x_m)}{\partial x_i} \right\}_{m, m} \end{aligned} \]

这样又好看又不好看(

但某篇论文里说这个形式不好用,所以我们可以用内向树形式重写这结构。取 \(\mathcal T = (V, E)\) 为 \(m+1\) 个点(编号为 \(0\sim m\))的完全图的一棵以 \(0\) 为根的内向生成树,并令 \(\bm x\) 为占位元列 \(x_1, x_2, \dots, x_m\),\(\bm f\) 为函数列 \(f_0(\bm x), f_1(\bm x), \dots, f_m(\bm x)\)。那么记

\[\frac{\partial \bm f}{\partial\mathcal T} = \prod_{v \in V} \left(\prod_{(u, v) \in E} \frac{\partial}{\partial x_u} \right) f_{v}(\bm x) \]

注意这里每条有向边的起点都不可能为 \(0\),所以占位元只有 \(m\) 个。同时如果 \(\text{deg }v = 0\),里层的 \(\prod\) 贡献为 \(1\),而其他时刻需要对 \(f_v\) 连续求导,而非分别求导后相乘。

\(\textbf{定理 2 } \text{(多元拉格朗日反演)}\)

对 \(m\) 个 \(m\) 元形式幂级数 \(f_1, f_2, \dots, f_m\),令其复合逆分别为 \(g_1, g_2, \dots, g_n\),满足对任意 \(k\) 有

\[f_k(\bm x) = x_k g_k(\bm f) \]

则对任意 \(m\) 元形式幂级数 \(h(\bm x)\),我们有

\[[\bm t^{\bm n}] h(\bm f(\bm t)) = \frac{1}{\prod_{i =1}^m n_i} [\bm x^{\bm n - \bm 1}] \sum_{\mathcal T} \frac{\partial(h, g_1^{n_1} , g_2^{n_2}, \dots, g_m^{n_m})}{\partial\mathcal T} \]

然后是练习环节。

\(\text{例 1 }\) erase(多元)

当 \(m = 1\) 时,我们有 \(f(x) = x g(f(x))\),这也就是说 \(f^{\langle - 1\rangle} (x) = \dfrac{x}{g(x)}\)。此时多元拉格朗日反演 \(1\) 退化为

\[\begin{aligned} [t^n] h(f(t)) & \ = [x^{n}] h(x) g(x)^n\left(1 - \frac{x}{g(x)}g'(x) \right) \\ & \ = [x^{n}] h(x) g(x)^{n} \left(\frac{x}{g(x)} \right)' g(x) \\ & \ = [x^{n}] h(x) \left(f^{\langle - 1\rangle} (x)\right)' \left(\frac{x}{f^{\langle - 1\rangle} (x)}\right)^{n + 1} \end{aligned}\]

这自然是另类拉格朗日反演的形式。

此时树只有一种形态,因此多元拉格朗日反演 \(2\) 退化为

\[\begin{aligned} [t^n] h(f(t)) & \ = \frac{1}{n} [x^{n - 1}] \frac{\partial(h, g^n)}{\partial\mathcal T} \\ & \ = \frac{1}{n} [x^{n - 1}] \left( \frac{\partial h(x)}{\partial x} \right) \left( g(x)^n \right) \\ & \ = \frac{1}{n} [x^{n - 1}] h'(x) \left( \frac{x}{f(x)} \right)^n \end{aligned}\]

这自然是拉格朗日反演的形式。

\(\text{例 2 }\) 变量数++

现在我们来到了二元生成函数,希望推导出一个普适的方程。

多元拉反 \(1\):

\[\begin{aligned} & [x^{n_1}y^{n_2}] h(g_1(x, y), g_2(x, y)) \\ = \ &[x^{n_1}y^{n_2}] h(x, y) g_1(x, y)^{n_1} g_2(x, y)^{n_2} \\ & \begin{bmatrix} 1 - \dfrac{x_1}{g_1(x, y)}\dfrac{\partial g_1(x,y)}{\partial x_1} & - \dfrac{x_1}{g_2(x, y)}\dfrac{\partial g_2(x,y)}{\partial x_1} \\ - \dfrac{x_2}{g_1(x, y)}\dfrac{\partial g_1(x,y)}{\partial x_2} & 1 - \dfrac{x_2}{g_2(x, y)}\dfrac{\partial g_2(x,y)}{\partial x_2} \end{bmatrix} \end{aligned}\]

多元拉反 \(2\):

\[\begin{aligned} & [x^{n_1}y^{n_2}] h(g_1(x, y), g_2(x, y)) \\ = \ & \frac{1}{n_1 n_2} [x^{n_1-1} y^{n_2-1}]\left( \left(\frac{\partial }{\partial x} \frac{\partial }{\partial y}h(x, y) \right) g_1(x, y)^{n_1} g_2(x, y)^{n_2} + \right. \\ & \left. \frac{\partial h(x, y)}{\partial x} \frac{\partial g_1(x, y)^{n_1}}{\partial y} g_2(x, y)^{n_2}+ \frac{\partial h(x, y)}{\partial y} \frac{\partial g_2(x, y)^{n_2}}{\partial x} g_1(x, y)^{n_1}\right) \end{aligned}\]

感觉很难再化简了啊!

\(\text{例 3 }\) P10324 洞察

我们已经有了

\[\left\{ \begin{aligned} & f_1(x, y) = x \exp(f_2) \\& f_2(x, y) = y \exp(f_1) \end{aligned} \right.\]

这就有

\[\left\{ \begin{aligned} & g_1(x, y) = \exp y \\& g_2(x, y) = \exp x \end{aligned} \right.\]

施多元拉反 \(1\) 有

\[\begin{aligned} & [x^{n}y^{n}] \exp (f_1 + f_2) \\ = \ & [x^{n}y^{n}] e^{x + y} \left(e^y\right)^n \left(e^y\right)^n \begin{bmatrix} 1 - \dfrac{x}{e^y}\dfrac{\partial e^y}{\partial x} & - \dfrac{x}{e^x}\dfrac{\partial e^x}{\partial x} \\ - \dfrac{y}{e^y}\dfrac{\partial e^y}{\partial y} & 1 - \dfrac{y}{e^x}\dfrac{\partial e^x}{\partial y} \end{bmatrix} \\ = \ & [x^{n}y^{n}] e^{(n + 1)(x + y)} \begin{bmatrix} 1 & -x \\ - y & 1 \end{bmatrix} \\ = \ & [x^{n}y^{n}] e^{(n + 1)(x + y)} (1-xy) \end{aligned}\]

施多元拉反 \(2\) 有

\[\begin{aligned} & [x^{n}y^{n}] \exp (f_1 + f_2) \\ = \ & \frac{1}{n^2} [x^{n-1} y^{n-1}]\left( \left( \frac{\partial}{\partial x} \frac{\partial}{\partial y} e^{x+y} \right) (e^y)^{n} (e^x)^{n} + \right. \\ & \left. \frac{\partial e^{x+y}}{\partial x} \frac{\partial (e^y)^{n}}{\partial y} (e^x)^{n}+ \frac{\partial e^{x+y}}{\partial y} \frac{\partial (e^x)^{n}}{\partial x} (e^y)^{n}\right) \\ = \ & \frac{1}{n^2} [x^{n-1} y^{n-1}]\left( e^{x + y} (e^y)^{n} (e^x)^{n} + e^{x + y} ne^{ny} (e^x)^{n}+ e^{x + y} ne^{nx} (e^y)^{n}\right) \\ = \ & \frac{1}{n^2} [x^{n-1} y^{n-1}]\left( e^{(n + 1)(x + y)} + n e^{(n + 1)(x + y)}+ n e^{(n + 1)(x + y)}\right) \\ = \ & \frac{1}{n^2} [x^{n-1} y^{n-1}] e^{(n + 1)(x + y)} \left(1 + 2n\right) \end{aligned}\]

这提取系数得到的结果与 \(1\) 相符。

其他例题待补(

标签:right,frac,23,bm,24.6,闲话,partial,aligned,left
From: https://www.cnblogs.com/joke3579/p/-/chitchat240623

相关文章

  • from 2024.6.23 to unkown
    中考?不做评价。Question1【AGC023E】给定一个长度为\(n\)的正整数序列\(A\),对于一个长度为\(n\)的全排列\(P\),记\(I(P)\)表示\(P\)的逆序对数量,求:\[\underset{\forall1\leqi\leqn,P_i\leqA_i}{\sum}I(P)\]\(1\leqn\leq2\times10^5,1\leqA_i\leqn\)......
  • 2024/6/23 本周总结
    DemoFusion:DemocratisingHigh-ResolutionImageGenerationWithNo$$$2024/5/11任意尺度超分生成一个\(K\)倍大小的图像,需要边长扩大为\(\sqrt{K}\),就是从潜在空间(latentsapce)\(\mathbb{R}^{c\timesh\timesw}\)到目标空间\(\mathbb{R}^{c\times{H}\timesW}\),其中\(......
  • Java 学习知识点汇集(2024.6)
    VSCode,run程序时,提示,错误:找不到或无法加载主类Exam_32猜测原因,目录中有中文字符?解决办法:**在Java中,final类不能作为父类被继承**。讯飞星火:在Java的LSP(LiskovSubstitutionPrinciple,里氏替换原则)中,如果一个类被设计为不可变的(immutable)或者已经完成的(complete),它应该......
  • Lightroom Classic 2023 for Mac(摄影后期图像编辑工具) v12.4版
    lightingClassic是Adobe公司推出的一款图像处理软件,是数字摄影后期制作的重要工具之一。与其他图像处理软件相比,LightroomClassic具有以下特点:LightroomClassic2023forMac(摄影后期图像编辑工具)软件地址高效的图像管理:LightroomClassic提供了强大的图像管理功能,可以......
  • 胡说八道(24.6.21)——认识通信(杂谈)
        昨天说了在无线电通信几个应用,虽然这些天说了一大推,但是这只是理论上的东西,没有实物是没用的,离深入了解它们还是相差甚远。我觉得人家对不同种通信的认知就非常好。人家写下来,咱们就跟了解了解,正所谓集思广益。继续看雷达的简单应用。        雷达的优点......
  • 胡说八道(24.6.17)——STM32以及通信杂谈
        之前的文章中咱们谈到了STM32的时钟,今天我们来联系实际,来看看内部时钟下和外部时钟下的两种不同时钟的电平翻转。本次终于有硬件了,是最基础的STM32F103C8T6。    首先是,内部时钟的配置操作。             系统的内部时钟是72MHz,由上图......
  • 胡说八道(24.6.18)——通信杂谈(科普知识)
        既聊完中国和西方的通信历史之后,咱们继续来看看与有线通信相对应的无线通信,至于有线通信线的类型这里就不多说,像电话线,光纤,电缆(用于有线电视信号传输、早期的计算机网络)等,这些都是有线通信的范围。今天,来看看无线通信。        有线传输通信解决了异地......
  • Leetcode 力扣 125. 验证回文串 (抖音号:708231408)
    如果在将所有大写字符转换为小写字符、并移除所有非字母数字字符之后,短语正着读和反着读都一样。则可以认为该短语是一个 回文串 。字母和数字都属于字母数字字符。给你一个字符串 s,如果它是 回文串 ,返回 true ;否则,返回 false 。示例1:输入:s="Aman,aplan,......
  • Leetcode 力扣 128. 最长连续序列 (抖音号:708231408)
    给定一个未排序的整数数组 nums ,找出数字连续的最长序列(不要求序列元素在原数组中连续)的长度。请你设计并实现时间复杂度为 O(n) 的算法解决此问题。示例1:输入:nums=[100,4,200,1,3,2]输出:4解释:最长数字连续序列是[1,2,3,4]。它的长度为4。示例2:输入:nums=......
  • 238. 除自身以外数组的乘积
    题目给你一个整数数组nums,返回数组answer,其中answer[i]等于nums中除nums[i]之外其余各元素的乘积。题目数据保证数组nums之中任意元素的全部前缀元素和后缀的乘积都在32位整数范围内。请不要使用除法,且在O(n)时间复杂度内完成此题。示例1:输入:nums=......