【学习笔记】gcd matrix
考察:基本行列式变换。
- 用一个数乘行列式的任一列(行),等于用这个数乘此行列式。
- 交换 A 的第 i 行和第 j 行,行列式的 值变号 。
- 把行列式的某行(列)的 k k k 倍加到另一行(列)上,行列式的值不变。
例子 1:
D = ∣ 3 2 2 2 2 3 2 2 2 2 3 2 2 2 2 3 ∣ 分析:各行元素之和为一定数,故累加到 同一行再提公因数。 D=\left| \begin{array}{cccc} 3 & 2 & 2 & 2 \\ 2 & 3 & 2 & 2\\ 2 & 2 & 3 & 2\\ 2 & 2 & 2 & 3 \\ \end{array} \right|分析:各行元素之和为一定数,故累加到\\同一行再提公因数。 D=∣ ∣3222232222322223∣ ∣分析:各行元素之和为一定数,故累加到同一行再提公因数。
原式= ∣ 9 9 9 9 2 3 2 2 2 2 3 2 2 2 2 3 ∣ = 9 ∣ 1 1 1 1 2 3 2 2 2 2 3 2 2 2 2 3 ∣ = 9 ∣ 1 1 1 1 0 1 0 0 0 0 1 0 0 0 0 1 ∣ = 9 ∗ 1 = 9 \left| \begin{array}{cccc} 9 & 9 & 9 & 9 \\ 2 & 3 & 2 & 2\\ 2 & 2 & 3 & 2\\ 2 & 2 & 2 & 3 \\ \end{array} \right|=9\left| \begin{array}{cccc} 1 & 1 & 1 & 1 \\ 2 & 3 & 2 & 2\\ 2 & 2 & 3 & 2\\ 2 & 2 & 2 & 3 \\ \end{array} \right|=9\left| \begin{array}{cccc} 1 & 1 & 1 & 1 \\ 0 & 1 & 0 & 0\\ 0 & 0 & 1 & 0\\ 0 & 0 & 0 & 1 \\ \end{array} \right|=9*1=9 ∣ ∣9222932292329223∣ ∣=9∣ ∣1222132212321223∣ ∣=9∣ ∣1000110010101001∣ ∣=9∗1=9
注意这里对系数的处理(利用对称性创造出均一系数),在很多行列式的问题(数论代数组合数学综合)中都有体现。
行列式求值的基本思想是消成上三角或下三角模型。此时行列式的值为 对角线数的乘积 。
例2:
已知 A i , j = gcd ( i , j ) A_{i,j} =\gcd(i,j) Ai,j=gcd(i,j) ,求这个行列式的值。
sol:
我们首先断言:
A i , j ′ = { ϕ ( i ) , i ∣ j 0 , i ∤ j A^\prime_{i,j}= \begin{cases} \phi(i),&i\mid j \\ 0,&i\nmid j \end{cases} Ai,j′={ϕ(i),0,i∣ji∤j
考虑这样的消元:在第 i i i 次消元中,把 i i i 的倍数行减去第 i i i 行。
我们考虑前 i − 1 i-1 i−1 行对 A i , j A_{i,j} Ai,j 的贡献。第 k k k 行对第 i i i 行有贡献当且仅当 k ∣ i k\mid i k∣i ,同时 A k , j ≠ 0 A_{k,j}\neq 0 Ak,j=0 当且仅当 k ∣ j k\mid j k∣j 。综上可得 k ∣ gcd ( i , j ) k|\gcd(i,j) k∣gcd(i,j) 。
那么 A i , j ′ = gcd ( i , j ) − ∑ d ∣ gcd ( i , j ) , d < i ϕ ( d ) A'_{i,j}=\gcd(i,j)-\sum_{d|\gcd(i,j),d<i}\phi(d) Ai,j′=gcd(i,j)−∑d∣gcd(i,j),d<iϕ(d)
我们知道 ∑ d ∣ n ϕ ( d ) = n \sum_{d|n}\phi(d)=n ∑d∣nϕ(d)=n (顺便捡一下数论的知识,虽然忘得差不多了 qwq) ,所以当 gcd ( i , j ) ≠ i \gcd(i,j)\neq i gcd(i,j)=i 时为 0,否则为 ϕ ( i ) \phi(i) ϕ(i) 。
同时我们证明了消元后一定是一个上三角矩阵。答案为 ∏ i = 1 n ϕ ( i ) \prod_{i=1}^n\phi(i) ∏i=1nϕ(i) 。
总结:本题中的断言是思维难点也是关键点,可以大大简化思维。同时借助数论中的 欧拉函数 完成了证明。
__EOF__

本文链接:https://www.cnblogs.com/cqbzly/p/17530217.html
关于博主:评论和私信会在第一时间回复。或者直接私信我。
版权声明:本博客所有文章除特别声明外,均采用 BY-NC-SA 许可协议。转载请注明出处!
声援博主:如果您觉得文章对您有帮助,可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力!
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
· DeepSeek 开源周回顾「GitHub 热点速览」