algebra 0.2.1
Sign up to get free protection for your applications and to get access to all the features.
- checksums.yaml +7 -0
- data/.gitignore +53 -0
- data/.travis.yml +7 -0
- data/CODE_OF_CONDUCT.md +49 -0
- data/Gemfile +7 -0
- data/Gemfile.lock +28 -0
- data/LICENSE.txt +21 -0
- data/README-ja.txt +60 -0
- data/README.md +61 -0
- data/Rakefile +43 -0
- data/algebra.gemspec +33 -0
- data/bin/console +14 -0
- data/bin/setup +8 -0
- data/doc-ja/README-ja.html +69 -0
- data/doc-ja/README-ja.rd +60 -0
- data/doc-ja/README.html +0 -0
- data/doc-ja/README.rd +0 -0
- data/doc-ja/algebra-ja.html +66 -0
- data/doc-ja/algebra-ja.rd +60 -0
- data/doc-ja/algebraic-equation-ja.html +36 -0
- data/doc-ja/algebraic-equation-ja.rd +26 -0
- data/doc-ja/algebraic-extension-field-ja.html +122 -0
- data/doc-ja/algebraic-extension-field-ja.rd +118 -0
- data/doc-ja/algebraic-parser-ja.html +87 -0
- data/doc-ja/algebraic-parser-ja.rd +95 -0
- data/doc-ja/changes.html +161 -0
- data/doc-ja/changes.rd +119 -0
- data/doc-ja/elementary-divisor-ja.html +56 -0
- data/doc-ja/elementary-divisor-ja.rd +45 -0
- data/doc-ja/euclidian-ring-ja.html +78 -0
- data/doc-ja/euclidian-ring-ja.rd +69 -0
- data/doc-ja/finite-group-ja.html +285 -0
- data/doc-ja/finite-group-ja.rd +276 -0
- data/doc-ja/finite-map-ja.html +224 -0
- data/doc-ja/finite-map-ja.rd +227 -0
- data/doc-ja/finite-set-ja.html +393 -0
- data/doc-ja/finite-set-ja.rd +408 -0
- data/doc-ja/index-ja.html +113 -0
- data/doc-ja/index-ja.rd +112 -0
- data/doc-ja/jordan-form-ja.html +106 -0
- data/doc-ja/jordan-form-ja.rd +99 -0
- data/doc-ja/localized-ring-ja.html +121 -0
- data/doc-ja/localized-ring-ja.rd +131 -0
- data/doc-ja/m-polynomial-ja.html +542 -0
- data/doc-ja/m-polynomial-ja.rd +596 -0
- data/doc-ja/matrix-algebra-ja.html +686 -0
- data/doc-ja/matrix-algebra-ja.rd +727 -0
- data/doc-ja/matrix-algebra-triplet-ja.html +124 -0
- data/doc-ja/matrix-algebra-triplet-ja.rd +129 -0
- data/doc-ja/permutation-group-ja.html +159 -0
- data/doc-ja/permutation-group-ja.rd +151 -0
- data/doc-ja/polynomial-converter-ja.html +64 -0
- data/doc-ja/polynomial-converter-ja.rd +54 -0
- data/doc-ja/polynomial-ja.html +446 -0
- data/doc-ja/polynomial-ja.rd +473 -0
- data/doc-ja/residue-class-ring-ja.html +112 -0
- data/doc-ja/residue-class-ring-ja.rd +113 -0
- data/doc-ja/sample-algebraic-equation01.rb.v.rd +11 -0
- data/doc-ja/sample-algebraic-equation02.rb.v.rd +11 -0
- data/doc-ja/sample-algebraic-root01.rb.v.rd +20 -0
- data/doc-ja/sample-algebraicfield01.rb.v.rd +26 -0
- data/doc-ja/sample-algebraicfield02.rb.v.rd +15 -0
- data/doc-ja/sample-cayleyhamilton01.rb.v.rd +14 -0
- data/doc-ja/sample-diagonalization01.rb.v.rd +48 -0
- data/doc-ja/sample-divmod01.rb.v.rd +17 -0
- data/doc-ja/sample-elementary-divisor01.rb.v.rd +48 -0
- data/doc-ja/sample-factorize01.rb.v.rd +11 -0
- data/doc-ja/sample-factorize02.rb.v.rd +13 -0
- data/doc-ja/sample-factorize03.rb.v.rd +14 -0
- data/doc-ja/sample-factorize04.rb.v.rd +14 -0
- data/doc-ja/sample-factorize05.rb.v.rd +53 -0
- data/doc-ja/sample-galois-group01.rb.v.rd +27 -0
- data/doc-ja/sample-gaussian-elimination01.rb.v.rd +19 -0
- data/doc-ja/sample-geometry01.rb.v.rd +33 -0
- data/doc-ja/sample-geometry02.rb.v.rd +34 -0
- data/doc-ja/sample-geometry03.rb.v.rd +45 -0
- data/doc-ja/sample-geometry04.rb.v.rd +33 -0
- data/doc-ja/sample-geometry07.rb.v.rd +62 -0
- data/doc-ja/sample-groebner01.rb.v.rd +13 -0
- data/doc-ja/sample-groebner02.rb.v.rd +18 -0
- data/doc-ja/sample-groebner03.rb.v.rd +19 -0
- data/doc-ja/sample-group01.rb.v.rd +21 -0
- data/doc-ja/sample-jordan-form01.rb.v.rd +56 -0
- data/doc-ja/sample-jordanform01.rb.v.rd +55 -0
- data/doc-ja/sample-lagrange-multiplier01.rb.v.rd +35 -0
- data/doc-ja/sample-m-factorize01.rb.v.rd +15 -0
- data/doc-ja/sample-m-factorize02.rb.v.rd +12 -0
- data/doc-ja/sample-m-polynomial01.rb.v.rd +10 -0
- data/doc-ja/sample-map01.rb.v.rd +9 -0
- data/doc-ja/sample-polynomial01.rb.v.rd +9 -0
- data/doc-ja/sample-polynomial02.rb.v.rd +11 -0
- data/doc-ja/sample-primefield01.rb.v.rd +13 -0
- data/doc-ja/sample-quotientfield01.rb.v.rd +11 -0
- data/doc-ja/sample-quotientfield02.rb.v.rd +18 -0
- data/doc-ja/sample-quotientfield03.rb.v.rd +16 -0
- data/doc-ja/sample-quotientfield04.rb.v.rd +16 -0
- data/doc-ja/sample-set01.rb.v.rd +18 -0
- data/doc-ja/sample-splitting-field01.rb.v.rd +19 -0
- data/doc-ja/samples-ja.html +885 -0
- data/doc-ja/samples-ja.rd +221 -0
- data/doc-ja/style.css +85 -0
- data/doc-ja/todo.html +20 -0
- data/doc-ja/todo.rd +9 -0
- data/lib/algebra.rb +36 -0
- data/lib/algebra/algebraic-equation.rb +56 -0
- data/lib/algebra/algebraic-extension-field.rb +135 -0
- data/lib/algebra/algebraic-parser.rb +160 -0
- data/lib/algebra/algebraic-system.rb +224 -0
- data/lib/algebra/annihilate.rb +52 -0
- data/lib/algebra/array-supplement.rb +53 -0
- data/lib/algebra/auto-require.rb +84 -0
- data/lib/algebra/chinese-rem-th.rb +135 -0
- data/lib/algebra/combinatorial.rb +145 -0
- data/lib/algebra/elementary-divisor.rb +193 -0
- data/lib/algebra/euclidian-ring.rb +161 -0
- data/lib/algebra/factors.rb +305 -0
- data/lib/algebra/finite-group.rb +374 -0
- data/lib/algebra/finite-map.rb +201 -0
- data/lib/algebra/finite-set.rb +456 -0
- data/lib/algebra/galois-group.rb +129 -0
- data/lib/algebra/gaussian-elimination.rb +385 -0
- data/lib/algebra/groebner-basis-coeff.rb +228 -0
- data/lib/algebra/groebner-basis.rb +197 -0
- data/lib/algebra/import-module-single-thread.rb +86 -0
- data/lib/algebra/import-module.rb +491 -0
- data/lib/algebra/jordan-form.rb +114 -0
- data/lib/algebra/linear-algebra.rb +143 -0
- data/lib/algebra/localized-ring.rb +294 -0
- data/lib/algebra/m-index.rb +282 -0
- data/lib/algebra/m-polynomial-factor-int.rb +186 -0
- data/lib/algebra/m-polynomial-factor-zp.rb +114 -0
- data/lib/algebra/m-polynomial-factor.rb +315 -0
- data/lib/algebra/m-polynomial-gcd.rb +40 -0
- data/lib/algebra/m-polynomial.rb +875 -0
- data/lib/algebra/matrix-algebra-triplet.rb +292 -0
- data/lib/algebra/matrix-algebra.rb +929 -0
- data/lib/algebra/numeric-supplement.rb +123 -0
- data/lib/algebra/permutation-group.rb +257 -0
- data/lib/algebra/polynomial-converter.rb +193 -0
- data/lib/algebra/polynomial-factor-alg.rb +148 -0
- data/lib/algebra/polynomial-factor-int.rb +252 -0
- data/lib/algebra/polynomial-factor-zp.rb +165 -0
- data/lib/algebra/polynomial-factor.rb +140 -0
- data/lib/algebra/polynomial.rb +592 -0
- data/lib/algebra/powers.rb +31 -0
- data/lib/algebra/prime-gen.rb +32 -0
- data/lib/algebra/rational.rb +81 -0
- data/lib/algebra/residue-class-ring.rb +219 -0
- data/lib/algebra/sets-system.rb +43 -0
- data/lib/algebra/splitting-field.rb +103 -0
- data/lib/algebra/version.rb +3 -0
- data/sample/Makefile +10 -0
- data/sample/do-testscripts.rb +100 -0
- data/sample/sample-algebraic-equation01.rb +6 -0
- data/sample/sample-algebraic-equation02.rb +6 -0
- data/sample/sample-algebraic-root01.rb +15 -0
- data/sample/sample-algebraicfield01.rb +21 -0
- data/sample/sample-algebraicfield02.rb +10 -0
- data/sample/sample-cayleyhamilton01.rb +9 -0
- data/sample/sample-diagonalization01.rb +43 -0
- data/sample/sample-divmod01.rb +12 -0
- data/sample/sample-elementary-divisor01.rb +43 -0
- data/sample/sample-factorize01.rb +6 -0
- data/sample/sample-factorize02.rb +8 -0
- data/sample/sample-factorize03.rb +9 -0
- data/sample/sample-factorize04.rb +9 -0
- data/sample/sample-factorize05.rb +48 -0
- data/sample/sample-galois-group01.rb +22 -0
- data/sample/sample-gaussian-elimination01.rb +14 -0
- data/sample/sample-geometry01.rb +28 -0
- data/sample/sample-geometry02.rb +29 -0
- data/sample/sample-geometry03.rb +40 -0
- data/sample/sample-geometry04.rb +28 -0
- data/sample/sample-geometry07.rb +58 -0
- data/sample/sample-groebner01.rb +8 -0
- data/sample/sample-groebner02.rb +13 -0
- data/sample/sample-groebner03.rb +14 -0
- data/sample/sample-group01.rb +16 -0
- data/sample/sample-jordan-form01.rb +51 -0
- data/sample/sample-lagrange-multiplier01.rb +30 -0
- data/sample/sample-m-factorize01.rb +10 -0
- data/sample/sample-m-factorize02.rb +7 -0
- data/sample/sample-m-polynomial01.rb +5 -0
- data/sample/sample-map01.rb +4 -0
- data/sample/sample-polynomial01.rb +4 -0
- data/sample/sample-polynomial02.rb +6 -0
- data/sample/sample-primefield01.rb +8 -0
- data/sample/sample-quotientfield01.rb +6 -0
- data/sample/sample-quotientfield02.rb +13 -0
- data/sample/sample-quotientfield03.rb +11 -0
- data/sample/sample-quotientfield04.rb +11 -0
- data/sample/sample-set01.rb +13 -0
- data/sample/sample-splitting-field01.rb +14 -0
- data/sample/test-00-cayley-hamilton.rb +76 -0
- data/sample/test-00-groebner-basis.rb +274 -0
- data/sample/test-00-polynomial-factor-alg.rb +89 -0
- data/sample/test-00.list +4 -0
- data/sample/time-trial.rb +65 -0
- metadata +373 -0
@@ -0,0 +1,112 @@
|
|
1
|
+
<?xml version="1.0" encoding="utf-8" ?>
|
2
|
+
<!DOCTYPE html
|
3
|
+
PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
|
4
|
+
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
|
5
|
+
<html xmlns="http://www.w3.org/1999/xhtml">
|
6
|
+
<head>
|
7
|
+
<meta http-equiv="Content-type" content="text/html; charset=utf-8" />
|
8
|
+
<title>doc-ja/residue-class-ring-ja.rd</title>
|
9
|
+
<link href="style.css" type="text/css" rel="stylesheet" />
|
10
|
+
</head>
|
11
|
+
<body>
|
12
|
+
<p>[<a href="index-ja.html">index-ja</a>] </p>
|
13
|
+
<h1><a name="label-0" id="label-0">Algebra::ResidueClassRing</a></h1><!-- RDLabel: "Algebra::ResidueClassRing" -->
|
14
|
+
<p><em>(剰余環クラス)</em></p>
|
15
|
+
<p>環からその1つの元を法とした剰余環を構成します。実際のクラスを生成するには環
|
16
|
+
と法とを指定して、クラスメソッド <a href="#label-7">::create</a> あるいは
|
17
|
+
関数 <a href="#label-5">Algebra.ResidueClassRing</a>() を用います。</p>
|
18
|
+
<h2><a name="label-1" id="label-1">ファイル名:</a></h2><!-- RDLabel: "ファイル名:" -->
|
19
|
+
<ul>
|
20
|
+
<li><var>residue-class-ring.rb</var></li>
|
21
|
+
</ul>
|
22
|
+
<h2><a name="label-2" id="label-2">スーパークラス:</a></h2><!-- RDLabel: "スーパークラス:" -->
|
23
|
+
<ul>
|
24
|
+
<li><var>Object</var></li>
|
25
|
+
</ul>
|
26
|
+
<h2><a name="label-3" id="label-3">インクルードしているモジュール:</a></h2><!-- RDLabel: "インクルードしているモジュール:" -->
|
27
|
+
<p>なし</p>
|
28
|
+
<h2><a name="label-4" id="label-4">関連する関数:</a></h2><!-- RDLabel: "関連する関数:" -->
|
29
|
+
<dl>
|
30
|
+
<dt><a name="label-5" id="label-5"><code>Algebra.ResidueClassRing(<var>ring</var>, <var>mod</var>)</code></a></dt><!-- RDLabel: "Algebra.ResidueClassRing" -->
|
31
|
+
<dd>
|
32
|
+
<a href="#label-7">::create</a>(ring, mod) と同じです。</dd>
|
33
|
+
</dl>
|
34
|
+
<h2><a name="label-6" id="label-6">クラスメソッド:</a></h2><!-- RDLabel: "クラスメソッド:" -->
|
35
|
+
<dl>
|
36
|
+
<dt><a name="label-7" id="label-7"><code>::create(<var>ring</var>, <var>mod</var>)</code></a></dt><!-- RDLabel: "::create" -->
|
37
|
+
<dd>
|
38
|
+
<p>クラス <var>ring</var> で表現されるを環とその環の元 <var>mod</var> から、
|
39
|
+
その元を法とした剰余環を表現するクラスを返します。</p>
|
40
|
+
<p>この戻り値は <a href="#label-0">Algebra::ResidueClassRing</a> クラスのサブクラスです。
|
41
|
+
このサブクラスにはクラスメソッドとして <var>::ground</var> と
|
42
|
+
<var>::modulus</var> と <code>[x]</code> が定義され、それぞれ、基礎環 <var>ring</var>、
|
43
|
+
法 <var>mod</var>、x を代表元とする剰余類を返します。</p>
|
44
|
+
<p>例: 多項式環を法 <code>x**2 + x + 1</code> で割る。</p>
|
45
|
+
<pre>require "rational"
|
46
|
+
require "polynomial"
|
47
|
+
require "residue-class-ring"
|
48
|
+
Px = Algebra.Polynomial(Rational, "x")
|
49
|
+
x = Px.var
|
50
|
+
F = Algebra.ResidueClassRing(Px, x**2 + x + 1)
|
51
|
+
p F[x + 1]**100 #=> -x - 1</pre>
|
52
|
+
<p><var>ring</var> が Integer である場合に限り、全ての逆数を予め計算して
|
53
|
+
保管します。また <code>0, 1, ... , mod-1</code> に対応する剰余類の配列を
|
54
|
+
<var>to_ary</var> で得ることができます。</p>
|
55
|
+
<p>例: modulo 7 の素体</p>
|
56
|
+
<pre>require "residue-class-ring"
|
57
|
+
F7 = Algebra::ResidueClassRing.create(Integer, 7)
|
58
|
+
a, b, c, d, e, f, g = F7
|
59
|
+
p [e + c, e - c, e * c, e * 2001, 3 + c, 1/c, 1/c * c]
|
60
|
+
#=> [6, 2, 1, 3, 5, 4, 1]
|
61
|
+
p( (1...7).collect{|i| F7[i]**6} )
|
62
|
+
#=> [1, 1, 1, 1, 1, 1]</pre></dd>
|
63
|
+
<dt><a name="label-8" id="label-8"><code>::[<var>x</var>]</code></a></dt><!-- RDLabel: "::[]" -->
|
64
|
+
<dd>
|
65
|
+
<var>x</var> で代表される剰余類を返します。</dd>
|
66
|
+
<dt><a name="label-9" id="label-9"><code>::zero</code></a></dt><!-- RDLabel: "::zero" -->
|
67
|
+
<dd>
|
68
|
+
零元を返します。</dd>
|
69
|
+
<dt><a name="label-10" id="label-10"><code>::unity</code></a></dt><!-- RDLabel: "::unity" -->
|
70
|
+
<dd>
|
71
|
+
単位元を返します。</dd>
|
72
|
+
</dl>
|
73
|
+
<h2><a name="label-11" id="label-11">メソッド:</a></h2><!-- RDLabel: "メソッド:" -->
|
74
|
+
<dl>
|
75
|
+
<dt><a name="label-12" id="label-12"><code>lift</code></a></dt><!-- RDLabel: "lift" -->
|
76
|
+
<dd>
|
77
|
+
剰余類の代表元を返します。</dd>
|
78
|
+
<dt><a name="label-13" id="label-13"><code>zero?</code></a></dt><!-- RDLabel: "zero?" -->
|
79
|
+
<dd>
|
80
|
+
零元であるとき真を返します。</dd>
|
81
|
+
<dt><a name="label-14" id="label-14"><code>zero</code></a></dt><!-- RDLabel: "zero" -->
|
82
|
+
<dd>
|
83
|
+
零元を返します。</dd>
|
84
|
+
<dt><a name="label-15" id="label-15"><code>unity</code></a></dt><!-- RDLabel: "unity" -->
|
85
|
+
<dd>
|
86
|
+
単位元を返します。</dd>
|
87
|
+
<dt><a name="label-16" id="label-16"><code>==(<var>other</var>)</code></a></dt><!-- RDLabel: "==" -->
|
88
|
+
<dd>
|
89
|
+
等しいとき真を返します。</dd>
|
90
|
+
<dt><a name="label-17" id="label-17"><code>+(<var>other</var>)</code></a></dt><!-- RDLabel: "+" -->
|
91
|
+
<dd>
|
92
|
+
和を計算します。</dd>
|
93
|
+
<dt><a name="label-18" id="label-18"><code>-(<var>other</var>)</code></a></dt><!-- RDLabel: "-" -->
|
94
|
+
<dd>
|
95
|
+
差を計算します。</dd>
|
96
|
+
<dt><a name="label-19" id="label-19"><code>*(<var>other</var>)</code></a></dt><!-- RDLabel: "*" -->
|
97
|
+
<dd>
|
98
|
+
積を計算します。</dd>
|
99
|
+
<dt><a name="label-20" id="label-20"><code>**(<var>n</var>)</code></a></dt><!-- RDLabel: "**" -->
|
100
|
+
<dd>
|
101
|
+
<var>n</var> 乗を計算します。</dd>
|
102
|
+
<dt><a name="label-21" id="label-21"><code>/(<var>other</var>)</code></a></dt><!-- RDLabel: "/" -->
|
103
|
+
<dd>
|
104
|
+
<a href="#label-22">inverse</a> を利用して商を計算します。</dd>
|
105
|
+
<dt><a name="label-22" id="label-22"><code>inverse</code></a></dt><!-- RDLabel: "inverse" -->
|
106
|
+
<dd>
|
107
|
+
基礎環がユークリッド環であることを仮定して、逆数を返します。
|
108
|
+
逆数が存在しない場合の値は nil です。</dd>
|
109
|
+
</dl>
|
110
|
+
|
111
|
+
</body>
|
112
|
+
</html>
|
@@ -0,0 +1,113 @@
|
|
1
|
+
=begin
|
2
|
+
[((<index-ja|URL:index-ja.html>))]
|
3
|
+
= Algebra::ResidueClassRing
|
4
|
+
((*(剰余環クラス)*))
|
5
|
+
|
6
|
+
環からその1つの元を法とした剰余環を構成します。実際のクラスを生成するには環
|
7
|
+
と法とを指定して、クラスメソッド ((<::create>)) あるいは
|
8
|
+
関数 ((<Algebra.ResidueClassRing>))() を用います。
|
9
|
+
|
10
|
+
|
11
|
+
== ファイル名:
|
12
|
+
* ((|residue-class-ring.rb|))
|
13
|
+
|
14
|
+
== スーパークラス:
|
15
|
+
|
16
|
+
* ((|Object|))
|
17
|
+
|
18
|
+
== インクルードしているモジュール:
|
19
|
+
|
20
|
+
なし
|
21
|
+
|
22
|
+
== 関連する関数:
|
23
|
+
|
24
|
+
--- Algebra.ResidueClassRing(ring, mod)
|
25
|
+
((<::create>))(ring, mod) と同じです。
|
26
|
+
|
27
|
+
== クラスメソッド:
|
28
|
+
|
29
|
+
--- ::create(ring, mod)
|
30
|
+
|
31
|
+
クラス ((|ring|)) で表現されるを環とその環の元 ((|mod|)) から、
|
32
|
+
その元を法とした剰余環を表現するクラスを返します。
|
33
|
+
|
34
|
+
この戻り値は ((<Algebra::ResidueClassRing>)) クラスのサブクラスです。
|
35
|
+
このサブクラスにはクラスメソッドとして ((|::ground|)) と
|
36
|
+
((|::modulus|)) と (({[x]})) が定義され、それぞれ、基礎環 ((|ring|))、
|
37
|
+
法 ((|mod|))、x を代表元とする剰余類を返します。
|
38
|
+
|
39
|
+
例: 多項式環を法 (({x**2 + x + 1})) で割る。
|
40
|
+
require "rational"
|
41
|
+
require "polynomial"
|
42
|
+
require "residue-class-ring"
|
43
|
+
Px = Algebra.Polynomial(Rational, "x")
|
44
|
+
x = Px.var
|
45
|
+
F = Algebra.ResidueClassRing(Px, x**2 + x + 1)
|
46
|
+
p F[x + 1]**100 #=> -x - 1
|
47
|
+
|
48
|
+
((|ring|)) が Integer である場合に限り、全ての逆数を予め計算して
|
49
|
+
保管します。また (({0, 1, ... , mod-1})) に対応する剰余類の配列を
|
50
|
+
((|to_ary|)) で得ることができます。
|
51
|
+
|
52
|
+
例: modulo 7 の素体
|
53
|
+
require "residue-class-ring"
|
54
|
+
F7 = Algebra::ResidueClassRing.create(Integer, 7)
|
55
|
+
a, b, c, d, e, f, g = F7
|
56
|
+
p [e + c, e - c, e * c, e * 2001, 3 + c, 1/c, 1/c * c]
|
57
|
+
#=> [6, 2, 1, 3, 5, 4, 1]
|
58
|
+
p( (1...7).collect{|i| F7[i]**6} )
|
59
|
+
#=> [1, 1, 1, 1, 1, 1]
|
60
|
+
|
61
|
+
--- ::[](x)
|
62
|
+
((|x|)) で代表される剰余類を返します。
|
63
|
+
|
64
|
+
#--- ::indeterminate(obj)
|
65
|
+
|
66
|
+
--- ::zero
|
67
|
+
零元を返します。
|
68
|
+
|
69
|
+
--- ::unity
|
70
|
+
単位元を返します。
|
71
|
+
|
72
|
+
|
73
|
+
== メソッド:
|
74
|
+
|
75
|
+
--- lift
|
76
|
+
剰余類の代表元を返します。
|
77
|
+
|
78
|
+
#--- monomial?
|
79
|
+
|
80
|
+
--- zero?
|
81
|
+
零元であるとき真を返します。
|
82
|
+
|
83
|
+
--- zero
|
84
|
+
零元を返します。
|
85
|
+
|
86
|
+
--- unity
|
87
|
+
単位元を返します。
|
88
|
+
|
89
|
+
--- ==(other)
|
90
|
+
等しいとき真を返します。
|
91
|
+
|
92
|
+
--- +(other)
|
93
|
+
和を計算します。
|
94
|
+
|
95
|
+
--- -(other)
|
96
|
+
差を計算します。
|
97
|
+
|
98
|
+
--- *(other)
|
99
|
+
積を計算します。
|
100
|
+
|
101
|
+
--- **(n)
|
102
|
+
((|n|)) 乗を計算します。
|
103
|
+
|
104
|
+
--- /(other)
|
105
|
+
((<inverse>)) を利用して商を計算します。
|
106
|
+
|
107
|
+
--- inverse
|
108
|
+
基礎環がユークリッド環であることを仮定して、逆数を返します。
|
109
|
+
逆数が存在しない場合の値は nil です。
|
110
|
+
|
111
|
+
#--- to_s
|
112
|
+
|
113
|
+
=end
|
@@ -0,0 +1,11 @@
|
|
1
|
+
=begin
|
2
|
+
# sample-algebraic-equation01.rb
|
3
|
+
|
4
|
+
require "algebra"
|
5
|
+
|
6
|
+
PQ = MPolynomial(Rational)
|
7
|
+
a, b, c = PQ.vars("abc")
|
8
|
+
p AlgebraicEquation.minimal_polynomial(a + b + c, a**2-2, b**2-3, c**2-5)
|
9
|
+
#=> x^8 - 40x^6 + 352x^4 - 960x^2 + 576
|
10
|
+
((<_|CONTENTS>))
|
11
|
+
=end
|
@@ -0,0 +1,11 @@
|
|
1
|
+
=begin
|
2
|
+
# sample-algebraic-equation02.rb
|
3
|
+
|
4
|
+
require "algebra"
|
5
|
+
|
6
|
+
PQ = MPolynomial(Rational)
|
7
|
+
a, b, c = PQ.vars("abc")
|
8
|
+
p AlgebraicEquation.minimal_polynomial(a + b + c, a**2-2, b**2-3, c**2-5)
|
9
|
+
#=> x^8 - 40x^6 + 352x^4 - 960x^2 + 576
|
10
|
+
((<_|CONTENTS>))
|
11
|
+
=end
|
@@ -0,0 +1,20 @@
|
|
1
|
+
=begin
|
2
|
+
# sample-algebraic-root01.rb
|
3
|
+
|
4
|
+
require "algebra"
|
5
|
+
|
6
|
+
R2, r2, r2_ = Root(Rational, 2) # r2 = sqrt(2), -sqrt(2)
|
7
|
+
p r2 #=> sqrt(2)
|
8
|
+
R3, r3, r3_ = Root(R2, 3) # r3 = sqrt(3), -sqrt(3)
|
9
|
+
p r3 #=> sqrt(3)
|
10
|
+
R6, r6, r6_ = Root(R3, 6) # R6 = R3, r6 = sqrt(6), -sqrt(6)
|
11
|
+
p r6 #=> -r2r3
|
12
|
+
p( (r6 + r2)*(r6 - r2) ) #=> 4
|
13
|
+
|
14
|
+
F, a, b = QuadraticExtensionField(Rational){|x| x**2 - x - 1}
|
15
|
+
# fibonacci numbers
|
16
|
+
(0..100).each do |n|
|
17
|
+
puts( (a**n - b**n)/(a - b) )
|
18
|
+
end
|
19
|
+
((<_|CONTENTS>))
|
20
|
+
=end
|
@@ -0,0 +1,26 @@
|
|
1
|
+
=begin
|
2
|
+
# sample-algebraicfield01.rb
|
3
|
+
|
4
|
+
require "algebra"
|
5
|
+
|
6
|
+
Px = Polynomial(Rational, "x")
|
7
|
+
x = Px.var
|
8
|
+
F = ResidueClassRing(Px, x**2 + x + 1)
|
9
|
+
x = F[x]
|
10
|
+
p( (x + 1)**100 )
|
11
|
+
#=> -x - 1
|
12
|
+
p( (x-1)** 3 / (x**2 - 1) )
|
13
|
+
#=> -3x - 3
|
14
|
+
|
15
|
+
G = Polynomial(F, "y")
|
16
|
+
y = G.var
|
17
|
+
p( (x + y + 1)** 7 )
|
18
|
+
#=> y^7 + (7x + 7)y^6 + 8xy^5 + 4y^4 + (4x + 4)y^3 + 5xy^2 + 7y + x + 1
|
19
|
+
|
20
|
+
H = ResidueClassRing(G, y**5 + x*y + 1)
|
21
|
+
y = H[y]
|
22
|
+
p( 1/(x + y + 1)**7 )
|
23
|
+
#=> (1798/3x + 1825/9)y^4 + (-74x + 5176/9)y^3 +
|
24
|
+
# (-6886/9x - 5917/9)y^2 + (1826/3x - 3101/9)y + 2146/9x + 4702/9
|
25
|
+
((<_|CONTENTS>))
|
26
|
+
=end
|
@@ -0,0 +1,15 @@
|
|
1
|
+
=begin
|
2
|
+
# sample-algebraicfield02.rb
|
3
|
+
|
4
|
+
require "algebra"
|
5
|
+
|
6
|
+
F = AlgebraicExtensionField(Rational, "x") {|x| x**2 + x + 1}
|
7
|
+
x = F.var
|
8
|
+
p( (x + 1)**100 )
|
9
|
+
p( (x-1)** 3 / (x**2 - 1) )
|
10
|
+
|
11
|
+
H = AlgebraicExtensionField(F, "y") {|y| y**5 + x*y + 1}
|
12
|
+
y = H.var
|
13
|
+
p( 1/(x + y + 1)**7 )
|
14
|
+
((<_|CONTENTS>))
|
15
|
+
=end
|
@@ -0,0 +1,14 @@
|
|
1
|
+
=begin
|
2
|
+
# sample-cayleyhamilton01.rb
|
3
|
+
|
4
|
+
require "algebra"
|
5
|
+
|
6
|
+
n = 4
|
7
|
+
R = MPolynomial(Integer)
|
8
|
+
MR = SquareMatrix(R, n)
|
9
|
+
m = MR.matrix{|i, j| R.var("x#{i}#{j}") }
|
10
|
+
Rx = Polynomial(R, "x")
|
11
|
+
ch = m.char_polynomial(Rx)
|
12
|
+
p ch.evaluate(m) #=> 0
|
13
|
+
((<_|CONTENTS>))
|
14
|
+
=end
|
@@ -0,0 +1,48 @@
|
|
1
|
+
=begin
|
2
|
+
# sample-diagonalization01.rb
|
3
|
+
|
4
|
+
require "algebra"
|
5
|
+
M = SquareMatrix(Rational, 3)
|
6
|
+
a = M[[1,-1,-1], [-1,1,-1], [2,1,-1]]
|
7
|
+
puts "A = "; a.display; puts
|
8
|
+
#A =
|
9
|
+
# 1, -1, -1
|
10
|
+
# -1, 1, -1
|
11
|
+
# 2, 1, -1
|
12
|
+
|
13
|
+
e = a.diagonalize
|
14
|
+
|
15
|
+
puts "Charactoristic Poly.: #{e.chpoly} => #{e.facts}";puts
|
16
|
+
#Charactoristic Poly.: t^3 - t^2 + t - 6 => (t - 2)(t^2 + t + 3)
|
17
|
+
|
18
|
+
puts "Algebraic Numbers:"
|
19
|
+
e.roots.each do |po, rs|
|
20
|
+
puts "#{rs.join(', ')} : roots of #{po} == 0"
|
21
|
+
end; puts
|
22
|
+
#Algebraic Numbers:
|
23
|
+
#a, -a - 1 : roots of t^2 + t + 3 == 0
|
24
|
+
|
25
|
+
puts "EigenSpaces: "
|
26
|
+
e.evalues.uniq.each do |ev|
|
27
|
+
puts "W_{#{ev}} = <#{e.espaces[ev].join(', ')}>"
|
28
|
+
end; puts
|
29
|
+
#EigenSpaces:
|
30
|
+
#W_{2} = <[4, -5, 1]>
|
31
|
+
#W_{a} = <[1/3a + 1/3, 1/3a + 1/3, 1]>
|
32
|
+
#W_{-a - 1} = <[-1/3a, -1/3a, 1]>
|
33
|
+
|
34
|
+
puts "Trans. Matrix:"
|
35
|
+
puts "P ="
|
36
|
+
e.tmatrix.display; puts
|
37
|
+
puts "P^-1 * A * P = "; (e.tmatrix.inverse * a * e.tmatrix).display; puts
|
38
|
+
#P =
|
39
|
+
# 4, 1/3a + 1/3, -1/3a
|
40
|
+
# -5, 1/3a + 1/3, -1/3a
|
41
|
+
# 1, 1, 1
|
42
|
+
#
|
43
|
+
#P^-1 * A * P =
|
44
|
+
# 2, 0, 0
|
45
|
+
# 0, a, 0
|
46
|
+
# 0, 0, -a - 1
|
47
|
+
((<_|CONTENTS>))
|
48
|
+
=end
|
@@ -0,0 +1,17 @@
|
|
1
|
+
=begin
|
2
|
+
# sample-divmod01.rb
|
3
|
+
|
4
|
+
require "algebra"
|
5
|
+
P = MPolynomial(Rational)
|
6
|
+
x, y, z = P.vars("xyz")
|
7
|
+
f = x**2*y + x*y**2 + y*2 + z**3
|
8
|
+
g = x*y-z**3
|
9
|
+
h = y*2-6*z
|
10
|
+
|
11
|
+
P.set_ord(:lex) # lex, grlex, grevlex
|
12
|
+
puts "(#{f}).divmod([#{g}, #{h}]) =>", "#{f.divmod(g, h).inspect}"
|
13
|
+
#=> (x^2y + xy^2 + 2y + z^3).divmod([xy - z^3, 2y - 6z]) =>
|
14
|
+
# [[x + y, 1/2z^3 + 1], xz^3 + 3z^4 + z^3 + 6z]
|
15
|
+
# = [[Quotient1,Quotient2], Remainder]
|
16
|
+
((<_|CONTENTS>))
|
17
|
+
=end
|
@@ -0,0 +1,48 @@
|
|
1
|
+
=begin
|
2
|
+
# sample-elementary-divisor01.rb
|
3
|
+
|
4
|
+
require "algebra"
|
5
|
+
|
6
|
+
M = SquareMatrix(Rational, 4)
|
7
|
+
a = M[
|
8
|
+
[2, 0, 0, 0],
|
9
|
+
[0, 2, 0, 0],
|
10
|
+
[0, 0, 2, 0],
|
11
|
+
[5, 0, 0, 2]
|
12
|
+
]
|
13
|
+
P = Polynomial(Rational, "x")
|
14
|
+
MP = SquareMatrix(P, 4)
|
15
|
+
|
16
|
+
ac = a._char_matrix(MP)
|
17
|
+
ac.display; puts #=>
|
18
|
+
#x - 2, 0, 0, 0
|
19
|
+
# 0, x - 2, 0, 0
|
20
|
+
# 0, 0, x - 2, 0
|
21
|
+
# -5, 0, 0, x - 2
|
22
|
+
|
23
|
+
p ac.elementary_divisor #=> [1, x - 2, x - 2, x^2 - 4x + 4]
|
24
|
+
|
25
|
+
require "algebra/matrix-algebra-triplet"
|
26
|
+
at = ac.to_triplet.e_diagonalize
|
27
|
+
|
28
|
+
at.body.display; puts #=>
|
29
|
+
# 1, 0, 0, 0
|
30
|
+
# 0, x - 2, 0, 0
|
31
|
+
# 0, 0, x - 2, 0
|
32
|
+
# 0, 0, 0, x^2 - 4x + 4
|
33
|
+
|
34
|
+
at.left.display; puts #=>
|
35
|
+
# 0, 0, 0, -1/5
|
36
|
+
# 0, 1, 0, 0
|
37
|
+
# 0, 0, 1, 0
|
38
|
+
# 5, 0, 0, x - 2
|
39
|
+
|
40
|
+
at.right.display; puts #=>
|
41
|
+
# 1, 0, 0, 1/5x - 2/5
|
42
|
+
# 0, 1, 0, 0
|
43
|
+
# 0, 0, 1, 0
|
44
|
+
# 0, 0, 0, 1
|
45
|
+
|
46
|
+
p at.left * ac * at.right == at.body #=> true
|
47
|
+
((<_|CONTENTS>))
|
48
|
+
=end
|