Towards High Performance Multivariate Factorization. Michael Monagan. This is joint work with Baris Tuncer.
|
|
- Karin Holt
- 6 years ago
- Views:
Transcription
1 Towards High Performance Multivariate Factorization Michael Monagan Center for Experimental and Constructive Mathematics Simon Fraser University British Columbia This is joint work with Baris Tuncer.
2 The Multivariate Diophantine Problem (MDP) Given A, B, C Z p [x 1, x 2,..., x n ] with gcd(a, B) = 1 solve σa + τb = C for σ, τ Z p [x 1,..., x n ] s.t. deg x1 σ < deg x1 B. Outline Context: factoring a Z[x 1,..., x n ] using Hensel lifting Wang s solution Kaltofen s sparse solution v. our sparse solution Eliminating multi-precision arithmetic Eliminating the error computation
3 Multivariate polynomial factorization Input a = fg Z[x 1,..., x n ] primitive and square-free. Output f, g. 1 Pick x 1 and factor LC(a) = h e1 1 he in Z[x 2,..., x n ]. 2 Pick α = α 2,..., α n Z such that LC(a)(α) 0 and gcd(a(α), a x 1 (α)) = 1 and h i (α) Z have distinct prime divisors. 3 Factor a(x 1, α) = f (x 1, α) g(x 1, α) over Z. 4 Distribute the factors of LC(a) on f (x 1, α) and g(x 1, α). 5 Pick a prime p > α i and L N s.t. p L > 2 f where f a and gcd(f (x 1, α), g(x 1, α)) = 1 in Z p [x 1 ]. Hensel lift x 2 then x 3 then... then x n doing arithmetic mod p L. P. Wang. An improved Multivariate Polynomial Factoring Algorithm. Mathematics of Computation, 32: AMS (1978) Algorithms for Computer Algebra, K.O. Geddes, S.R. Czapor, G. Labahn. Kluwer, 1992.
4 Wang s Multivariate Hensel Lifting (MHL) Input a Z[x 1,..., x n], α = (α 2,... α n), (f 0, g 0) Z[x 1] s.t. (i) a(x 1, α) f 0g 0 = 0 and (ii) gcd(f 0, g 0) = 1 in Z p[x 1]. Output f, g satifying a = fg or FAIL. 1 If n = 1 output (f 0, g 0). 2 Lift x 2,..., x n 1: (f 0, g 0) := MHL( a(x n = α n), α, (f 0, g 0) ) Idea: f = df k=0 f k(x n α n) k and g = dg k=0 g k(x n α n) k.
5 Wang s Multivariate Hensel Lifting (MHL) Input a Z[x 1,..., x n], α = (α 2,... α n), (f 0, g 0) Z[x 1] s.t. (i) a(x 1, α) f 0g 0 = 0 and (ii) gcd(f 0, g 0) = 1 in Z p[x 1]. Output f, g satifying a = fg or FAIL. 1 If n = 1 output (f 0, g 0). 2 Lift x 2,..., x n 1: (f 0, g 0) := MHL( a(x n = α n), α, (f 0, g 0) ) Idea: f = df k=0 f k(x n α n) k and g = dg k=0 g k(x n α n) k. 3 Initialize: (f, g) := (f 0, g 0); error := a fg. 4 For k = 1,..., deg x1 a while error 0 do T k := Taylor coeff( error, (x n α n) k ) If T k 0 then Solve f k g 0 + g k f 0 = T k for f k, g k Z p[x 1,..., x n 1]. Set f := f + f k (x n α n) k and g := g + g k (x n α n) k Set error := a fg 5 If error = 0 output (f, g) else output FAIL.
6 Wang s Multivariate Diophantine Lifting (MDP) Input A, B, C Z p[x 1,..., x n], α = α 2,... α n ) Output σ, τ Z p[x 1,..., x n] satisfying σa + τb = C 1 If n = 1 solve σa + τb = C using the Euclidean algorithm in Z p[x 1]. Let σ = df k=0 σ k(x n α n) k and τ = dg k=0 τ k(x n α n) k. 2 (σ 0, τ 0) := MultiDioLift( A(x n = α n), B(x n α n), C(x n = α n), α )
7 Wang s Multivariate Diophantine Lifting (MDP) Input A, B, C Z p[x 1,..., x n], α = α 2,... α n ) Output σ, τ Z p[x 1,..., x n] satisfying σa + τb = C 1 If n = 1 solve σa + τb = C using the Euclidean algorithm in Z p[x 1]. Let σ = df k=0 σ k(x n α n) k and τ = dg k=0 τ k(x n α n) k. 2 (σ 0, τ 0) := MultiDioLift( A(x n = α n), B(x n α n), C(x n = α n), α ) 3 Initialize: (σ, τ) := (σ 0, τ 0) error := C σa τb 4 For k = 1, 2,... while error 0 do T k := Taylor coeff( error, (x n α n) k ) If T k 0 then σ k, τ k := MultiDioLift( σ 0, τ 0, T k, α ) σ, τ := σ + σ k (x n α n) k, τ + τ k (x n α n) k ) error := error σ k (x n α n) k A τ k (x n α n) k B 5 output (σ, τ).
8 Wang s Multivariate Diophantine Lifting (MDP) Input A, B, C Z p[x 1,..., x n], α = α 2,... α n ) Output σ, τ Z p[x 1,..., x n] satisfying σa + τb = C 1 If n = 1 solve σa + τb = C using the Euclidean algorithm in Z p[x 1]. Let σ = df k=0 σ k(x n α n) k and τ = dg k=0 τ k(x n α n) k. 2 (σ 0, τ 0) := MultiDioLift( A(x n = α n), B(x n α n), C(x n = α n), α ) 3 Initialize: (σ, τ) := (σ 0, τ 0) error := C σa τb 4 For k = 1, 2,... while error 0 do T k := Taylor coeff( error, (x n α n) k ) If T k 0 then σ k, τ k := MultiDioLift( σ 0, τ 0, T k, α ) σ, τ := σ + σ k (x n α n) k, τ + τ k (x n α n) k ) error := error σ k (x n α n) k A τ k (x n α n) k B 5 output (σ, τ). Let M(n) count calls to the Euclidean algorithm and d max(deg xi f, deg xi g). Then M(1) = 1, M(n) dm(n 1) = M(n) d n 1.
9 Recall that f = df i=0 f i(x n α n) i and g = dg i=0 g i(x n α n) i. Solve the MDP f k g 0 + g k f 0 = T k for f k, g k Z p[x 1,..., x n 1]. Interpolate x 2,..., x n 1 from images in Z p[x 1] using sparse interpolation? We tried Zippel s variable at a time sparse interpolation from 1979.
10 Recall that f = df i=0 f i(x n α n) i and g = dg i=0 g i(x n α n) i. Solve the MDP f k g 0 + g k f 0 = T k for f k, g k Z p[x 1,..., x n 1]. Interpolate x 2,..., x n 1 from images in Z p[x 1] using sparse interpolation? We tried Zippel s variable at a time sparse interpolation from Theorem [MT] The Strong Sparse Hensel Assumption If α n is random then, with high probability supp(f 0) supp(f 1) supp(f df ) and supp(g 0) supp(g 1) supp(g dg ). Example (n = 8, #f = 10, 000, df = 16, α = 3) f i = 9877, 7043, 4932, 3374, 2310, 1545, 1001, 654, 418, 245, 141, 81, 34, 13, 5, 2, 1.
11 Recall that f = df i=0 f i(x n α n) i and g = dg i=0 g i(x n α n) i. Solve the MDP f k g 0 + g k f 0 = T k for f k, g k Z p[x 1,..., x n 1]. Interpolate x 2,..., x n 1 from images in Z p[x 1] using sparse interpolation? We tried Zippel s variable at a time sparse interpolation from Theorem [MT] The Strong Sparse Hensel Assumption If α n is random then, with high probability supp(f 0) supp(f 1) supp(f df ) and supp(g 0) supp(g 1) supp(g dg ). Example (n = 8, #f = 10, 000, df = 16, α = 3) f i = 9877, 7043, 4932, 3374, 2310, 1545, 1001, 654, 418, 245, 141, 81, 34, 13, 5, 2, 1. Kaltofen s Sparse Hensel Lifting. Let supp(f 0) = {M 1,..., M r ) and supp(g 0) = {N 1,..., N s} in (x 1,..., x n 1). Equate the coefficients c i of ( r ) ( s ) c i M i g 0 + c r+i N i f 0 = T k in Z p[x 1,..., x n 1] i=1 in monomials in (x 1,..., x n 1) and solve the (r + s) (r + s) linear system. i=1
12 MTSHL : Monagan and Tuncer s Sparse Hensel Lifting Solve f k g 0 + g k f 0 = T k for f k, g k Z p [x 1,..., x n 1 ]. Let f k 1 = df i=1 a i(x 2,..., x n 1)x1 i and supp(a i ) = {M i1,..., M iri }, and g k 1 = dg i=1 b i(x 2,..., x n 1)x1 i and supp(b i ) = {N i1,..., N isi }. Assume f k = ( df ri ) i=0 j=1 a ijm ij x1 i and g k = ( dg si ) i=0 j=1 b ijn ij x1. i
13 MTSHL : Monagan and Tuncer s Sparse Hensel Lifting Solve f k g 0 + g k f 0 = T k for f k, g k Z p [x 1,..., x n 1 ]. Let f k 1 = df i=1 a i(x 2,..., x n 1)x1 i and supp(a i ) = {M i1,..., M iri }, and g k 1 = dg i=1 b i(x 2,..., x n 1)x1 i and supp(b i ) = {N i1,..., N isi }. Assume f k = ( df ri ) i=0 j=1 a ijm ij x1 i and g k = ( dg si ) i=0 j=1 b ijn ij x1. i Let t = max(#a i, #b i ) #a + #b. Let β = (β 2,..., β n 1) Z n 2 p. Solve σ j g 0(β j ) + τ j f 0(β j ) = T k (β j ) for σ j, τ j Z p[x 1]. for 1 j t. Solve df + dg Vandermonde systems {coeff (f k (β j ), x i n) = coeff (σ j, x i n) for 1 j t} {coeff (g k (β j ), x i n) = coeff (τ, x i n) for 1 j t} Cost of MTSHL: number of arithmetic operations in Z p where d = max deg xi a (n 1)d O( #f #g + #a +(#a + #f + #g) t }{{}}{{}}{{} error Taylor eval + d 2 t }{{} UniDio + }{{} dt 2 ). Vandermonde
14 MTSHL with division f k g 0 + g k f 0 = T k = g k = (T k f k g 0 )/f 0. Let f k 1 = df i=1 a i(x 2,..., x n 1)x1 i and supp(a i ) = {M i1,..., M iri }, Assume f k = df i=0 ( ri j=1 a ijm ij ) x i 1. Let t = max(#a i, #b i ) and β = (β 2,..., β n 1) Z n 2 p. Solve σ j g 0(β j ) + τ j f 0(β j ) = T k (β j ) for σ j, τ j Z p[x 1]. for 1 j t. Solve df +dg Vandermonde systems obtained by equating coefficients of x 1 f k (β j ) = σ j and g k (β j ) = τ j for 1 j t Finally compute g k := (T k f k g 0)/f 0 in Z p[x 1,..., x n 1]. Cost of MTSHL: number of arithmetic operations in Z p where d = max deg xi a (n 1)d O( #f #g + #a +(#a + #f +#g) t }{{}}{{}}{{} error Taylor eval + d 2 t }{{} UniDio + }{{} dt 2 + #f #g ). }{{} Vandermonde g k
15 n/d/t Wang (MDP) Kaltofen (MDP) MTSHL (MDP) 4/35/ (11.95) 1.75 (1.18) 1.51 (0.24) 5/35/ (86.28) 3.75 (2.57) 1.16 (0.36) 7/35/ (797.0) 5.04 (4.08) 1.58 (0.59) 9/35/ (4449.4) 8.13 (6.22) 2.94 (0.56) 4/35/ (26.48) (635.1) (0.82) 5/35/ (402.5) (1899.6) 26.0 (4.86) 7/35/ (3842.2) (2305.5) 43.1 (6.84) 9/35/500 > (3805.9) 79.6 (9.71) > M := product( x i, i = 1..n ) ; > f := x d 1 + M randpoly( [x 1,..., x n], terms = T, degree = d ) ;
16 Chapter 6 Newton s Iteration and the Hensel Construction Run the entire Hensel lifting mod p L where p L > 2 max( a, f, g ). Why? To avoid the expression swell when solving σg 0(α) + τf 0(α) = T k (α) in Z[x 1] in MDP. Lemma [Gelfond] Let a Z[x 1,..., x n] with d i = deg xi a. If f a then f e d 1+d d n a where e =
17 Let p be a prime and let f = hf k=0 f kp k and g = hg k=0 g kp k. Example f = 2x 1 + (5 + 0 p + 2p 2 )x 2 + (7 + p)x 3 = (2x 1 + 5x 2 + 7x 3) + 3x 3 p + 2x }{{}}{{} 2 p }{{} f 0 f 1 f 2 2. supp(f 0) = {x 1, x 2, x 3} supp(f 1) = {x 3} supp(f 2) = {x 2}.
18 Let p be a prime and let f = hf k=0 f kp k and g = hg k=0 g kp k. Example f = 2x 1 + (5 + 0 p + 2p 2 )x 2 + (7 + p)x 3 = (2x 1 + 5x 2 + 7x 3) + 3x 3 p + 2x }{{}}{{} 2 p }{{} f 0 f 1 f 2 2. supp(f 0) = {x 1, x 2, x 3} supp(f 1) = {x 3} supp(f 2) = {x 2}. Theorem If p is chosen at random from [2 b 1, 2 b ] for b sufficiently large then supp(f 0) supp(f 1) supp(f 2)... supp(f hf ) and supp(g 0) supp(g 1) supp(g 2)... supp(g hg ) whp. Idea: Run the entire Hensel lifting modulo a machine prime p to obtain f 0, g 0 Z p[x 1,..., x n] satisfying a f 0g 0 mod p = 0. Now do a p-adic lift on f 0, g 0 to get f, g.
19 Multivariate p-adic Lifting Input a Z[x 1,..., x n], (f 0, g 0) Z p[x 1,..., x n] satisfying a f 0g 0 mod p = 0 prime p and a lifting bound L satisfying p L > f where f a. Output f, g Z[x 1,..., x n] satifying a = fg or FAIL. 1 Initialize: (f, g) := (f 0, g 0); error := a fg; 2 For k = 1,..., L while error 0 do T k := ( error p k ) mod p. If T k 0 then Solve the MDP f k g 0 + g k f 0 = T k for f k, g k Z p[x 1,..., x n]. Set f k := mods(σ, p) and g k := mods(τ, p). Set f := f + f k p k ; g := g + g k p k ; error := a fg. 3 If error = 0 output (f, g) else output FAIL. The MDP is solved in Z p[x 1,..., x n] where p is a machine prime!
20 n/d/t m L MTSHL mod p L MTSHL mod p + lift 5/10/ (5.10) /10/ (6.10) /10/ (7.596) /10/ (12.82) /10/ (15.18) /10/ (17.22) Table: CPU timings (seconds) for p adic lifting for p = > f := x d 1 + randpoly( [x 1,..., x n], terms = T, degree = d, coeffs = rand(p m ) ) ;
21 Eliminating the error multiplication a fg Let f (k) = k 1 i=0 f i(x n α n) i and g (k) = k 1 i=0 g i(x n α n) i. Main Idea: error(β j ) = a(β j ) f (k) (β j ) g (k) (β j ) in Z p[x 1, x n]. Compute a(β j ) for 1 j t before main loop. Compute f k 1 (β j ) and g k 1 (β j ) for 1 j t. Let f (k) (β j ) = k 1 i=0 a i(x n α n) i and g (k) (β j ) = k 1 i=0 b i(x n α n) i. Then T k (β j ) = coeff (a(β j ), (x n α) k ) k 1 i=1 }{{} b i(x 1)c k i (x 1). }{{} O(d 2 ) (n 1)d O( #f #g + #a }{{} + }{{}}{{} d 2 t + (#a/d + #f + #g) t }{{} error Taylor Taylor eval We eliminated the the error computation in Z p[x 1,..., x n 1]! O(d 2 ) + d 2 t }{{} UniDio + }{{} dt 2 ). Vandermonde
22 Factoring the determinants of Cyclic and Toeplitz matrices Let C n denote the n n cyclic matrix and let T n denote the n n symmetric Toeplitz matrix below. x 1 x 2... x n 1 x n x 1 x 2 x n 1 x n x n x 1... x n 2 x n 1 x 2 x 1 x n 2 x n 1 C n =..... T n = x 3 x 4... x 1 x 2 x n 1 x n 2 x 1 x 2 x 2 x 3... x n x 1 x n x n 1 x 2 x 1 The determinants of C n and T n are homogeneous polynomials in n variables. Example det T 4 = ( x 1 2 x 1x 2 x 1x 4 x x 2x 3 + x 2x 4 x 3 2 ) ( x1 2 + x 1x 2 + x 1x 4 x x 2x 3 + x 2x 4 x 3 2 ) det C 4 = (x 4 + x 3 + x 1 + x 2) (x 4 x 3 x 1 + x 2) ( x x 1x 3 + x x 2x 4 + x x 4 2 )
23 n #d n #factors Maple (MDP) MTSHL Magma Singular , % , % , % , % , % , % , % , % , % Table: Factorization timings (seconds) for det T n evaluated at x n = 1 Notes: Intel Core i GHz, 16 gigs RAM. Maple 17: kernelopts(numcpus=1), Magma , Singular 3 1 6,
24 n #d n #fmax Maple (MDP) MTSHL Magma Singular % % % % % % % ?? % > % Table: Factorization timings (seconds) for det C n evaluated at x n = 1 Notes:?? = I cannot compute det(c n) nor read in det(c n) nor it s factors.
Towards High Performance Multivariate Factorization. Michael Monagan. This is joint work with Baris Tuncer.
Towards High Performance Multivariate Factorization Michael Monagan Center for Experimental and Constructive Mathematics Simon Fraser University British Columbia This is joint work with Baris Tuncer. To
More informationComputing with polynomials: Hensel constructions
Course Polynomials: Their Power and How to Use Them, JASS 07 Computing with polynomials: Hensel constructions Lukas Bulwahn March 28, 2007 Abstract To solve GCD calculations and factorization of polynomials
More informationComputing Characteristic Polynomials of Matrices of Structured Polynomials
Computing Characteristic Polynomials of Matrices of Structured Polynomials Marshall Law and Michael Monagan Department of Mathematics Simon Fraser University Burnaby, British Columbia, Canada mylaw@sfu.ca
More informationAn Efficient Algorithm for Computing Parametric Multivariate Polynomial GCD
An Efficient Algorithm for Computing Parametric Multivariate Polynomial GCD Dong Lu Key Laboratory of Mathematics Mechanization Academy of Mathematics and Systems Science, CAS Joint work with Deepak Kapur,
More informationAn Approach to Hensel s Lemma
Irish Math. Soc. Bulletin 47 (2001), 15 21 15 An Approach to Hensel s Lemma gary mcguire Abstract. Hensel s Lemma is an important tool in many ways. One application is in factoring polynomials over Z.
More informationSPARSE POLYNOMIAL INTERPOLATION AND THE FAST EUCLIDEAN ALGORITHM
SPARSE POLYNOMIAL INTERPOLATION AND THE FAST EUCLIDEAN ALGORITHM by Soo H. Go B.Math., University of Waterloo, 2006 a Thesis submitted in partial fulfillment of the requirements for the degree of Master
More informationIn-place Arithmetic for Univariate Polynomials over an Algebraic Number Field
In-place Arithmetic for Univariate Polynomials over an Algebraic Number Field Seyed Mohammad Mahdi Javadi 1, Michael Monagan 2 1 School of Computing Science, Simon Fraser University, Burnaby, B.C., V5A
More informationA fast parallel sparse polynomial GCD algorithm.
A fast parallel sparse polynomial GCD algorithm. Jiaxiong Hu and Michael Monagan Department of Mathematics, Simon Fraser University Abstract We present a parallel GCD algorithm for sparse multivariate
More informationOn Factorization of Multivariate Polynomials over Algebraic Number and Function Fields
On Factorization of Multivariate Polynomials over Algebraic Number and Function Fields ABSTRACT Seyed Mohammad Mahdi Javadi School of Computing Science Simon Fraser University Burnaby, B.C. Canada. sjavadi@cecm.sfu.ca.
More informationExact Arithmetic on a Computer
Exact Arithmetic on a Computer Symbolic Computation and Computer Algebra William J. Turner Department of Mathematics & Computer Science Wabash College Crawfordsville, IN 47933 Tuesday 21 September 2010
More informationWhat s the best data structure for multivariate polynomials in a world of 64 bit multicore computers?
What s the best data structure for multivariate polynomials in a world of 64 bit multicore computers? Michael Monagan Center for Experimental and Constructive Mathematics Simon Fraser University British
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 informationResolving zero-divisors using Hensel lifting
Resolving zero-divisors using Hensel lifting John Kluesner and Michael Monagan Department of Mathematics, Simon Fraser University Burnaby, British Columbia, V5A-1S6, Canada jkluesne@sfu.ca mmonagan@sfu.ca
More informationSparse Polynomial Multiplication and Division in Maple 14
Sparse Polynomial Multiplication and Division in Maple 4 Michael Monagan and Roman Pearce Department of Mathematics, Simon Fraser University Burnaby B.C. V5A S6, Canada October 5, 9 Abstract We report
More informationProbabilistic Algorithms for Computing Resultants
Probabilistic Algorithms for Computing Resultants Michael Monagan Centre for Experimental and Constructive Mathematics Simon Fraser University, Burnaby, B.C. V5A 1S6, CANADA mmonagan@cecm.sfu.ca ABSTRACT
More informationAlgorithms for the Non-monic Case of the Sparse Modular GCD Algorithm
Algorithms for the Non-monic Case of the Sparse Modular GCD Algorithm Jennifer de Kleine Department of Mathematics Simon Fraser University Burnaby, B.C. Canada. dekleine@cecm.sfu.ca. Michael Monagan Department
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 informationPolynomial multiplication and division using heap.
Polynomial multiplication and division using heap. Michael Monagan and Roman Pearce Department of Mathematics, Simon Fraser University. Abstract We report on new code for sparse multivariate polynomial
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 informationParallel Sparse Polynomial Interpolation over Finite Fields
Parallel Sparse Polynomial Interpolation over Finite Fields Seyed Mohammad Mahdi Javadi School of Computing Science, Simon Fraser University, Burnaby, B.C. Canada Michael Monagan Department of Mathematics,
More informationComputing over Z, Q, K[X]
Computing over Z, Q, K[X] Clément PERNET M2-MIA Calcul Exact Outline Introduction Chinese Remainder Theorem Rational reconstruction Problem Statement Algorithms Applications Dense CRT codes Extension to
More informationLecture Notes. Advanced Discrete Structures COT S
Lecture Notes Advanced Discrete Structures COT 4115.001 S15 2015-01-13 Recap Divisibility Prime Number Theorem Euclid s Lemma Fundamental Theorem of Arithmetic Euclidean Algorithm Basic Notions - Section
More informationSection III.6. Factorization in Polynomial Rings
III.6. Factorization in Polynomial Rings 1 Section III.6. Factorization in Polynomial Rings Note. We push several of the results in Section III.3 (such as divisibility, irreducibility, and unique factorization)
More informationLocal properties of plane algebraic curves
Chapter 7 Local properties of plane algebraic curves Throughout this chapter let K be an algebraically closed field of characteristic zero, and as usual let A (K) be embedded into P (K) by identifying
More informationOptimizing and Parallelizing Brown s Modular GCD Algorithm
Optimizing and Parallelizing Brown s Modular GCD Algorithm Matthew Gibson, Michael Monagan October 7, 2014 1 Introduction Consider the multivariate polynomial problem over the integers; that is, Gcd(A,
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 information1. Algebra 1.5. Polynomial Rings
1. ALGEBRA 19 1. Algebra 1.5. Polynomial Rings Lemma 1.5.1 Let R and S be rings with identity element. If R > 1 and S > 1, then R S contains zero divisors. Proof. The two elements (1, 0) and (0, 1) are
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 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 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 informationObjective Type Questions
DISTANCE EDUCATION, UNIVERSITY OF CALICUT NUMBER THEORY AND LINEARALGEBRA Objective Type Questions Shyama M.P. Assistant Professor Department of Mathematics Malabar Christian College, Calicut 7/3/2014
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 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 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 informationLecture Notes on Computer Algebra
Lecture Notes on Computer Algebra Ziming Li Abstract These notes record seven lectures given in the computer algebra course in the fall of 2004. The theory of subresultants is not required for the final
More informationAlgorithms for exact (dense) linear algebra
Algorithms for exact (dense) linear algebra Gilles Villard CNRS, Laboratoire LIP ENS Lyon Montagnac-Montpezat, June 3, 2005 Introduction Problem: Study of complexity estimates for basic problems in exact
More informationCHAPTER 14. Ideals and Factor Rings
CHAPTER 14 Ideals and Factor Rings Ideals Definition (Ideal). A subring A of a ring R is called a (two-sided) ideal of R if for every r 2 R and every a 2 A, ra 2 A and ar 2 A. Note. (1) A absorbs elements
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 informationWORKING WITH MULTIVARIATE POLYNOMIALS IN MAPLE
WORKING WITH MULTIVARIATE POLYNOMIALS IN MAPLE JEFFREY B. FARR AND ROMAN PEARCE Abstract. We comment on the implementation of various algorithms in multivariate polynomial theory. Specifically, we describe
More informationFactoring univariate polynomials over the rationals
Factoring univariate polynomials over the rationals Tommy Hofmann TU Kaiserslautern November 21, 2017 Tommy Hofmann Factoring polynomials over the rationals November 21, 2017 1 / 31 Factoring univariate
More informationAn introduction to the algorithmic of p-adic numbers
An introduction to the algorithmic of p-adic numbers David Lubicz 1 1 Universté de Rennes 1, Campus de Beaulieu, 35042 Rennes Cedex, France Outline Introduction 1 Introduction 2 3 4 5 6 7 8 When do we
More informationPolynomial Rings. i=0. i=0. n+m. i=0. k=0
Polynomial Rings 1. Definitions and Basic Properties For convenience, the ring will always be a commutative ring with identity. Basic Properties The polynomial ring R[x] in the indeterminate x with coefficients
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 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 informationCOMP239: Mathematics for Computer Science II. Prof. Chadi Assi EV7.635
COMP239: Mathematics for Computer Science II Prof. Chadi Assi assi@ciise.concordia.ca EV7.635 The Euclidean Algorithm The Euclidean Algorithm Finding the GCD of two numbers using prime factorization is
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 informationArithmetic properties of lacunary sums of binomial coefficients
Arithmetic properties of lacunary sums of binomial coefficients Tamás Mathematics Department Occidental College 29th Journées Arithmétiques JA2015, July 6-10, 2015 Arithmetic properties of lacunary sums
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 informationOutput-sensitive algorithms for sumset and sparse polynomial multiplication
Output-sensitive algorithms for sumset and sparse polynomial multiplication Andrew Arnold Cheriton School of Computer Science University of Waterloo Waterloo, Ontario, Canada Daniel S. Roche Computer Science
More informationECEN 604: Channel Coding for Communications
ECEN 604: Channel Coding for Communications Lecture: Introduction to Cyclic Codes Henry D. Pfister Department of Electrical and Computer Engineering Texas A&M University ECEN 604: Channel Coding for Communications
More informationAlgebraic structures I
MTH5100 Assignment 1-10 Algebraic structures I For handing in on various dates January March 2011 1 FUNCTIONS. Say which of the following rules successfully define functions, giving reasons. For each one
More informationGroups, Rings, and Finite Fields. Andreas Klappenecker. September 12, 2002
Background on Groups, Rings, and Finite Fields Andreas Klappenecker September 12, 2002 A thorough understanding of the Agrawal, Kayal, and Saxena primality test requires some tools from algebra and elementary
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 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 informationChuck Garner, Ph.D. May 25, 2009 / Georgia ARML Practice
Some Chuck, Ph.D. Department of Mathematics Rockdale Magnet School for Science Technology May 25, 2009 / Georgia ARML Practice Outline 1 2 3 4 Outline 1 2 3 4 Warm-Up Problem Problem Find all positive
More informationBasic Algorithms in Number Theory
Basic Algorithms in Number Theory Algorithmic Complexity... 1 Basic Algorithms in Number Theory Francesco Pappalardi Discrete Logs, Modular Square Roots & Euclidean Algorithm. July 20 th 2010 Basic Algorithms
More informationMath 312/ AMS 351 (Fall 17) Sample Questions for Final
Math 312/ AMS 351 (Fall 17) Sample Questions for Final 1. Solve the system of equations 2x 1 mod 3 x 2 mod 7 x 7 mod 8 First note that the inverse of 2 is 2 mod 3. Thus, the first equation becomes (multiply
More informationGRÖBNER BASES AND POLYNOMIAL EQUATIONS. 1. Introduction and preliminaries on Gróbner bases
GRÖBNER BASES AND POLYNOMIAL EQUATIONS J. K. VERMA 1. Introduction and preliminaries on Gróbner bases Let S = k[x 1, x 2,..., x n ] denote a polynomial ring over a field k where x 1, x 2,..., x n are indeterminates.
More informationInverting integer and polynomial matrices. Jo60. Arne Storjohann University of Waterloo
Inverting integer and polynomial matrices Jo60 Arne Storjohann University of Waterloo Integer Matrix Inverse Input: An n n matrix A filled with entries of size d digits. Output: The matrix inverse A 1.
More informationSolving Sparse Rational Linear Systems. Pascal Giorgi. University of Waterloo (Canada) / University of Perpignan (France) joint work with
Solving Sparse Rational Linear Systems Pascal Giorgi University of Waterloo (Canada) / University of Perpignan (France) joint work with A. Storjohann, M. Giesbrecht (University of Waterloo), W. Eberly
More informationRational Univariate Reduction via Toric Resultants
Rational Univariate Reduction via Toric Resultants Koji Ouchi 1,2 John Keyser 1 Department of Computer Science, 3112 Texas A&M University, College Station, TX 77843-3112, USA Abstract We describe algorithms
More informationPOLY : A new polynomial data structure for Maple 17 that improves parallel speedup.
: A new polynomial data structure for Maple 7 that improves parallel speedup. Department of Mathematics, Simon Fraser University British Columbia, CANADA Parallel Computer Algebra Applications ACA 22,
More informationReal Solving on Algebraic Systems of Small Dimension
Real Solving on Algebraic Systems of Small Dimension Master s Thesis Presentation Dimitrios I. Diochnos University of Athens March 8, 2007 D. I. Diochnos (Univ. of Athens, µ Q λ ) Real Solving on Bivariate
More informationA Fast Las Vegas Algorithm for Computing the Smith Normal Form of a Polynomial Matrix
A Fast Las Vegas Algorithm for Computing the Smith Normal Form of a Polynomial Matrix Arne Storjohann and George Labahn Department of Computer Science University of Waterloo Waterloo, Ontario, Canada N2L
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 informationAdaptive decoding for dense and sparse evaluation/interpolation codes
Adaptive decoding for dense and sparse evaluation/interpolation codes Clément PERNET INRIA/LIG-MOAIS, Grenoble Université joint work with M. Comer, E. Kaltofen, J-L. Roch, T. Roche Séminaire Calcul formel
More informationCover Page. The handle holds various files of this Leiden University dissertation
Cover Page The handle http://hdl.handle.net/1887/32076 holds various files of this Leiden University dissertation Author: Junjiang Liu Title: On p-adic decomposable form inequalities Issue Date: 2015-03-05
More information2a 2 4ac), provided there is an element r in our
MTH 310002 Test II Review Spring 2012 Absractions versus examples The purpose of abstraction is to reduce ideas to their essentials, uncluttered by the details of a specific situation Our lectures built
More informationIn-place Arithmetic for Polynomials over Zn
In-place Arithmetic for Polynomials over Zn Michael Monagan Institut ffir Wissenschaftliches Rechnen ETH-Zentrum, CH-8092 Zfirich, Switzerland monagan~inf.ethz.ch Abstract. We present space and time efficient
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 informationNew algebraic decoding method for the (41, 21,9) quadratic residue code
New algebraic decoding method for the (41, 21,9) quadratic residue code Mohammed M. Al-Ashker a, Ramez Al.Shorbassi b a Department of Mathematics Islamic University of Gaza, Palestine b Ministry of education,
More informationPolynomials over UFD s
Polynomials over UFD s Let R be a UFD and let K be the field of fractions of R. Our goal is to compare arithmetic in the rings R[x] and K[x]. We introduce the following notion. Definition 1. A non-constant
More informationInteger Multiplication
Integer Multiplication in almost linear time Martin Fürer CSE 588 Department of Computer Science and Engineering Pennsylvania State University 1/24/08 Karatsuba algebraic Split each of the two factors
More informationChapter 4 Finite Fields
Chapter 4 Finite Fields Introduction will now introduce finite fields of increasing importance in cryptography AES, Elliptic Curve, IDEA, Public Key concern operations on numbers what constitutes a number
More informationGiac and GeoGebra: improved Gröbner basis computations
Giac and GeoGebra: improved Gröbner basis computations Z. Kovács, B. Parisse JKU Linz, University of Grenoble I November 25, 2013 Two parts talk 1 GeoGebra (Z. Kovács) 2 (B. Parisse) History of used CAS
More informationMath 120 HW 9 Solutions
Math 120 HW 9 Solutions June 8, 2018 Question 1 Write down a ring homomorphism (no proof required) f from R = Z[ 11] = {a + b 11 a, b Z} to S = Z/35Z. The main difficulty is to find an element x Z/35Z
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 informationQuestion 1: The graphs of y = p(x) are given in following figure, for some Polynomials p(x). Find the number of zeroes of p(x), in each case.
Class X - NCERT Maths EXERCISE NO:.1 Question 1: The graphs of y = p(x) are given in following figure, for some Polynomials p(x). Find the number of zeroes of p(x), in each case. (i) (ii) (iii) (iv) (v)
More information6.S897 Algebra and Computation February 27, Lecture 6
6.S897 Algebra and Computation February 7, 01 Lecture 6 Lecturer: Madhu Sudan Scribe: Mohmammad Bavarian 1 Overview Last lecture we saw how to use FFT to multiply f, g R[x] in nearly linear time. We also
More informationHandout - Algebra Review
Algebraic Geometry Instructor: Mohamed Omar Handout - Algebra Review Sept 9 Math 176 Today will be a thorough review of the algebra prerequisites we will need throughout this course. Get through as much
More informationElliptic Curve Discrete Logarithm Problem
Elliptic Curve Discrete Logarithm Problem Vanessa VITSE Université de Versailles Saint-Quentin, Laboratoire PRISM October 19, 2009 Vanessa VITSE (UVSQ) Elliptic Curve Discrete Logarithm Problem October
More informationDiscrete logarithms: Recent progress (and open problems)
Discrete logarithms: Recent progress (and open problems) CryptoExperts Chaire de Cryptologie de la Fondation de l UPMC LIP6 February 25 th, 2014 Discrete logarithms Given a multiplicative group G with
More informationPolynomials, Ideals, and Gröbner Bases
Polynomials, Ideals, and Gröbner Bases Notes by Bernd Sturmfels for the lecture on April 10, 2018, in the IMPRS Ringvorlesung Introduction to Nonlinear Algebra We fix a field K. Some examples of fields
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 informationReal Solving on Bivariate Systems with Sturm Sequences and SLV Maple TM library
Real Solving on Bivariate Systems with Sturm Sequences and SLV Maple TM library Dimitris Diochnos University of Illinois at Chicago Dept. of Mathematics, Statistics, and Computer Science September 27,
More informationMATH 433 Applied Algebra Lecture 4: Modular arithmetic (continued). Linear congruences.
MATH 433 Applied Algebra Lecture 4: Modular arithmetic (continued). Linear congruences. Congruences Let n be a postive integer. The integers a and b are called congruent modulo n if they have the same
More informationAn 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 informationbe any ring homomorphism and let s S be any element of S. Then there is a unique ring homomorphism
21. Polynomial rings Let us now turn out attention to determining the prime elements of a polynomial ring, where the coefficient ring is a field. We already know that such a polynomial ring is a UFD. Therefore
More informationPolynomial factorization and curve decomposition algorithms
UNIVERSITÁ DEGLI STUDI DI TORINO FACOLTÁ DI SCIENZE MATEMATICHE FISICHE E NATURALI UNIVERSITÉ DE NICE SOPHIA ANTIPOLIS FACULTÉ DES SCIENCES DE NICE THÈSE DE DOCTORAT pour l obtention du titre de Dottore
More informationDiscrete Mathematics and Probability Theory Spring 2016 Rao and Walrand Discussion 6B Solution
CS 70 Discrete Mathematics and Probability Theory Spring 016 Rao and Walrand Discussion 6B Solution 1. GCD of Polynomials Let A(x) and B(x) be polynomials (with coefficients in R or GF(m)). We say that
More informationLinear Algebra III Lecture 11
Linear Algebra III Lecture 11 Xi Chen 1 1 University of Alberta February 13, 2015 Outline Minimal Polynomial 1 Minimal Polynomial Minimal Polynomial The minimal polynomial f (x) of a square matrix A is
More informationPolynomial Review Problems
Polynomial Review Problems 1. Find polynomial function formulas that could fit each of these graphs. Remember that you will need to determine the value of the leading coefficient. The point (0,-3) is on
More informationA Modular Algorithm for Computing Greatest Common Right Divisors of Ore Polynomials 1)
MM Research Preprints, 69 88 No. 15, April 1997. Beijing 69 A Modular Algorithm for Computing Greatest Common Right Divisors of Ore Polynomials 1) Ziming Li 2) and István Nemes 3) Abstract. This paper
More informationMath 2070BC Term 2 Weeks 1 13 Lecture Notes
Math 2070BC 2017 18 Term 2 Weeks 1 13 Lecture Notes Keywords: group operation multiplication associative identity element inverse commutative abelian group Special Linear Group order infinite order cyclic
More informationAlgebra Review 2. 1 Fields. A field is an extension of the concept of a group.
Algebra Review 2 1 Fields A field is an extension of the concept of a group. Definition 1. A field (F, +,, 0 F, 1 F ) is a set F together with two binary operations (+, ) on F such that the following conditions
More informationax b mod m. has a solution if and only if d b. In this case, there is one solution, call it x 0, to the equation and there are d solutions x m d
10. Linear congruences In general we are going to be interested in the problem of solving polynomial equations modulo an integer m. Following Gauss, we can work in the ring Z m and find all solutions to
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 informationPrimary Decomposition
Primary Decomposition p. Primary Decomposition Gerhard Pfister pfister@mathematik.uni-kl.de Departement of Mathematics University of Kaiserslautern Primary Decomposition p. Primary Decomposition:References
More informationModern Computer Algebra
Modern Computer Algebra JOACHIM VON ZUR GATHEN and JURGEN GERHARD Universitat Paderborn CAMBRIDGE UNIVERSITY PRESS Contents Introduction 1 1 Cyclohexane, cryptography, codes, and computer algebra 9 1.1
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 information