Ruppert matrix as subresultant mapping

Size: px
Start display at page:

Download "Ruppert matrix as subresultant mapping"

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)

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 information

Math 4310 Solutions to homework 7 Due 10/27/16

Math 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 information

g(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.

g(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 information

Polynomials. Chapter 4

Polynomials. 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 information

CHAPTER I. Rings. Definition A ring R is a set with two binary operations, addition + and

CHAPTER 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 information

IRREDUCIBILITY OF POLYNOMIALS MODULO p VIA NEWTON POLYTOPES. 1. Introduction

IRREDUCIBILITY 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 information

Computing Approximate GCD of Univariate Polynomials by Structure Total Least Norm 1)

Computing 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 information

Mathematical Olympiad Training Polynomials

Mathematical 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 information

Public-key Cryptography: Theory and Practice

Public-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 information

Fast 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 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 information

MATH 431 PART 2: POLYNOMIAL RINGS AND FACTORIZATION

MATH 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 information

D-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 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 information

MTH310 EXAM 2 REVIEW

MTH310 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 information

The Sylvester Resultant

The 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 information

Solutions of exercise sheet 11

Solutions 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 information

0 Sets and Induction. Sets

0 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 information

Simultaneous Linear, and Non-linear Congruences

Simultaneous 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 information

Chapter 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: 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 information

Structured Matrix Methods Computing the Greatest Common Divisor of Polynomials

Structured 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 information

The Mathematica Journal Symbolic-Numeric Algebra for Polynomials

The 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 information

Math 4320 Final Exam

Math 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 information

Chapter 4. Greatest common divisors of polynomials. 4.1 Polynomial remainder sequences

Chapter 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 information

CS 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 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 information

x 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

x 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 information

Lattice reduction of polynomial matrices

Lattice 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 =

ϕ : 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 information

NOTES ON DIOPHANTINE APPROXIMATION

NOTES 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 information

CHAPTER 10: POLYNOMIALS (DRAFT)

CHAPTER 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 information

Parity of the Number of Irreducible Factors for Composite Polynomials

Parity 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 information

Lecture Notes Math 371: Algebra (Fall 2006) by Nathanael Leedom Ackerman

Lecture 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 information

Multivariate Statistical Analysis

Multivariate 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 information

Tensor Decompositions for Machine Learning. G. Roeder 1. UBC Machine Learning Reading Group, June University of British Columbia

Tensor 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 information

Lecture 6. Numerical methods. Approximation of functions

Lecture 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 information

Explicit Criterion to Determine the Number of Positive Roots of a Polynomial 1)

Explicit 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 information

Lecture 7: Polynomial rings

Lecture 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 information

Course 2316 Sample Paper 1

Course 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 information

Chinese Remainder Theorem

Chinese 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 information

computing an irreducible decomposition of A B looking for solutions where the Jacobian drops rank

computing 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 information

MATH FINAL EXAM REVIEW HINTS

MATH 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 information

Basic Algorithms in Number Theory

Basic 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 information

Chapter 5. Linear Algebra

Chapter 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 information

Fast Polynomial Multiplication

Fast 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 information

Galois Fields and Hardware Design

Galois 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 information

Chapter 4. Remember: F will always stand for a field.

Chapter 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 information

arxiv: v3 [math.ac] 11 May 2016

arxiv: 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 information

Finite Fields. Mike Reiter

Finite 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 information

Math 547, Exam 2 Information.

Math 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 information

Groebner 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 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 information

Sylvester Matrix and GCD for Several Univariate Polynomials

Sylvester 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 information

Lecture 02 Linear Algebra Basics

Lecture 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 information

Class Notes; Week 7, 2/26/2016

Class 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 information

Further linear algebra. Chapter II. Polynomials.

Further 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 information

Computer Algebra: General Principles

Computer 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 information

Linear and Bilinear Algebra (2WF04) Jan Draisma

Linear 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 information

Outline. MSRI-UP 2009 Coding Theory Seminar, Week 2. The definition. Link to polynomials

Outline. 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 information

Example Linear Algebra Competency Test

Example 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 information

12x + 18y = 30? ax + by = m

12x + 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 information

Section 33 Finite fields

Section 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 information

Construction of latin squares of prime order

Construction 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 information

Homework 8 Solutions to Selected Problems

Homework 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 information

Factorization in Polynomial Rings

Factorization 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 information

Conjugacy classes of torsion in GL_n(Z)

Conjugacy 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 information

Linear Algebra (Review) Volker Tresp 2017

Linear 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 information

RINGS: SUMMARY OF MATERIAL

RINGS: 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 information

Chapter 5: The Integers

Chapter 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 information

Finite Fields and Error-Correcting Codes

Finite 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 information

Chapter 8. P-adic numbers. 8.1 Absolute values

Chapter 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 information

Contents Regularization and Matrix Computation in Numerical Polynomial Algebra

Contents 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 information

ACI-matrices all of whose completions have the same rank

ACI-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 information

4 Powers of an Element; Cyclic Groups

4 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 information

a = 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

a = 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 information

The Fundamental Theorem of Arithmetic

The 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 information

Integration of Rational Functions by Partial Fractions

Integration 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 information

Linear Algebra in Actuarial Science: Slides to the lecture

Linear 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 information

Algebraic function fields

Algebraic 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 information

5. Some theorems on continuous functions

5. 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 information

8 Further theory of function limits and continuity

8 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 information

Section IV.23. Factorizations of Polynomials over a Field

Section 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 information

Review of Some Concepts from Linear Algebra: Part 2

Review 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 information

Cylindrical Algebraic Decomposition in Coq

Cylindrical 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 information

Math 121 Homework 3 Solutions

Math 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 information

The 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 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 information

Linear Algebra (Review) Volker Tresp 2018

Linear 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 information

Finite Fields. SOLUTIONS Network Coding - Prof. Frank H.P. Fitzek

Finite 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 information

Introduction to finite fields

Introduction 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 information

1 Introduction. 2 Determining what the J i blocks look like. December 6, 2006

1 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 information

Solving Approximate GCD of Multivariate Polynomials By Maple/Matlab/C Combination

Solving 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 information

Integration of Rational Functions by Partial Fractions

Integration 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 information

Algebra Exam Topics. Updated August 2017

Algebra 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 information

Markov Chains, Stochastic Processes, and Matrix Decompositions

Markov 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 information

Chapter 3. Rings. The basic commutative rings in mathematics are the integers Z, the. Examples

Chapter 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 information

Coding Theory ( Mathematical Background I)

Coding 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 information

Lesson 2.1: Quadratic Functions

Lesson 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 information

When is n a member of a Pythagorean Triple?

When 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 information

Forms as sums of powers of lower degree forms

Forms 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 information

Problem 1A. Use residues to compute. dx x

Problem 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 information

Chap 3. Linear Algebra

Chap 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 information

Solutions 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 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 information

Math 473: Practice Problems for Test 1, Fall 2011, SOLUTIONS

Math 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 information

1 Positive definiteness and semidefiniteness

1 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