May 2005 - mathfun@mailman.xmission.com

[math-fun] Re: Packing circles into a finite rectangle
by James Propp 10 Feb '06

10 Feb '06

>> If you have circles of radii 1/1, 1/2, 1/3, ..., 1/n for all n >= 1, >> then their total area is finite, although the sum of their radii is >> infinite. >> >> Can one fit these circles into a finite-sized rectangle? > > yes, they will fit inside a 1 x (11/6) rectangle for example. put each > circle inside a square, and then pack the squares like this: > http://www.stetson.edu/~efriedma/harmonic/ It is apparently still unknown whether one can pack squares of sides-lengths 1/1, 1/2, 1/3, etc. into a 1-by-(pi squared over six) rectangle. See Clive Tooth's web-page www.pisquaredoversix.force9.co.uk/Tiling.htm for an example of how one might get started. I asked Clive for background on this tiling, and he kindly forwarded to me copies of some relevant sci.math postings from last year (one by Clive himself and one by David Cantrell): ************************************************************************** From: "The Last Danish Pastry" <TheLastDanishPas...(a)yahoo.com> Newsgroups: sci.math Subject: Packing squares into a 1 x Zeta[2] rectangle Date: Sun, 14 Mar 2004 13:41:34 -0000 Since the sum of the reciprocals of the squares of the positive integers is pi^2/6, the question arises as to whether squares with sides 1, 1/2, 1/3, etc can be packed into a rectangle of size 1 by pi^2/6. A picture of such a packing appears at http://www.pisquaredoversix.force9.co.uk/Tiling.htm I know of no proof that such a packing is possible or impossible. A program which I wrote has recently packed the first one million such squares into such a rectangle. I wrote the program in Mathematica and it uses exact arithmetic throughout. The algorithm is simple: The program maintains a list of rectangles. Each member of this list is itself a list containing exactly two members: the short side of a rectangle followed by the long side of the rectangle. The program does not store the position or orientation of any rectangle. The list is stored in ascending short side order. Initially the list contains one entry: {1, pi^2/6}. The program inserts the squares in decreasing size order, starting with the 1x1 square. At each stage, to insert the square with side 1/n, the program finds the rectangle with the shortest short side which will accommodate the square. [In the case that there are two or more such rectangles with equal short sides, the program will pick the youngest one.] Suppose that this rectangle is a x b, where 1/n<=a<b. The program deletes this rectangle from the list and inserts two new rectangles (b-1/n) x a and (a-1/n) x 1/n into the list. - If any side of one of these new rectangles is zero, the rectangle is not inserted. - Each new rectangle is stored in {short side, long side} form. Thus, the list of rectangles grows by one per inserted square, except in the rare cases in which a rectangle of exactly the right width is found. A trivial optimisation, which took me a long time to notice, is that if the number of squares to be inserted is decided at the beginning, as N, say, then any rectangles narrower than 1/N can be discarded as soon as they appear. This has a dramatic effect on the memory requirements and running time of the program. In the case of N=10^6 the maximum length of the rectangle list is only 7,493 which occurs just after inserting the square with side 1/55,205. Thereafter the length of the list slowly decreases. After inserting the last square, with side 10^-6, there are just 7 rectangles in the list. Their dimensions are, roughly: 0.000001000264 x 0.000371508565 [1,000,000] 0.000001000289 x 0.000006561855 [ 999,736] 0.000001000332 x 0.000011950546 [ 999,710] 0.000001000424 x 0.000099911008 [ 999,668] 0.000001001705 x 0.000013366392 [ 999,575] 0.000001001771 x 0.000555285230 [ 998,297] 0.000621399972 x 0.000621402936 [ 998,232] The number in square brackets is the square insertion step on which the rectangle was created. Note that each rectangle, except the last, has a short side just exceeding 10^-6. The last one is huge in comparison to the others. Also, it is almost square, having an aspect ratio of about 1.00000477 . The exact size of that last rectangle may be expressed as: (pi^2/6 - ShortNumerator/ShortDenominator) by (LongNumerator/LongDenominator) Where ShortNumerator = 88148577597184641314023698687804198322622003067950 90277185502203375757868045576236778623499141581835 91095283100445795053089253938670874859361703421694 93589311823604254188777249571387340885945805010588 59158565482204748239750452173829861295626045336259 37611823762149119593932203343599960000593198902456 62151237628335467123582313067030461952701063806541 09667253352708582423851542412389977650553557102776 48867766805593783542952432276691372915179943637567 30921614687225494242879103415386639772227105415075 60643644416702161056703911777435451470868967261015 03285333687384626370723064421360623232628962182637 33785897530162588640506103091960735379749803817168 27337105343328553249474140253583305664128829433056 13496001065570446619534699711595669897241621797339 93388113729020710706312912651983440135427698863738 63455691221381048718834049233446561869585869886892 11066449759482338440176987851999633382542943944743 83052271691048660471857179765488027759900589960305 37378583137068870187987573969424873832345344118539 33368635075465276924658419905001027235726511505469 90353609271920928926254763071708700239940565565014 69168650106503612867362486635886531701170797485510 99375741498187773039269448405524049571692907463933 52292512770766034352574830725344567781634199947344 39128864912738319438009565992988317045941602472518 84296172734535520954249457165155850357846464984168 55299222999724423825677947429849212042460805178788 23214473791347512926915739195451188103089711423995 69891250713418835019305494967439261535807066257765 86799633983386818438909781178910531233498247023179 22489232680161587870238609995903489290522620579734 27913072409961919211050524749783993630732833119534 51221464170732847284504244076797981093088600479209 053071348915420741889 ShortDenominator = 53608160645434703140761211994317971700783074368336 23325975705041657962199994733358382994458318582660 77598233503372093023408990096331259105144227503324 99931830157705570560424971647310855286927378535483 51519419895212616165032657829920140448566211296935 08925318960134998789890442610040864611429495599467 55671385275409484022080990210675967224987014208598 03335236373663242927780013944648864149878730486125 39041613975393755631755935261748017063478349180906 68435262002209072815068783227098957825899390123056 87116154057047656878878907753207321889684985791605 16008019019744195226005265347943729214648741668979 90732236259531674523664568538104488968273268732774 74645971065348455713992347183400863276778739619943 89625588842652615699820699583835957017913827580111 39619239553717364061586421301662144411944880262864 13395969932087668263688142360819354641427777701996 11575083489100468980371536759296335301136691198138 39559100156995546068054947710363960181181289364408 47219039882502818687162480694856587042340603392188 93192097190331493671006154660882637398429894947140 66724834302759980279336075917538263518035481261929 09959571761983813160895335069812335369906345916646 57307263315174319702690110108769119880417800559541 00846379467033278426493445187132240090150811316237 11870955392306082427560405688265993364931620972062 42335013731582665984021398103261291575096750373611 78793366106368150944101270721191063173365978409196 32419016511795819355629216137822167151147383280072 17692953976640188298210375720290857531988737370969 01207985676688854219832216404101906724620701738590 65777947161752685898547108814129100143613045915161 23374032620084188029394058765229087587568694264890 44233796349785411756140522085279648064146173069741 478636801148014080000 LongNumerator = 86047691660315035868968439680618750650385496662337 69341618423224868557176707904056350747739644607152 19341301092244142985694427567208621603825424715808 50112490209730250303812592410205852634552809617406 54471220118615832220663814127702644186714156802514 89272173167125109946591599082250813475432136844437 68434150604827750301142594110541508119955806775639 48463821703879298880841441517861092382325647927555 17115984413353561811152746794319549286695056542231 75370181749913246697604501180515750468588947850814 96706860842334628047016110265861987131951269186913 55908484624794372227266755207694315503681976948427 59732973255837108405693520922709400140190847636749 60091554508265275883766083172672994896381234279724 35485615485280999872894822261731570552709443192476 07685741096374165452217373319152585412622236248087 75856424918727148070254116048744362672852782317829 99207549345144244124571110563211160343448412194752 79110755652853486307738886342254777707077572587992 29405091182393418319601369763056120887765700618802 27139753549245993550447831136020976532760818138519 98216552674127086312240961717818956038535226187300 43026597993187644795910836481706823332909020130603 09729716420449811214331443400182398534278484763571 27904525989212662382850630094162630962359551947273 97311922973034069310043168935737831265241199089047 46559558404697050977201086017144902884602248144108 71005578441670661613107026657746377160408235952122 70087369832576028781964223879066392055590320740847 64872057694761838719127407266757016871496087930551 58594303661605591550655979493954945142285417155556 97061447617516531128930865552526142619566789459050 16802730874016347300046435875463493515247956092289 80788952041557673741641666771939979618817508420910 11050430185421809507580791626321726210797393934779 33362836600904869583181 LongDenominator = 13847326202197605511604348541676281864027795818043 31494335574789145450920249049910013265993162156505 60831474662540178414423726822350848233301461067293 21368995360010230519749558542277057071473577106811 23548536664652068929533780687547853336157369555086 86435599910456774645119837370840414802536199558214 07527909448749766750406409444125057733978666254122 19353807926168688488715863115960118447262678159252 66974533884810419996115934738184899912380787401311 77856867483197701324817440576351765023155329523060 68640715686158929369832278778842403496825357945180 51999476544062632945516973292520225337991834072383 82846228280041100692245330583478338241014202245803 30323143236417275386997321382053499979484896390576 45426647084624716067636294896758072966317096343334 28169063530634898773783098553798071542939121086910 08006721058240228461124891602875033925435245224172 90909661066024272174911721070106393261965722571136 18608587173321485889447514408382277645342459166650 66976600733716756193252940889415454560334136529577 70915746535414564177067690689656587269581203051464 54295311030944948219808305598301111195229732416082 42817235468158091436079088616794858652804473035517 82817817155337872004025893330437425894039012377193 76447539452009335637045807403996380650616246802471 05691575600545667349590304669247493320414626833379 52612693846867728391114410877249463229225980473996 60123182871902607011087292407939606651288648107916 20645836880002779839862309998364888537462475216562 83515669632032306317521726424095159019732485843596 26389117354057755490411586373128011034072899097872 57238596906896860163242623358201011919041838989031 36072177994374094882330699298066695784371692394667 11940033790850742395135563324255474227834715815498 61014376098413906103604397447388393432391565219123 065068758213935671227536000 -- Clive Tooth http://www.clivetooth.dk ************************************************************************** Subject: Re: Packing squares into a 1 x Zeta[2] rectangle From: David W. Cantrell <DWCantr...(a)sigmaxi.org> Date: 18 Mar 2004 16:31:52 GMT Newsgroups: sci.math "The Last Danish Pastry" <TheLastDanishPas...(a)yahoo.com> wrote: > Since the sum of the reciprocals of the squares of the positive integers > is pi^2/6, the question arises as to whether squares with sides 1, 1/2, > 1/3, etc can be packed into a rectangle of size 1 by pi^2/6. > > A picture of such a packing appears at > http://www.pisquaredoversix.force9.co.uk/Tiling.htm > > I know of no proof that such a packing is possible or impossible. [snip] It would be good if people interested in this thread were aware of some of the history of this problem. To that end: _Unsolved Problems in Geometry_ by Croft, Falconer and Guy (Springer-Verlag, 1991) "D5. Packing unequal rectangles and squares in a square." (pp. 112-113) mentions a question similar to the one considered in this thread. D5 caused me, in 1996, to consider several packing problems of this sort. I found simple algorithms, which I implemented in Mathematica, which seemed to pack rectangles of decreasing size "perfectly". The algorithms worked for classic problems, such as packing the rectangles 1/n by 1/(n+1) into the unit square, as well as for interesting new packing problems. I did not publish my results then, hoping first to prove that the algorithms were adequate for packing all of the rectangles. In November 1999, spurred by Clive Tooth, I posted to sci.math a brief description of what I had done. See both parts of "Packing Rectangles of Decreasing Size" at <http://mathforum.org/discuss/sci.math/a/t/241167>. (Readers may also be interested in looking at the packing algorithm which Clive had used back then.) Unbeknownst to me in 1999, an algorithm for such packings had appeared previously in the literature: An Algorithm for Packing Squares, Marc M. Paulhus, _J. Combin. Theory Ser. A_ 82 (1997) 147-157. His algorithm is somewhat different from any of the algorithms which I had used (but I wonder if perhaps it is the same as Clive's current algorithm). For packing squares of side lengths 1/2, 1/3, 1/4,... into a rectangle 1/2 by 2(pi^2/6 - 1), Paulhus gives two rules: "Rule 1. Always place the next square in the corner of the smallest width rectangle into which it will fit." "Rule 2. After placing a square into the corner of a rectangle, always cut the remaining area into two rectangular pieces by cutting from the free corner of the square tot the longer side of the original rectangle." For older references, see those cited by Paulhus. The most recent articles known to me concerning such packings -- but again, I may not be up-to-date with the literature -- are the following: Perfect Square Packings, Adam Chalcraft, _J. Combin. Theory Ser. A_ 92 (2000) 158-172. Compactness Theorems for Geometric Packings, Greg Martin, <http://arxiv.org/PS_cache/math/pdf/0005/0005054.pdf>. Perfect Packings of Squares Using the Stack-Pack Strategy, Johan Wastlund, <www.mai.liu.se/~jowas/dcg.ps <http://www.mai.liu.se/~jowas/dcg.ps>>. Thoughtful comments are welcome. Regards, David W. Cantrell ************************************************************************** (Clive tells me that there's more discussion of this problem in sci.math, though I haven't checked.) Jim Propp

4 3

[math-fun] Found in the IHP cellar! 2/2
by Christian Boyer 21 Nov '05

21 Nov '05

---- FAUQUEMBERGUE computations And I have the chance to have discovered, again in the cellar of the Institut Henri Poincaré, the original handwritten computations on Mersenne numbers M101, M103, and M109 as sent by Fauquembergue to André Gérardin. http://www.utm.edu/research/primes/mersenne/LukeMirror/lit/lit_008s.htm These computations done in 1913-1914 were... false, as found by Robinson 40 years later: http://www.utm.edu/research/primes/mersenne/LukeMirror/lit/lit_024s.htm Within the huge number of handwritten 0s and 1s, the game would be to locate the first line where the error occurs! I can send .JPG photos to interested people. Christian.

1 1

RE: [math-fun] FW: Beautiful Zeta(3) formulas!
by Meeussen Wouter (bkarnd) 30 May '05

30 May '05

Harrumph indeed! Simply replacing n!^p /(n+w)!^p by 1/Product[n+k,{k,1,w}] makes it easy: Sum[n!^p/(n+w)!^p,{n,0,oo}] == HypergeometricPFQRegularized[{1,1 ..p+1 copies}, {1 + w, 1 + w,.. p copies}, 1] In[1]:= Table[Sum[1/(Product[n + k, {k, 1, w}])^p, {n, 0, \[Infinity]}] // FullSimplify, {w, 1, 7, 2}, {p, 3, 7, 2}] // InputForm // Timing Out[1]={2.781 Second, { {Zeta[3], Zeta[5], Zeta[7]} , {29/32 - (3*Zeta[3])/ 4, (-3*(-791 + 520*Zeta[3] + 160*Zeta[5]))/ 512, (-3*(-35731 + 21448*Zeta[3] + 8288*Zeta[5] + 1344*Zeta[7]))/ 4096} , {(5*(-2077 + 1728*Zeta[3]))/ 1327104, (-1186062547 + 781704000*Zeta[3] + 237634560*Zeta[5])/ 330225942528, (5*(-1404612185513 + 848813853888*Zeta[3] + 322757084160*Zeta[5] + 49203044352*Zeta[7]))/ 4565043429507072} , {(7*(389467 - 324000*Zeta[3]))/ 503884800000, (-7*(-784046040073 + 516756240000*Zeta[3] + 157075200000*Zeta[5]))/ 125382662553600000000, (156827579578950685957 - 94800261054346848000*Zeta[3] - 36021058152710400000*Zeta[5] - 5475321411840000000*Zeta[7])/584985350410076160000000000000} } } W. -----Original Message----- From: R. William Gosper [mailto:rwg@osots.com] Sent: maandag 30 mei 2005 9:10 To: math-fun(a)mailman.xmission.com Subject: Re: [math-fun] FW: Beautiful Zeta(3) formulas! rcs>Well, cute anyway. I assume these aren't news to some of our readers. I should hope not. >-----Original Message----- From: Number Theory List [mailto:NMBRTHRY@LISTSERV.NODAK.EDU] On Behalf Of Troy Kessler Sent: Friday, May 27, 2005 6:46 AM To: NMBRTHRY(a)LISTSERV.NODAK.EDU Subject: Beautiful Zeta(3) formulas! 29/32 -3*Zeta[3]/4 = sum n=0 to oo (n)!^3/(n+3)!^3 (mathematica can do) 5*(-2077+1728Zeta[3])/1327104 = sum n=0 to oo (n)!^3/(n+5)!^3 Mathematica gives an answer that is 5 lines long for this second sum. Harrumph. (c154) closedform(sum(n!^3/(n+5)!^3,n,1,inf)) Time= 1131 msec. 5 zeta(3) 1298221 (d154) --------- - --------- 768 165888000 (c156) closedform(sum(n!^3/(n+7)!^3,n,1,inf)) Time= 2163 msec. 935108917 7 zeta(3) (d156) --------------- - --------- 172832486400000 1555200 =============================== This email is confidential and intended solely for the use of the individual to whom it is addressed. If you are not the intended recipient, be advised that you have received this email in error and that any use, dissemination, forwarding, printing, or copying of this email is strictly prohibited. You are explicitly requested to notify the sender of this email that the intended recipient was not reached.

1 0

[math-fun] FW: Beautiful Zeta(3) formulas!
by Schroeppel, Richard 30 May '05

30 May '05

Well, cute anyway. I assume these aren't news to some of our readers. Rich -----Original Message----- From: Number Theory List [mailto:NMBRTHRY@LISTSERV.NODAK.EDU] On Behalf Of Troy Kessler Sent: Friday, May 27, 2005 6:46 AM To: NMBRTHRY(a)LISTSERV.NODAK.EDU Subject: Beautiful Zeta(3) formulas! 29/32 -3*Zeta[3]/4 = sum n=0 to oo (n)!^3/(n+3)!^3 (mathematica can do) 5*(-2077+1728Zeta[3])/1327104 = sum n=0 to oo (n)!^3/(n+5)!^3 Mathematica gives an answer that is 5 lines long for this second sum. It can be shown that if 3 or 5 is replaced with a odd number the result is rational + rational *Zeta[3]. (Use Zeilberger in A=B)

2 1

[math-fun] number theory question
by wouter meeussen 28 May '05

28 May '05

since 1/Product[(1 - 1/Prime[k]^w), {k, \[Infinity]}] equals Zeta[w], I wondered what this product becomes when specialising for either primes p with 4|(p+1) or primes q with 4|(q-1). Numerics at w=2 point to 1.690 and 1.335 aboutish. W.

2 1

[math-fun] another Roots-of-unity puzzle
by Scott Huddleston 27 May '05

27 May '05

Dan Asimov <dasimov(a)earthlink.net> posed the following Roots-of-unity puzzle: > For each n >= 2, find f(n) = the size of the largest > subset of the n complex nth roots of unity that > does not contain k equally spaced points for any > k >= 2. > > (I.e., which does not contain a congruent copy of > the k complex kth roots of unity for any k >= 2.) I'll pose a modification -- find the subsets from Dan's puzzle whose root sum has minimum absolute value. Can someone prove or disprove that this minimum absolute value > 0 ? - Scott

1 0

[math-fun] Roots-of-unity puzzle
by dasimov＠earthlink.net 27 May '05

27 May '05

For each n >= 2, find f(n) = the size of the largest subset of the n complex nth roots of unity that does not contain k equally spaced points for any k >= 2. (I.e., which does not contain a congruent copy of the k complex kth roots of unity for any k >= 2.) --Dan

1 0

[math-fun] dgf's and binary sequences
by Mike Stay 27 May '05

27 May '05

Where's the best place to find out what's known about Dirichlet generating functions for binary sequences? -- Mike Stay staym(a)clear.net.nz http://www.cs.auckland.ac.nz/~msta039

1 0

RE: [math-fun] fun with A := solid_angle(a,b,c)
by dasimov＠earthlink.net 26 May '05

26 May '05

This paragraph should have had the highlighted words: << Back to K: As M approaches K smoothly (at least in the C^3 topology) the Gauss map is no longer well-defined at edges or vertices of M. But its image can be replaced by the set of outward unit normals to all planes P in space that (assuming K convex) that *intersect K and* keep K inside of one of the closed half-spaces defined by P. The area of the spherical image of all the unit normals must still be 2pi*X(K) = 4*pi. >> P.S. (This also works for nonconvex K but requires more words.)

1 0

RE: [math-fun] fun with A := solid_angle(a,b,c)
by dasimov＠earthlink.net 26 May '05

26 May '05

Rich writes: << There's a plane geometry theorem that the sum of the vertex angles of a polygon is pi*(#vertices-2). Is there a 3D analog? >> In short, the external spherical angles at all vertices of a polyhedron K must add up to 4*pi. (As in 2D this can be thought of as a limiting case of the embedded Gauss-Bonnet theorem for smooth surfaces in 3D (resp. smooth curves in 2D). Approximate a polyhedron K by an (internal) smooth surface M. M will have a Gauss map G: M -> S^2 -- G(x) = the point of S^2 having its outward unit normal = the same vector as the outward unit normal of M at x. Let JDG: M -> R be the Jacobian determinant, i.e. the factor by which G magnifies area at each point of. Then Gauss-Bonnet tells us that the (integral over M of JDG dA) = 2pi*X(M), where X(M) = its Euler characteristic = 2 for K a topological sphere. Back to K: As M approaches K smoothly (at least in the C^3 topology) the Gauss map is no longer well-defined at edges or vertices of M. But its image can be replaced by the set of outward unit normals to all planes P in space that (assuming K convex) keep K inside of one of the closed half-spaces defined by P. The area of the spherical image of all the unit normals must still be 2pi*X(K) = 4*pi. This is the sum over all vertices of the external spherical angle. It's a nice little exercise to calculate the formula for the spherical angle at a vertex where 3 edges meet, whose unit edge vectors (i.e., vectors leaving the vertex in the directions of the edges it's in) are say u,v,w. --Dan

1 0