数学
文章平均质量分 92
山登绝顶我为峰 3(^v^)3
只有平凡理想的域(✪ω✪)
展开
-
快速乘法技巧:Karatsuba, Toom, Good, Schonhage, Strassen, Nussbaumer
[Ber01] 详细地调研了上百篇有关各类乘法算法的文献,以环同构的视角解释这些算法。原创 2023-11-15 17:40:52 · 341 阅读 · 0 评论 -
分园域的一些结论
代数整数、分园数域、分园环、理想分解原创 2023-11-02 23:56:35 · 240 阅读 · 0 评论 -
加权自动机:在 Semirings 上建模
文章目录:幺半群(Monoids)半环(Semirings)形式幂级数广义矩阵加权自动机经典自动机原创 2023-10-10 23:43:36 · 106 阅读 · 0 评论 -
分式理想 & 对偶群 & 对偶空间
分式理想、对偶群、对偶空间原创 2023-09-28 16:33:18 · 732 阅读 · 0 评论 -
Branching Program(5-PBP)
将一个**布尔电路转化为 5-PBP** 的方法如下:1. 在 Output-Gate 的位置开始递归,将识别 $A$ 的深度 $d$ 的布尔电路,转化为两个深度 $d-1$ 的分别识别 $A_1,A_2$ 的布尔电路以及一个 AND-Gate2. 接着对深度 $d-1$ 的两个布尔电路分别构造出 5-PBP,将它们按照 AND-Gate 的串联程序组合起来,如果两个 5-PBP 的输出 $\sigma_1=\sigma_2=\sigma$ 则使用引理换一个 $\sigma_2=\tau \neq \原创 2023-04-07 16:25:32 · 666 阅读 · 4 评论 -
有限域上的所有不可约多项式
之前写过,已知分圆多项式都是上的。根据剩余类判别法,存在一些素域,使得分圆多项式在其上也是不可约的。那么,给定一个有限域,哪些分圆多项式是不可约的?进一步的,所有的不可约多项式是哪些?原创 2022-10-12 21:59:00 · 2042 阅读 · 2 评论 -
剩余类环上可逆矩阵的计数
令m∈N+m \in \mathbb N^+m∈N+是任意正整数,利用数论基本定理,可以写作m=p1e1p2e2⋯psesm = p_1^{e_1}p_2^{e_2} \cdots p_s^{e_s}m=p1e1p2e2⋯pses其中pi,i=1,⋯ ,sp_i,i=1,\cdots,spi,i=1,⋯,s都是素数,ei,i=1,⋯ ,se_i,i=1,\cdots,sei,i=1,⋯,s是正整数。下面我们考虑一般线性群GLn(Zm)GL_n(\mathbb Z_m)GLn(Zm原创 2022-09-14 10:29:10 · 874 阅读 · 0 评论 -
函子(Functor)
给定两个范畴C,D\mathcal C,\mathcal DC,D,函子(functor)是一个函数T:C→DT:\mathcal C \to \mathcal DT:C→D,满足如果对象A∈obj(C)A \in obj(\mathcal C)A∈obj(C),那么T(A)∈obj(D)T(A) \in obj(\mathcal D)T(A)∈obj(D),函子将对象映射到对象。如果态射f:A→A′f: A \to A'f:A→A′在范畴C\mathcal CC里,那么T(f):T(A)→T(A′)T(f原创 2022-06-19 21:42:56 · 393 阅读 · 0 评论 -
范畴(Category)
在”朴素集合论“中,“集合”的定义十分模糊,它使用“a collection of things”这一直觉来定义集合。如果我们定义“所有自包含集合”的集合AAA,那么集合AAA是否包含自己?如果AAA包含自己,那么AAA成为AAA的一个成员(member);如果AAA不包含自己,那么AAA就不是AAA的成员;两者都是自恰的,这是不可判定问题。虽然它没有对朴素集合论造成毁灭性打击,但我们接下来考虑下“不自包含集合”。罗素集合(Russell set):一个集合SSS,它不包含自己,即S∉SS \notin S原创 2022-06-17 23:12:32 · 680 阅读 · 0 评论 -
几何定理的自动化证明(Automatic Geometric Theorem Proving)
对于几何定理(geometric theorems)的证明,可以通过引入笛卡尔坐标系,然后将前提(hypotheses)和结论(conclusions)描述为关于命题中点的坐标的多项式方程。令A,B,C,D,E,FA,B,C,D,E,FA,B,C,D,E,F都是实数空间中的点,这些基本的几何性质,都可以用关于点坐标的多项式来描述。我们说一个几何定理是可接受的(admissible),如果它的前提和结论都可以被翻译为多项式方程。给定任意一个几何定理,总有一些点的坐标是任意的(arbitrary),也总有另一些原创 2022-06-15 15:48:00 · 1262 阅读 · 0 评论 -
代数几何:簇的态射
令V⊆km,W⊆knV \subseteq k^m,W \subseteq k^nV⊆km,W⊆kn,一个函数ϕ:V→W\phi:V \to Wϕ:V→W叫做多项式映射(polynomial mapping、regular mapping),如果存在多项式f1,⋯ ,fn∈k[x1,⋯ ,xm]f_1,\cdots,f_n \in k[x_1,\cdots,x_m]f1,⋯,fn∈k[x1,⋯,xm],使得对于任意的a∈V\pmb a \in Vaaa∈V都有ϕ(a1,⋯ ,am)=(f1(a1,原创 2022-06-15 15:45:37 · 474 阅读 · 0 评论 -
代数几何:簇的分解、理想的分解
利用理想的升链条件(ACC),以及理想和簇的对应关系,可以容易得到如下条件。The Descending Chain Condition(DCC):任意的仿射簇的递减序列V1⊇V2⊇V3⊇⋯V_1 \supseteq V_2 \supseteq V_3 \supseteq \cdotsV1⊇V2⊇V3⊇⋯最终都会稳定(stabilize),即存在一个正整数NNN,使得VN=VN+1=⋯V_N = V_{N+1} = \cdotsVN=VN+1=⋯令V⊆knV \subseteq k^nV⊆原创 2022-06-10 15:48:21 · 474 阅读 · 0 评论 -
代数几何:Zariski Closures、不可约簇-素理想
# Zariski Closures任给一个集合$S \in k^n$,无论它是否是仿射簇,$$\pmb I(S) = \{ f \in k[x_1,\cdots,x_n]: f(a)=0,\forall a \in S \}$$它总是一个环$k[x_1,\cdots,x_n]$上的根理想,于是根据 ideal–variety correspondence,$\pmb V(\pmb I(S))$是簇。实际上,仿射簇$\pmb V(\pmb I(S))$是**包含集合$S$的最小簇**(sma原创 2022-06-09 22:06:59 · 1002 阅读 · 0 评论 -
代数几何:仿射簇-理想的对应关系,理想的和、积、交
令$V$是任意的仿射簇,如果对于某正整数$m$,有$f^m \in \pmb I(V)$,那么$f \in \pmb I(V)$。因为任意点$a \in V$,有$(f(a))^m = 0$,如果$f(a) \neq 0$,那么矛盾。**根理想**:理想$I$是**根的**(radical),它满足如果$\exist m\ge 1,f^m \in I$则导致$f \in I$。一个簇$V$,$I(V)$是根理想。**理想的根**:令$I \subseteq k[x_1,\cdots,x_n]$是理想原创 2022-06-09 15:50:46 · 1273 阅读 · 0 评论 -
代数几何:消元(Elimination)、扩展(Extension)
参考文献:Ideals, Varieties, and Algorithms (4th ed.) [Cox, Little & O’Shea 2015-06-14]前置文章:仿射簇 和 Groebner基域kkk,给定理想I=⊆k[x1,⋯ ,xn]I= \subseteq k[x_1,\cdots,x_n]I=⊆k[x1,⋯,xn],它的lll次消元理想( l-th elimination ideal)定义为Il=I∩k[xl......原创 2022-06-08 20:40:30 · 615 阅读 · 0 评论 -
如何判断多项式是否不可约
令环RRR上的nnn次多项式为:a(x)=∑i=0naixi∈R[x]a(x) = \sum_{i=0}^n a_i x^i \in R[x]a(x)=i=0∑naixi∈R[x]复数域上复数域是的代数封闭的。因此,复数域上的多项式a(x)a(x)a(x)都可以写成:a(x)=∏i=1n(x−αi)a(x) = \prod_{i=1}^{n}(x-\alpha_i)a(x)=i=1∏n(x−αi)因此,a(x)a(x)a(x)在复域上不可约 ⟺ n=1\iff n=1⟺n=1原创 2022-05-27 23:02:36 · 7055 阅读 · 0 评论 -
度量空间(metric space)
度量空间定义度量空间(metric space)是二元组(M,d)(M,d)(M,d),其中MMM是非空集合,度量(metric)是d:M×M→Rd:M \times M \to Rd:M×M→R是实值函数,它有如下性质:正定性(positive definiteness):∀x,y∈M, d(x,y)≥0\forall x,y \in M,\,\, d(x,y) \ge 0∀x,y∈M,d(x,y)≥0,并且d(x,y)=0 ⟺ x=yd(x,y)=0 \iff x=yd(x,y)=0⟺x=原创 2022-05-24 21:28:30 · 2299 阅读 · 0 评论 -
特殊卷积的多项式环 以及 Walsh变换
特殊卷积Convolution:hk=∑i+j=kmod naibjh_k = \sum_{i+j=k \mod n} a_i b_jhk=i+j=kmodn∑aibjXOR-Convolution:hk=∑i⊕j=kaibjh_k = \sum_{i \oplus j=k} a_i b_jhk=i⊕j=k∑aibjOR-Convolution:hk=∑i∨j=kaibjh_k = \sum_{i \vee j=k} a_i b_jhk=i∨j=k∑原创 2022-05-12 19:42:21 · 359 阅读 · 0 评论 -
代数结构:群、环、域、模、线性空间、格
群定义设GGG是非空集合,代数结构(G,∘)(G,\circ)(G,∘)叫做群,如果GGG上的二元运算∘\circ∘满足结合律:(a∘b)∘c=a∘(b∘c)(a \circ b) \circ c = a \circ (b \circ c)(a∘b)∘c=a∘(b∘c)单位元:e∘g=g∘e=ge \circ g = g \circ e = ge∘g=g∘e=g逆元:g∘g−1=g−1∘g=eg \circ g^{-1} = g^{-1} \circ g = eg∘g−1=g−1∘g=e环原创 2022-05-03 23:33:04 · 2651 阅读 · 0 评论 -
代数结构:模(Module)
# 模(Modules)### 定义环$R$,$R-$模是一个**加群**(交换加法群)$(M,+)$,伴随有**环作用**(action,也叫作数乘)$$R \times M \rightarrow M\\(r,u) \mapsto ru\\$$这个环作用满足:1. 右分配律:$r(u+v)=ru+rv$2. 左分配律:$(r+s)v=rv+sv$3. 对$R$的结合律:$r(sv)=(rs)v$4. 稳定性:$1_Rv=v$其中$r,s \in R$,$u,v \in M原创 2022-05-03 22:21:00 · 2288 阅读 · 0 评论 -
奇异值分解(SVD)与 主成分分析(PCA)
线性映射线性映射同构于矩阵乘。线性空间的一组基α=(α1,⋯ ,αn)\alpha=(\alpha_1, \cdots, \alpha_n)α=(α1,⋯,αn),一个点的坐标为x=(x1,⋯ ,xn)x=(x_1,\cdots,x_n)x=(x1,⋯,xn),那么点可以记做两者内积α⋅x=∑i=1nxiαi\alpha \cdot x = \sum_{i=1}^{n} x_i\alpha_iα⋅x=∑i=1nxiαi。注意,基不一定是向量,可以是任何线性无关的对象(比如,三角函数)。线性原创 2022-04-26 22:22:51 · 1024 阅读 · 0 评论 -
计算复杂性理论
参考:OI Wiki - 计算理论基础问题字母表(alphabet):有限非空集合Σ\SigmaΣ,其中的元素被称为符号(symbol)语言(language):字母表中元素组成的串的子集,L∈Σ∗L \in \Sigma^*L∈Σ∗;任何语言都可以转化为比特串,一般令Σ={0,1}\Sigma=\{0,1\}Σ={0,1}判定问题:判定一个串xxx是否属于一个语言LLL,只能回答Yes/No。f:Σ∗→{0,1}f:\Sigma^* \rightarrow \{0,1\}f:Σ∗→{0,1}功能原创 2022-04-17 13:49:48 · 637 阅读 · 0 评论 -
代数几何:Affine Variety 和 Groebner Basis
本文参考书籍:Ideals, Varieties, and Algorithms (4th ed.) [Cox, Little & O’Shea 2015-06-14]给定域kkk,令f1,⋯ ,fsf_1,\cdots,f_sf1,⋯,fs是多元多项式环k[x1,⋯ ,xn]k[x_1,\cdots,x_n]k[x1,⋯,xn]上的元素,定义仿射簇:V(f1,⋯ ,fs):={(a1,⋯ ,an)∈kn∣fi(a1,⋯ ,an)=0,∀1≤i≤s}V(f_1,\cdots,f_s):=\{(a原创 2022-04-06 16:29:38 · 1266 阅读 · 0 评论 -
快速Euclidean算法
扩展欧几里得算法给定a,ba,ba,b,计算(x,y,d)(x,y,d)(x,y,d),使满足:ax+by=d, d=gcd(a,b)ax+by=d,\, d=gcd(a,b)ax+by=d,d=gcd(a,b)易知:a⋅1+b⋅0=aa⋅0+b⋅1=b\begin{aligned}a \cdot 1 + b \cdot 0 = a\\a \cdot 0 + b \cdot 1 = b\\\end{aligned}a⋅1+b⋅0=aa⋅0+b⋅1=b以它为初值,令:a⋅x1+b⋅y1=原创 2021-12-24 21:44:20 · 438 阅读 · 4 评论 -
深入理解NTT
深入理解NTT中国剩余定理 (CRT):如果f,g∈R[x]f,g \in R[x]f,g∈R[x]且gcd(f,g)=1gcd(f,g)=1gcd(f,g)=1,那么存在环同构ϕ:R[x]/(f(x)g(x))≅R[x]/(f(x))×R[x]/(g(x))ϕ(h)=(hmod f, hmod g)\begin{aligned}\phi:R[x]/(f(x)g(x)) &\cong R[x]/(f(x)) \times R[x]/(g(x))\\\phi(h) &= (h原创 2021-11-21 11:07:40 · 1307 阅读 · 4 评论 -
数论变换(NTT)及其变体
CRT令RRR是基环 (base ring),那么R[x]R[x]R[x]是一元多项式环。如果f,g∈R[x]f,g \in R[x]f,g∈R[x]且gcd(f,g)=1gcd(f,g)=1gcd(f,g)=1,那么存在环同构:ϕ:R[x]/(f(x)g(x))≅R[x]/(f(x))×R[x]/(g(x))ϕ(h)=(hmod f, hmod g)\begin{aligned}\phi:R[x]/(f(x)g(x)) &\cong R[x]/(f(x)) \times R[x]/(原创 2021-11-19 21:56:17 · 1229 阅读 · 0 评论 -
快速傅里叶变换(FFT):蝶形算法(CT蝴蝶、GS蝴蝶)
DFT离散傅里叶变换和相应的逆变换为:Xk=∑l=0n−1xlwnkl, k=0,1,⋯ ,n−1xl=1n∑k=0n−1Xkwn−kl, l=0,1,⋯ ,n−1\begin{aligned}X_k = \sum_{l=0}^{n-1} x_l w_n^{kl},\,\,k=0,1,\cdots,n-1 \\x_l = \dfrac{1}{n} \sum_{k=0}^{n-1} X_k w_n^{-kl},\,\,l=0,1,\cdots,n-1 \\\end{aligned}Xk原创 2021-11-18 23:06:34 · 3103 阅读 · 1 评论 -
线性算法及应用(4)
正交内积映射:(u,v)∈V×V↦<u,v>∈R(u,v) \in V \times V \mapsto <u,v> \in R(u,v)∈V×V↦<u,v>∈R,满足对称性:<u,v>=<v,u><u,v>=<v,u><u,v>=<v,u>非负性:<u,u>≥0<u,u> \ge 0<u,u>≥0正定性:<u,u>=0 ⟺ u=0&l原创 2021-11-17 14:12:35 · 152 阅读 · 0 评论 -
线性算法及应用(3)
n阶差分方程n阶齐次线性差分方程yk+n+a1yk+n−1+⋯+anyk=0y_{k+n}+a_1y_{k+n-1}+\cdots+a_ny_{k}=0yk+n+a1yk+n−1+⋯+anyk=0,可以写成一阶差分方程xk+1=Axk, A∈Rn×nx_{k+1}=Ax_{k},\,A \in R^{n \times n}xk+1=Axk,A∈Rn×n对于xk+1=Axkx_{k+1}=Ax_{k}xk+1=Axk,如果AAA拥有n个线性无关的特征向量v1,⋯ ,vnv_1,\原创 2021-11-16 16:03:31 · 562 阅读 · 0 评论 -
线性算法及应用(2)
基础对于A=[a1 ⋯ an]∈Rm×nA=[a_1\,\,\cdots\,\,a_n] \in R^{m \times n}A=[a1⋯an]∈Rm×n,列空间(Column Space):Col(A):=Span{a1,⋯ ,an}={Ax∣x∈Rn}⊆RmCol(A) := Span\{a_1,\cdots,a_n\} = \{Ax|x \in R^n\} \sube R^mCol(A):=Span{a1,⋯,an}={Ax∣x∈Rn}⊆Rm零空间(Null Space)原创 2021-11-15 14:56:39 · 617 阅读 · 0 评论 -
线性算法及应用(1)
有助于快速计算的公式当det A≠0det\,A \not = 0detA=0时,2阶方阵AAA的逆矩阵的通项表达式为:A=[abcd]A−1=1ad−bc[d−b−ca]\begin{aligned} A = \begin{bmatrix} a & b\\ c & d\\ \end{bmatrix}\\ A^{-1} = \dfrac{1}{ad-bc} \begin{bmatrix} d & -原创 2021-11-13 21:17:28 · 707 阅读 · 0 评论 -
(高维)离散傅里叶变换DFT 和 快速傅里叶变换FFT
DFT 令p1,p2,...,pbp_1,p_2,...,p_bp1,p2,...,pb是b个整数,令 θj=e−2πi/pj\theta_j = e^{-2 \pi i/p_j}θj=e−2πi/pj 是pjp_jpj阶单位根,其中 i2+1=0i^2+1=0i2+1=0 定义G:=ZP1×Zp2×⋯×ZpbG:= Z_{P_1} \times Z_{p_2} \times \dots \times Z_{p_b}G:=ZP1×Zp2×⋯×Zpb,容易验证G是拥有p1p原创 2021-10-23 16:20:50 · 806 阅读 · 1 评论 -
Galois理论
共轭域将域K上所有代数元都添加到K中得到域KacK^{ac}Kac,称为K的代数闭包(algebraic closure)。域KacK^{ac}Kac是代数封闭的 (LLL的任何扩域中元素α\alphaα在L上代数,则α∈L\alpha \in Lα∈L;即LLL不存在更大代数扩域)。代数封闭域一定是无限域。有理数域QQQ的有限扩域 (自然是代数扩域) KKK叫做代数数域 (简称数域)。代数扩张L/KL/KL/K,有:Lac=KavL^{ac}=K^{av}Lac=Kav;数域KKK的代数闭包Ka原创 2021-10-30 22:16:18 · 980 阅读 · 0 评论