arxiv: v2 [math.co] 27 Nov 2008
|
|
- Imogene Hunter
- 6 years ago
- Views:
Transcription
1 Stablty results for unquely determned sets from two drectons n dscrete tomography Brgt van Dalen arxv: v [mathco] 7 Nov 008 Mathematsch Insttuut, Unverstet Leden, Nels Bohrweg 1, 333 CA Leden, The Netherlands Abstract In ths paper we prove several new stablty results for the reconstructon of bnary mages from two projectons We consder an orgnal mage that s unquely determned by ts projectons and possble reconstructons from slghtly dfferent projectons We show that for a gven dfference n the projectons, the reconstructon can only be dsjont from the orgnal mage f the sze of the mage s not too large We also prove an upper bound for the sze of the mage gven the error n the projectons and the sze of the ntersecton between the mage and the reconstructon Key words: Dscrete tomography, Stablty, Image reconstructon, Bnary mage 1 Introducton Dscrete tomography s concerned wth problems such as reconstructng bnary mages on a lattce from gven projectons n lattce drectons [6] Each pont of a bnary mage has a value equal to zero or one The lne sum of a lne through the mage s the sum of the values of the ponts on ths lne The projecton of the mage n a certan lattce drecton conssts of all the lne sums of the lnes through the mage n ths drecton Several problems related to the reconstructon of bnary mages from two or more projectons have been descrbed n the lterature [6,7] Already n 1957, Ryser gave an algorthm to reconstruct bnary mages from ther horzontal Emal address: dalen@mathledenunvnl (Brgt van Dalen) Preprnt submtted to Elsever 6 January 01
2 and vertcal projectons and charactersed the set of projectons that correspond to a unque bnary mage [11] For any set of drectons, t s possble to construct mages that are not unquely determned by ther projectons n those drectons [6, Theorem 31] The problem of decdng whether an mage s unquely determned by ts projectons and the problem of reconstructng t are NP-hard for any set of more than two drectons [] Asde from varous nterestng theoretcal problems, dscrete tomography also has applcatons n a wde range of felds The most mportant are electron mcroscopy [8] and medcal magng [5,13], but there are also applcatons n nuclear scence [9,10] and varous other felds [1,15] An nterestng problem n dscrete tomography s the stablty of reconstructons Even f an mage s unquely determned by ts projectons, a very small error n the projectons may lead to a completely dfferent reconstructon [1,3] Alpers et al [1,] showed that n the case of two drectons a total error of at most n the projectons can only cause a small dfference n the reconstructon They also proved a lower bound on the error f the reconstructon s dsjont from the orgnal mage In ths paper we mprove ths bound, and we resolve the open problem of stablty wth a projecton error greater than Notaton and statement of the problems Let F 1 and F be two fnte subsets of Z wth characterstc functons χ 1 and χ (That s, χ h (x, y) = 1 f and only f (x, y) F h, h {1, }) For Z, we defne row as the set {(x, y) Z : x = } We call the ndex of the row For j Z, we defne column j as the set {(x, y) Z : y = j} We call j the ndex of the column Followng matrx notaton, we use row numbers that ncrease when gong downwards and column numbers that ncrease when gong to the rght The row sum r (h) s the number of elements of F h n row, that s r (h) = j Z χ h (, j) The column sum c (h) j of F h s the number of elements of F h n column j, that s c (h) j = Z χ h (, j) We refer to both row and column sums as the lne sums of F h Throughout ths paper, we assume that F 1 s unquely determned by ts row and column sums Such sets were studed by, among others, Ryser [11] and Wang [1] Let a be the number of rows and b the number of columns that
3 Fg 1 A unquely determned set wth the assumed row and column orderng contan elements of F 1 We renumber the rows and columns such that we have r (1) 1 r (1) r (1) a > 0, c (1) 1 c (1) c (1) b > 0, and such that all elements of F are contaned n rows and columns wth postve ndces By [1, Theorem 3] we have the followng property of F 1 (see Fgure 1): n row the elements of F 1 are precsely the ponts (, 1), (, ),, (, r (1) ), n column j the elements of F 1 are precsely the ponts (1, j), (, j),, (c (1) j, j) We wll refer to ths property as the trangular shape of F 1 Everywhere except n Secton 6 we assume that F 1 = F Note that we do not assume F to be unquely determned As F 1 and F are dfferent and F 1 s unquely determned by ts lne sums, F cannot have exactly the same lne sums as F 1 Defne the dfference or error n the lne sums as c (1) j c () j + r (1) r () j 1 1 As n general t s t + s mod, the above expresson s congruent to j 1 ( (1) c j + c () j ) + 1 ( (1) r + r () ) F1 + F 0 mod, hence the error n the lne sums s always even We wll denote t by α, where α s a postve nteger For notatonal convenence, we wll often wrte p for F 1 F We consder two problems concernng stablty Problem 1 Suppose F 1 F = How large can F 1 be n terms of α? 3
4 Alpers et al [, Theorem 9] proved that F 1 α They also showed that there s no constant c such that F 1 cα for all F 1 and F In Secton of ths paper we wll prove the new bound F 1 α(1 + log α) and show that ths bound s asymptotcally sharp Problem How small can F 1 F be n terms of F 1 and α, or, equvalently, how large can F 1 be n terms of F 1 F and α? Alpers ([1, Theorem 5118]) showed n the case α = 1 that F 1 F F F Ths bound s sharp: f F 1 = 1 n(n + 1) for some postve nteger n, then there exsts an example for whch equalty holds A smlar result s stated n [, Theorem 19] Whle [1,] only deal wth the case α = 1, we wll gve stablty results for general α In Secton 5 we wll gve two dfferent upper bounds for F 1 The bounds have dfferent asymptotc behavour Wrtng p for F 1 F, the second bound reduces to F 1 p p + 1 n case α = 1, whch s equvalent to p F 1 F 1 Hence the second new bound can be vewed as a generalsaton of Alpers bound The frst new bound s dfferent and better n the case that α s very large In Secton 6 we wll generalse the results to the case F 1 F 3 Starcases Alpers ntroduced the noton of a starcase to characterse F 1 F n the case α = 1 We wll use a slghtly dfferent defnton and then show that for general α the symmetrc dfference F 1 F conssts of α starcases Defnton 3 A set of ponts (p 1, p,, p n ) n Z s called a starcase f the followng two condtons are satsfed: for each wth 1 n 1 one of the ponts p and p +1 s an element of F 1 \F and the other s an element of F \F 1 ;
5 Fg A starcase The set F 1 conssts of the whte and the black-and-whte ponts, whle F conssts of the black and the black-and-whte ponts The starcase s ndcated by the dashed lne segments ether for all the ponts p and p +1 are n the same column and the ponts p +1 and p + are n the same row, or for all the ponts p and p +1 are n the same row and the ponts p +1 and p + are n the same column Ths defnton s dfferent from [1,] n the followng way Frstly, the number of ponts does not need to be even Secondly, the ponts p 1 and p n can both be ether n F 1 \F or n F \F 1 So ths defnton s slghtly more general than the one used n [1,] for the case α = 1 Consder a pont p F 1 \F of a starcase (p 1, p,, p n ) Assume p 1 s n the same column as p and p +1 s n the same row as p Because of the trangular shape of F 1, the row ndex of p 1 must be larger than the row ndex of p, and the column ndex of p +1 must be larger than the column ndex of p Therefore, the starcase looks lke a real-world starcase (see Fgure ) From now on, we assume for all starcases that p 1 s the pont wth the largest row ndex and the smallest column ndex, whle p n s the pont wth the smallest row ndex and the largest column ndex We say that the starcase begns wth p 1 and ends wth p n Lemma Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums, and F 1 = F Let α be defned as n Secton Then the set F 1 F s the dsjont unon of α starcases PROOF We wll construct the starcases one by one and delete them from F 1 F For a subset A of F 1 F, defne ρ (A) = {j Z : (, j) A F 1 } {j Z : (, j) A F }, Z, σ j (A) = { Z : (, j) A F 1 } { Z : (, j) A F }, j Z, τ(a) = ρ (A) + j σ j (A) We have α = τ(f 1 F ) 5
6 Assume that the rows and columns are ordered as n Secton Because of the trangular shape of F 1, for any pont (, j) F 1 \F and any pont (k, l) F \F 1 we then have k > or l > j Suppose we have deleted some starcases and are now left wth a non-empty subset A of F 1 F Let (p 1, p,, p n ) be a starcase of maxmal length that s contaned n A Let (x 1, y 1 ) and (x n, y n ) be the coordnates of the ponts p 1 and p n respectvely Each of those two ponts can be ether n A F 1 or n A F, so there are four dfferent cases (If n = 1, so p 1 and p n are the same pont, then there are only two cases) We consder two cases; the other two are smlar Frst suppose p 1 A F 1 and p n A F If (x, y 1 ) s a pont of A F n the same column as p 1, then x > x 1, so we can extend the starcase by addng ths pont That contradcts the maxmal length of the starcase So there are no ponts of A F n column y 1 Therefore σ y1 (A) > 0 Smlarly, snce p n A F, there are no ponts of A F 1 n the same column as p n Therefore σ yn (A) < 0 All rows and all columns that contan ponts of the starcase, except columns y 1 and y n, contan exactly two ponts of the starcase, one n A F 1 and one n A F Let A = A\{p 1, p,, p n } Then ρ (A ) = ρ (A) for all, and σ j (A ) = σ j (A) for all j y 1, y n Furthermore, σ y1 (A ) = σ y1 (A) 1 and σ yn (A ) = σ yn (A) + 1 Snce σ y1 (A) > 0 and σ yn (A) < 0, ths gves τ(a ) = τ(a) Now consder the case p 1 A F 1 and p n A F 1 As above, we have σ y1 (A) > 0 Suppose (x n, y) s a pont of A F n the same row as p n Then y > y n, so we can extend the starcase by addng ths pont That contradcts the maxmal length of the starcase So there are no ponts of A F n row x n Therefore ρ xn (A) > 0 All rows and all columns that contan ponts of the starcase, except column y 1 and row x n, contan exactly two ponts of the starcase, one n A F 1 and one n A F Let A = A\{p 1, p,, p n } Then ρ (A ) = ρ (A) for all x n, and σ j (A ) = σ j (A) for all j y 1 Furthermore, σ y1 (A ) = σ y1 (A) 1 and ρ xn (A ) = ρ xn (A) 1 Snce σ y1 (A) > 0 and ρ xn (A) > 0, ths gves τ(a ) = τ(a) We can contnue deletng starcases n ths way untl all ponts of F 1 F have been deleted Snce τ(a) 0 for all subsets A F 1 F, ths must happen after deletng exactly α starcases Remark 5 Some remarks about the above lemma and ts proof 6
7 () The α starcases from the prevous lemma have α endponts n total (where we count the same pont twce n case of a starcase consstng of one pont) Each endpont contrbutes a dfference of 1 to the lne sums n one row or column Snce all these dfferences must add up to α, they cannot cancel each other () A starcase consstng of more than one pont can be splt nto two or more starcases So t may be possble to wrte F 1 F as the dsjont unon of more than α starcases However, n that case some of the contrbutons of the endponts of starcases to the dfference n the lne sums cancel each other On the other hand, t s mpossble to decompose F 1 F nto fewer than α starcases () The endponts of a starcase can be n F 1 \F or F \F 1 For a starcase T of whch the two endponts are n dfferent sets, we have T F 1 = T F For a starcase T of whch the two endponts are n the same set, we have T F 1 = 1 + T F or T F = 1 + T F 1 Snce F 1 \F = F \F 1, the number of starcases wth two endponts n F 1 \F must be equal to the number of starcases wth two endponts n F \F 1 Ths mples that of the α endponts, exactly α are n the set F 1 \F and α are n the set F \F 1 Consder a decomposton of F 1 F as n the proof of Lemma We wll now show that for our purposes we may assume that all these starcases begn wth a pont p 1 F 1 \F and end wth a pont p n F \F 1 Suppose there s a starcase begnnng wth a pont (x, y) F \F 1 Then there also exsts a starcase endng wth a pont (x, y ) F 1 \F : otherwse more than half of the α endponts would be n F \F 1, whch s a contradcton to Remark 5() Because of Remark 5() we must have r (1) x < r () x and r (1) x > r() x Let y be such that (x, y ) F 1 F Delete the pont (x, y) from F and add the pont (x, y ) to F Then r x () decreases by 1 and r () x ncreases by 1, so the dfference n the row sums decreases by Meanwhle, the dfference n the column sums ncreases by at most So α does not ncrease, whle F 1, F and F 1 F do not change So the new stuaton s just as good or better than the old one The starcase that began wth (x, y) n the old stuaton now begns wth a pont of F 1 \F The pont that we added becomes the new endpont of the starcase that prevously ended wth (x, y ) Therefore, n our nvestgatons we may assume that all starcases begn wth a pont of F 1 \F and end wth a pont of F \F 1 Ths s an mportant assumpton that we wll use n the proofs throughout the paper An mmedate consequence of t s that r (1) = r () for all The only dfference between correspondng lne sums occurs n the columns 7
8 A new bound for the dsjont case Usng the concept of starcases, we can prove a new bound for Problem 1 Theorem 6 Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums, F 1 = F, and F 1 F = Let α be defned as n Secton Then α α F 1 PROOF Assume that the rows and columns are ordered as n Secton Let a be the number of rows and b the number of columns that contan elements of F 1 Let (k, l) F 1 Then all the ponts n the rectangle {(, j) : 1 k, 1 j l} are elements of F 1 Snce F 1 and F are dsjont, none of the ponts n ths rectangle s an element of F, and all the ponts belong to F 1 F So all of the kl ponts must belong to dfferent starcases, whch mples α kl For all wth 1 a we have (, r (1) ) F 1, hence r (1) an nteger, we have F 1 = a Snce (a, 1) F 1, we have a α, so r (1) α α F 1 a α α Snce r(1) must be Corollary 7 Let F 1, F and α be defned as n Theorem 6 Then F 1 α(1 + log α) PROOF We have α α F 1 α α ( 1 α ) α x dx = α (1 + log α) 8
9 Fg 3 The constructon from Example 8 wth m = 3 The followng example shows that the upper bound cannot even be mproved by a factor 1 log 07 Example 8 (taken from [1]) Let m 1 be an nteger We construct sets F 1 and F as follows (see also Fgure 3) Row 1: (1, j) F 1 for 1 j m, (1, j) F for m + 1 j m+1 Let 0 l m 1 Row, where l + 1 l+1 : (, j) F 1 for 1 j m l 1, (, j) F for m l j m l The constructon s almost completely symmetrcal: f (, j) F 1, then (j, ) F 1 ; and f (, j) F wth > 1, then (j, ) F Snce t s clear from the constructon that each row contans exactly as many ponts of F 1 as ponts of F, we conclude that each column j wth j m contans exactly as many ponts of F 1 as ponts of F as well The only dfference n the lne sums occurs n the frst column (whch has m ponts of F 1 and none of F ) and n columns m + 1 up to m+1 (each of whch contans one pont of F and none of F 1 ) So we have Furthermore, F 1 = m + m 1 l=0 α = m Hence for ths famly of examples t holds that l m l 1 = m + m m 1 F 1 = α + 1 α log α, whch s very close to the bound we proved n Corollary 7 9
10 5 Two bounds for general α In case F 1 and F are not dsjont, we can use an approach very smlar to Secton n order to derve a bound for Problem Theorem 9 Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums, and F 1 = F Let α be defned as n Secton, and let p = F 1 F Then F 1 α+p α + p PROOF Assume that the rows and columns are ordered as n Secton Let (k, l) F 1 Then all the ponts n the rectangle {(, j) : 1 k, 1 j l} are elements of F 1 At most p of the ponts n ths rectangle are elements of F, so at least kl p ponts belong to F 1 F None of the ponts n the rectangle s an element of F \F 1, so all of the kl p ponts of F 1 F n the rectangle must belong to dfferent starcases, whch mples α + p kl For all wth 1 a we have (, r (1) ) F 1, hence r (1) an nteger, we have F 1 = a r (1) Snce (a, 1) F 1, we have a α + p, so F 1 α+p a α + p α + p α+p Snce r (1) Corollary 10 Let F 1, F, α and p be defned as n Theorem 9 Then F 1 (α + p)(1 + log(α + p)) must be PROOF Analogous to the proof of Corollary 7 The followng example shows that the upper bound cannot even be mproved 1 p+1 by a factor 07, provded that α > log(p + 1) log log 1 Example 11 10
11 Fg The constructon from Example 11 wth k = 3 and m = Let k and m be ntegers satsfyng k and m k We construct sets F 1 and F as follows (see also Fgures and 5) Row 1: (1, j) F 1 F for 1 j k 1, (1, j) F 1 for k j m k 1 + 1, (1, j) F for m k 1 + j m+1 k k 1 + Let 0 l k Row, where l + 1 l+1 : (, 1) F 1 F, (, j) F 1 for j m l 1 k l + 1, (, j) F for m l 1 k l + j m l k l Let k 1 l m k Row, where l + 1 l+1 : (, j) F 1 for 1 j m l 1, (, j) F for m l j m l Let m k + 1 l m 1 Row, where l l m+k 1 + l+1 l m+k + 1: (, j) F 1 for 1 j m l 1, (, j) F for m l j m l The constructon s almost symmetrcal: f (, j) F 1, then (j, ) F 1 ; f (, j) F 1 F, then (j, ) F 1 F ; and f (, j) F wth > 1, then (j, ) F Snce t s clear from the constructon that each row contans exactly as many ponts of F 1 as ponts of F, we conclude that each column j wth j m k contans exactly as many ponts of F 1 as ponts of F as well The only dfference n the lne sums occurs n the frst column (whch has m k ponts of F 1 and only k 1 of F ) and n columns m k 1 + up to m+1 k k 1 + (each of whch contans one pont of F and none of F 1 ) So we have 11
12 Fg 5 The constructon from Example 11 wth k = and m = α = 1 ( ( m k 1 + 1) k 1 + ( m+1 k k 1 + ) ( m k 1 + 1) ) = m k + 1 It s easy to see that p = F 1 F = k 1 Now we count the number of elements of F 1 Row 1 contans m k elements of F 1 Let 0 l k Rows l + 1 up to l+1 together contan l ( m l 1 k l + 1) = m 1 k + l elements of F 1 Let k 1 l m k Rows l + 1 up to l+1 together contan l m l 1 = m 1 elements of F 1 Let m k + 1 l m 1 Rows l l m+k 1 + up to l+1 l m+k + 1 together contan ( l l m+k 1 )( m l 1 ) = m 1 k elements of F 1 Hence the number of elements of F 1 s F 1 = m k (k 1)( m 1 k ) + k l l=0 +(m k + ) m 1 + (k 1)( m 1 k ) = m + m m 1 + k 1 k k 1 For ths famly of examples we now have F 1 = α + p + α + p log (α + p) + p + 1 p + 1 log (p + 1) 1
13 We wll now prove another bound, whch s better f p = F 1 F s large compared to α Let u be an nteger such that u = F 1 F We wll frst derve an upper bound on u n terms of a, b and α Then we wll derve a lower bound on F 1 n terms of a, b and α By combnng these two, we fnd an upper bound on u n terms of α and p Lemma 1 Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums, and F 1 = F Let α, a and b be defned as n Secton Defne u as u = F 1 F Then we have u α (a + b)(a + b + α 1) PROOF Decompose F 1 F nto α starcases as n Lemma, and let T be the set consstng of these starcases Let T T be a starcase and a + 1 a postve nteger Consder the elements of T F n rows, + 1,, a If such elements exst, then let w (T ) be the largest column ndex that occurs among these elements If there are no elements of T F n those rows, then let w (T ) be equal to the smallest column ndex of an element of T F 1 (no longer restrcted to rows,, a) We have w (T ) 1 Defne W = T T w (T ) Let d be the number of elements of F 1 \F n row Let y 1 < < y d be the column ndces of the elements of F 1 \F n row, and let y 1 < < y d be the column ndces of the elements of F \F 1 n row Let T T be the set of starcases wth elements n row The elements n F \F 1 of these starcases are n columns y 1, y,, y d, hence the set {w (T ) : T T } s equal to the set {y 1, y,, y d } The elements n F 1 \F are n columns y 1, y,, y d and are ether the frst element of a starcase or correspond to an element of F \F 1 n the same column but n a row wth ndex at least + 1 In ether case, for a starcase T T we have w +1 (T ) = y j for some j Hence the set {w +1 (T ) : T T } s equal to the set {y 1, y,, y d } We have d d w +1 (T ) = y j (y d j + 1) = d y d 1 T T j=1 j=1 (d 1)d, and Hence d d w (T ) = y j (y d + j) = d y d + 1 T T j=1 j=1 (d + 1)d 13
14 W = W +1 + (w (T ) w +1 (T )) T T W (d + 1)d + 1 (d 1)d = W +1 + d Snce W a+1 α, we fnd W 1 α + d d a We may assume that f (x, y) s the endpont of a starcase, then (x, y ) s an element of F 1 F for 1 y < y (e there are no gaps between the endponts and other elements of F 1 F on the same row) After all, by movng the endpont of a starcase to another empty poston on the same row, the error n the columns can only become smaller (f the new poston of the endpont happens to be n the same column as the frst pont of another starcase, n whch case the two starcases fuse together to one) but not larger, and u, a and b do not change So on the other hand, as W 1 s the sum of the column ndces of the endponts of the starcases, we have W 1 (b + 1) + (b + ) + + (b + α) = αb + 1 α(α + 1) We conclude a α + d αb + 1 α(α + 1) Note that a d = u By the Cauchy-Schwarz nequalty, we have so From ths t follows that or, equvalently, By symmetry we also have ( a ) ( a ) ( a ) d 1 d = u, a d u a αb + 1 u α(α + 1) α + a, u αab + 1 α(α 1)a u αab + 1 α(α 1)b 1
15 Hence u αab + 1 α(α 1)(a + b) Usng that ab a+b, we fnd ( ) (a + b) u (α 1)(a + b) α + = α (a + b)(a + b + α 1) Lemma 13 Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums, and F 1 = F Let α, a and b be defned as n Secton Then we have F 1 (a + b) (α + 1) PROOF Wthout loss of generalty, we may assume that all rows and columns that contan elements of F 1 also contan at least one pont F 1 F : f a row or column does not contan any ponts of F 1 F, we may delete t By dong so, F 1 F does not change, whle F 1 becomes smaller, so the stuaton becomes better Frst consder the case r (1) +1 < r (1) α for some We wll show that ths s mpossble If a column does not contan an element of F \F 1, then by the assumpton above t contans an element of F 1 \F, whch must then be the frst pont of a starcase Consder all ponts of F \F 1 and all frst ponts of starcases n columns r , r +1 +,, r Snce these are more than α columns, at least two of those ponts must belong to the same starcase On the other hand, f (x, y) F 1 \F s the frst pont of a starcase wth r +1 < y r, then we have x, so the second pont (x, y ) n the starcase, whch s n F \F 1, must satsfy x and therefore y > r So the second pont cannot also be n one of the columns r , r +1 +,, r If two ponts of F \F 1 n columns r , r +1 +,, r belong to the same starcase, then they must be connected by a pont of F 1 \F n the same columns However, by a smlar argument ths forces the next pont to be outsde the mentoned columns, whle we assumed that t was n those columns We conclude that t s mpossble for row sums of two consecutve rows to dffer by more than α By the same argument, column sums of two consecutve columns cannot dffer by more than α Hence we have r (1) +1 r (1) α for all, and c (1) j+1 c (1) j α for all j 15
16 Fg 6 The number of ponts of F 1 (ndcated by small black dots) s equal to the grey area We now have r (1) b α, r (1) 3 b α, and so on Also, c (1) a α, c (1) 3 a α, and so on Usng ths, we can derve a lower bound on F 1 for fxed a and b Consder Fgure 6 The ponts of F 1 are ndcated by black dots The number of ponts s equal to the grey area n the pcture, whch conssts of all 1 1-squares wth a pont of F 1 n the upper left corner We can estmate ths area from below by drawng a lne wth slope α through the pont (a + 1, 1) and a lne wth slope 1 through the pont (b + 1, 1); the α area closed n by these two lnes and the two axes s less than or equal to the number of ponts of F 1 For α = 1 those lnes do not have a pont of ntersecton Under the assumpton we made at the begnnng of ths proof, we must n ths case have a = b and the number of ponts of F 1 s equal to a(a + 1) a α + 1 = (a + b) (α + 1), so n ths case we are done In order to compute the area for α we swtch to the usual coordnates n R, see Fgure 7 The equaton of the frst lne s y = αx a, and the equaton of the second lne s y = 1 x 1 b We fnd that the pont of ntersecton s α α gven by ( aα b (x, y) = α 1, bα + a ) α 1 The area of the grey part of Fgure 7 s equal to We now have 1 a aα b α b bα a α 1 = a α + b α ab (α 1) F 1 α(a + b ) ab (α 1) α (a+b) (a+b) (α 1) = (a + b) (α + 1) Theorem 1 Let F 1 and F be fnte subsets of Z such that 16
17 (0, 0) (b, 0) y = 1 α x 1 α b (0, a) y = αx a Fg 7 Computng the area bounded by the two lnes and the two axes F 1 s unquely determned by ts row and column sums, and F 1 = F Let α be defned as n Secton, and let p = F 1 F Wrte β = α(α + 1) Then ( α F 1 p + β + β(α 1) + (α + 1)p + β + α 1 ) (α 1) α 16 PROOF Wrte s = a + b for convenence of notaton From Lemma 1 we derve ( α u s + α 1 ) We substtute F 1 = u + p n Lemma 13 and use the above bound for u: ( α s + α 1 ) s + p F 1 (α + 1) Solvng for s, we fnd s α(α + 1) + α(α 1) + (α + 1)p + α(α + 1) = β + β(α + 1) + (α + 1)p + β Fnally we substtute ths n Lemma 1: u ( α β + β(α 1) + (α + 1)p + β + α 1 Ths, together wth F 1 = u + p, yelds the clamed result ) (α 1) α 16 Remark 15 By a straghtforward generalsaton of [, Proposton 13 and Lemma 16], we fnd a bound very smlar to the one n Theorem 1: F 1 p + (α + 1)(α 1 ) + (α + 1) p + (α 1) 17
18 N α N α Fg 8 The constructon from Example 16 wth N = and α = 3 Theorem 1 says that F 1 s asymptotcally bounded by p + α p + α The next example shows that F 1 can be asymptotcally as large as p + αp + α Example 16 Let N be a postve nteger We construct F 1 and F wth total dfference n the lne sums equal to α as follows (see also Fgure 8) Let (, j) F 1 F for 1 N, 1 j N Furthermore, for 1 N: Let (, j), (j, ) F 1 F for N + 1 j N + (N )α Let (, j), (j, ) F 1 for N + (N )α + 1 j N + (N + 1)α Let (, j), (j, ) F for N + (N + 1)α + 1 j N + (N + )α Fnally, for 1 t α, let (, j) F wth = N + t and j = N + α + 1 t The only dfferences n the lne sums occur n the frst column (a dfference of α) and n columns N + Nα + 1 up to N + Nα + α (a dfference of 1 n each column) We have p = N + 1 N(N 1)α = N + N α Nα, and F 1 = N + 1 N(N + 1)α = N + N α + Nα 18
19 From the frst equalty we derve N = α (α + 1) + p α α (α + 1) Hence F 1 = p + Nα = p + α α α p α α (α + 1) 6 Generalsaton to unequal szes Untl now, we have assumed that F 1 = F However, we can easly generalse all the results to the case F 1 F Suppose F 1 > F Then there must be a row wth r (1) > r () Let j > b be such that (, j) F and defne F 3 = F {(, j)} We have r (3) = r () + 1, so the error n row has decreased by one, whle the error n column j has ncreased by one In ths way, we can keep addng ponts untl F together wth the extra ponts s just as large as F 1, whle the total dfference n the lne sums s stll α Note that p = F 1 F and F 1 have not changed durng ths process, so the results from Theorem 1 and Corollary 10 are stll vald n exactly the same form Suppose on the other hand that F 1 < F Then there must be a row wth r (1) < r () Let j be such that (, j) F \F 1 and defne F 3 = F \{(, j)} The error n row has now decreased by one, whle the error n column j has at most ncreased by one, so the total error n the lne sums has not ncreased We can keep deletng ponts of F untl there are exactly F 1 ponts left, whle the total dfference n the lne sums s at most α By usng F 1 F = ( F 1 p), we can state the results from Theorem 1 and Corollary 10 n a more symmetrc way, not dependng on the sze of F 1 Theorem 17 Let F 1 and F be fnte subsets of Z such that F 1 s unquely determned by ts row and column sums Let α be defned as n Secton, and let p = F 1 F Wrte β = α(α + 1) Then (1) F 1 F α + (α + p) log(α + p) () F 1 F α ( β + β(α 1) + (α + 1)p + β + α 1 ) (α 1) α 19
20 References [1] A Alpers, Instablty and stablty n dscrete tomography, PhD thess, Technsche Unverstät München, ISBN X, Shaker Verlag, Aachen (003) [] A Alpers, S Brunett, Stablty results for the reconstructon of bnary pctures from two projectons, Image and Vson Computng 5 (007) [3] A Alpers, P Grtzmann, L Thorens, Stablty and nstablty n dscrete tomography, Lectures Notes n Computer Scence 3: Dgtal and Image Geometry (001) [] RJ Gardner, P Grtzmann, D Prangenberg, On the computatonal complexty of reconstructng lattce sets from ther X-rays, Dscrete Mathematcs 0 (1999) 5-71 [5] GT Herman, A Kuba, Dscrete tomography n medcal magng, Proceedngs of IEEE 91 (003) [6] GT Herman, A Kuba, edtors, Dscrete Tomography: Foundatons, Algorthms and Applcatons, Brkhäuser, Boston (1999) [7] GT Herman, A Kuba, edtors, Advances n Dscrete Tomography and Its Applcatons, Brkhäuser, Boston (007) [8] JR Jnschek, KJ Batenburg, HA Calderon, R Klaas, V Radmlovc and C Kselowsk, 3-D reconstructon of the atomc postons n a smulated gold nanocrystal based on dscrete tomography, Ultramcroscopy (007), do:101016/jultramc [9] A Kuba, L Rodek, Z Kss, L Ruskó, A Nagy, M Balaskó, Dscrete tomography n neutron radography, Nucl Instrum Methods Phys Res, Sect A 5 (005) [10] JC Palacos, LC Longora, J Santos, RT Perry, A PC-based dscrete tomography magng software system for assayng radoactve waste contaners, Nucl Instrum Methods Phys Res, Sect A 508 (003) [11] HJ Ryser, Combnatoral propertes of matrces of zeros and ones, Can J Math 9 (1957) [1] B Schllnger, Proposed combnaton of CAD data and dscrete tomography for the detecton of cokng and lubrcants n turbne blades or engnes, Electronc Notes n Dscrete Mathematcs 0 (005) [13] H Slump, JJ Gerbrands, A network flow approach to reconstructon of the left ventrcle from two projectons, Computer graphcs and mage processng 18 (198) [1] YR Wang, Characterzaton of bnary patterns and ther projectons, IEEE Trans Comput (1975)
21 [15] L Wang, JY Park, Y Fu, Representaton of real partcles for DEM smulaton usng X-ray tomography, Constructon and Buldng Materals 1 (007)
Affine transformations and convexity
Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationAppendix B. Criterion of Riemann-Stieltjes Integrability
Appendx B. Crteron of Remann-Steltes Integrablty Ths note s complementary to [R, Ch. 6] and [T, Sec. 3.5]. The man result of ths note s Theorem B.3, whch provdes the necessary and suffcent condtons for
More informationarxiv: v1 [math.co] 1 Mar 2014
Unon-ntersectng set systems Gyula O.H. Katona and Dánel T. Nagy March 4, 014 arxv:1403.0088v1 [math.co] 1 Mar 014 Abstract Three ntersecton theorems are proved. Frst, we determne the sze of the largest
More informationThe lower and upper bounds on Perron root of nonnegative irreducible matrices
Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College
More informationFACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP
C O L L O Q U I U M M A T H E M A T I C U M VOL. 80 1999 NO. 1 FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP BY FLORIAN K A I N R A T H (GRAZ) Abstract. Let H be a Krull monod wth nfnte class
More informationFinding Primitive Roots Pseudo-Deterministically
Electronc Colloquum on Computatonal Complexty, Report No 207 (205) Fndng Prmtve Roots Pseudo-Determnstcally Ofer Grossman December 22, 205 Abstract Pseudo-determnstc algorthms are randomzed search algorthms
More informationProblem Set 9 Solutions
Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem
More informationLinear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space.
Lnear, affne, and convex sets and hulls In the sequel, unless otherwse specfed, X wll denote a real vector space. Lnes and segments. Gven two ponts x, y X, we defne xy = {x + t(y x) : t R} = {(1 t)x +
More informationDiscrete tomography with two directions
Discrete tomography with two directions Proefschrift ter verkrijging van de graad van Doctor aan de Universiteit Leiden, op gezag van Rector Magnificus prof.mr. P.F. van der Heijden, volgens besluit van
More information2.3 Nilpotent endomorphisms
s a block dagonal matrx, wth A Mat dm U (C) In fact, we can assume that B = B 1 B k, wth B an ordered bass of U, and that A = [f U ] B, where f U : U U s the restrcton of f to U 40 23 Nlpotent endomorphsms
More informationComplete subgraphs in multipartite graphs
Complete subgraphs n multpartte graphs FLORIAN PFENDER Unverstät Rostock, Insttut für Mathematk D-18057 Rostock, Germany Floran.Pfender@un-rostock.de Abstract Turán s Theorem states that every graph G
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationFirst day August 1, Problems and Solutions
FOURTH INTERNATIONAL COMPETITION FOR UNIVERSITY STUDENTS IN MATHEMATICS July 30 August 4, 997, Plovdv, BULGARIA Frst day August, 997 Problems and Solutons Problem. Let {ε n } n= be a sequence of postve
More informationSelf-complementing permutations of k-uniform hypergraphs
Dscrete Mathematcs Theoretcal Computer Scence DMTCS vol. 11:1, 2009, 117 124 Self-complementng permutatons of k-unform hypergraphs Artur Szymańsk A. Paweł Wojda Faculty of Appled Mathematcs, AGH Unversty
More information20. Mon, Oct. 13 What we have done so far corresponds roughly to Chapters 2 & 3 of Lee. Now we turn to Chapter 4. The first idea is connectedness.
20. Mon, Oct. 13 What we have done so far corresponds roughly to Chapters 2 & 3 of Lee. Now we turn to Chapter 4. The frst dea s connectedness. Essentally, we want to say that a space cannot be decomposed
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationRandić Energy and Randić Estrada Index of a Graph
EUROPEAN JOURNAL OF PURE AND APPLIED MATHEMATICS Vol. 5, No., 202, 88-96 ISSN 307-5543 www.ejpam.com SPECIAL ISSUE FOR THE INTERNATIONAL CONFERENCE ON APPLIED ANALYSIS AND ALGEBRA 29 JUNE -02JULY 20, ISTANBUL
More informationAnti-van der Waerden numbers of 3-term arithmetic progressions.
Ant-van der Waerden numbers of 3-term arthmetc progressons. Zhanar Berkkyzy, Alex Schulte, and Mchael Young Aprl 24, 2016 Abstract The ant-van der Waerden number, denoted by aw([n], k), s the smallest
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
More informationGraph Reconstruction by Permutations
Graph Reconstructon by Permutatons Perre Ille and Wllam Kocay* Insttut de Mathémathques de Lumny CNRS UMR 6206 163 avenue de Lumny, Case 907 13288 Marselle Cedex 9, France e-mal: lle@ml.unv-mrs.fr Computer
More informationCOMPLEX NUMBERS AND QUADRATIC EQUATIONS
COMPLEX NUMBERS AND QUADRATIC EQUATIONS INTRODUCTION We know that x 0 for all x R e the square of a real number (whether postve, negatve or ero) s non-negatve Hence the equatons x, x, x + 7 0 etc are not
More informationA new construction of 3-separable matrices via an improved decoding of Macula s construction
Dscrete Optmzaton 5 008 700 704 Contents lsts avalable at ScenceDrect Dscrete Optmzaton journal homepage: wwwelsevercom/locate/dsopt A new constructon of 3-separable matrces va an mproved decodng of Macula
More informationLecture 5 Decoding Binary BCH Codes
Lecture 5 Decodng Bnary BCH Codes In ths class, we wll ntroduce dfferent methods for decodng BCH codes 51 Decodng the [15, 7, 5] 2 -BCH Code Consder the [15, 7, 5] 2 -code C we ntroduced n the last lecture
More informationA CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS
Journal of Mathematcal Scences: Advances and Applcatons Volume 25, 2014, Pages 1-12 A CHARACTERIZATION OF ADDITIVE DERIVATIONS ON VON NEUMANN ALGEBRAS JIA JI, WEN ZHANG and XIAOFEI QI Department of Mathematcs
More informationBOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS
BOUNDEDNESS OF THE IESZ TANSFOM WITH MATIX A WEIGHTS Introducton Let L = L ( n, be the functon space wth norm (ˆ f L = f(x C dx d < For a d d matrx valued functon W : wth W (x postve sem-defnte for all
More informationVolume 18 Figure 1. Notation 1. Notation 2. Observation 1. Remark 1. Remark 2. Remark 3. Remark 4. Remark 5. Remark 6. Theorem A [2]. Theorem B [2].
Bulletn of Mathematcal Scences and Applcatons Submtted: 016-04-07 ISSN: 78-9634, Vol. 18, pp 1-10 Revsed: 016-09-08 do:10.1805/www.scpress.com/bmsa.18.1 Accepted: 016-10-13 017 ScPress Ltd., Swtzerland
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationFoundations of Arithmetic
Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an
More informationLecture 20: Lift and Project, SDP Duality. Today we will study the Lift and Project method. Then we will prove the SDP duality theorem.
prnceton u. sp 02 cos 598B: algorthms and complexty Lecture 20: Lft and Project, SDP Dualty Lecturer: Sanjeev Arora Scrbe:Yury Makarychev Today we wll study the Lft and Project method. Then we wll prove
More informationEdge Isoperimetric Inequalities
November 7, 2005 Ross M. Rchardson Edge Isopermetrc Inequaltes 1 Four Questons Recall that n the last lecture we looked at the problem of sopermetrc nequaltes n the hypercube, Q n. Our noton of boundary
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More information10-801: Advanced Optimization and Randomized Methods Lecture 2: Convex functions (Jan 15, 2014)
0-80: Advanced Optmzaton and Randomzed Methods Lecture : Convex functons (Jan 5, 04) Lecturer: Suvrt Sra Addr: Carnege Mellon Unversty, Sprng 04 Scrbes: Avnava Dubey, Ahmed Hefny Dsclamer: These notes
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
More informationOn the size of quotient of two subsets of positive integers.
arxv:1706.04101v1 [math.nt] 13 Jun 2017 On the sze of quotent of two subsets of postve ntegers. Yur Shtenkov Abstract We obtan non-trval lower bound for the set A/A, where A s a subset of the nterval [1,
More informationTheoretical Computer Science. Reconstruction of convex lattice sets from tomographic projections in quartic time
Theoretcal Computer Scence 406 (2008) 55 62 Contents lsts avalable at ScenceDrect Theoretcal Computer Scence journal homepage: www.elsever.com/locate/tcs Reconstructon of convex lattce sets from tomographc
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationLecture 3. Ax x i a i. i i
18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest
More information= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.
Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and
More informationMin Cut, Fast Cut, Polynomial Identities
Randomzed Algorthms, Summer 016 Mn Cut, Fast Cut, Polynomal Identtes Instructor: Thomas Kesselhem and Kurt Mehlhorn 1 Mn Cuts n Graphs Lecture (5 pages) Throughout ths secton, G = (V, E) s a mult-graph.
More informationAn efficient algorithm for multivariate Maclaurin Newton transformation
Annales UMCS Informatca AI VIII, 2 2008) 5 14 DOI: 10.2478/v10065-008-0020-6 An effcent algorthm for multvarate Maclaurn Newton transformaton Joanna Kapusta Insttute of Mathematcs and Computer Scence,
More information5 The Rational Canonical Form
5 The Ratonal Canoncal Form Here p s a monc rreducble factor of the mnmum polynomal m T and s not necessarly of degree one Let F p denote the feld constructed earler n the course, consstng of all matrces
More informationAssortment Optimization under MNL
Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.
More informationMaximizing the number of nonnegative subsets
Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum
More informationA 2D Bounded Linear Program (H,c) 2D Linear Programming
A 2D Bounded Lnear Program (H,c) h 3 v h 8 h 5 c h 4 h h 6 h 7 h 2 2D Lnear Programmng C s a polygonal regon, the ntersecton of n halfplanes. (H, c) s nfeasble, as C s empty. Feasble regon C s unbounded
More informationSubset Topological Spaces and Kakutani s Theorem
MOD Natural Neutrosophc Subset Topologcal Spaces and Kakutan s Theorem W. B. Vasantha Kandasamy lanthenral K Florentn Smarandache 1 Copyrght 1 by EuropaNova ASBL and the Authors Ths book can be ordered
More informationMATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS
MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS These are nformal notes whch cover some of the materal whch s not n the course book. The man purpose s to gve a number of nontrval examples
More informationWeek 2. This week, we covered operations on sets and cardinality.
Week 2 Ths week, we covered operatons on sets and cardnalty. Defnton 0.1 (Correspondence). A correspondence between two sets A and B s a set S contaned n A B = {(a, b) a A, b B}. A correspondence from
More informationLecture 2: Gram-Schmidt Vectors and the LLL Algorithm
NYU, Fall 2016 Lattces Mn Course Lecture 2: Gram-Schmdt Vectors and the LLL Algorthm Lecturer: Noah Stephens-Davdowtz 2.1 The Shortest Vector Problem In our last lecture, we consdered short solutons to
More informationSupplement: Proofs and Technical Details for The Solution Path of the Generalized Lasso
Supplement: Proofs and Techncal Detals for The Soluton Path of the Generalzed Lasso Ryan J. Tbshran Jonathan Taylor In ths document we gve supplementary detals to the paper The Soluton Path of the Generalzed
More informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
More informationTHE WEIGHTED WEAK TYPE INEQUALITY FOR THE STRONG MAXIMAL FUNCTION
THE WEIGHTED WEAK TYPE INEQUALITY FO THE STONG MAXIMAL FUNCTION THEMIS MITSIS Abstract. We prove the natural Fefferman-Sten weak type nequalty for the strong maxmal functon n the plane, under the assumpton
More informationFINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN
FINITELY-GENERTED MODULES OVER PRINCIPL IDEL DOMIN EMMNUEL KOWLSKI Throughout ths note, s a prncpal deal doman. We recall the classfcaton theorem: Theorem 1. Let M be a fntely-generated -module. (1) There
More informationFormulas for the Determinant
page 224 224 CHAPTER 3 Determnants e t te t e 2t 38 A = e t 2te t e 2t e t te t 2e 2t 39 If 123 A = 345, 456 compute the matrx product A adj(a) What can you conclude about det(a)? For Problems 40 43, use
More informationON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES. 1. Introduction
ON THE EXTENDED HAAGERUP TENSOR PRODUCT IN OPERATOR SPACES TAKASHI ITOH AND MASARU NAGISA Abstract We descrbe the Haagerup tensor product l h l and the extended Haagerup tensor product l eh l n terms of
More informationCSCE 790S Background Results
CSCE 790S Background Results Stephen A. Fenner September 8, 011 Abstract These results are background to the course CSCE 790S/CSCE 790B, Quantum Computaton and Informaton (Sprng 007 and Fall 011). Each
More informationa b a In case b 0, a being divisible by b is the same as to say that
Secton 6.2 Dvsblty among the ntegers An nteger a ε s dvsble by b ε f there s an nteger c ε such that a = bc. Note that s dvsble by any nteger b, snce = b. On the other hand, a s dvsble by only f a = :
More informationA combinatorial problem associated with nonograms
A combnatoral problem assocated wth nonograms Jessca Benton Ron Snow Nolan Wallach March 21, 2005 1 Introducton. Ths work was motvated by a queston posed by the second named author to the frst named author
More informationTHE CHVÁTAL-ERDŐS CONDITION AND 2-FACTORS WITH A SPECIFIED NUMBER OF COMPONENTS
Dscussones Mathematcae Graph Theory 27 (2007) 401 407 THE CHVÁTAL-ERDŐS CONDITION AND 2-FACTORS WITH A SPECIFIED NUMBER OF COMPONENTS Guantao Chen Department of Mathematcs and Statstcs Georga State Unversty,
More informationSalmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to
More informationExercise Solutions to Real Analysis
xercse Solutons to Real Analyss Note: References refer to H. L. Royden, Real Analyss xersze 1. Gven any set A any ɛ > 0, there s an open set O such that A O m O m A + ɛ. Soluton 1. If m A =, then there
More informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationSTEINHAUS PROPERTY IN BANACH LATTICES
DEPARTMENT OF MATHEMATICS TECHNICAL REPORT STEINHAUS PROPERTY IN BANACH LATTICES DAMIAN KUBIAK AND DAVID TIDWELL SPRING 2015 No. 2015-1 TENNESSEE TECHNOLOGICAL UNIVERSITY Cookevlle, TN 38505 STEINHAUS
More informationNP-Completeness : Proofs
NP-Completeness : Proofs Proof Methods A method to show a decson problem Π NP-complete s as follows. (1) Show Π NP. (2) Choose an NP-complete problem Π. (3) Show Π Π. A method to show an optmzaton problem
More informationProblem Do any of the following determine homomorphisms from GL n (C) to GL n (C)?
Homework 8 solutons. Problem 16.1. Whch of the followng defne homomomorphsms from C\{0} to C\{0}? Answer. a) f 1 : z z Yes, f 1 s a homomorphsm. We have that z s the complex conjugate of z. If z 1,z 2
More informationModule 9. Lecture 6. Duality in Assignment Problems
Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept
More informationPerron Vectors of an Irreducible Nonnegative Interval Matrix
Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of
More informationOn a direct solver for linear least squares problems
ISSN 2066-6594 Ann. Acad. Rom. Sc. Ser. Math. Appl. Vol. 8, No. 2/2016 On a drect solver for lnear least squares problems Constantn Popa Abstract The Null Space (NS) algorthm s a drect solver for lnear
More informationModule 2. Random Processes. Version 2 ECE IIT, Kharagpur
Module Random Processes Lesson 6 Functons of Random Varables After readng ths lesson, ou wll learn about cdf of functon of a random varable. Formula for determnng the pdf of a random varable. Let, X be
More information8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS
SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars
More informationYong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )
Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More informationTHE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens
THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of
More informationThe Second Eigenvalue of Planar Graphs
Spectral Graph Theory Lecture 20 The Second Egenvalue of Planar Graphs Danel A. Spelman November 11, 2015 Dsclamer These notes are not necessarly an accurate representaton of what happened n class. The
More informationMAT 578 Functional Analysis
MAT 578 Functonal Analyss John Qugg Fall 2008 Locally convex spaces revsed September 6, 2008 Ths secton establshes the fundamental propertes of locally convex spaces. Acknowledgment: although I wrote these
More informationEXPANSIVE MAPPINGS. by W. R. Utz
Volume 3, 978 Pages 6 http://topology.auburn.edu/tp/ EXPANSIVE MAPPINGS by W. R. Utz Topology Proceedngs Web: http://topology.auburn.edu/tp/ Mal: Topology Proceedngs Department of Mathematcs & Statstcs
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationRepresentation theory and quantum mechanics tutorial Representation theory and quantum conservation laws
Representaton theory and quantum mechancs tutoral Representaton theory and quantum conservaton laws Justn Campbell August 1, 2017 1 Generaltes on representaton theory 1.1 Let G GL m (R) be a real algebrac
More informationEuropean Journal of Combinatorics
European Journal of Combnatorcs 0 (009) 480 489 Contents lsts avalable at ScenceDrect European Journal of Combnatorcs journal homepage: www.elsever.com/locate/ejc Tlngs n Lee metrc P. Horak 1 Unversty
More informationExample: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,
The greatest common dvsor of two ntegers a and b (not both zero) s the largest nteger whch s a common factor of both a and b. We denote ths number by gcd(a, b), or smply (a, b) when there s no confuson
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationProblem Solving in Math (Math 43900) Fall 2013
Problem Solvng n Math (Math 43900) Fall 2013 Week four (September 17) solutons Instructor: Davd Galvn 1. Let a and b be two nteger for whch a b s dvsble by 3. Prove that a 3 b 3 s dvsble by 9. Soluton:
More informationk(k 1)(k 2)(p 2) 6(p d.
BLOCK-TRANSITIVE 3-DESIGNS WITH AFFINE AUTOMORPHISM GROUP Greg Gamble Let X = (Z p d where p s an odd prme and d N, and let B X, B = k. Then t was shown by Praeger that the set B = {B g g AGL d (p} s the
More informationRandom Partitions of Samples
Random Parttons of Samples Klaus Th. Hess Insttut für Mathematsche Stochastk Technsche Unverstät Dresden Abstract In the present paper we construct a decomposton of a sample nto a fnte number of subsamples
More informationSpectral Graph Theory and its Applications September 16, Lecture 5
Spectral Graph Theory and ts Applcatons September 16, 2004 Lecturer: Danel A. Spelman Lecture 5 5.1 Introducton In ths lecture, we wll prove the followng theorem: Theorem 5.1.1. Let G be a planar graph
More informationGeneral viscosity iterative method for a sequence of quasi-nonexpansive mappings
Avalable onlne at www.tjnsa.com J. Nonlnear Sc. Appl. 9 (2016), 5672 5682 Research Artcle General vscosty teratve method for a sequence of quas-nonexpansve mappngs Cuje Zhang, Ynan Wang College of Scence,
More informationGenericity of Critical Types
Genercty of Crtcal Types Y-Chun Chen Alfredo D Tllo Eduardo Fangold Syang Xong September 2008 Abstract Ely and Pesk 2008 offers an nsghtful characterzaton of crtcal types: a type s crtcal f and only f
More informationREGULAR POSITIVE TERNARY QUADRATIC FORMS. 1. Introduction
REGULAR POSITIVE TERNARY QUADRATIC FORMS BYEONG-KWEON OH Abstract. A postve defnte quadratc form f s sad to be regular f t globally represents all ntegers that are represented by the genus of f. In 997
More informationThe path of ants Dragos Crisan, Andrei Petridean, 11 th grade. Colegiul National "Emil Racovita", Cluj-Napoca
Ths artcle s wrtten by students. It may nclude omssons and mperfectons, whch were dentfed and reported as mnutely as possble by our revewers n the edtoral notes. The path of ants 07-08 Students names and
More informationHMMT February 2016 February 20, 2016
HMMT February 016 February 0, 016 Combnatorcs 1. For postve ntegers n, let S n be the set of ntegers x such that n dstnct lnes, no three concurrent, can dvde a plane nto x regons (for example, S = {3,
More information12 MATH 101A: ALGEBRA I, PART C: MULTILINEAR ALGEBRA. 4. Tensor product
12 MATH 101A: ALGEBRA I, PART C: MULTILINEAR ALGEBRA Here s an outlne of what I dd: (1) categorcal defnton (2) constructon (3) lst of basc propertes (4) dstrbutve property (5) rght exactness (6) localzaton
More informationPolynomials. 1 More properties of polynomials
Polynomals 1 More propertes of polynomals Recall that, for R a commutatve rng wth unty (as wth all rngs n ths course unless otherwse noted), we defne R[x] to be the set of expressons n =0 a x, where a
More informationCalculation of time complexity (3%)
Problem 1. (30%) Calculaton of tme complexty (3%) Gven n ctes, usng exhaust search to see every result takes O(n!). Calculaton of tme needed to solve the problem (2%) 40 ctes:40! dfferent tours 40 add
More informationLine Drawing and Clipping Week 1, Lecture 2
CS 43 Computer Graphcs I Lne Drawng and Clppng Week, Lecture 2 Davd Breen, Wllam Regl and Maxm Peysakhov Geometrc and Intellgent Computng Laboratory Department of Computer Scence Drexel Unversty http://gcl.mcs.drexel.edu
More informationModule 14: THE INTEGRAL Exploring Calculus
Module 14: THE INTEGRAL Explorng Calculus Part I Approxmatons and the Defnte Integral It was known n the 1600s before the calculus was developed that the area of an rregularly shaped regon could be approxmated
More information