[seqfan] Re: Puzzle

Thu Jan 22 14:33:48 CET 2015

It seems that a real (or complex) number r is in David Wilson's set if
and only if there exist positive integers N and k_1, ..., k_N such
that r is a zero of the following polynomial of degree N-1 with
integer coefficients:
\sum_{i=1}^N a_i X^{i-1}
where each a_i is the following sum of i-fold products of distinct k_j's:
the summands are the products k_{j_1}...k_{j_i} such that 1 <= j_1 <
... < j_i <= N and the difference between two consecutive j_m's is
odd, and this product is given the sign (-1)^{j_1}.

I hope I didn't make too many mistakes while quickly computing this.
My method was to multiply matrices of the form (0,1;1,k_i r).

Benoit

On Thu, Jan 22, 2015 at 2:59 AM,  <hv at crypt.org> wrote:
> Very nice. :)
>
> If I counted right, that's the sequence [
>   18, 3, 201, 3, 3, 1, 39, 117, 1, 3, 1, 25, 24, 2, 13, 5, 3, 1, 1, 5, 1
> ].
>
> When trying to find such a solution, I guessed that the best chance lay in
> ascending always to a rational with numerator divisible by 9 but not by 27,
> but I see that your progression frequently switches with them divisible
> only by 3, and in one case not divisible by 3 at all.
>
> Hugo
>
> Matthijs Coster <seqfan at matcos.nl> wrote:
> :However there is a solution for r = 7/9:
> :I jump immediately from t to 1/(r+t)
> :
> :1, 16/9, 23/9, 10/3, 37/9, 44/9, 17/3, 58/9, 65/9, 8, 79/9, 86/9, 31/3,
> :100/9, 107/9, 38/3, 121/9, 128/9, 1/15, 38/45, 73/45, 5/12, 43/36,
> :71/36, 11/4, 127/36, 155/36, 61/12, 211/36, 239/36, 89/12, 295/36,
> :323/36, 39/4, 379/36, 407/36, 145/12, 463/36, 491/36, 173/12, 547/36,
> :575/36, 67/4, 631/36, 659/36, 229/12, 715/36, 743/36, 257/12, 799/36,
> :827/36, 95/4, 883/36, 911/36, 313/12, 967/36, 995/36, 341/12, 1051/36,
> :1079/36, 123/4, 1135/36, 1163/36, 397/12, 1219/36, 1247/36, 425/12,
> :1303/36, 1331/36, 151/4, 1387/36, 1415/36, 481/12, 1471/36, 1499/36,
> :509/12, 1555/36, 1583/36, 179/4, 1639/36, 1667/36, 565/12, 1723/36,
> :1751/36, 593/12, 1807/36, 1835/36, 207/4, 1891/36, 1919/36, 649/12,
> :1975/36, 2003/36, 677/12, 2059/36, 2087/36, 235/4, 2143/36, 2171/36,
> :733/12, 2227/36, 2255/36, 761/12, 2311/36, 2339/36, 263/4, 2395/36,
> :2423/36, 817/12, 2479/36, 2507/36, 845/12, 2563/36, 2591/36, 291/4,
> :2647/36, 2675/36, 901/12, 2731/36, 2759/36, 929/12, 2815/36, 2843/36,
> :319/4, 2899/36, 2927/36, 985/12, 2983/36, 3011/36, 1013/12, 3067/36,
> :3095/36, 347/4, 3151/36, 3179/36, 1069/12, 3235/36, 3263/36, 1097/12,
> :3319/36, 3347/36, 375/4, 3403/36, 3431/36, 1153/12, 3487/36, 3515/36,
> :1181/12, 3571/36, 3599/36, 403/4, 3655/36, 3683/36, 1237/12, 3739/36,
> :3767/36, 1265/12, 3823/36, 3851/36, 431/4, 3907/36, 3935/36, 1321/12,
> :3991/36, 4019/36, 1349/12, 4075/36, 4103/36, 459/4, 4159/36, 4187/36,
> :1405/12, 4243/36, 4271/36, 1433/12, 4327/36, 4355/36, 487/4, 4411/36,
> :4439/36, 1489/12, 4495/36, 4523/36, 1517/12, 4579/36, 4607/36, 515/4,
> :4663/36, 4691/36, 1573/12, 4747/36, 4775/36, 1601/12, 4831/36, 4859/36,
> :543/4, 4915/36, 4943/36, 1657/12, 4999/36, 5027/36, 1685/12, 5083/36,
> :5111/36, 571/4, 5167/36, 5195/36, 1741/12, 5251/36, 5279/36, 1769/12,
> :5335/36, 5363/36, 599/4, 5419/36, 5447/36, 1825/12, 5503/36, 5531/36,
> :1853/12, 5587/36, 5615/36, 4/627, 1475/1881, 2938/1881, 209/489,
> :1768/1467, 2909/1467, 163/450, 50/57, 283/171, 416/171, 61/19, 682/171,
> :815/171, 316/57, 1081/171, 1214/171, 449/57, 1480/171, 1613/171, 194/19,
> :1879/171, 2012/171, 715/57, 2278/171, 2411/171, 848/57, 2677/171,
> :2810/171, 327/19, 3076/171, 3209/171, 1114/57, 3475/171, 3608/171,
> :1247/57, 3874/171, 4007/171, 460/19, 4273/171, 4406/171, 1513/57,
> :4672/171, 4805/171, 1646/57, 5071/171, 5204/171, 19/593, 4322/5337,
> :8473/5337, 4208/1779, 16775/5337, 20926/5337, 8359/1779, 29228/5337,
> :33379/5337, 4170/593, 41681/5337, 45832/5337, 16661/1779, 54134/5337,
> :58285/5337, 20812/1779, 66587/5337, 70738/5337, 8321/593, 79040/5337,
> :83191/5337, 29114/1779, 91493/5337, 95644/5337, 33265/1779, 103946/5337,
> :108097/5337, 12472/593, 116399/5337, 120550/5337, 41567/1779,
> :128852/5337, 133003/5337, 45718/1779, 141305/5337, 145456/5337,
> :16623/593, 153758/5337, 157909/5337, 54020/1779, 166211/5337,
> :170362/5337, 58171/1779, 178664/5337, 182815/5337, 20774/593,
> :191117/5337, 195268/5337, 66473/1779, 203570/5337, 207721/5337,
> :70624/1779, 216023/5337, 220174/5337, 24925/593, 228476/5337,
> :232627/5337, 78926/1779, 240929/5337, 245080/5337, 83077/1779,
> :253382/5337, 257533/5337, 29076/593, 265835/5337, 269986/5337,
> :91379/1779, 278288/5337, 282439/5337, 95530/1779, 290741/5337,
> :294892/5337, 33227/593, 303194/5337, 307345/5337, 103832/1779,
> :315647/5337, 319798/5337, 107983/1779, 328100/5337, 332251/5337,
> :37378/593, 340553/5337, 344704/5337, 116285/1779, 353006/5337,
> :357157/5337, 120436/1779, 365459/5337, 369610/5337, 41529/593,
> :377912/5337, 382063/5337, 128738/1779, 390365/5337, 394516/5337,
> :132889/1779, 402818/5337, 406969/5337, 45680/593, 415271/5337,
> :419422/5337, 141191/1779, 427724/5337, 431875/5337, 145342/1779,
> :440177/5337, 444328/5337, 49831/593, 452630/5337, 456781/5337,
> :153644/1779, 465083/5337, 469234/5337, 157795/1779, 477536/5337,
> :481687/5337, 593/53982, 5998/4731, 29033/14193, 40072/14193, 1577/5679,
> :631/666, 383/222, 1667/666, 2185/666, 901/222, 3221/666, 3739/666,
> :473/74, 4775/666, 5293/666, 1937/222, 6329/666, 6847/666, 2455/222,
> :7883/666, 8401/666, 991/74, 9437/666, 9955/666, 3491/222, 10991/666,
> :11509/666, 4009/222, 12545/666, 13063/666, 74/1509, 3743/4527,
> :7264/4527, 3595/1509, 14306/4527, 17827/4527, 2372/503, 24869/4527,
> :28390/4527, 10637/1509, 35432/4527, 38953/4527, 14158/1509, 45995/4527,
> :49516/4527, 5893/503, 56558/4527, 60079/4527, 21200/1509, 67121/4527,
> :70642/4527, 24721/1509, 77684/4527, 81205/4527, 503/9414, 7825/9414,
> :1046/1683, 785/561, 3664/1683, 4973/1683, 698/187, 7591/1683, 8900/1683,
> :3403/561, 11518/1683, 12827/1683, 4712/561, 15445/1683, 16754/1683,
> :187/2007, 1748/2007, 1103/669, 4870/2007, 6431/2007, 223/888, 2741/2664,
> :4813/2664, 296/765, 85/99, 11/18, 25/18, 13/6, 53/18, 67/18, 2/9, 1
> :
> :Code in sage:
> :
> :def qsort(A):
> :   if A == []:
> :     return []
> :
> :   i = randrange(len(A))
> :   A_lo = qsort([x for x in (A[:i] + A[i + 1 :]) if x < A[i]])
> :   A_hi = qsort([x for x in (A[:i] + A[i + 1 :]) if x >= A[i]])
> :
> :   return A_lo + [A[i]] + A_hi
> :
> :d = 7
> :k = 9
> :z = d/k
> :f = 1
> :print f
> :S = [[d+k, 1, 1-z]]
> :G = []
> :h = 0
> :while (len(S) > 0) & (not(h == 1)):
> :   G.append(S[0])
> :   h = S[0][2]
> :   del(S[0])
> :   h1 = 1/h-z
> :   h2 = h-z
> :   if (h1 > 0) & (not(([h1.numerator()+h1.denominator(), h, h1] in S) |
> :([h1.numerator()+h1.denominator(), h, h1] in G))):
> :     S.append([h1.numerator()+h1.denominator(), h, h1])
> :   if (h2 > 0) & (not(([h2.numerator()+h2.denominator(), h, h2] in S) |
> :([h2.numerator()+h2.denominator(), h, h2] in G))):
> :     S.append([h2.numerator()+h2.denominator(), h, h2])
> :   S = qsort(S)
> :print
> :G = qsort(G)
> :LG = len(G)
> :if h == 1:
> :   W = []
> :   print "Solved!"
> :   #backtrack
> :   while not(h == 1-z):
> :     i = 0
> :     while not(G[i][2] == h):
> :       i = i+1
> :     #print G[i]
> :     h = G[i][1]
> :     W.append(G[i][2])
> :     del(G[i])
> :   print W
> :print "ready"
> :
> :Matthijs Coster
> :
> :
> :hv at crypt.org schreef op 21-1-2015 om 15:56:
> :> When n == 1 (mod m), m/n solves [ n(n-1)/m, (n-1)/m ].
> :> 3/(6n+2) solves [ n(6n+2), 2n ].
> :> 3/(6n+5) solves [ (4n+3)(6n+5), (n+1)(6n+4), 2n+1 ].
> :>
> :> That's enough to cover all 1/n, 2/n, 3/n.
> :>
> :> I suspect that all rationals r < 1 are solvable, and that it should be
> :> provable. But I've had no luck with 7/9 yet, either.
> :>
> :> Hugo
> :>
> :> Frank Adams-Watters <franktaw at netscape.net> wrote:
> :> :1/(N-1) - 1/(N*(N-1)) = 1/N, so take a = N*(N-1), b=N-1.
> :> :
> :> :Franklin T. Adams-Watters
> :> :
> :> :-----Original Message-----
> :> :From: M. F. Hasler <oeis at hasler.fr>
> :> :To: Sequence Fanatics Discussion list <seqfan at list.seqfan.eu>
> :> :Sent: Tue, Jan 20, 2015 9:38 pm
> :> :Subject: [seqfan] Re: Puzzle
> :> :
> :> :
> :> :Hugo makes a good point :
> :> :the sequence of operations can be coded as a list of integers (m(k),
> :> :k=1..n),
> :> :which mean "add m(k) times r, then take the reciprocal"
> :> :(the last reciprocal being "useless" if it ends in 1).
> :> :
> :> :e.g., for r=1/sqrt(2), one has that operations [1,4,1]
> :> : yield the values [ 1/(1+r)=2-2r ; 1/(2+2r) = 1-r ; 1 ]
> :> :
> :> :for r=1/sqrt(3), one has that [1,3,1] yields the
> :> : values [ 1/(1+r)=3/2-r*3/2 ; 1/(3/2+r*3/2) = 1-r ; 1 ]
> :> :
> :> :E.g., for r = 1/N, one has that [N(N-1), N-1] yields
> :> : values  [ 1/(1+N-1) = 1/N ; 1/(1/N + (N-1)/N) = 1 ]
> :> :
> :> :(Not any 1/N is of the form (a-b)/ab, though, is it?)
> :> :Maximilian
> :> :
> :> :
> :> :> -----Original Message-----
> :> :> From: hv <hv at crypt.org>
> :> :> To: Sequence Fanatics Discussion list <seqfan at list.seqfan.eu>
> :> :> Sent: Tue, Jan 20, 2015 8:24 pm
> :> :> Subject: [seqfan] Re: Puzzle
> :> :>
> :> :>
> :> :> It is never useful to use the reciprocal twice in a row, nor to start
> :> :or
> :> :> end the sequence with it, so we're talking about a sequence [a, b, c,
> :> :... ]
> :> :> of additions, with a single reciprocal between each set of additions.
> :> :>
> :> :> When the sequence is length 1 (ie with no reciprocals), the only
> :> :solution
> :> :> is r = 0; when the sequence is length 2, we have r = 0 or r = (a - b)
> :> :/ ab,
> :> :> with a and b any positive integers.
> :> :>
> :> :> At length 3, we get an ugly quadratic, which simplifies to the same
> :> :case
> :> :> when a + c = b: abcr^2 + (ab - bc)r + (a + c - b) = 0, if I have it
> :> :right.
> :> :> I guess it'll only get uglier for longer sequences.
> :> :>
> :> :> The rationals satisfying (a - b) / ab with a, b in Z+ seem like quite
> :> :> an interesting subset of Q to characterize, though.
> :> :>
> :> :> I'm not sure "no irrational" is correct though, I think for example
> :> :> that 1/sqrt(3) works via the sequence [add, rec, add, add, add, rec,
> :> :add].
> :> :>
> :> :> Hugo
> :> :
> :> :_______________________________________________
> :> :
> :> :Seqfan Mailing list - http://list.seqfan.eu/
> :> :
> :> :
> :> :
> :> :
> :> :_______________________________________________
> :> :
> :> :Seqfan Mailing list - http://list.seqfan.eu/
> :>
> :> _______________________________________________
> :>
> :> Seqfan Mailing list - http://list.seqfan.eu/
> :
> :
> :_______________________________________________
> :
> :Seqfan Mailing list - http://list.seqfan.eu/
>
> _______________________________________________
>
> Seqfan Mailing list - http://list.seqfan.eu/