June 2013 - mathfun@mailman.xmission.com

[math-fun] Dean Hickerson
by James Propp 21 Oct '19

21 Oct '19

Anyone have a current email address for Dean Hickerson? "dean at math daht ucdavis daht edu" and "dean at ucdmath daht ucdavis daht edu" aren't working anymore. Thanks, Jim Propp

5 5

[math-fun] ellipsoid from girths: multivariate series reversion
by Bill Gosper 25 Jun '19

25 Jun '19

Corey notes that ellipsoid girths can't be {large, tiny, tiny} and challenges us to find the "ellipsoid girth inequality" analogous to the triangle inequality. What are the dimensions (semiaxes a,b, and c) of say, the ellipsoidoidal rock by the road near the Foothill College fire station, given the girths gab, gac, gbc; a>b>c? (Assuming you can get the damned thing off your tape measure once you've rolled it on.) In terms of a, b, c, the girths are {gab -> 4 a EllipticE[1 - b^2/a^2], gac -> 4 a EllipticE[1 - c^2/a^2], gbc -> 4 b EllipticE[1 - c^2/b^2]} This is a routine problem for FindRoot, but it would be nice to see how a, b, c depend on the girths in the low eccentricity case, say, by reverting the series expansions of the girths. Julian, Corey, and I were surprised by the absence of multivariate reversion from Macsyma and Mathematica, and set about building our own. After several hours of cruel hoodwinking by Mma's Series facility, we emerged having only solved the ellipsoid case. Surprisingly(?), it has at least eight (2^3) solutions, which may help to explain the shortage of multivariate reverters. To fourth order, {a -> (gab/(2 Pi))(1 - eac + eac^2 + ebc + (eac ebc)/2 - (3 ebc^2)/2 + 1/8 (-8 eac^3 - eac^2 ebc - 9 eac ebc^2 + 18 ebc^3) + 1/16 (16 eac^4 + 2 eac^3 ebc - 3 eac^2 ebc^2 + 42 eac ebc^3 - 57 ebc^4)), a -> (gac/(2 Pi))(1 + ebc + 3/2 (eac ebc - ebc^2) + 3/8 (eac^2 ebc - 7 eac ebc^2 + 6 ebc^3) - 3/16 (7 eac^2 ebc^2 - 26 eac ebc^3 + 19 ebc^4))} {b -> (gab/(2 Pi))(1 + eac - (3 eac^2)/2 - ebc + (eac ebc)/2 + ebc^2 + 1/8 (18 eac^3 - 9 eac^2 ebc - eac ebc^2 - 8 ebc^3) + 1/16 (-57 eac^4 + 42 eac^3 ebc - 3 eac^2 ebc^2 + 2 eac ebc^3 + 16 ebc^4)), b -> (gac/(2 Pi))(1 + 2 eac - eac^2/2 - ebc - (eac ebc)/2 + ebc^2 + 1/8 (6 eac^3 - 5 eac^2 ebc + 7 eac ebc^2 - 8 ebc^3) + 1/16 (-21 eac^4 + 24 eac^3 ebc - 5 eac^2 ebc^2 - 14 eac ebc^3 + 16 ebc^4))} {c -> (gab/(2 Pi))(1 - eac + eac^2 + eac^4 - ebc - (eac ebc)/2 - (eac^3 ebc)/ 8 + ebc^2 - (5 eac^2 ebc^2)/16 - (eac ebc^3)/8 + ebc^4 + 1/8 (-8 eac^3 + eac^2 ebc + eac ebc^2 - 8 ebc^3)), c -> (gac/(2 Pi))(1 - ebc - (3 eac ebc)/2 - (3 eac^2 ebc)/8 + ebc^2 + (9 eac ebc^2)/8 - (3 eac^2 ebc^2)/16 - ebc^3 - ( 9 eac ebc^3)/8 + ebc^4)} where {eac -> gab/gac - 1, ebc -> gab/gbc - 1}, which are small for low eccentricity. E.g., for {a -> 5., b -> 4., c -> 3.}, the girth formulæ give {gab -> 28.3617, gac -> 25.527, gbc -> 22.1035}. Plugging these back into the reversions: {{a -> 4.97828, a -> 4.98561}, {b -> 4.00689, b -> 4.00371}, {c -> 3.0062, c -> 3.00353}}. --rwg

3 5

[math-fun] muffin problem -- clusters and bounds
by Huddleston, Scott 01 Dec '17

01 Dec '17

I enclose examples of some of the bounds in my 28 Feb. post, and a derivation of the most involved example. S(12/7) = 3/7 = r/4 = (3-r)/3 Muffins: 1/7 * { [3,4]*12 } Persons: 1/7 * { [4*3]*4, [3*4]*3 } S(7/4) = 5/12 = (3-r)/3 = (2r-1)/6 Muffins: 1/12 * { [5,7]*6, [6*2] } Persons: 1/12 * { [5*3,6]*2, [7*3]*2 } S(40/23) = 19/46 = (2r-1)/6 = (12-5r)/8 Muffins: 1/46 * { [19,27]*24, [20,26]*12, [23*2]*4 } Persons: 1/46 * { [19*3,23]*8, [20*4]*3, [27*2,26]*12 } S(228/131) = 54/131 = (12-5r)/8 = (7r-6)/15 Muffins: 1/131 * { [54,77]*120, [55,76]*24, [57,74]*60, [65,66]*24 } Persons: 1/131 * { [54*3,66]*24, [54*2,55,65]*24, [57*4]*15, [77*2,74]*60, [76*3]*8 } S(449/258) = 319/774 = (7r-6)/15 = (73-32r)/42 Muffins: 1/774 * { [319,455]*210, [320,454]*30, [325,449]*42, [336,438]*30, [337,437]*90, [384,398]*42, [387*2]*5 } Persons: 1/774 * { [319*3,390]*42, [319*2,325,384]*42, [320*3,387]*10, [336,337*3]*30, [455*2,437]*90, [455,454,438]*30, [449*3]*14 } Here's a derivation of 449/258 example. 449 muffins -> 898 pieces, used in 124 4-parts with 496 pieces and 134 3-parts with 402 pieces. We split muffins between a 3-part and a 4-part when possible, but 4-parts contain 94 more pieces than 3-parts. We divide these 94 pieces (47 muffins) as evenly as possible among the 124 4-parts. 47 pairs of 4-parts contain both halves of a muffin, which I'll call an A^1 cluster and write A^1 = [a1*3,m]*2 where m denotes the two muffin halves and the other pieces are denoted a1, indicating a set of pieces of possibly different sizes that occure in A^1 clusters. 34 4-parts contain only muffin pieces whose mates appear in 3-parts. These I call A^0 clusters and write A^0 = [a0*4] The 449/258 example uses 30 A^0 clusters ([a0*4])*30 (120 a0 pieces, 30 4-parts) 47 A^1 clusters ([a1*3,m]*2)*47 (282 a1 pieces, 94 4-parts) Note that the A^1 cluster by itself sets an upper bound on S -- the cluster is 2 persons (size 2r), it cntains 1 complete muffin (size 1, pieces m*2), and it has 6 unmatched muffin pieces (a1*6), giving S(449/258) <= avg(a1) = (2r-1)/6. A^0 clusters also bound S by S(449/258) <= avg(a0) = r/4, but the A^1 cluster bound is tighter. Note that this (2r-1)/6 bound holds for examples 7/4 and 40/23. The muffin mates of a0 pieces (resp. a1 pieces) are z0 pieces (resp. z1 pieces). We distribute z0 and z1 pieces as evenly as possible across the 134 3-parts (402 pieces). With 282 z1 pieces, at least 14 3-parts must consist of only z1 pieces. These are Z^0 clusters, Z^0 = [z1*3] = [z10*3] where I refine piece labels so that z1 pieces appearing in Z^0 clusters are labelled z10. I create general Z^k clusters by finding minimal sets of persons that contain complete a0/z0 muffins. Maximally even distribution of the 120 z0 pieces places one per 3-part. Then a Z^1 cluster, closed under a0/z0 muffins, is Z^1 = [a0*4],[z0,z1*2]*4 = A^0,[z0,z11*2]*4 where z1 pieces appearing in Z^1 clusters are relabelled z11. Z^0 clusters bound S by S(449/258) <= avg(a10) = 1-avg(z10) = 1-r/3 = (3-r)/3 and Z^1 clusters provide an even tighter bound S(449/258) <= avg(a11) = 1-avg(z11) = 1-(5r-4)/8 = (12-5r)/8. The 449/258 example uses 14 Z^0 clusters = ([z10*3])*14 ( 42 z10 pieces, 14 3-parts) 30 Z^1 clusters = (A^0,[z0,z11*2]*4)*30 (240 z11 pieces, 120 3-parts) As before, we refine the labelling of a1 pieces into a10 and a11 pieces, according to z10 or z11 complements. The next level of clustering is constructed by taking closure under a10/z10 muffins. The 282 a1 pieces in the 47 A^1 clusters are divided into 240 a11 pieces and 42 a10 pieces. Thus at least 5 A^1 clusters must contain only a11 pieces, and with maximally even distribution we have 42 A^1 clusters with 1 a10 piece and 5 a11 pieces. Continuing to refine pieces according to their appearance in A^10 or A^11 clusters, we have A^10 = A^1 = {a110*6, ...} (6 a110 pieces and some other pieces) A^11 = Z^0, (A^1)*3 = (a111*15, ...) with 449/258 consisting of 5 A^10 clusters ( 5 A^1 clusters, 30 a110 pieces) 14 A^11 clusters (42 A^1 clusters, 252 a111 pieces) A^11 clusters provide bound S(449/258) <= avg(a111) = (7r-6)/15 The 449/258 example lies on the intersection of this bound and the next level cluster bound. The next cluster level takes a110/z110 muffin closure and induces clusters Z^10 = Z^1 Z^11 = A^10, (Z^1)*6 with 449/258 consisting of 0 Z^10 clusters 5 Z^11 clusters (i.e., only Z^11 clusters) and bound S(449/258) <= avg(a1111) = 1-avg(z1111) = (73-32r)/42. I conclude with some computer generated additional information on these bounds. If anyone sort of understands this post and would like to use my program as is (i.e., no manual and documentation is mostly my posts), drop me a note and I'll send the C++ source (300 lines). B<4>(r) = (1r-0)/4; S'(r) <= B<4>(r) on (3/2,2/1] tight bound on (3/2,8/5] and at r = (4+8k)/(2+5k) (3/2,8/5] contains mediant(3/2,8/5) = 11/7 S'( (3/2,8/5] ) = (3/8,2/5]; S'(11/7) = 11/28 B<>(r) = (1-0r)/2 B<4>(r) = B<>(r) at r = 2/1 = 2.0000000 S'(2/1) = 1/2 = 0.50000000 B<4,1>(r) = (3-1r)/3; S'(r) <= B<4,1>(r) on [12/7,2/1) tight bound on [9/5,2/1) and at r = (12+9k)/(7+5k) [9/5,2/1) contains mediant(9/5,2/1) = 11/6 S'( [9/5,2/1) ) = [2/5,1/3); S'(11/6) = 7/18 B<4>(r) = (1r-0)/4 B<4,1>(r) = B<4>(r) at r = 12/7 = 1.7142857 S'(12/7) = 3/7 = 0.42857143 B<4,1,1>(r) = (2r-1)/6; S'(r) <= B<4,1,1>(r) on (12/7,7/4] tight bound on (12/7,59/34] and at r = (21+59k)/(12+34k) (12/7,59/34] contains mediant(12/7,7/4) = 19/11 S'( (12/7,59/34] ) = (17/42,7/17]; S'(19/11) = 9/22 B<4,1>(r) = (3-1r)/3 B<4,1,1>(r) = B<4,1>(r) at r = 7/4 = 1.7500000 S'(7/4) = 5/12 = 0.41666667 B<4,1,1,1>(r) = (12-5r)/8; S'(r) <= B<4,1,1,1>(r) on [40/23,7/4) tight bound on [148/85,7/4) and at r = (80+148k)/(46+85k) [148/85,7/4) contains mediant(47/27,7/4) = 54/31 S'( [148/85,7/4) ) = [7/17,13/32); S'(54/31) = 51/124 B<4,1,1>(r) = (2r-1)/6 B<4,1,1,1>(r) = B<4,1,1>(r) at r = 40/23 = 1.7391304 S'(40/23) = 19/46 = 0.41304348 B<4,1,1,1,1>(r) = (7r-6)/15; S'(r) <= B<4,1,1,1,1>(r) on (40/23,228/131] tight bound on (40/23,1119/643] and at r = (228+1119k)/(131+643k) (40/23,1119/643] contains mediant(40/23,47/27) = 87/50 S'( (40/23,1119/643] ) = (142/345,265/643]; S'(87/50) = 103/250 B<4,1,1,1>(r) = (12-5r)/8 B<4,1,1,1,1>(r) = B<4,1,1,1>(r) at r = 228/131 = 1.7404580 S'(228/131) = 54/131 = 0.41221374 B<4,1,1,1,1,1>(r) = (73-32r)/42; S'(r) <= B<4,1,1,1,1,1>(r) on [449/258,228/131) tight bound on [9349/5372,228/131) and at r = (1347+9349k)/(774+5372k) [9349/5372,228/131) contains mediant(134/77,47/27) = 181/104 S'( [9349/5372,228/131) ) = [1107/2686,2267/5502); S'(181/104) = 75/182 B<4,1,1,1,1>(r) = (7r-6)/15 B<4,1,1,1,1,1>(r) = B<4,1,1,1,1>(r) at r = 449/258 = 1.7403101 S'(449/258) = 319/774 = 0.41214470 B<4,1,1,1,1,1,1>(r) = (61r-60)/112; S'(r) <= B<4,1,1,1,1,1,1>(r) on (449/258,764/439] tight bound on (449/258,149516/85913] and at r = (10696+149516k)/(6146+85913k) (449/258,149516/85913] contains mediant(449/258,764/439) = 1213/697 S'( (449/258,149516/85913] ) = (11909/28896,35408/85913]; S'(1213/697) = 32173/78064 B<4,1,1,1,1,1>(r) = (73-32r)/42 B<4,1,1,1,1,1,1>(r) = B<4,1,1,1,1,1>(r) at r = 764/439 = 1.7403189 S'(764/439) = 2533/6146 = 0.41213798 B<4,1,1,1,1,1,1,1>(r) = (1101-487r)/615; S'(r) <= B<4,1,1,1,1,1,1,1>(r) on [160212/92059,764/439) tight bound on [1313469/754729,764/439) and at r = (160212+6567345k)/(92059+3773645k) [1313469/754729,764/439) contains mediant(19549/11233,764/439) = 20313/11672 S'( [1313469/754729,764/439) ) = [1555262/3773645,111271/269985); S'(20313/11672) = 986147/2392760 B<4,1,1,1,1,1,1>(r) = (61r-60)/112 B<4,1,1,1,1,1,1,1>(r) = B<4,1,1,1,1,1,1>(r) at r = 160212/92059 = 1.7403187 S'(160212/92059) = 37941/92059 = 0.41213787

2 1

[math-fun] A problem about finite sets
by Victor Miller 09 Sep '16

09 Sep '16

Say that a finite set of integers, A, has *property P* is there are two non-empty disjoint subsets S and T of A, such that sum(S) = sum(T) (where sum(X) is the sum of the elements of X). 1) Show that every 10 element subset of {1,...,100} has property P (this is fairly easy). 2) Is it true that every 9 element subset of {1,...,100} has property P? Victor

8 9

Re: [math-fun] Raghavendra's algorithm for solving linear systems over finite fields
by Henry Baker 06 Jun '16

06 Jun '16

I can see applications in distributed/parallel implementations. You can do most of the checking & probably a lot of the randomization in parallel. At 09:02 AM 8/15/2012, Warren Smith wrote: >>Preliminary draft: >> http://www.eecs.berkeley.edu/~prasad/linsystems.pdf > >--I don't understand why anybody would want to use this algorithm >for any purpose. It's slower than everything else I can think of and also >randomized i.e. can fail. > >However... perhaps it might be useful if changed to be about finite >RINGS not fields. In rings, you cannot divide, hence gaussian elimination is >not usable. Hence this algorithm then has less competition. >(Even then I'm dubious, but then it'd have at least a hope.)

8 9

[math-fun] xkcd
by rcs＠xmission.com 17 Dec '15

17 Dec '15

Today's XKCD cartoon has a nice collection of approximations. http://xkcd.com/1047/ There's the ancient expression 6 pi^5 for the mass ratio of the proton to the electron, and the more accurate (e^8 - 10)/phi. Robert Munafo & REIS are at least partially responsible. Rich

5 4

[math-fun] Julian pimped my pump
by Bill Gosper 17 Dec '15

17 Dec '15

To run it <http://gosper.org/pumper_2.mp4> continuously, download with Miro, QuickTime-View-Loop. The pumped shapes split lengthwise resemble boats, so a grand entrance to a math museum would be a chain of these boats in an annular canal, passing through four rotors acting as an (approximate) air|water lock, with another four for the exit. Oh, the insurance. --rwg No manual entry for overboard

5 8

[math-fun] Terminology: "Characteristic Generating Function"?
by Bill Gosper 22 Oct '15

22 Oct '15

"Indicator Generating Function"? "Predicate Generating Function"? What is the name for the G.F. that encodes an integer sequence in the exponents with coefficients in {0,1} ? E.g., the OGF for the squares is Sum[n^2*x^n, {n, ∞}] == -((x*(1 + x))/(-1 + x)^3), but the thing I'm looking for is Sum[x^n^2,{n,0,∞}]==(1/2)*(1 + EllipticTheta[3, 0, x]). (Years ago I sent some of you the GF for expressions for k as the sum of n (optionally distinct) members of such a set, in terms of the GF predicating membership in the set. It turned out to be a superspecial case of a formula in Generatingfunctionology that was so general I didn't even recognize it as the same problem.) The reason I ask is that Neil just came up with 1/2 (-1+EllipticTheta[2,0,x]/x^(1/4)+EllipticTheta[3,0,x]) predicating the quartersquares<http://en.wikipedia.org/wiki/Multiplication_algorithm#Quarter_square_multip…> A002620 <https://oeis.org/A002620>, from which he got (1/(2 (-1 + x)^2))((-1 + x) x^(3/4) EllipticTheta[2, 0, x] + x (3 - x + (-1 + x) EllipticTheta[3, 0, x])) generating what we think is the number of volumes you can pump with n of the recently discussed rotors, arranged optimally. --rwg

4 3

Re: [math-fun] Another(?) stick length illusion
by Dan Asimov 25 Aug '15

25 Aug '15

Huh? At least to the nearest pixel, both lines are the same length. They are corresponding diagonals of congruent regular pentagons. I find this to be an impressive illusion. --Dan << << To me the stick that is in fact longer looks like it is in fact longer. If they're supposed to be the same length I think your image needs a bit of editing. >> On Thu, Aug 19, 2010 at 4:31 PM, Bill Gosper <billgosper(a)gmail.com> wrote: > While hacking the pentagon fill, tutor Julian found > www.tweedledum.com/rwg/6pents.png, which seems to work backwards: > The stick larger than its context appears smaller than the > stick smaller than its context. >> _____________________________________________________________________ "It don't mean a thing if it ain't got that certain je ne sais quoi." --Peter Schickele

5 4

[math-fun] the 19th singular value
by Bill Gosper 16 Jul '15

16 Jul '15

ETA(%E^-(2*SQRT(19)*%PI)) = ((3*SQRT(57)+1)^(1/3)-8/(3*SQRT(57)+1)^(1/3))^(1/8)*(GAMMA(1/38)*GAMMA(5/38)*GAMMA(9/38)*GAMMA(11/38)*GAMMA(7/19)*GAMMA(17/38)*GAMMA(23/38)*GAMMA(25/38)/(38*GAMMA(13/38)*GAMMA(8/19)))^(1/4)/(2*2^(53/152)*%PI) where Corey's denester neatly whacks the 8th root: (-(8/(1 + 3*Sqrt[57])^(1/3)) + (1 + 3*Sqrt[57])^(1/3))^(1/8) == -(2^(3/8)/3) + ((-1 + 3*Sqrt[57])^(1/3)*(4*(-5 + Sqrt[57]) + (11 + Sqrt[57])*(-1 + 3*Sqrt[57])^(1/3)))/(96*2^(5/8)) . It also denests LambdaStar[19], but the result is badly undersimplified. Fullsimplify thuggishly reverts to the original Root object, even when given a ComplexityFunction that enormously penalizes _Root ! But it is possible to guess the general form from the denester mess and determine coefficients with PSLQ: LambdaStar[19]==(((-Sqrt[57] + Sqrt[19] + 7*Sqrt[3] + 11)*(3*Sqrt[57] + 1)^(1/3)/ 16) - ((Sqrt[57] + Sqrt[19] - 7*Sqrt[3] + 11)/(2*(3*Sqrt[57] + 1)^(1/3))) - ((Sqrt[19] - 1)/2))/(3*Sqrt[2]) I was unable to Google whether this denesting is known. Is there a table of singular values more extensive than MathWorld's? --rwg Joerg>Let's not forget Broadhurst's very neat expressions: Message-ID: <200803140250.m2E2oHi4002315(a)philter.princeton.idaccr.org> Date: Thu, 13 Mar 2008 22:50:17 -0400 From: David Broadhurst <D.Broadhurst(a)open.ac.uk> To: NMBRTHRY(a)LISTSERV.NODAK.EDU Reply-To: David Broadhurst <D.Broadhurst(a)open.ac.uk> Subject: Singular value for Euler's numerus idoneus N=1365 default(realprecision,1050); a_1 = 550*sqrt(13) + 318*sqrt(39) + 750*sqrt(7) + 433*sqrt(21); a_2 = 1986 + 1145*sqrt(3) + 208*sqrt(91) + 120*sqrt(273); b = ((a_1-a_2)^2 - 2)/2; k_1365 = sqrt(1/2 - b^2/4*sqrt(4-b^2) - (2-b^2)/4*sqrt(1-b^2)); --------------- The denester also reduced the depth of k_1365, but made a proverbial Very Large Expression, more than Mathematica could further chew.

2 2