Testwiki:Reference desk/Archives/Mathematics/2016 June 7

From testwiki
Revision as of 02:53, 14 June 2016 by imported>Scsbot (edited by robot: archiving June 7)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Template:Error:not substituted

{| width = "100%"

|- ! colspan="3" align="center" | Mathematics desk |- ! width="20%" align="left" | < June 6 ! width="25%" align="center"|<< May | June | Jul >> ! width="20%" align="right" |Current desk > |}

Welcome to the Wikipedia Mathematics Reference Desk Archives
The page you are currently viewing is a transcluded archive page. While you can leave answers for any questions shown below, please ask new questions on one of the current reference desk pages.


June 7

Algorithm to enumerate (n^2+n) sets of size n from n^2 items

I am seeking an algorithm to enumerate (n^2+n) sets of size n from n^2 items where every pair of items are in common in one and only one set.

Firstly, I believe it is possible. Secondly I have an algorithm that works for prime n: Arrange the items in a square. Read the columns for the "+n" sets. Then (for the n^2 sets) form n bunches of n sets, each bunch starting with an item from the left hand column and taking one item from each subsequent column using a different "drop per column" for each bunch. I'm using 'bunch' as a non-mathsy word instead of the more natural word "group" which might have math meaning.

Example

a b c
d e f
g h i

Columns adg, beh, cfi;
Drop 0 bunch abc, def, ghi;
Drop 1 bunch aei, dhc, gbf;
Drop 2 bunch ahf, dbi, gec.

This algorithm fails when n is non-prime. I can manually enumerate 20 sets for n=4, but can't extend the method to 6, 8 or 9 yet.

Any suggestions please? -- SGBailey (talk) 06:14, 7 June 2016 (UTC)

Would that be the Steiner system S(2,n,n2)? I think this may be slightly useful to get started. 5.11 seems to imply that the sets you wish to generate might not be possible for non-prime power n. —  crh 23  (Talk) 14:10, 7 June 2016 (UTC)
As you say, it looks as though Steiner is what I want. I (believe I) can do S(2, prime, prime^2) and S(2,4,16). My first stumbling block would be S(2,6,36). Your refs need careful reading. Whether 5.11 says it is impossible I need to study. Thanks - this will keep me going for a week or too. -- SGBailey (talk) 16:27, 7 June 2016 (UTC)
You might try [1], section 4 in particular. It mentions the Thirty-six officers problem in regard to S(2,6,36). --RDBury (talk) 08:18, 8 June 2016 (UTC)

Elegant Solution to Sharp Inequality

How could we prove, without the aid of a calculator, that 125512>12 ? Moving the negative term to the right hand side, and then exponentiating, is —for painfully obvious reasons— unfeasible. Perhaps some clever manipulation might show the way out of this impasse, but I fail to see how... — 82.137.54.252 (talk) 07:51, 7 June 2016 (UTC)

Why do you think there is an easy way? 1255121.643751.14353=0.50022 That is so close to 0.5 that it doesn't leave much room for approximations, and you seem to reject both calculating the root or expanding a long exponential, so I don't see where you have any path forward. Dragons flight (talk) 08:10, 7 June 2016 (UTC)
Given my username and the section title, I feel somewhat obligated to attempt to reply to this, but I don't see a way either. The difference is so close to 0.5 that any approximations that would be simple enough to use without a calculator would not be accurate enough to prove the inequality. Double sharp (talk) 09:01, 7 June 2016 (UTC)
(Huh. So the word "inequation" is apparently a thing. I confess that I do not remember having heard it.) Double sharp (talk) 13:38, 7 June 2016 (UTC)
Hm. inequation says that it is a statement of inequality (which is a relation). But equation says that one of those is an equality. This is a troublesome state of affairs... SemanticMantis (talk) 14:27, 7 June 2016 (UTC)
I've tweaked the opening sentence of equation accordingly. Loraof (talk) 15:17, 7 June 2016 (UTC)
The usual term for these sorts of sentences is "inequality", not "inequation". I smell the odor of MathWorld here, which is one of the two references at inequation, the other being a book called "The A to Z of Mathematics" or some such, another name that does not inspire high confidence (sounds like a tertiary ref).
Wikipedia should not be used for these sorts of language-reform efforts. My opinion is that we should rework the articles to refer to the sentences primarily as "inequalities", while mentioning somewhere that the term "inequation" also exists. --Trovatore (talk) 16:01, 7 June 2016 (UTC)
I have a hazy recollection from long ago, an undergrad math professor using these words a bit differently to distinguish assertions from assignations, different from the notion in our articles of distinguishing statements from relations. E.g. 4=2+2 is an assertion (and a statement of equality) but 4=x is more like an assignment; it is an equation by fiat only. Does that ring a bell with anyone else? My first thought when seeing 'inequation' was that maybe somebody was taking this distinction to inequality as well, e.g. maybe 4>2 is considered an inequation while 4>x is considered an inequality. I cannot support this idea with any references at present though, so maybe my hazy recollections are playing tricks on me :) Finally, I'll add that cursory googling seems to indicate that 'inequation' has more use in BrEng (e.g. BBC [2]), so there may be an WP:ENGVAR issue conflated with a semantic one. SemanticMantis (talk) 18:29, 7 June 2016 (UTC)
I think you need to look up assignation :-) I doubt that's the issue. 4=x is a fine equation in general, albeit one containing a free variable. Assignments are generally distinguished by context (for example, "let x=4"), and there's no obvious corresponding concept in inequalities (you can say "assume x<4", but that's an assumption, which is different from an assignment, though you might be able to dispense with assignments by replacing them with assumptions). However the ENGVAR thing is a possibility; I would be interested to hear more evidence on that. --Trovatore (talk) 18:40, 7 June 2016 (UTC)
It came out as slip of the tongue, but then I did look it up, and decided it was fine, as the second definition in my OSX copy of NOAD is pretty much what we mean when we talk about assignments for variables, and this is also supported by sense 2 of your your wiktionary link. In case you're curious, I haveAssignation: "the allocation or attribution of someone or something as belonging to something." Pretty much a synonym to the second sense of Assignment: "the attribution of someone or something as belonging." So maybe you could benefit by looking in to Muphry's_law ;)
Anyway, I do know how to talk about free variables and assumptions and assignments, I just thought I'd share a thought I had about a possible distinction. SemanticMantis (talk) 19:01, 7 June 2016 (UTC)
  • Does this work?
It can be shown by multiplication and division that 263/160 is an underestimate of 125
It can be shown by multiplication and division that 183/160 is an overestimate of 512
263/160183/160=80/160=0.5 exactly.
Hence, (125a)(512+b)=0.5; where both a and b are positive.
125a512b=0.5
125512=0.5+a+b; where both a and b are positive.
125512>0.5.

--NorwegianBlue talk 21:38, 8 June 2016 (UTC)

I found two fractions with smaller numerators and denominators for which the same logic holds: 166/101 and 143/87. --NorwegianBlue talk 19:41, 9 June 2016 (UTC)
Maybe you meant to subtract 1/2 from the smaller of your two fractions? (Your first example has the smallest max denominator for any pair of rational numbers differing by at least 1/2 in the relevant interval.) --JBL (talk) 20:23, 9 June 2016 (UTC)
My point was that 166/101 and 143/87 both are underestimates of 125
166/1011/2 and 143/871/2 are both overestimates of 512
Since the numerator of 143/87 is odd, we will need to double the denominator when subtracting 1/2, i.e. 199/174.
If there is anything wrong with my attempt at showing the feasibility of proving Sharp's inequality by manual (albeit tedious) calculation, please let me know! --NorwegianBlue talk 07:51, 10 June 2016 (UTC)
No, nothing wrong, I didn't understand that you meant one number to stand for a pair. (Of course both your examples produce one fraction with denominator larger than 160.) --JBL (talk) 19:53, 10 June 2016 (UTC)

It is done in integer arithmetic like this. The number z=512125 satisfies a polynomial equation with integer coefficients, which are determined by eliminating x and y from the equations

z=yx
y12=5
x5=12.

The elimination goes like this. First eliminate y by substituting y=x+z into y12=5

(x+z)12=5.

Expanding the binomial

i=012(12i)xiz12i=5.

Split the sum and substitute x5=12.

i=04(12i)xiz12i+i=59(12i)12xi5z12i+i=1012(12i)122xi10z12i=5.

Change indexes

i=04(12i)xiz12i+i=04(12i+5)12xiz125i+i=02(12i+10)122xiz1210i=5.

obtaining the degree 4 equation

i=04Aixi=0

where

A0=(120)z12+(125)12z7+(1210)122z25
A1=(121)z11+(126)12z6+(1211)122z
A2=(122)z10+(127)12z5+(1212)122
A3=(123)z9+(128)12z4
A4=(124)z8+(129)12z3

A second equation of degree 4 is obtained by multiplying by x

i=04Aixi+1=0

and using x5=12

i=03Aixi+1+12A4=0.

So now x5 is eliminated! The next step is to eliminate x4

i=03Bixi=0
i=03Cixi=0

where

B0=A0A312A42
B1=A1A3A0A4
B2=A2A3A1A4
B3=A32A2A4.

and

C0=B3A0
C1=B3A1A4B0
C2=B3A2A4B1
C3=B3A3A4B2

Then eliminate x3 :

i=02Dixi=0
i=02Eixi=0

where

D0=B0C3B3C0
D1=B1C3B3C1
D2=B2C3B3C2

and

E0=C0D2
E1=C3D0C1D2
E2=C3D1C2D2

Then eliminate x2 :

i=01Fixi=0
i=01Gixi=0

where

F0=D0E2D2E0
F1=D1E2D2E1

and

G0=E0F1
G1=E2F0E1F1

Finally eliminate x :

F0G1F1G0=0

This polynomial equation in z has the root z=512125 Bo Jacoby (talk) 14:18, 10 June 2016 (UTC).

Thank you, Bo Jacoby, that's fascinating—I was wondering how it could be done! Two comments: (1) B4 appears in the definitions of all Ci but is not itself defined. (2) What is the degree of this equation in z? It looks very high, but I keep getting lost when I try to count the degrees (and I don't know the degree of B4). Loraof (talk) 18:03, 11 June 2016 (UTC)

You are absolutely right. There is no B4. Thanks! I'll correct it. The C's are computed from this.

C0+C1x+C2x2+C3x3
=B3(A0+A1x+A2x2+A3x3+A4x4)
A4x(B0+B1x+B2x2+B3x3)

The degree of the equation is 5*12=60. Bo Jacoby (talk) 16:24, 13 June 2016 (UTC).

Also, is there a way to find out which root is the one we are looking for? Suppose we plot it against z, and find that there is a value below 1/2 and another value above 1/2 that both satisfy the equation. Does our desired value stand out in some way? Loraof (talk) 18:09, 11 June 2016 (UTC)

The equation y12=5 has 12 complex roots, forming a regular 12-gon, and the equation x5=12 has 5 complex roots forming a regular 5-gon, so of all the 60 differences z=yx only one is close to zero.

The equation

0=F0G1F1G0

is computed by backwards substitution. Using

G0=E0F1
G1=E2F0E1F1

get

0=F0(E2F0E1F1)F1(E0F1)

simplify

0=E2F02E1F0F1+E0F12

Using

F0=D0E2D2E0
F1=D1E2D2E1

get

0=E2(D0E2D2E0)2E1(D0E2D2E0)(D1E2D2E1)+E0(D1E2D2E1)2

simplify

0=+D02E23+D22E02E22D0D2E0E22D0D1E1E22+D1D2E0E1E2+D0D2E12E2D22E0E12+D12E0E22+D22E0E122D1D2E0E1E2

and so on. It isn't easy, but it is without round-off errors. Bo Jacoby (talk) 16:24, 13 June 2016 (UTC).