Short multipliers for the extended gcd problem
|
|
- Jeffrey Warner
- 5 years ago
- Views:
Transcription
1 Short multipliers for the extended gcd problem Keith Matthews Abstract For given non zero integers s 1,, s m, the problem of finding integers a 1,, a m satisfying s = gcd (s 1,, s m ) = a 1 s a m s m, with a a m minimal, is thought to be computationally hard In this paper, we present an algorithm which takes as its starting point the recent LLL based algorithm of Havas, Majewski and Matthews and which often finds a shorter vector (a 1,, a m ) 1 Introduction Let s 1,, s m be integers and s = gcd (s 1,, s m ) In a recent paper [Havas, Majewski, Matthews 1998], the author and his collaborators used variants of the LLL algorithm to find multiplier vectors (a 1,, a m ) of small Euclidean length X = (a a m) 1/ such that s = a 1 s a m s m In each case a unimodular m m matrix P is produced such that P [s 1,, s m ] t = [0,, 0, s] t Rows p 1,, p of P constitute a basis of short vectors for the (m 1) dimensional lattice Λ formed by the vectors X = (a 1,, a m ) with a 1,, a m Z, satisfying a 1 s a m s m = 0 In particular, every such X can be expressed uniquely as an integer linear combination X = z 1 p z p In addition, p m, the last row of P, is a short multiplier vector and the general multiplier vector p is given by p = p m + x 1 p x p, where x 1,, x Z The matrix P has further properties: If the Gram Schmidt basis corresponding to rows p 1,, p m is denoted by p 1,, p m, where p 1 = p 1, k 1 p k = p k µ kj p j, µ kj = p k p j p j, (1) p j j=1 1
2 then (a) µ kj 1/ for 1 j < k m, (b) p k p k (α µ k k 1)p k 1 p k 1 for k m 1 () (Here 1/ < α 1) In what follows we assume α = 1, so that () becomes From the equations a multiplier vector p may be written as where p k p k (1 µ k k 1)p k 1 p k 1 () k 1 p 1 = p 1, p k = p k + µ kj p j, () p = p m + k=1 y k = x k + j=1 x k p k = p m + i=k+1 The orthogonality of p 1,, p m then implies where p = p m + k=1 k=1 y k p k, (5) µ ik x i + µ mk (6) x k p k = p m + Q(x 1,, x ) = B (x + µ m, ) and B k = p k for k = 1,, m k=1 y k p k = B m + Q(x 1,, x ), (7) + B m (x m + µ,m x + µ m,m ) + B 1 (x 1 + µ,1 x + + µ,1 x + µ m,1 ) (8)
3 The equation P [s 1,, s m ] t = [0,, 0, s] t and the fact that [s 1,, s m ] is orthogonal to each of p 1,, p together imply Hence p m = s (s s s 1,, s m ) (9) m B m = p m = (10) s s m From equations () and the fact that the determinant of the orthogonal matrix whose rows are p 1,, p m, is equal to det P = ±1, we also have Consequently B 1 B m = det(p i p j) = (det P ) = 1 s = (det Λ) = B 1 B = (s s m)/s (11) Also () becomes for k m 1 B k (1 µ k k 1)B k 1, (1) The motivation for the algorithm Before we give our algorithm for generating possibly shorter multipliers than p m, we give some background The minimum value 0 of the quadratic expression in equation (8) occurs at the point (ρ 1,, ρ ) Q, where ρ = µ m, ρ k = ( i=k+1 µ ik ρ i + µ mk ), 1 k < m 1 (1) It is then an easy exercise in determinants to show that for 1 k m 1 ρ k = k /, (1) where k is the (m 1) (m 1) determinant formed from the Gram determinant = det(p i p j ), by replacing the k th column by (p m p 1 ),, (p m p )
4 Consequently k is an integer We remark that in practice all ρ k tend to be small We have also observed that each shortest multiplier is always associated via (5) with a point (x 1,, x ) Z in the vicinity of (ρ 1,, ρ ) Let D be the set of points (z 1,, z ) Q satisfying z k + i=k+1 µ ik z i + µ mk < 1, k = m 1,, 1 (15) Then (ρ 1,, ρ ) D Also (0,, 0) D Definition We say Property G holds if for each shortest multiplier p, the x 1,, x of equation (5) satisfy (x 1,, x ) D Remarks 1 In view of (15), if the integer vector (x 1,, x ) D, there are at most two possibilities for each x k So if property G holds for a given m tuple (s 1,, s m ), then the number N of shortest multipliers is at most Given that a shortest multiplier p satisfies p p m, (7) and (8) show that property G holds if B k > p m for 1 k m 1 This is the case in Example 1 below, but does not always hold, as Example shows The examples where s i = for 1 i m 1, s m = m if m is odd, but s m = m 1 if m is even, produce values ( ) ( and ) m, respectively for N We prove below that property G holds when m = The least value of m for which it fails to hold appears to be 11, as in Example 5 Failures occur extremely rarely 5 We remark that in contrast to our situation, [Rosser 19] gave an example of the quadratic expression ( 7x 1y) + ( 7x y), which assumes its minimum value for integers x and y at the point ( 11, 0), whereas the minimum value 0 for real numbers x, y occurs at ( /, 80/) The algorithm We construct a sequence of multiplier vectors X K, K = m 1,, 1, which correspond to points (x 1,, x ) in D, close to (ρ 1,, ρ ), as follows Define x = 0,, x K+1 = 0 Then define x K,, x 1 Z recursively by:
5 (i) 0 if µ mk = 0 x K = 1 if µ mk < 0 1 if µ mk > 0; (ii) for 1 k < K, x k = σ k, where σ k = i=k+1 µ ik x i + µ mk (16) and θ is the nearest integer symbol, with θ = θ 1, if θ is a non negative half integer, but θ + 1 if θ is a negative half integer We also let X 0 = p m Remark For m =, a perusal of the list of shortest multipliers in the Appendix reveals that our algorithm will produce all the shortest multipliers For m = the algorithm appears to deliver at least one shortest multiplier For 5 m 10, whenever the algorithm fails to deliver a shortest multiplier, the excess length squared is always observed to be 1, as in Example In example 5, the excess is The case m = Lemma Let (x 1,, x ) Z correspond to a shortest multiplier p via equation (7) Then x + µ m ( ( ) 1 m ) Proof The inequality p p m implies Q(x 1,, x ) Q(0,, 0) Then equation (8) gives (x + µ m ) B µ mb + + µ m1b 1 (x + µ m ) µ m + µ B m mm + + µ m1 B ( B1 B ) m 5
6 Now () gives B k (1 µ kk 1 )B k 1 B k 1 Hence (x + µ m ) µ m + µ mm ( + + µ m1 ) ) m ( + + = 1 ( ( ) ) 1 1 = ( ) m ( ) m Corollary Property G holds if m = Proof Assume m = and that (x 1, x ) Z defines a minimum point for Q(x 1, x ) = (x + µ ) B + (x 1 + µ 1 x + µ 1 ) B 1 Then from the Lemma, we have x + µ 7/1 (17) Also the inequalities x (x + µ ) 0 and Q(x 1, x ) Q(0, 0) give x 1 + σ 1 = x 1 + µ 1 x + µ 1 µ 1 (18) Remark From the Corollary, it follows that N if m = In fact N if m = For if N =, we see from (17) and (18) that µ = ɛ 1 = ±1, µ 1 = ɛ = ±1, µ 1 = 0 Then ɛ 1 p 1 1 = (det P ) = det (p i p j ) = p p ɛ p ɛ 1 p 1 ɛ p p which gives = p 1 p ( p p 1 p p ) This leads to a contradiction, as the p i p j are integers The example (s 1, s, s ) = (1,, 9) shows that the bound N = is attained Here the quadratic expression in (8) is Q(x 1, x ) = 591(x ) + 6(x 1 7 x ) The unimodular matrix P = (ρ 1, ρ ) = ( 185, Also X, X 1, X 0 are given by 591 ) , 6
7 K (x 1, x ) X K X K (0, 1) (,, 6) 61 1 ( 1, 0) (,, 6) 61 0 (0, 0) (6, 0, 5) 61 and are the shortest multiplier vectors 5 Numerical results Example 1 Take s 1, s, s to be, 6, 9 The unimodular matrix P = The quadratic expression in (8) and the ρ k of (1) are given by ( ) Q(x 1, x ) = 1 x 6 (, x1 6 x 1 1) 6 (ρ 1, ρ ) = ( 90, 6 1 1) Also X, X 1, X 0 are given by 5 K (x 1, x ) X K X K (1, 1) (1, 1, 1) 1 (1, 0) (1,, 1) 6 0 (0, 0) (, 0, 1) 5 The shortest multiplier is X = p + p 1 + p Example Take s 1,, s 5 to be 10, 51, 10, 177, 07 The unimodular matrix P = 6 The quadratic expression in (8) is Q(x 1,, x ) = (ρ 1, ρ, ρ, ρ ) = x 6885 «+ 770 x x x x , , , x « x 1 10 x + 1 «10 ««+ 10 x 1 Also X,, X 0 are given by 7
8 K (x 1, x, x, x ) X K X K (0, 1, 0, 1) (1, 1,,, 0) 15 (0, 0, 1, 0) (0,, 1, 0, 1) 18 (0, 1, 0, 0) (, 0, 1,, ) 18 1 ( 1, 0, 0, 0) ( 1, 1,,, 0) 0 (0, 0, 0, 0) (, 0,, 1, 0) 1 The shortest multiplier is p = p 5 + p = [0, 1,,, ], with p = 1 Property G holds here Example (Example 7 of [Havas, Majewski, Matthews 1998]) Take s 1,, s 10 to be 7686, , 1119, , , 07775, 1179, 1675, , The unimodular matrix P = Then X 9,, X 0 are given by K (x 1,, x 9 ) X K X K 9 ( 1, 1, 1, 0, 1, 1, 0, 0, 1) (, 1, 1,, 1,,,,, ) 6 8 (0, 1, 0, 0, 1, 1, 0, 1, 0) (, 0,,,, 1, 1, 1, 1, 5) 6 7 (0, 0, 0, 0, 1, 1, 1, 0, 0) ( 1,, 1, 0, 0,,,,, ) 1 6 (0, 1, 0, 0, 1, 1, 0, 0, 0) ( 1,, 1, 0, 5, 0, 1, 0,, 1) 7 5 (0, 1, 1, 1, 1, 0, 0, 0, 0) (,, 1, 1,, 1, 5, 1,, 0) 55 (0, 1, 0, 1, 0, 0, 0, 0, 0) ( 1,, 1,,, 1, 1, 1,, 0) 50 (0, 0, 1, 0, 0, 0, 0, 0, 0) (, 0, 1,,, 0,, 0, 1, ) 51 (0, 1, 0, 0, 0, 0, 0, 0, 0) ( 1, 1, 1, 5,,, 0, 1,, 1) 59 1 ( 1, 0, 0, 0, 0, 0, 0, 0, 0) (1, 0,,, 1,,, 1, 1, 0) 5 0 (0, 0, 0, 0, 0, 0, 0, 0, 0) ( 1, 0, 1,, 1,,,,, ) Truncated to decimals, ρ = ( 01, 080, 0, 00, 076, 07, 05, 05, 00) X 9 is the shortest multiplier vector Property G holds here Example Take s 1,, s 0 to be 55, , 56, , 6565, , , , , , 65656, , , 555, , 56668, , , , , , , , , , , , , 8 7 5
9 , , , , , , , , 678, , , Here LLL delivers a p 0 with p 0 = 0 and our algorithm gives a multiplier X 7 with X 7 = 18 There are shortest multipliers, lengthsquared 1: p 0 p p p 5 p 6 + p 8 p 1 + p 15 + p 16 p 17 + p 18 p 1 + p + p = ( 1, 1, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1), p 0 p 1 p + p + p 5 + p 6 + p 8 p 9 + p 11 + p 1 + p 1 + p 16 p 17 + p 18 p 0 + p 1 + p + p 7 p 8 = (0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 1, 1, 0, 1, 0, 0, 1) Truncated to decimals, ρ = ( 01, 16, 007, 191, 07, 01, 06, 061, 09, 06, 0, 00, 185, 066, 009,\ 11, 167, 075, 06, 06, 019, 068, 019, 00, 015, 09, 0, 0, 000, 000, 000, 000,\ 000, 000, 000, 000, 000, 000, 000) Property G holds here Example 5 Take s 1,, s 11 to be , , , 07189, , , , , , , 1819 The shortest X K is X[7] = p 11 + p p 7 = (,,,,, 0,, 0,, 1, 1) with X 7 = 81 The shortest multiplier is p = p 11 + p p p 5 p 10 = (, 1, 1,,, 0,, 6, 1,, 0), with p = 79 Property G does not hold, as x 10 = 1 and σ 10 = / Example 6 The following random examples illustrate the improved multipliers X K that are produced by the algorithm in section The shortest multiplier vectors are unknown here 6 Appendix m p m X K X 96 = X 17 = X 80 = X 96 = 9 In this section we present a complete classification of the possible shortest multipliers when m =, based on inequalities (17) and (18): x + µ 7 1 and x 1 + µ 1 x + µ 1 µ 1 9
10 1 µ = 0, µ 1 < 1/: p µ = 0, µ 1 = 1/: p and p p 1 µ = 0, µ 1 = 1/: p and p + p 1 0 < µ 1/, 0 µ 1 < 1/: (i) µ 1 0: p or p p (ii) µ 1 < 0: p or p p 1 p 5 0 < µ 1/, µ 1 = 1/: (Here p = p p 1 ) (i) µ 1 > 0: p and p p 1, or p p (ii) µ 1 < 0: p and p p 1, or p p 1 p (iii) µ 1 = 0: p and p p 1 Note: µ < 1/ here For if µ = 1/, we have shortest multipliers: 6 0 < µ 1/, 1/ < µ 1 0: (i) µ 1 > 0: p or p + p 1 p (ii) µ 1 0: p or p p p, p p 1, p p, p p 1 p 7 0 < µ 1/, µ 1 = 1/: (Here p = p + p 1 ) (i) µ 1 > 0: p and p + p 1, or p + p 1 p (ii) µ 1 < 0: p and p + p 1, or p p (iii) µ 1 = 0: p and p + p 1 (µ < 1/) 8 1/ µ < 0, 1/ < µ 1 0: (i) µ 1 0: p or p + p (ii) µ 1 < 0: p or p + p 1 + p 9 1/ µ < 0, µ 1 = 1/: (Here p = p + p 1 ) (i) µ 1 > 0: p and p + p 1, or p + p 10
11 (ii) µ 1 < 0: p and p + p 1, or p + p 1 + p (iii) µ 1 = 0: p and p + p 1 ( 1/ < µ ) 10 1/ µ < 0, 0 µ 1 < 1/: (i) µ 1 0: p or p + p (ii) µ 1 > 0: p or p p 1 + p 11 1/ µ < 0, µ 1 = 1/: (Here p = p p 1 ) (i) µ 1 > 0: p and p p 1, or p p 1 + p (ii) µ 1 < 0: p and p p 1, or p + p (iii) µ 1 = 0: p and p p 1 ( 1/ < µ ) References [Havas, Majewski, Matthews 1998] G Havas, BS Majewski, KR Mat thews, Extended gcd and Hermite normal form algorithms via lattice reduction, Experimental Mathematics 7 No (1998) [Rosser 19] JB Rosser, A generalization of the Euclidean algorithm to several dimensions, Duke Math 9 (19) K R Matthews, 11
Extended gcd and Hermite normal form algorithms via lattice basis reduction
Extended gcd and Hermite normal form algorithms via lattice basis reduction George Havas School of Information Technology The University of Queensland Queensland 4072, Australia URL: http://www.it.uq.edu.au/personal/havas/
More informationComplexity of the Havas, Majewski, Matthews LLL. Mathematisch Instituut, Universiteit Utrecht. P.O. Box
J. Symbolic Computation (2000) 11, 1{000 Complexity of the Havas, Majewski, Matthews LLL Hermite Normal Form algorithm WILBERD VAN DER KALLEN Mathematisch Instituut, Universiteit Utrecht P.O. Box 80.010
More informationMinimal multipliers for consecutive Fibonacci numbers
ACTA ARITHMETICA LXXV3 (1996) Minimal multipliers for consecutive Fibonacci numbers by K R Matthews (Brisbane) 1 Introduction The Fibonacci and Lucas numbers F n, L n (see [2]) are defined by Since F 1
More informationReduction of Smith Normal Form Transformation Matrices
Reduction of Smith Normal Form Transformation Matrices G. Jäger, Kiel Abstract Smith normal form computations are important in group theory, module theory and number theory. We consider the transformation
More informationHomework #2 Solutions Due: September 5, for all n N n 3 = n2 (n + 1) 2 4
Do the following exercises from the text: Chapter (Section 3):, 1, 17(a)-(b), 3 Prove that 1 3 + 3 + + n 3 n (n + 1) for all n N Proof The proof is by induction on n For n N, let S(n) be the statement
More informationLattice Basis Reduction and the LLL Algorithm
Lattice Basis Reduction and the LLL Algorithm Curtis Bright May 21, 2009 1 2 Point Lattices A point lattice is a discrete additive subgroup of R n. A basis for a lattice L R n is a set of linearly independent
More informationCSE 206A: Lattice Algorithms and Applications Spring Basis Reduction. Instructor: Daniele Micciancio
CSE 206A: Lattice Algorithms and Applications Spring 2014 Basis Reduction Instructor: Daniele Micciancio UCSD CSE No efficient algorithm is known to find the shortest vector in a lattice (in arbitrary
More informationA New Algorithm and Refined Bounds for Extended Gcd Computation
A New Algorithm and Refined Bounds for Extended Gcd Computation David Ford* and George Havas** Department of Computer Science, Concordia University, Montrfial, Qufibec, Canada H3G 1M8 and Department of
More informationMath 324 Summer 2012 Elementary Number Theory Notes on Mathematical Induction
Math 4 Summer 01 Elementary Number Theory Notes on Mathematical Induction Principle of Mathematical Induction Recall the following axiom for the set of integers. Well-Ordering Axiom for the Integers If
More informationShort solutions for a linear Diophantine equation
Lecturas Matemáticas Volumen 27 (2006), páginas 5 16 Short solutions for a linear Diophantine equation H. Esmaeili Bu-Ali Sina University, Hamedan, Iran Abstract. It is known that finding the shortest
More information1 Shortest Vector Problem
Lattices in Cryptography University of Michigan, Fall 25 Lecture 2 SVP, Gram-Schmidt, LLL Instructor: Chris Peikert Scribe: Hank Carter Shortest Vector Problem Last time we defined the minimum distance
More informationCSE 206A: Lattice Algorithms and Applications Winter The dual lattice. Instructor: Daniele Micciancio
CSE 206A: Lattice Algorithms and Applications Winter 2016 The dual lattice Instructor: Daniele Micciancio UCSD CSE 1 Dual Lattice and Dual Basis Definition 1 The dual of a lattice Λ is the set ˆΛ of all
More informationNeville s Method. MATH 375 Numerical Analysis. J. Robert Buchanan. Fall Department of Mathematics. J. Robert Buchanan Neville s Method
Neville s Method MATH 375 Numerical Analysis J. Robert Buchanan Department of Mathematics Fall 2013 Motivation We have learned how to approximate a function using Lagrange polynomials and how to estimate
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 informationProblem Set 5 Solutions
Problem Set 5 Solutions Section 4.. Use mathematical induction to prove each of the following: a) For each natural number n with n, n > + n. Let P n) be the statement n > + n. The base case, P ), is true
More informationIrreducible Polynomials over Finite Fields
Chapter 4 Irreducible Polynomials over Finite Fields 4.1 Construction of Finite Fields As we will see, modular arithmetic aids in testing the irreducibility of polynomials and even in completely factoring
More information47-831: Advanced Integer Programming Lecturer: Amitabh Basu Lecture 2 Date: 03/18/2010
47-831: Advanced Integer Programming Lecturer: Amitabh Basu Lecture Date: 03/18/010 We saw in the previous lecture that a lattice Λ can have many bases. In fact, if Λ is a lattice of a subspace L with
More informationMATH 324 Summer 2011 Elementary Number Theory. Notes on Mathematical Induction. Recall the following axiom for the set of integers.
MATH 4 Summer 011 Elementary Number Theory Notes on Mathematical Induction Principle of Mathematical Induction Recall the following axiom for the set of integers. Well-Ordering Axiom for the Integers If
More informationThe Shortest Vector Problem (Lattice Reduction Algorithms)
The Shortest Vector Problem (Lattice Reduction Algorithms) Approximation Algorithms by V. Vazirani, Chapter 27 - Problem statement, general discussion - Lattices: brief introduction - The Gauss algorithm
More informationMathematical Optimisation, Chpt 2: Linear Equations and inequalities
Mathematical Optimisation, Chpt 2: Linear Equations and inequalities Peter J.C. Dickinson p.j.c.dickinson@utwente.nl http://dickinson.website version: 12/02/18 Monday 5th February 2018 Peter J.C. Dickinson
More informationDS-GA 1002 Lecture notes 0 Fall Linear Algebra. These notes provide a review of basic concepts in linear algebra.
DS-GA 1002 Lecture notes 0 Fall 2016 Linear Algebra These notes provide a review of basic concepts in linear algebra. 1 Vector spaces You are no doubt familiar with vectors in R 2 or R 3, i.e. [ ] 1.1
More informationCHAPTER 8: EXPLORING R
CHAPTER 8: EXPLORING R LECTURE NOTES FOR MATH 378 (CSUSM, SPRING 2009). WAYNE AITKEN In the previous chapter we discussed the need for a complete ordered field. The field Q is not complete, so we constructed
More informationThus f is continuous at x 0. Matthew Straughn Math 402 Homework 6
Matthew Straughn Math 402 Homework 6 Homework 6 (p. 452) 14.3.3, 14.3.4, 14.3.5, 14.3.8 (p. 455) 14.4.3* (p. 458) 14.5.3 (p. 460) 14.6.1 (p. 472) 14.7.2* Lemma 1. If (f (n) ) converges uniformly to some
More informationExtended GCD and Hermite Normal Form Algorithms via Lattice Basis Reduction
Extended GCD and Hermite Normal Form Algorithms via Lattice Basis Reduction George Havas, Bohdan S. Majewski, and Keith R. Matthews CONTENTS 1. Introduction. The LLL Algorithm. Analysis of Algorithm. Algorithm
More informationThe Abundancy index of divisors of odd perfect numbers Part III
Notes on Number Theory Discrete Mathematics Print ISSN 1310 5132, Online ISSN 2367 8275 Vol 23, 2017, No 3, 53 59 The Abundancy index of divisors of odd perfect numbers Part III Jose Arnaldo B Dris Department
More information1: Introduction to Lattices
CSE 206A: Lattice Algorithms and Applications Winter 2012 Instructor: Daniele Micciancio 1: Introduction to Lattices UCSD CSE Lattices are regular arrangements of points in Euclidean space. The simplest
More information1 2 3 style total. Circle the correct answer; no explanation is required. Each problem in this section counts 5 points.
1 2 3 style total Math 415 Examination 3 Please print your name: Answer Key 1 True/false Circle the correct answer; no explanation is required. Each problem in this section counts 5 points. 1. The rings
More informationMathematical Induction
Mathematical Induction MAT231 Transition to Higher Mathematics Fall 2014 MAT231 (Transition to Higher Math) Mathematical Induction Fall 2014 1 / 21 Outline 1 Mathematical Induction 2 Strong Mathematical
More informationCSC 2414 Lattices in Computer Science September 27, Lecture 4. An Efficient Algorithm for Integer Programming in constant dimensions
CSC 2414 Lattices in Computer Science September 27, 2011 Lecture 4 Lecturer: Vinod Vaikuntanathan Scribe: Wesley George Topics covered this lecture: SV P CV P Approximating CVP: Babai s Nearest Plane Algorithm
More informationNote that r = 0 gives the simple principle of induction. Also it can be shown that the principle of strong induction follows from simple induction.
Proof by mathematical induction using a strong hypothesis Occasionally a proof by mathematical induction is made easier by using a strong hypothesis: To show P(n) [a statement form that depends on variable
More informationBALANCED INTEGER SOLUTIONS OF LINEAR EQUATIONS
BALANCED INTEGER SOLUTIONS OF LINEAR EQUATIONS KONSTANTINOS A. DRAZIOTIS Abstract. We use lattice based methods in order to get an integer solution of the linear equation a x + +a nx n = a 0, which satisfies
More informationA strongly polynomial algorithm for linear systems having a binary solution
A strongly polynomial algorithm for linear systems having a binary solution Sergei Chubanov Institute of Information Systems at the University of Siegen, Germany e-mail: sergei.chubanov@uni-siegen.de 7th
More informationBest approximation in the 2-norm
Best approximation in the 2-norm Department of Mathematical Sciences, NTNU september 26th 2012 Vector space A real vector space V is a set with a 0 element and three operations: Addition: x, y V then x
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationProof Terminology. Technique #1: Direct Proof. Learning objectives. Proof Techniques (Rosen, Sections ) Direct Proof:
Proof Terminology Proof Techniques (Rosen, Sections 1.7 1.8) TOPICS Direct Proofs Proof by Contrapositive Proof by Contradiction Proof by Cases Theorem: statement that can be shown to be true Proof: a
More informationRATIONAL REALIZATION OF MAXIMUM EIGENVALUE MULTIPLICITY OF SYMMETRIC TREE SIGN PATTERNS. February 6, 2006
RATIONAL REALIZATION OF MAXIMUM EIGENVALUE MULTIPLICITY OF SYMMETRIC TREE SIGN PATTERNS ATOSHI CHOWDHURY, LESLIE HOGBEN, JUDE MELANCON, AND RANA MIKKELSON February 6, 006 Abstract. A sign pattern is a
More informationMath Circle Beginners Group February 28, 2016 Euclid and Prime Numbers Solutions
Math Circle Beginners Group February 28, 2016 Euclid and Prime Numbers Solutions Warm-up Problems 1. What is a prime number? Give an example of an even prime number and an odd prime number. A prime number
More information12. Perturbed Matrices
MAT334 : Applied Linear Algebra Mike Newman, winter 208 2. Perturbed Matrices motivation We want to solve a system Ax = b in a context where A and b are not known exactly. There might be experimental errors,
More informationa s 1.3 Matrix Multiplication. Know how to multiply two matrices and be able to write down the formula
Syllabus for Math 308, Paul Smith Book: Kolman-Hill Chapter 1. Linear Equations and Matrices 1.1 Systems of Linear Equations Definition of a linear equation and a solution to a linear equations. Meaning
More informationHow to Factor N 1 and N 2 When p 1 = p 2 mod 2 t
How to Factor N 1 and N 2 When p 1 = p 2 mod 2 t Kaoru Kurosawa and Takuma Ueda Ibaraki University, Japan Abstract. Let N 1 = p 1q 1 and N 2 = p 2q 2 be two different RSA moduli. Suppose that p 1 = p 2
More informationCSE 206A: Lattice Algorithms and Applications Spring Minkowski s theorem. Instructor: Daniele Micciancio
CSE 206A: Lattice Algorithms and Applications Spring 2014 Minkowski s theorem Instructor: Daniele Micciancio UCSD CSE There are many important quantities associated to a lattice. Some of them, like the
More informationA Subrecursive Refinement of FTA. of the Fundamental Theorem of Algebra
A Subrecursive Refinement of the Fundamental Theorem of Algebra P. Peshev D. Skordev Faculty of Mathematics and Computer Science University of Sofia Computability in Europe, 2006 Outline Introduction 1
More informationSolutions to Second In-Class Exam: Math 401 Section 0201, Professor Levermore Friday, 16 April 2010
Solutions to Second In-Class Exam: Math 4 Section, Professor Levermore Friday, 6 April No notes, books, or electrontics You must show your reasoning for full credit Good luck! [] Consider the polynomials
More informationThe Geometric Distribution
MATH 382 The Geometric Distribution Dr. Neal, WKU Suppose we have a fixed probability p of having a success on any single attempt, where p > 0. We continue to make independent attempts until we succeed.
More informationTopics in Basis Reduction and Integer Programming
Topics in Basis Reduction and Integer Programming Mustafa Kemal Tural A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial fulfillment of the requirements
More informationInvertible Matrices over Idempotent Semirings
Chamchuri Journal of Mathematics Volume 1(2009) Number 2, 55 61 http://www.math.sc.chula.ac.th/cjm Invertible Matrices over Idempotent Semirings W. Mora, A. Wasanawichit and Y. Kemprasit Received 28 Sep
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 informationPYTHAGOREAN DESCENT KEITH CONRAD
PYTHAGOREAN DESCENT KEITH CONRAD 1. Introduction A Pythagorean triple is a triple of positive integers (a, b, c) such that a 2 + b 2 = c 2. Examples include (3, 4, 5), (5, 12, 13), and (6, 8, 10). A Pythagorean
More informationLEGENDRE S THEOREM, LEGRANGE S DESCENT
LEGENDRE S THEOREM, LEGRANGE S DESCENT SUPPLEMENT FOR MATH 370: NUMBER THEORY Abstract. Legendre gave simple necessary and sufficient conditions for the solvablility of the diophantine equation ax 2 +
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More informationComposition of Axial Functions of Product of Finite Sets
Composition of Axial Functions of Product of Finite Sets Krzysztof P lotka Department of Mathematics University of Scranton Scranton, PA 18510, USA plotkak2@scranton.edu September 14, 2004 Abstract To
More informationPYTHAGOREAN DESCENT KEITH CONRAD
PYTHAGOREAN DESCENT KEITH CONRAD 1. Introduction A Pythagorean triple is a triple of positive integers (a, b, c) such that a 2 + b 2 = c 2. Examples include (3, 4, 5), (5, 12, 13), and (6, 8, 10). A Pythagorean
More informationMath 443 Differential Geometry Spring Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook.
Math 443 Differential Geometry Spring 2013 Handout 3: Bilinear and Quadratic Forms This handout should be read just before Chapter 4 of the textbook. Endomorphisms of a Vector Space This handout discusses
More informationDefinition 6.1 (p.277) A positive integer n is prime when n > 1 and the only positive divisors are 1 and n. Alternatively
6 Prime Numbers Part VI of PJE 6.1 Fundamental Results Definition 6.1 (p.277) A positive integer n is prime when n > 1 and the only positive divisors are 1 and n. Alternatively D (p) = { p 1 1 p}. Otherwise
More informationHomework #2 solutions Due: June 15, 2012
All of the following exercises are based on the material in the handout on integers found on the class website. 1. Find d = gcd(475, 385) and express it as a linear combination of 475 and 385. That is
More informationShort Vectors of Planar Lattices via Continued Fractions. Friedrich Eisenbrand Max-Planck-Institut für Informatik
Short Vectors of Planar Lattices via Continued Fractions Friedrich Eisenbrand Max-Planck-Institut für Informatik Outline Lattices Definition of planar integral lattices Shortest Vectors Applications The
More informationAn LLL-Reduction Algorithm with Quasi-linear Time Complexity 1
An LLL-Reduction Algorithm with Quasi-linear Time Complexity 1 Andrew Novocin, Damien Stehlé, and Gilles Villard CNRS, ENS de Lyon, INRIA, UCBL, U. Lyon Laboratoire LIP 46 Allée d Italie, 69364 Lyon Cedex
More informationMath 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.
Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,
More informationBricklaying and the Hermite Normal Form
Bricklaying and the Hermite Normal Form William J. Gilbert Pure Mathematics Department University of Waterloo Waterloo, Ontario Canada N2L 3G1 AMS Classifications: Primary: 15A36 Appeared in the American
More informationNotes on Systems of Linear Congruences
MATH 324 Summer 2012 Elementary Number Theory Notes on Systems of Linear Congruences In this note we will discuss systems of linear congruences where the moduli are all different. Definition. Given the
More informationLinear Models Review
Linear Models Review Vectors in IR n will be written as ordered n-tuples which are understood to be column vectors, or n 1 matrices. A vector variable will be indicted with bold face, and the prime sign
More informationPrimitive sets in a lattice
Primitive sets in a lattice Spyros. S. Magliveras Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431, U.S.A spyros@fau.unl.edu Tran van Trung Institute for Experimental
More informationJournal of Integer Sequences, Vol. 5 (2002), Article
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 5 (2002), Article 02.2.7 ON SHANKS ALGORITHM FOR COMPUTING THE CONTINUED FRACTION OF log b a TERENCE JACKSON 1 AND KEITH MATTHEWS 2 1 Department of Mathematics
More informationShort proofs of the universality of certain diagonal quadratic forms
Arch. Math. 91 (008), 44 48 c 008 Birkhäuser Verlag Basel/Switzerland 0003/889X/010044-5, published online 008-06-5 DOI 10.1007/s00013-008-637-5 Archiv der Mathematik Short proofs of the universality of
More informationLooking back at lattice-based cryptanalysis
September 2009 Lattices A lattice is a discrete subgroup of R n Equivalently, set of integral linear combinations: α 1 b1 + + α n bm with m n Lattice reduction Lattice reduction looks for a good basis
More informationProof by Contradiction
Proof by Contradiction MAT231 Transition to Higher Mathematics Fall 2014 MAT231 (Transition to Higher Math) Proof by Contradiction Fall 2014 1 / 12 Outline 1 Proving Statements with Contradiction 2 Proving
More informationMATRIX ALGEBRA. or x = (x 1,..., x n ) R n. y 1 y 2. x 2. x m. y m. y = cos θ 1 = x 1 L x. sin θ 1 = x 2. cos θ 2 = y 1 L y.
as Basics Vectors MATRIX ALGEBRA An array of n real numbers x, x,, x n is called a vector and it is written x = x x n or x = x,, x n R n prime operation=transposing a column to a row Basic vector operations
More informationPlotkin s Bound in Codes Equipped with the Euclidean Weight Function
Tamsui Oxford Journal of Mathematical Sciences 5() (009) 07-4 Aletheia University Plotkin s Bound in Codes Equipped with the Euclidean Weight Function Sapna Jain Department of Mathematics, University of
More informationTopics in Graph Theory
Topics in Graph Theory September 4, 2018 1 Preliminaries A graph is a system G = (V, E) consisting of a set V of vertices and a set E (disjoint from V ) of edges, together with an incidence function End
More informationj=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.
Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. Let u = [u
More informationNote on shortest and nearest lattice vectors
Note on shortest and nearest lattice vectors Martin Henk Fachbereich Mathematik, Sekr. 6-1 Technische Universität Berlin Straße des 17. Juni 136 D-10623 Berlin Germany henk@math.tu-berlin.de We show that
More informationLattice Reductions over Euclidean Rings with Applications to Cryptanalysis
IMACC 2017 December 12 14, 2017 Lattice Reductions over Euclidean Rings with Applications to Cryptanalysis Taechan Kim and Changmin Lee NTT Secure Platform Laboratories, Japan and Seoul National University,
More informationA Course in Computational Algebraic Number Theory
Henri Cohen 2008 AGI-Information Management Consultants May be used for personal purporses only or by libraries associated to dandelon.com network. A Course in Computational Algebraic Number Theory Springer
More informationMath Circle Beginners Group February 28, 2016 Euclid and Prime Numbers
Math Circle Beginners Group February 28, 2016 Euclid and Prime Numbers Warm-up Problems 1. What is a prime number? Give an example of an even prime number and an odd prime number. (a) Circle the prime
More informationSection 19 Integral domains
Section 19 Integral domains Instructor: Yifan Yang Spring 2007 Observation and motivation There are rings in which ab = 0 implies a = 0 or b = 0 For examples, Z, Q, R, C, and Z[x] are all such rings There
More informationDiscrete Math, Second Problem Set (June 24)
Discrete Math, Second Problem Set (June 24) REU 2003 Instructor: Laszlo Babai Scribe: D Jeremy Copeland 1 Number Theory Remark 11 For an arithmetic progression, a 0, a 1 = a 0 +d, a 2 = a 0 +2d, to have
More informationAlgorithms for Picture Analysis. Lecture 07: Metrics. Axioms of a Metric
Axioms of a Metric Picture analysis always assumes that pictures are defined in coordinates, and we apply the Euclidean metric as the golden standard for distance (or derived, such as area) measurements.
More information12x + 18y = 50. 2x + v = 12. (x, v) = (6 + k, 2k), k Z.
Math 3, Fall 010 Assignment 3 Solutions Exercise 1. Find all the integral solutions of the following linear diophantine equations. Be sure to justify your answers. (i) 3x + y = 7. (ii) 1x + 18y = 50. (iii)
More informationLattices. A Lattice is a discrete subgroup of the additive group of n-dimensional space R n.
Lattices A Lattice is a discrete subgroup of the additive group of n-dimensional space R n. Lattices have many uses in cryptography. They may be used to define cryptosystems and to break other ciphers.
More information10.1 Radical Expressions and Functions Math 51 Professor Busken
0. Radical Expressions and Functions Math 5 Professor Busken Objectives Find square roots without a calculator Simplify expressions of the form n a n Evaluate radical functions and find the domain of radical
More informationHomework 4, 5, 6 Solutions. > 0, and so a n 0 = n + 1 n = ( n+1 n)( n+1+ n) 1 if n is odd 1/n if n is even diverges.
2..2(a) lim a n = 0. Homework 4, 5, 6 Solutions Proof. Let ɛ > 0. Then for n n = 2+ 2ɛ we have 2n 3 4+ ɛ 3 > ɛ > 0, so 0 < 2n 3 < ɛ, and thus a n 0 = 2n 3 < ɛ. 2..2(g) lim ( n + n) = 0. Proof. Let ɛ >
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 informationMath 109 HW 9 Solutions
Math 109 HW 9 Solutions Problems IV 18. Solve the linear diophantine equation 6m + 10n + 15p = 1 Solution: Let y = 10n + 15p. Since (10, 15) is 5, we must have that y = 5x for some integer x, and (as we
More informationHard Instances of Lattice Problems
Hard Instances of Lattice Problems Average Case - Worst Case Connections Christos Litsas 28 June 2012 Outline Abstract Lattices The Random Class Worst-Case - Average-Case Connection Abstract Christos Litsas
More informationPOSITIVE DEFINITE n-regular QUADRATIC FORMS
POSITIVE DEFINITE n-regular QUADRATIC FORMS BYEONG-KWEON OH Abstract. A positive definite integral quadratic form f is called n- regular if f represents every quadratic form of rank n that is represented
More informationDefinite Quadratic Forms over F q [x]
Definite Quadratic Forms over F q [x] Larry J. Gerstein Department of Mathematics University of California Santa Barbara, CA 93106 E-mail: gerstein@math.ucsb.edu Version: September 30, 2002 ABSTRACT. Let
More informationUSING SHANKS BABY-STEP GIANT-STEP METHOD TO SOLVE THE GENERALIZED PELL EQUATION x 2 Dy 2 = N. Copyright 2009 by John P. Robertson. 1.
USING SHANKS BABY-STEP GIANT-STEP METHOD TO SOLVE THE GENERALIZED PELL EQUATION x 2 Dy 2 = N Abstract. For D > 0 not a square, and N 0, the continued fraction algorithm can be used to solve the generalized
More informationMathematical Induction
Mathematical Induction MAT30 Discrete Mathematics Fall 018 MAT30 (Discrete Math) Mathematical Induction Fall 018 1 / 19 Outline 1 Mathematical Induction Strong Mathematical Induction MAT30 (Discrete Math)
More informationSydney University Mathematical Society Problems Competition Solutions.
Sydney University Mathematical Society Problems Competition 005 Solutions 1 Suppose that we look at the set X n of strings of 0 s and 1 s of length n Given a string ɛ = (ɛ 1,, ɛ n ) X n, we are allowed
More information6 Permutations Very little of this section comes from PJE.
6 Permutations Very little of this section comes from PJE Definition A permutation (p147 of a set A is a bijection ρ : A A Notation If A = {a b c } and ρ is a permutation on A we can express the action
More information198 VOLUME 46/47, NUMBER 3
LAWRENCE SOMER Abstract. Rotkiewicz has shown that there exist Fibonacci pseudoprimes having the forms p(p + 2), p(2p 1), and p(2p + 3), where all the terms in the products are odd primes. Assuming Dickson
More informationLecture 5. Symmetric Shearer s Lemma
Stanfor University Spring 208 Math 233: Non-constructive methos in combinatorics Instructor: Jan Vonrák Lecture ate: January 23, 208 Original scribe: Erik Bates Lecture 5 Symmetric Shearer s Lemma Here
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 informationMATH 51H Section 4. October 16, Recall what it means for a function between metric spaces to be continuous:
MATH 51H Section 4 October 16, 2015 1 Continuity Recall what it means for a function between metric spaces to be continuous: Definition. Let (X, d X ), (Y, d Y ) be metric spaces. A function f : X Y is
More informationNumber Theory Proof Portfolio
Number Theory Proof Portfolio Jordan Rock May 12, 2015 This portfolio is a collection of Number Theory proofs and problems done by Jordan Rock in the Spring of 2014. The problems are organized first by
More informationDynamic Programming: Shortest Paths and DFA to Reg Exps
CS 374: Algorithms & Models of Computation, Spring 207 Dynamic Programming: Shortest Paths and DFA to Reg Exps Lecture 8 March 28, 207 Chandra Chekuri (UIUC) CS374 Spring 207 / 56 Part I Shortest Paths
More informationHOMEWORK 11 MATH 4753
HOMEWORK 11 MATH 4753 Recall that R = Z[x]/(x N 1) where N > 1. For p > 1 any modulus (not necessarily prime), R p = (Z/pZ)[x]/(x N 1). We do not assume p, q are prime below unless otherwise stated. Question
More informationMath 61CM - Solutions to homework 6
Math 61CM - Solutions to homework 6 Cédric De Groote November 5 th, 2018 Problem 1: (i) Give an example of a metric space X such that not all Cauchy sequences in X are convergent. (ii) Let X be a metric
More informationIn N we can do addition, but in order to do subtraction we need to extend N to the integers
Chapter The Real Numbers.. Some Preliminaries Discussion: The Irrationality of 2. We begin with the natural numbers N = {, 2, 3, }. In N we can do addition, but in order to do subtraction we need to extend
More informationOn the second smallest prime non-residue
On the second smallest prime non-residue Kevin J. McGown 1 Department of Mathematics, University of California, San Diego, 9500 Gilman Drive, La Jolla, CA 92093 Abstract Let χ be a non-principal Dirichlet
More information