Ruppert matrix as subresultant mapping
|
|
- Marsha King
- 5 years ago
- Views:
Transcription
1 Ruppert matrix as subresultant mapping Kosaku Nagasaka Kobe University JAPAN This presentation is powered by Mathematica : 29 : Ruppert matrix as subresultant mapping Prev Next
2 2 CASC2007slideshow.nb Outline of this talk How to compute polynomial GCDs Sylvester matrix and Subresultant mapping Ruppert matrix and factoring polynomials Ruppert matrix for GCD of two polynomials Ruppert matrix for GCD of several polynomials Conclusion 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
3 CASC2007slideshow.nb 3 How to compute polynomial GCDs Euclidean Algorithm Eculidean algorithm for polynomials f x and gx over reals. function GCD(f(x), g(x)) if g(x) = 0 return f(x) else return GCD( g(x), remainder of f(x) by g(x) ) QR factoring of Sylvester matrix The usual subresultant mapping is important. See the next slide... Lots of other powerful approximate algorithms 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
4 4 CASC2007slideshow.nb QR factoring of Sylvester matrix Sylvester matrix of the following f x and gx : f x f n x n f n1 x n1 f 1 x f 0, gx g m x m g m1 x m1 g 1 x g 0. f n f n1 f 1 f f n f n1 f 1 f 0 0 Syl f, g 0 0 f n f n1 f 1 f 0 g m g m1 g 1 g g m g m1 g 1 g g m g m1 g 1 g 0 Well known relation between GCD and QR factoring of Sylvester matrix. The last non-zero row vector of R where Syl f, g Q R, is the coefficient vector of GCD of f x and gx. (See: Laidacker, M.A., 1969 and others) 09 : 29 : Ruppert matrix as subresultant mapping Prev Next.
5 CASC2007slideshow.nb 5 Subresultant mapping Subresultant mapping S r of f 0 x and f 1 x S r : n 1 r1 n0 r1 n0 n 1 r1 u 0, u 1 u 1 f 0 u 0 f 1 where f 0 x f 0,n0 x n 0 f0,1 x f 0,0 f 1 x f 1,n1 x n 1 f1,1 x f 1,0 Sylvester subresultant matrix S r f 0, f 1 S r f 0, f 1 C n0 r f 1 C n1 r f 0 where C k p is a k-th convolution matrix of px and S 0 f 0, f 1 is the transpose of Sylvester matrix of f 0 x and f 1 x. 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
6 6 CASC2007slideshow.nb Convolution matrix k-th convolution matrix of polynomial px. p n p n1 p n p n1 0 C k p p 0 p n 0 0 p 0 p n1 p n 0 p n1 p p 0 where px p n x n p 1 x p 0. Subresultant mapping and Sylvester matrix n1 r1 n0 r1 n0 n 1 r1 S r : u 0, u 1 u 1 f 0 u 0 f 1 S r f 0, f 1 C n0 r f 1 C n1 r f 0, S 0 f 0, f 1 t Syl f 0, f 1 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
7 CASC2007slideshow.nb 7 Subresultant mapping and GCD Well known fact (see Rupprecht, D. (1999) and so on) (1) If r is the largest integer that the subresultant mapping S r is not injective, we can compute the GCD of f 0 x and f 1 x from the right null vector of S r f 0, f 1. (2) The dimension of the null space is the degree of the polynomial GCD. 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
8 8 CASC2007slideshow.nb Ruppert criterion Reducibility of bivariate polynomials (Ruppert, W.M., 1999) f x, y is absolutely irreducible the following differential equation does not have any non-trivial solutions gx, y and hx, y. g f f h f y g y h x f x 0, g, h x, y with degree constraints: deg x g deg x f 1, deg y g deg y f, deg x h deg x f, deg y h deg y f 2 Newton polytope ver. by Gao, S. and Rodrigues, V.M. (2003) Multivariate version by May, J.P. (2005) 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
9 CASC2007slideshow.nb 9 Ruppert matrix Rf The differential equation w.r.t. gx, y and hx, y a linear equation w.r.t. unknown coefficients. g f y f g y f h x h f x 0, g, h x, y, deg x g deg x f 1, deg y g deg y f, deg x h deg x f, deg y h deg y f 2 R f t coefficient vectors of g, h SVD of Ruppert matrix: irreducibility radius. Ruppert matrix depends term orders: lexicographic in this talk. 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
10 10 CASC2007slideshow.nb Factoring with Ruppert and Gao s matrix Known Algorithms using Corollary 1 Corollary 1 For a given f x, y x, y that is square-free over y, the dimension (over ) of the null space of R f is equal to "(the number of absolutely irreducible factors of f x, y over ) 1". Correction: Corollary 1 from John May s thesis. and others by Kaltofen, Gao and so on. Kaltofen, E., May, J., Yang, Z., and Zhi, L. Approximate factorization of multivariate polynomials using singular value decomposition. Manuscript, 22 pages. Submitted, Jan : 29 : Ruppert matrix as subresultant mapping Prev Next
11 CASC2007slideshow.nb 11 Purpose of this talk GCD and Factorization Primitive operations in Symbolic (and Numeric) computations. Common Properties Can be computed by some matrix decompositions. Sylvester and Ruppert matices Is there any relation between them? 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
12 12 CASC2007slideshow.nb How to link them f x, y f 0 x f 1 xy f x, y is reducible if f 0 x and f 1 x have a non-trivial GCD. f x, y is irreducible if they don t have any non-trivial GCD. Coprimeness of f 0 x and f 1 x can be tested by a rank deficiency of their Sylvester matrix. Irreducibility of f x, y f 0 x f 1 xy can be tested by a rank deficiency of its Ruppert matrix. 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
13 CASC2007slideshow.nb 13 Simple Result Lemma 1 For any polynomials f 0 x and f 1 x, the Sylvester matrix of f 0 x and f 1 x and the Ruppert matrix of f x, y f 0 x f 1 xy have the same information for computing their GCD, with the Ruppert s original differential equation and degree constraints. The proof is in the proceedings. 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
14 14 CASC2007slideshow.nb Simple Result (Sylvester and Ruppert) fa5 fa4 fa3 fa2 fa1 fa fa5 fa4 fa3 fa2 fa1 fa fa5 fa4 fa3 fa2 fa1 fa fa5 fa4 fa3 fa2 fa1 fa fa5 fa4 fa3 fa2 fa1 fa0 fb5 fb4 fb3 fb2 fb1 fb fb5 fb4 fb3 fb2 fb1 fb fb5 fb4 fb3 fb2 fb1 fb fb5 fb4 fb3 fb2 fb1 fb fb5 fb4 fb3 fb2 fb1 fb0 0 fa5 0 fa4 0 fa3 0 fa2 0 fa1 0 fa fb5 0 fb4 0 fb3 0 fb2 0 fb1 0 fb fa5 0 fa4 0 fa3 0 fa2 0 fa1 0 fa fb5 0 fb4 0 fb3 0 fb2 0 fb1 0 fb fa5 0 fa4 0 fa3 0 fa2 0 fa1 0 fa fb5 0 fb4 0 fb3 0 fb2 0 fb1 0 fb fa5 0 fa4 0 fa3 0 fa2 0 fa1 0 fa fb5 0 fb4 0 fb3 0 fb2 0 fb1 0 fb fa5 0 fa4 0 fa3 0 fa2 0 fa1 0 fa fb5 0 fb4 0 fb3 0 fb2 0 fb1 0 fb0 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
15 CASC2007slideshow.nb 15 Alternative Result 1 Theorem 1 The polynomial GCD of f 0 x and f 1 x can be computed by Singular Value Decomposition (SVD) of Ruppert matrix of f x, y f 0 x f 1 xy with the J.P.May s differential equation and degree constraints, if f x, y is square-free over y. Difference between Ruppert and May s f g y f h x h f x 0, g deg x f 1, deg y g deg y f, h deg x f, deg y h deg y f 2 The proof is in the proceedings. 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
16 16 CASC2007slideshow.nb Alternative Result 1(Ruppert with May) fb4 fb fb3 0 2 fb fb2 fb3 fb4 3 fb fb1 2 fb2 0 2 fb4 4 fb fb0 3 fb1 fb2 fb3 3 fb4 5 fb fb0 2 fb1 0 2 fb3 4 fb fb0 fb1 fb2 3 fb fb0 0 2 fb fb0 fb fa4 fa fa fb fa3 0 2 fa fa4 fa5 0 0 fb4 fb fa2 fa3 fa4 3 fa5 0 0 fa3 fa4 fa5 0 fb3 fb4 fb5 0 4 fa1 2 fa2 0 2 fa4 4 fa5 0 fa2 fa3 fa4 fa5 fb2 fb3 fb4 fb5 5 fa0 3 fa1 fa2 fa3 3 fa4 5 fa5 fa1 fa2 fa3 fa4 fb1 fb2 fb3 fb4 0 4 fa0 2 fa1 0 2 fa3 4 fa4 fa0 fa1 fa2 fa3 fb0 fb1 fb2 fb fa0 fa1 fa2 3 fa3 0 fa0 fa1 fa2 0 fb0 fb1 fb fa0 0 2 fa2 0 0 fa0 fa1 0 0 fb0 fb fa0 fa fa fb0 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
17 CASC2007slideshow.nb 17 Alternative Result 2 Theorem 2 The polynomial GCD of f 0 x and f 1 x can be computed by QR factoring of the transpose of the last 3n 0 rows of their Ruppert matrix R f R f 0 x f 1 xy with J.P.May s equation. The last non-zero row vector of the triangular matrix is the coefficient vector of their polynomial GCD. The proof is in the proceedings. 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
18 18 CASC2007slideshow.nb Alternative Result 2(Last 3n 0 rows) fb4 fb fb3 0 2 fb fb2 fb3 fb4 3 fb fb1 2 fb2 0 2 fb4 4 fb fb0 3 fb1 fb2 fb3 3 fb4 5 fb fb0 2 fb1 0 2 fb3 4 fb fb0 fb1 fb2 3 fb fb0 0 2 fb fb0 fb fa4 fa fa fb fa3 0 2 fa fa4 fa5 0 0 fb4 fb fa2 fa3 fa4 3 fa5 0 0 fa3 fa4 fa5 0 fb3 fb4 fb5 0 4 fa1 2 fa2 0 2 fa4 4 fa5 0 fa2 fa3 fa4 fa5 fb2 fb3 fb4 fb5 5 fa0 3 fa1 fa2 fa3 3 fa4 5 fa5 fa1 fa2 fa3 fa4 fb1 fb2 fb3 fb4 0 4 fa0 2 fa1 0 2 fa3 4 fa4 fa0 fa1 fa2 fa3 fb0 fb1 fb2 fb fa0 fa1 fa2 3 fa3 0 fa0 fa1 fa2 0 fb0 fb1 fb fa0 0 2 fa2 0 0 fa0 fa1 0 0 fb0 fb fa0 fa fa fb0 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
19 CASC2007slideshow.nb 19 Generalized Sylvester matrix Generalized Subresultant mapping S r of f 0 x,, f k x. k k ni r1 n0 n i r1 i 0 i 1 S r : u 0 u 1 f 0 u 0 f 1, n i deg x f i x u k u k f 0 u 0 f k Generalized Sylvester matrix S r f 0,, f k. S r f 0,, f k C n0 r f 1 C n1 r f 0 (see Rupprecht, D. (1999) and so on) C n0 r f 2 C n2 r f 0 C n0 r f k C nk r f 0 If r is the largest integer that the subresultant mapping S r is not injective, we can compute the GCD of f 0 x,, f k x from the right null vector of S r f 0,, f k. 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
20 20 CASC2007slideshow.nb How to link them f x, f 0 x f 1 xy 1 f k xy k f x, is reducible if f 0 x,, f k x have a non-trivial GCD. f x, is irreducible if they don t have any non-trivial GCD. Coprimeness of f 0 x,, f k x can be tested by a rank deficiency of their Sylvester matrix. Irreducibility of f x, f 0 x k i 1 by a rank deficiency of its Ruppert matrix. f i xy i can be tested 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
21 CASC2007slideshow.nb 21 Relation for several polynomials Postponed as a future work: J.P.May s differential equation becomes f 0 g i g 0 f i Λ i f i f 0 f 0 f i 0 f 1 g i g 1 f i Λ i f i f 1 f 1 f i 0 i 1,, k f k g i g k f i Λ i f i f k f k f i 0 This means that u 0 x g 0 x Λ i f 0 u i x g i x Λ i f i i 1,, k will be a part of proof of several polynomial version of theorem. 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
22 22 CASC2007slideshow.nb Conclusion Sylvester and Ruppert matrices have some relations for computing GCDs. Their relations lead to no algorithm which computes GCDs efficiently. The author hopes that the relations in this talk will make some progress. Thank you 09 : 29 : Ruppert matrix as subresultant mapping Prev Next
An Algorithm for Approximate Factorization of Bivariate Polynomials 1)
MM Research Preprints, 402 408 MMRC, AMSS, Academia Sinica No. 22, December 2003 An Algorithm for Approximate Factorization of Bivariate Polynomials 1) Zhengfeng Yang and Lihong Zhi 2) Abstract. In this
More informationMath 4310 Solutions to homework 7 Due 10/27/16
Math 4310 Solutions to homework 7 Due 10/27/16 1. Find the gcd of x 3 + x 2 + x + 1 and x 5 + 2x 3 + x 2 + x + 1 in Rx. Use the Euclidean algorithm: x 5 + 2x 3 + x 2 + x + 1 = (x 3 + x 2 + x + 1)(x 2 x
More informationg(x) = 1 1 x = 1 + x + x2 + x 3 + is not a polynomial, since it doesn t have finite degree. g(x) is an example of a power series.
6 Polynomial Rings We introduce a class of rings called the polynomial rings, describing computation, factorization and divisibility in such rings For the case where the coefficients come from an integral
More informationPolynomials. Chapter 4
Chapter 4 Polynomials In this Chapter we shall see that everything we did with integers in the last Chapter we can also do with polynomials. Fix a field F (e.g. F = Q, R, C or Z/(p) for a prime p). Notation
More informationCHAPTER I. Rings. Definition A ring R is a set with two binary operations, addition + and
CHAPTER I Rings 1.1 Definitions and Examples Definition 1.1.1. A ring R is a set with two binary operations, addition + and multiplication satisfying the following conditions for all a, b, c in R : (i)
More informationIRREDUCIBILITY OF POLYNOMIALS MODULO p VIA NEWTON POLYTOPES. 1. Introduction
IRREDUCIBILITY OF POLYNOMIALS MODULO p VIA NEWTON POLYTOPES SHUHONG GAO AND VIRGÍNIA M. RODRIGUES Abstract. Ostrowski established in 1919 that an absolutely irreducible integral polynomial remains absolutely
More informationComputing Approximate GCD of Univariate Polynomials by Structure Total Least Norm 1)
MM Research Preprints, 375 387 MMRC, AMSS, Academia Sinica No. 24, December 2004 375 Computing Approximate GCD of Univariate Polynomials by Structure Total Least Norm 1) Lihong Zhi and Zhengfeng Yang Key
More informationMathematical Olympiad Training Polynomials
Mathematical Olympiad Training Polynomials Definition A polynomial over a ring R(Z, Q, R, C) in x is an expression of the form p(x) = a n x n + a n 1 x n 1 + + a 1 x + a 0, a i R, for 0 i n. If a n 0,
More informationPublic-key Cryptography: Theory and Practice
Public-key Cryptography Theory and Practice Department of Computer Science and Engineering Indian Institute of Technology Kharagpur Chapter 2: Mathematical Concepts Divisibility Congruence Quadratic Residues
More informationFast algorithms for factoring polynomials A selection. Erich Kaltofen North Carolina State University google->kaltofen google->han lu
Fast algorithms for factoring polynomials A selection Erich Kaltofen North Carolina State University google->kaltofen google->han lu Overview of my work 1 Theorem. Factorization in K[x] is undecidable
More informationMATH 431 PART 2: POLYNOMIAL RINGS AND FACTORIZATION
MATH 431 PART 2: POLYNOMIAL RINGS AND FACTORIZATION 1. Polynomial rings (review) Definition 1. A polynomial f(x) with coefficients in a ring R is n f(x) = a i x i = a 0 + a 1 x + a 2 x 2 + + a n x n i=0
More informationD-MATH Algebra I HS18 Prof. Rahul Pandharipande. Solution 6. Unique Factorization Domains
D-MATH Algebra I HS18 Prof. Rahul Pandharipande Solution 6 Unique Factorization Domains 1. Let R be a UFD. Let that a, b R be coprime elements (that is, gcd(a, b) R ) and c R. Suppose that a c and b c.
More informationMTH310 EXAM 2 REVIEW
MTH310 EXAM 2 REVIEW SA LI 4.1 Polynomial Arithmetic and the Division Algorithm A. Polynomial Arithmetic *Polynomial Rings If R is a ring, then there exists a ring T containing an element x that is not
More informationThe Sylvester Resultant
Lecture 10 The Sylvester Resultant We want to compute intersections of algebraic curves F and G. Let F and G be the vanishing sets of f(x,y) and g(x, y), respectively. Algebraically, we are interested
More informationSolutions of exercise sheet 11
D-MATH Algebra I HS 14 Prof Emmanuel Kowalski Solutions of exercise sheet 11 The content of the marked exercises (*) should be known for the exam 1 For the following values of α C, find the minimal polynomial
More information0 Sets and Induction. Sets
0 Sets and Induction Sets A set is an unordered collection of objects, called elements or members of the set. A set is said to contain its elements. We write a A to denote that a is an element of the set
More informationSimultaneous Linear, and Non-linear Congruences
Simultaneous Linear, and Non-linear Congruences CIS002-2 Computational Alegrba and Number Theory David Goodwin david.goodwin@perisic.com 09:00, Friday 18 th November 2011 Outline 1 Polynomials 2 Linear
More informationChapter 7 Polynomial Functions. Factoring Review. We will talk about 3 Types: ALWAYS FACTOR OUT FIRST! Ex 2: Factor x x + 64
Chapter 7 Polynomial Functions Factoring Review We will talk about 3 Types: 1. 2. 3. ALWAYS FACTOR OUT FIRST! Ex 1: Factor x 2 + 5x + 6 Ex 2: Factor x 2 + 16x + 64 Ex 3: Factor 4x 2 + 6x 18 Ex 4: Factor
More informationStructured Matrix Methods Computing the Greatest Common Divisor of Polynomials
Spec Matrices 2017; 5:202 224 Research Article Open Access Dimitrios Christou, Marilena Mitrouli*, and Dimitrios Triantafyllou Structured Matrix Methods Computing the Greatest Common Divisor of Polynomials
More informationThe Mathematica Journal Symbolic-Numeric Algebra for Polynomials
This article has not been updated for Mathematica 8. The Mathematica Journal Symbolic-Numeric Algebra for Polynomials Kosaku Nagasaka Symbolic-Numeric Algebra for Polynomials (SNAP) is a prototype package
More informationMath 4320 Final Exam
Math 4320 Final Exam 2:00pm 4:30pm, Friday 18th May 2012 Symmetry, as wide or as narrow as you may define its meaning, is one idea by which man through the ages has tried to comprehend and create order,
More informationChapter 4. Greatest common divisors of polynomials. 4.1 Polynomial remainder sequences
Chapter 4 Greatest common divisors of polynomials 4.1 Polynomial remainder sequences If K is a field, then K[x] is a Euclidean domain, so gcd(f, g) for f, g K[x] can be computed by the Euclidean algorithm.
More informationCS 829 Polynomial systems: geometry and algorithms Lecture 3: Euclid, resultant and 2 2 systems Éric Schost
CS 829 Polynomial systems: geometry and algorithms Lecture 3: Euclid, resultant and 2 2 systems Éric Schost eschost@uwo.ca Summary In this lecture, we start actual computations (as opposed to Lectures
More informationx 3 2x = (x 2) (x 2 2x + 1) + (x 2) x 2 2x + 1 = (x 4) (x + 2) + 9 (x + 2) = ( 1 9 x ) (9) + 0
1. (a) i. State and prove Wilson's Theorem. ii. Show that, if p is a prime number congruent to 1 modulo 4, then there exists a solution to the congruence x 2 1 mod p. (b) i. Let p(x), q(x) be polynomials
More informationLattice reduction of polynomial matrices
Lattice reduction of polynomial matrices Arne Storjohann David R. Cheriton School of Computer Science University of Waterloo Presented at the SIAM conference on Applied Algebraic Geometry at the Colorado
More informationϕ : Z F : ϕ(t) = t 1 =
1. Finite Fields The first examples of finite fields are quotient fields of the ring of integers Z: let t > 1 and define Z /t = Z/(tZ) to be the ring of congruence classes of integers modulo t: in practical
More informationNOTES ON DIOPHANTINE APPROXIMATION
NOTES ON DIOPHANTINE APPROXIMATION Jan-Hendrik Evertse January 29, 200 9 p-adic Numbers Literature: N. Koblitz, p-adic Numbers, p-adic Analysis, and Zeta-Functions, 2nd edition, Graduate Texts in Mathematics
More informationCHAPTER 10: POLYNOMIALS (DRAFT)
CHAPTER 10: POLYNOMIALS (DRAFT) LECTURE NOTES FOR MATH 378 (CSUSM, SPRING 2009). WAYNE AITKEN The material in this chapter is fairly informal. Unlike earlier chapters, no attempt is made to rigorously
More informationParity of the Number of Irreducible Factors for Composite Polynomials
Parity of the Number of Irreducible Factors for Composite Polynomials Ryul Kim Wolfram Koepf Abstract Various results on parity of the number of irreducible factors of given polynomials over finite fields
More informationLecture Notes Math 371: Algebra (Fall 2006) by Nathanael Leedom Ackerman
Lecture Notes Math 371: Algebra (Fall 2006) by Nathanael Leedom Ackerman October 17, 2006 TALK SLOWLY AND WRITE NEATLY!! 1 0.1 Factorization 0.1.1 Factorization of Integers and Polynomials Now we are going
More informationMultivariate Statistical Analysis
Multivariate Statistical Analysis Fall 2011 C. L. Williams, Ph.D. Lecture 4 for Applied Multivariate Analysis Outline 1 Eigen values and eigen vectors Characteristic equation Some properties of eigendecompositions
More informationTensor Decompositions for Machine Learning. G. Roeder 1. UBC Machine Learning Reading Group, June University of British Columbia
Network Feature s Decompositions for Machine Learning 1 1 Department of Computer Science University of British Columbia UBC Machine Learning Group, June 15 2016 1/30 Contact information Network Feature
More informationLecture 6. Numerical methods. Approximation of functions
Lecture 6 Numerical methods Approximation of functions Lecture 6 OUTLINE 1. Approximation and interpolation 2. Least-square method basis functions design matrix residual weighted least squares normal equation
More informationExplicit Criterion to Determine the Number of Positive Roots of a Polynomial 1)
MM Research Preprints, 134 145 No. 15, April 1997. Beijing Explicit Criterion to Determine the Number of Positive Roots of a Polynomial 1) Lu Yang and Bican Xia 2) Abstract. In a recent article, a complete
More informationLecture 7: Polynomial rings
Lecture 7: Polynomial rings Rajat Mittal IIT Kanpur You have seen polynomials many a times till now. The purpose of this lecture is to give a formal treatment to constructing polynomials and the rules
More informationCourse 2316 Sample Paper 1
Course 2316 Sample Paper 1 Timothy Murphy April 19, 2015 Attempt 5 questions. All carry the same mark. 1. State and prove the Fundamental Theorem of Arithmetic (for N). Prove that there are an infinity
More informationChinese Remainder Theorem
Chinese Remainder Theorem Theorem Let R be a Euclidean domain with m 1, m 2,..., m k R. If gcd(m i, m j ) = 1 for 1 i < j k then m = m 1 m 2 m k = lcm(m 1, m 2,..., m k ) and R/m = R/m 1 R/m 2 R/m k ;
More informationcomputing an irreducible decomposition of A B looking for solutions where the Jacobian drops rank
Diagonal Homotopies 1 Intersecting Solution Sets computing an irreducible decomposition of A B 2 The Singular Locus looking for solutions where the Jacobian drops rank 3 Solving Systems restricted to an
More informationMATH FINAL EXAM REVIEW HINTS
MATH 109 - FINAL EXAM REVIEW HINTS Answer: Answer: 1. Cardinality (1) Let a < b be two real numbers and define f : (0, 1) (a, b) by f(t) = (1 t)a + tb. (a) Prove that f is a bijection. (b) Prove that any
More informationBasic Algorithms in Number Theory
Basic Algorithms in Number Theory Algorithmic Complexity... 1 Basic Algorithms in Number Theory Francesco Pappalardi #2 - Discrete Logs, Modular Square Roots, Polynomials, Hensel s Lemma & Chinese Remainder
More informationChapter 5. Linear Algebra
Chapter 5 Linear Algebra The exalted position held by linear algebra is based upon the subject s ubiquitous utility and ease of application. The basic theory is developed here in full generality, i.e.,
More informationFast Polynomial Multiplication
Fast Polynomial Multiplication Marc Moreno Maza CS 9652, October 4, 2017 Plan Primitive roots of unity The discrete Fourier transform Convolution of polynomials The fast Fourier transform Fast convolution
More informationGalois Fields and Hardware Design
Galois Fields and Hardware Design Construction of Galois Fields, Basic Properties, Uniqueness, Containment, Closure, Polynomial Functions over Galois Fields Priyank Kalla Associate Professor Electrical
More informationChapter 4. Remember: F will always stand for a field.
Chapter 4 Remember: F will always stand for a field. 4.1 10. Take f(x) = x F [x]. Could there be a polynomial g(x) F [x] such that f(x)g(x) = 1 F? Could f(x) be a unit? 19. Compare with Problem #21(c).
More informationarxiv: v3 [math.ac] 11 May 2016
GPGCD: An iterative method for calculating approximate GCD of univariate polynomials arxiv:12070630v3 [mathac] 11 May 2016 Abstract Akira Terui Faculty of Pure and Applied Sciences University of Tsukuba
More informationFinite Fields. Mike Reiter
1 Finite Fields Mike Reiter reiter@cs.unc.edu Based on Chapter 4 of: W. Stallings. Cryptography and Network Security, Principles and Practices. 3 rd Edition, 2003. Groups 2 A group G, is a set G of elements
More informationMath 547, Exam 2 Information.
Math 547, Exam 2 Information. 3/19/10, LC 303B, 10:10-11:00. Exam 2 will be based on: Homework and textbook sections covered by lectures 2/3-3/5. (see http://www.math.sc.edu/ boylan/sccourses/547sp10/547.html)
More informationGroebner Bases, Toric Ideals and Integer Programming: An Application to Economics. Tan Tran Junior Major-Economics& Mathematics
Groebner Bases, Toric Ideals and Integer Programming: An Application to Economics Tan Tran Junior Major-Economics& Mathematics History Groebner bases were developed by Buchberger in 1965, who later named
More informationSylvester Matrix and GCD for Several Univariate Polynomials
Sylvester Matrix and GCD for Several Univariate Polynomials Manuela Wiesinger-Widi Doctoral Program Computational Mathematics Johannes Kepler University Linz 4040 Linz, Austria manuela.wiesinger@dk-compmath.jku.at
More informationLecture 02 Linear Algebra Basics
Introduction to Computational Data Analysis CX4240, 2019 Spring Lecture 02 Linear Algebra Basics Chao Zhang College of Computing Georgia Tech These slides are based on slides from Le Song and Andres Mendez-Vazquez.
More informationClass Notes; Week 7, 2/26/2016
Class Notes; Week 7, 2/26/2016 Day 18 This Time Section 3.3 Isomorphism and Homomorphism [0], [2], [4] in Z 6 + 0 4 2 0 0 4 2 4 4 2 0 2 2 0 4 * 0 4 2 0 0 0 0 4 0 4 2 2 0 2 4 So {[0], [2], [4]} is a subring.
More informationFurther linear algebra. Chapter II. Polynomials.
Further linear algebra. Chapter II. Polynomials. Andrei Yafaev 1 Definitions. In this chapter we consider a field k. Recall that examples of felds include Q, R, C, F p where p is prime. A polynomial is
More informationComputer Algebra: General Principles
Computer Algebra: General Principles For article on related subject see SYMBOL MANIPULATION. Computer algebra is a branch of scientific computation. There are several characteristic features that distinguish
More informationLinear and Bilinear Algebra (2WF04) Jan Draisma
Linear and Bilinear Algebra (2WF04) Jan Draisma CHAPTER 3 The minimal polynomial and nilpotent maps 3.1. Minimal polynomial Throughout this chapter, V is a finite-dimensional vector space of dimension
More informationOutline. MSRI-UP 2009 Coding Theory Seminar, Week 2. The definition. Link to polynomials
Outline MSRI-UP 2009 Coding Theory Seminar, Week 2 John B. Little Department of Mathematics and Computer Science College of the Holy Cross Cyclic Codes Polynomial Algebra More on cyclic codes Finite fields
More informationExample Linear Algebra Competency Test
Example Linear Algebra Competency Test The 4 questions below are a combination of True or False, multiple choice, fill in the blank, and computations involving matrices and vectors. In the latter case,
More information12x + 18y = 30? ax + by = m
Math 2201, Further Linear Algebra: a practical summary. February, 2009 There are just a few themes that were covered in the course. I. Algebra of integers and polynomials. II. Structure theory of one endomorphism.
More informationSection 33 Finite fields
Section 33 Finite fields Instructor: Yifan Yang Spring 2007 Review Corollary (23.6) Let G be a finite subgroup of the multiplicative group of nonzero elements in a field F, then G is cyclic. Theorem (27.19)
More informationConstruction of latin squares of prime order
Construction of latin squares of prime order Theorem. If p is prime, then there exist p 1 MOLS of order p. Construction: The elements in the latin square will be the elements of Z p, the integers modulo
More informationHomework 8 Solutions to Selected Problems
Homework 8 Solutions to Selected Problems June 7, 01 1 Chapter 17, Problem Let f(x D[x] and suppose f(x is reducible in D[x]. That is, there exist polynomials g(x and h(x in D[x] such that g(x and h(x
More informationFactorization in Polynomial Rings
Factorization in Polynomial Rings Throughout these notes, F denotes a field. 1 Long division with remainder We begin with some basic definitions. Definition 1.1. Let f, g F [x]. We say that f divides g,
More informationConjugacy classes of torsion in GL_n(Z)
Electronic Journal of Linear Algebra Volume 30 Volume 30 (2015) Article 32 2015 Conjugacy classes of torsion in GL_n(Z) Qingjie Yang Renmin University of China yangqj@ruceducn Follow this and additional
More informationLinear Algebra (Review) Volker Tresp 2017
Linear Algebra (Review) Volker Tresp 2017 1 Vectors k is a scalar (a number) c is a column vector. Thus in two dimensions, c = ( c1 c 2 ) (Advanced: More precisely, a vector is defined in a vector space.
More informationRINGS: SUMMARY OF MATERIAL
RINGS: SUMMARY OF MATERIAL BRIAN OSSERMAN This is a summary of terms used and main results proved in the subject of rings, from Chapters 11-13 of Artin. Definitions not included here may be considered
More informationChapter 5: The Integers
c Dr Oksana Shatalov, Fall 2014 1 Chapter 5: The Integers 5.1: Axioms and Basic Properties Operations on the set of integers, Z: addition and multiplication with the following properties: A1. Addition
More informationFinite Fields and Error-Correcting Codes
Lecture Notes in Mathematics Finite Fields and Error-Correcting Codes Karl-Gustav Andersson (Lund University) (version 1.013-16 September 2015) Translated from Swedish by Sigmundur Gudmundsson Contents
More informationChapter 8. P-adic numbers. 8.1 Absolute values
Chapter 8 P-adic numbers Literature: N. Koblitz, p-adic Numbers, p-adic Analysis, and Zeta-Functions, 2nd edition, Graduate Texts in Mathematics 58, Springer Verlag 1984, corrected 2nd printing 1996, Chap.
More informationContents Regularization and Matrix Computation in Numerical Polynomial Algebra
Contents 1 Regularization and Matrix Computation in Numerical Polynomial Algebra.................................................... 1 Zhonggang Zeng 1.1 Notation and preliminaries..................................
More informationACI-matrices all of whose completions have the same rank
ACI-matrices all of whose completions have the same rank Zejun Huang, Xingzhi Zhan Department of Mathematics East China Normal University Shanghai 200241, China Abstract We characterize the ACI-matrices
More information4 Powers of an Element; Cyclic Groups
4 Powers of an Element; Cyclic Groups Notation When considering an abstract group (G, ), we will often simplify notation as follows x y will be expressed as xy (x y) z will be expressed as xyz x (y z)
More informationa = a i 2 i a = All such series are automatically convergent with respect to the standard norm, but note that this representation is not unique: i<0
p-adic Numbers K. Sutner v0.4 1 Modular Arithmetic rings integral domains integers gcd, extended Euclidean algorithm factorization modular numbers add Lemma 1.1 (Chinese Remainder Theorem) Let a b. Then
More informationThe Fundamental Theorem of Arithmetic
Chapter 1 The Fundamental Theorem of Arithmetic 1.1 Primes Definition 1.1. We say that p N is prime if it has just two factors in N, 1 and p itself. Number theory might be described as the study of the
More informationIntegration of Rational Functions by Partial Fractions
Title Integration of Rational Functions by MATH 1700 MATH 1700 1 / 11 Readings Readings Readings: Section 7.4 MATH 1700 2 / 11 Rational functions A rational function is one of the form where P and Q are
More informationLinear Algebra in Actuarial Science: Slides to the lecture
Linear Algebra in Actuarial Science: Slides to the lecture Fall Semester 2010/2011 Linear Algebra is a Tool-Box Linear Equation Systems Discretization of differential equations: solving linear equations
More informationAlgebraic function fields
Algebraic function fields 1 Places Definition An algebraic function field F/K of one variable over K is an extension field F K such that F is a finite algebraic extension of K(x) for some element x F which
More information5. Some theorems on continuous functions
5. Some theorems on continuous functions The results of section 3 were largely concerned with continuity of functions at a single point (usually called x 0 ). In this section, we present some consequences
More information8 Further theory of function limits and continuity
8 Further theory of function limits and continuity 8.1 Algebra of limits and sandwich theorem for real-valued function limits The following results give versions of the algebra of limits and sandwich theorem
More informationSection IV.23. Factorizations of Polynomials over a Field
IV.23 Factorizations of Polynomials 1 Section IV.23. Factorizations of Polynomials over a Field Note. Our experience with classical algebra tells us that finding the zeros of a polynomial is equivalent
More informationReview of Some Concepts from Linear Algebra: Part 2
Review of Some Concepts from Linear Algebra: Part 2 Department of Mathematics Boise State University January 16, 2019 Math 566 Linear Algebra Review: Part 2 January 16, 2019 1 / 22 Vector spaces A set
More informationCylindrical Algebraic Decomposition in Coq
Cylindrical Algebraic Decomposition in Coq MAP 2010 - Logroño 13-16 November 2010 Assia Mahboubi INRIA Microsoft Research Joint Centre (France) INRIA Saclay Île-de-France École Polytechnique, Palaiseau
More informationMath 121 Homework 3 Solutions
Math 121 Homework 3 Solutions Problem 13.4 #6. Let K 1 and K 2 be finite extensions of F in the field K, and assume that both are splitting fields over F. (a) Prove that their composite K 1 K 2 is a splitting
More informationThe Seven Dwarfs of Symbolic Computation and the Discovery of Reduced Symbolic Models. Erich Kaltofen North Carolina State University google->kaltofen
The Seven Dwarfs of Symbolic Computation and the Discovery of Reduced Symbolic Models Erich Kaltofen North Carolina State University google->kaltofen Outline 2 The 7 Dwarfs of Symbolic Computation Discovery
More informationLinear Algebra (Review) Volker Tresp 2018
Linear Algebra (Review) Volker Tresp 2018 1 Vectors k, M, N are scalars A one-dimensional array c is a column vector. Thus in two dimensions, ( ) c1 c = c 2 c i is the i-th component of c c T = (c 1, c
More informationFinite Fields. SOLUTIONS Network Coding - Prof. Frank H.P. Fitzek
Finite Fields In practice most finite field applications e.g. cryptography and error correcting codes utilizes a specific type of finite fields, namely the binary extension fields. The following exercises
More informationIntroduction to finite fields
Chapter 7 Introduction to finite fields This chapter provides an introduction to several kinds of abstract algebraic structures, particularly groups, fields, and polynomials. Our primary interest is in
More information1 Introduction. 2 Determining what the J i blocks look like. December 6, 2006
Jordan Canonical Forms December 6, 2006 1 Introduction We know that not every n n matrix A can be diagonalized However, it turns out that we can always put matrices A into something called Jordan Canonical
More informationSolving Approximate GCD of Multivariate Polynomials By Maple/Matlab/C Combination
Solving Approximate GCD of Multivariate Polynomials By Maple/Matlab/C Combination Kai Li Lihong Zhi Matu-Tarow Noda Department of Computer Science Ehime University, Japan {likai,lzhi,noda}@hpc.cs.ehime-u.ac.jp
More informationIntegration of Rational Functions by Partial Fractions
Title Integration of Rational Functions by Partial Fractions MATH 1700 December 6, 2016 MATH 1700 Partial Fractions December 6, 2016 1 / 11 Readings Readings Readings: Section 7.4 MATH 1700 Partial Fractions
More informationAlgebra Exam Topics. Updated August 2017
Algebra Exam Topics Updated August 2017 Starting Fall 2017, the Masters Algebra Exam will have 14 questions. Of these students will answer the first 8 questions from Topics 1, 2, and 3. They then have
More informationMarkov Chains, Stochastic Processes, and Matrix Decompositions
Markov Chains, Stochastic Processes, and Matrix Decompositions 5 May 2014 Outline 1 Markov Chains Outline 1 Markov Chains 2 Introduction Perron-Frobenius Matrix Decompositions and Markov Chains Spectral
More informationChapter 3. Rings. The basic commutative rings in mathematics are the integers Z, the. Examples
Chapter 3 Rings Rings are additive abelian groups with a second operation called multiplication. The connection between the two operations is provided by the distributive law. Assuming the results of Chapter
More informationCoding Theory ( Mathematical Background I)
N.L.Manev, Lectures on Coding Theory (Maths I) p. 1/18 Coding Theory ( Mathematical Background I) Lector: Nikolai L. Manev Institute of Mathematics and Informatics, Sofia, Bulgaria N.L.Manev, Lectures
More informationLesson 2.1: Quadratic Functions
Quadratic Functions: Lesson 2.1: Quadratic Functions Standard form (vertex form) of a quadratic function: Vertex: (h, k) Algebraically: *Use completing the square to convert a quadratic equation into standard
More informationWhen is n a member of a Pythagorean Triple?
Abstract When is n a member of a Pythagorean Triple? Dominic and Alfred Vella ' Theorem is perhaps the best known result in the whole of mathematics and yet many things remain unknown (or perhaps just
More informationForms as sums of powers of lower degree forms
Bruce Reznick University of Illinois at Urbana-Champaign SIAM Conference on Applied Algebraic Geometry Algebraic Geometry of Tensor Decompositions Fort Collins, Colorado August 2, 2013 Let H d (C n ) denote
More informationProblem 1A. Use residues to compute. dx x
Problem 1A. A non-empty metric space X is said to be connected if it is not the union of two non-empty disjoint open subsets, and is said to be path-connected if for every two points a, b there is a continuous
More informationChap 3. Linear Algebra
Chap 3. Linear Algebra Outlines 1. Introduction 2. Basis, Representation, and Orthonormalization 3. Linear Algebraic Equations 4. Similarity Transformation 5. Diagonal Form and Jordan Form 6. Functions
More informationSolutions or answers to Final exam in Error Control Coding, October 24, G eqv = ( 1+D, 1+D + D 2)
Solutions or answers to Final exam in Error Control Coding, October, Solution to Problem a) G(D) = ( +D, +D + D ) b) The rate R =/ and ν i = ν = m =. c) Yes, since gcd ( +D, +D + D ) =+D + D D j. d) An
More informationMath 473: Practice Problems for Test 1, Fall 2011, SOLUTIONS
Math 473: Practice Problems for Test 1, Fall 011, SOLUTIONS Show your work: 1. (a) Compute the Taylor polynomials P n (x) for f(x) = sin x and x 0 = 0. Solution: Compute f(x) = sin x, f (x) = cos x, f
More information1 Positive definiteness and semidefiniteness
Positive definiteness and semidefiniteness Zdeněk Dvořák May 9, 205 For integers a, b, and c, let D(a, b, c) be the diagonal matrix with + for i =,..., a, D i,i = for i = a +,..., a + b,. 0 for i = a +
More information