[math-fun] Divisibility sequences in OEIS.
Apologies for multiple copies of this rather long message, but I'm keen to find a few potential editors for the following scheme. If you are able to help, send to me personally, not to bother the rest of the lists. I can then send you attachments of even longer files, and we can divvy up the work so that there is only enough duplication to ensure checking. Here is the preface to a file: This is a list, in OEIS numerical order, of divisibility sequences. The suggestion is that these should be given the codeword ``div'' (or ``divi'' if codewords all have four letters). Where this is not already done, it should be remarked that they are divisibility sequences. Many are values of Chebyshev polynomials: for these the polynomial should be quoted. Those which satisfy recurrence relations should include the recurrence and the generating function. A `closed' formula is also always possible. Divisibility sequences satisfying recurrence relations of order 4 or higher may contain prime divisors with more than one rank of apparition. At least the first few examples of this should be noted. Many higher order sequences factor into ones of lower order: this should be noted. Check list for div seqs in (and not yet in) OEIS: 1. Is it stated that it's a divisibility seq ? Should there be a keyword, say `divi' ? 2. Check the offset. The ``right'' offset is $a(0)=0$, but there may be a good reason for having two entries with different offsets. For example, the Fibonacci sequence would normally be given as $a(0)=0$, 1, 1, 2, 3, 5, \ldots, but it often manifests itself as $a(0)=1$, 1, 2, 3, 5, \ldots, e.g., the number of ways of packing $n$ dominoes in a $2\times n$ box. 3. If a divisibility sequence is multiplied or divided through by an integer, it remains a divisibility sequence, so, in order to be sure of locating all manifestations, if a (the) ``natural'' occurrence of the sequence has $a(1)=k$, then the division by $k$, starting $a(0)=0$, $a(1)=1$ should also be given, and often the versions with $a(1)=2$, 3, etc., up to a point, should be given as well. 4. Is there a recurrence relation ? Note that, paradoxically, there may be more than one, and of different orders, because of factorization. 5. Is the generating function given ? 6. Is there a ``closed'' formula of type $A\alpha^n+B\beta^n+C\gamma^n+\cdots$ ? 7. Does the sequence factor into recurring sequences of lower order (which are not necessarily divisibility sequences) over the integers ? over other number fields ? 8. Are the terms values of one or more Chebyshev polynomials ? 9. Note that every second term forms a divisibility sequence which may be worthy of mention in its own right. Similarly for every third term, every fourth term, \ldots . 10. Numbertheoretic properties, e.g., the rank(s) of apparition of primes, may be worthy of note. In sequences of order $2^k$ a prime may have as many as $k$ ranks of apparition. Second order linear recurrences are often numerators or denominators of convergents to continued fractions of quadratic irrationals. Many 2nd order sequences are related to solutions of Bramahgupta-Bhaskara-Pell equations 11. Cross-references to other sequences, especially those arising from 2., 3., 6., 8., \& 9.\ above. 12. Does the sequence arise ``in real life'' as the number of points on an elliptic curve; or the number of perfect matchings in a family of graphs; or the numer of spanning trees (sometimes called the `complexity') thereof, or in other ways? Note that perfect matchings are also called dimer problems or domino tilings or packings. 13. If it's a second order sequence, then I believe that $a_n=F_{rn}/s$ for some $(r,s)$, where $F_n$ is the $n$\,th Fibonacci number. If so, then $r$ and $s$ should be given. 14. Other things that I've missed. For a start, compare and contrast the descriptions of A010892, A000027, A001906, A001353, A004254, A001109, A004187, A001090, A018913, A004189, A004190, A004191, A078362, A007655, A078364, A077412, A078366, A049660, A078368, A075843, A092449, A077421, A097778, A077423, A097780, A097309, A097781, A097311, A097782, A097313, and note that the varied properties that are mentioned are in fact shared by ALL of them, though rarely are more than one or two mentioned. When, as with the above examples, a sequence belongs to a family, it would be good to give the family a name, and say which member it is. The companion sequence of sequences is A000045, A000129, A006190, A001076, A052918, A005668, A054413, A041025, A099371, A041041, A049666, A041061, A140455, A041085, A154597, A041113, missing, A041145, missing, A041181, missing, A041221, missing, A041265, missing, A041313, missing, A041365, A049667, A041421, missing, A041481, missing, A041545, missing, A041613, missing, A041685, missing, A041761, missing, A041841, missing, A041925, missing, A042013, missing, A042105, missing, A042201, missing, A042301, missing, A042405, missing, A042513, missing, A042625, missing, A042741, [perhaps itself a more interesting sequence than many that are submitted?] The major interest of divisibility sequences is the factorization of the terms. It would be good to have this, perhaps on a separate page. In this sense, all of the 'Cunningham project' should be accessible via OEIS. Thankyou for your patience, if you've got as far as this. R.
I'm told that, assuming that it is possible to pick a positive integer at random, that one can derive some interesting (and false) consequences, such as 3/8 = 1/2. Does anyone know of such? Thanks, Bill C.
I don't know about the false consequences, but, for example Persi Diaconis' Ph.D. thesis dealt with probability measures on the natural numbers (once one has a probability measure on the natural numbers, it's easy to turn it into one for all the integers). One can start with trying to define a set A of natural numbers as measurable if the limit lim_{n-->infinity} A(n)/n exists, where A(n) means the number of elements of A which are <=n. When the limit exists this is also the measure. Unfortunately, lots of interesting sets are not measurable -- e.g. the set of natural numbers whose first decimal digit is 1. However, one can enlarge the collection of measurable sets A to those where the following limit exists: lim_{n-->infinity} sum_{k in A, k<=n} 1/k/sum_{k <=n} 1/k One can show that any set measurable the first way ("natural density") is measurable in the second ("harmonic density") and has the same measure. The first digit set above does have harmonic measure of log_{10} 2. There are other ways of enlarging the collection of measurable sets, for example zeta density lim{s --> 1 (from above)) (s-1) sum_{n in A} 1/n^s One can show that any set that has harmonic density or natural density also has zeta density (and is the same) Victor On Fri, Dec 11, 2009 at 3:05 PM, Cordwell, William R <wrcordw@sandia.gov>wrote:
I'm told that, assuming that it is possible to pick a positive integer at random, that one can derive some interesting (and false) consequences, such as 3/8 = 1/2.
Does anyone know of such?
Thanks, Bill C.
_______________________________________________ math-fun mailing list math-fun@mailman.xmission.com http://mailman.xmission.com/cgi-bin/mailman/listinfo/math-fun
This is leading to the question: What is the most slowly converging series? --Gene ________________________________ From: victor miller <victorsmiller@gmail.com> To: math-fun <math-fun@mailman.xmission.com> Sent: Fri, December 11, 2009 12:44:05 PM Subject: Re: [math-fun] Random Numbers I don't know about the false consequences, but, for example Persi Diaconis' Ph.D. thesis dealt with probability measures on the natural numbers (once one has a probability measure on the natural numbers, it's easy to turn it into one for all the integers). One can start with trying to define a set A of natural numbers as measurable if the limit lim_{n-->infinity} A(n)/n exists, where A(n) means the number of elements of A which are <=n. When the limit exists this is also the measure. Unfortunately, lots of interesting sets are not measurable -- e.g. the set of natural numbers whose first decimal digit is 1. However, one can enlarge the collection of measurable sets A to those where the following limit exists: lim_{n-->infinity} sum_{k in A, k<=n} 1/k/sum_{k <=n} 1/k One can show that any set measurable the first way ("natural density") is measurable in the second ("harmonic density") and has the same measure. The first digit set above does have harmonic measure of log_{10} 2. There are other ways of enlarging the collection of measurable sets, for example zeta density lim{s --> 1 (from above)) (s-1) sum_{n in A} 1/n^s One can show that any set that has harmonic density or natural density also has zeta density (and is the same) Victor On Fri, Dec 11, 2009 at 3:05 PM, Cordwell, William R <wrcordw@sandia.gov>wrote:
I'm told that, assuming that it is possible to pick a positive integer at random, that one can derive some interesting (and false) consequences, such as 3/8 = 1/2.
Does anyone know of such?
Thanks, Bill C.
I don't think that there is one. For example you could have a(n) = n * log(n) * log log(n) * .... * log^k (n)^{1+epsilon} where log^k means the k-fold iterated logarithm. Then sum_n 1/a(n) converges (but really really slowly). Victor On Fri, Dec 11, 2009 at 5:29 PM, Eugene Salamin <gene_salamin@yahoo.com>wrote:
This is leading to the question: What is the most slowly converging series?
--Gene
________________________________ From: victor miller <victorsmiller@gmail.com> To: math-fun <math-fun@mailman.xmission.com> Sent: Fri, December 11, 2009 12:44:05 PM Subject: Re: [math-fun] Random Numbers
I don't know about the false consequences, but, for example Persi Diaconis' Ph.D. thesis dealt with probability measures on the natural numbers (once one has a probability measure on the natural numbers, it's easy to turn it into one for all the integers). One can start with trying to define a set A of natural numbers as measurable if the limit lim_{n-->infinity} A(n)/n exists, where A(n) means the number of elements of A which are <=n. When the limit exists this is also the measure. Unfortunately, lots of interesting sets are not measurable -- e.g. the set of natural numbers whose first decimal digit is 1. However, one can enlarge the collection of measurable sets A to those where the following limit exists:
lim_{n-->infinity} sum_{k in A, k<=n} 1/k/sum_{k <=n} 1/k
One can show that any set measurable the first way ("natural density") is measurable in the second ("harmonic density") and has the same measure. The first digit set above does have harmonic measure of log_{10} 2. There are other ways of enlarging the collection of measurable sets, for example zeta density
lim{s --> 1 (from above)) (s-1) sum_{n in A} 1/n^s
One can show that any set that has harmonic density or natural density also has zeta density (and is the same)
Victor
On Fri, Dec 11, 2009 at 3:05 PM, Cordwell, William R <wrcordw@sandia.gov
wrote:
I'm told that, assuming that it is possible to pick a positive integer at random, that one can derive some interesting (and false) consequences, such as 3/8 = 1/2.
Does anyone know of such?
Thanks, Bill C.
_______________________________________________ math-fun mailing list math-fun@mailman.xmission.com http://mailman.xmission.com/cgi-bin/mailman/listinfo/math-fun
NJAS> Here is the preface to a file:
This is a list, in OEIS numerical order, of divisibility sequences.
I confess to having looked up only a few of the 60-odd sequences, but it may come as news to some that there are elliptic divisibility sequences of polynomials. E.g. s[1]:s[2]:1,s[3]:-1,s[4]:x,s[n+4] :=(s[n+3]*s[n+1] + s[n+2]^2)/s[n] gives 1 1 2 1 3 - 1 4 x 5 x + 1 2 6 x - x - 1 3 7 - (x + x + 1) 8 - x (3 x + 2) 5 4 2 9 x - x + 3 x + 3 x + 1 5 4 3 2 10 - (x + 1) (x - x + 3 x + 2 x - 2 x - 1) 7 6 5 4 3 2 11 - (x - 2 x + 3 x + 9 x + 5 x + 3 x + 3 x + 1) 2 6 4 3 2 12 - x (x - x - 1) (x + 2 x + 5 x + 9 x + 9 x + 3) . . . where s[<prime>] is irreducible, etc. --rwg
AhA!! Well done Bill! At last someone is begining to see the light. But one doesn't need to be so esoteric as to go into elliptic divisibility (= Somos??) sequences. You can get all second order linear recurrences (which are divisibility sequences if you take a(0)=0) by tipping up Omar Khayyam's triangle (see A011973 in OEIS): (0) 1 1 1 1 1 2 1 3 1 1 4 3 1 5 6 1 1 6 10 4 1 7 15 10 1 1 8 21 20 5 1 9 28 35 15 1 1 10 36 56 35 6 1 11 45 84 70 21 1 .. .. .. .. .. .. .. and then get a TWO-parameter family by multiplying the columns from right to left by a^0, a^1, a^2, ... and the left-falling diagonals downwards by b^0, b^1, b^2, ..., giving the sequence, which I'll call S(a,b), of polynomials 0, 1, a, a^2 + b, a^3 + 2ab, a^4 + 3a^2* b + b^2, a^5 + 4a^3*b + b^2, a^6 + 5a^4*b + 6a^2*b^2 + b^3, a^7 + 6a^5*b + 10a^3*b^2 + 4b^3, ... which factor into algebraic & primitive parts, just like Bill's. I believe that S(2x,-1) are (one sort of) Chebyshev polynomials. The lists that I gave earlier were for S(k,1) and S(k,-1). All the sequences for all a & b have almost all of the properties mentioned in my checklist, and people are beginning to add more. Every second term of the sequence is S(a^2+2b,-b^2), every third term is S(a^3+3ab,b^3), etc. Note that S(-a,b) = (-1)^(n+1) * S(a,b). Everyone knows that S(1,1) is the Fibs, but it may not be so well known that S(4,1) is F(3n)/2, for example. Can anyone state the theorem behind item 13 in the checklist? Just in case there's anyone rueing having deleted my earlier message, I'd be happy to resend to them personally. I don't want to bother the whole list, but I would like to enlist help from a few keen types to do the enormous amount of spring cleaning that needs to be done. As another example, here is a list of the rather motley collection of present titles of S(k,-2): ...... k=-3 missing (-1)^(n+1) * S(3,-2) k=-2 A108520 see k=2 below. k=-1 A001607 a(n)=-a(n-1)-2a(n-2) k=0 missing 0,1,0,-2,0,4,0,-8,0,16,0,... k=1 A107920 Lucas & Lehmer numbers ... k=2 missing, though is (-1)^(n+1) * A108520, which is Expansion of 1/(1+2x+2x^2) k=3) A000225 Mersenne numbers, repeated as k=3) A168604 No. of ways of partitioning the multiset ... k=4 A007070 a(n)=4a(n-1)-2a(n-2) k=5 A107839 a(n)=5a(n-1)-2a(n-2),a(0)=1,a(1)=5 (though these initial values are ``wrong'') [cf. A005824 and A109165 which each interweave this seq with another. Also A159289, which has curious initial values.] k=6 A154244 a(n)=((3+sqrt7)^n-(3-sqrt7)^n)/2*sqrt7 k=7 missing from here on. ........... Thankyou again for your patience, if you got this far. R. On Sun, 13 Dec 2009, rwg@sdf.lonestar.org wrote:
NJAS> Here is the preface to a file:
This is a list, in OEIS numerical order, of divisibility sequences.
I confess to having looked up only a few of the 60-odd sequences, but it may come as news to some that there are elliptic divisibility sequences of polynomials. E.g. s[1]:s[2]:1,s[3]:-1,s[4]:x,s[n+4] :=(s[n+3]*s[n+1] + s[n+2]^2)/s[n] gives 1 1 2 1 3 - 1 4 x 5 x + 1 2 6 x - x - 1 3 7 - (x + x + 1) 8 - x (3 x + 2) 5 4 2 9 x - x + 3 x + 3 x + 1 5 4 3 2 10 - (x + 1) (x - x + 3 x + 2 x - 2 x - 1) 7 6 5 4 3 2 11 - (x - 2 x + 3 x + 9 x + 5 x + 3 x + 3 x + 1) 2 6 4 3 2 12 - x (x - x - 1) (x + 2 x + 5 x + 9 x + 9 x + 3) . . . where s[<prime>] is irreducible, etc. --rwg _______________________________________________ math-fun mailing list math-fun@mailman.xmission.com http://mailman.xmission.com/cgi-bin/mailman/listinfo/math-fun
This is my second attempt to reply from laughably flaky Squirrelmail. Sorry for any repetitions.
AhA!! Well done Bill! At last someone is begining to see the light. But one doesn't need to be so esoteric as to go into elliptic divisibility (= Somos??) sequences. You can get all second order linear recurrences (which are divisibility sequences if you take a(0)=0) by tipping up Omar Khayyam's triangle (see A011973 in OEIS): (0) 1 1 1 1 1 2 1 3 1 1 4 3 1 5 6 1 1 6 10 4 1 7 15 10 1 1 8 21 20 5 1 9 28 35 15 1 1 10 36 56 35 6 1 11 45 84 70 21 1 .. .. .. .. .. .. ..
and then get a TWO-parameter family by multiplying the columns from right to left by a^0, a^1, a^2, ... and the left-falling diagonals downwards by b^0, b^1, b^2, ..., giving the sequence, which I'll call S(a,b), of polynomials
0, 1, a, a^2 + b, a^3 + 2ab, a^4 + 3a^2* b + b^2, a^5 + 4a^3*b + b^2, a^6 + 5a^4*b + 6a^2*b^2 + b^3, a^7 + 6a^5*b + 10a^3*b^2 + 4b^3, ...
which factor into algebraic & primitive parts, just like Bill's. I believe that S(2x,-1) are (one sort of) Chebyshev polynomials.
The lists that I gave earlier were for S(k,1) and S(k,-1). All the sequences for all a & b have almost all of the properties mentioned in my checklist, and people are beginning to add more.
Every second term of the sequence is S(a^2+2b,-b^2), every third term is S(a^3+3ab,b^3), etc.
Note that S(-a,b) = (-1)^(n+1) * S(a,b).
Everyone knows that S(1,1) is the Fibs, but it may not be so well known that S(4,1) is F(3n)/2, for example. Can anyone state the theorem behind item 13 in the checklist?
Just in case there's anyone rueing having deleted my earlier message, I'd be happy to resend to them personally. I don't want to bother the whole list, but I would like to enlist help from a few keen types to do the enormous amount of spring cleaning that needs to be done.
As another example, here is a list of the rather motley collection of present titles of S(k,-2):
...... k=-3 missing (-1)^(n+1) * S(3,-2) k=-2 A108520 see k=2 below. k=-1 A001607 a(n)=-a(n-1)-2a(n-2) k=0 missing 0,1,0,-2,0,4,0,-8,0,16,0,... k=1 A107920 Lucas & Lehmer numbers ... k=2 missing, though is (-1)^(n+1) * A108520, which is Expansion of 1/(1+2x+2x^2) k=3) A000225 Mersenne numbers, repeated as k=3) A168604 No. of ways of partitioning the multiset ... k=4 A007070 a(n)=4a(n-1)-2a(n-2) k=5 A107839 a(n)=5a(n-1)-2a(n-2),a(0)=1,a(1)=5 (though these initial values are ``wrong'') [cf. A005824 and A109165 which each interweave this seq with another. Also A159289, which has curious initial values.] k=6 A154244 a(n)=((3+sqrt7)^n-(3-sqrt7)^n)/2*sqrt7 k=7 missing from here on. ...........
Thankyou again for your patience, if you got this far. R.
On Sun, 13 Dec 2009, rwg@sdf.lonestar.org wrote:
NJAS> Here is the preface to a file:
This is a list, in OEIS numerical order, of divisibility sequences.
I confess to having looked up only a few of the 60-odd sequences, but it may come as news to some that there are elliptic divisibility sequences of polynomials. E.g. s[1]:s[2]:1,s[3]:-1,s[4]:x,s[n+4] :=(s[n+3]*s[n+1] + s[n+2]^2)/s[n] gives 1 1 2 1 3 - 1 4 x 5 x + 1 2 6 x - x - 1 3 7 - (x + x + 1) 8 - x (3 x + 2) 5 4 2 9 x - x + 3 x + 3 x + 1 5 4 3 2 10 - (x + 1) (x - x + 3 x + 2 x - 2 x - 1) 7 6 5 4 3 2 11 - (x - 2 x + 3 x + 9 x + 5 x + 3 x + 3 x + 1) 2 6 4 3 2 12 - x (x - x - 1) (x + 2 x + 5 x + 9 x + 9 x + 3) . . . where s[<prime>] is irreducible, etc. --rwg
(Sorry about ascribing your words to Neil.) Note that these polynomials satisfy a stronger condition than gcd(s[a],s[b]) = s[gcd(a,b)], namely a|n -> s[a]|s[n]. E.g., 2 5 4 2 18 (x - x - 1) (x - x + 3 x + 3 x + 1) 13 12 11 10 9 8 7 6 5 (x + x + 7 x + 19 x + 25 x + 78 x + 133 x + 108 x + 79 x 4 3 2 + 65 x + 24 x - 6 x - 6 x - 1) 22 21 20 19 18 17 16 19 - (x - 5 x + 15 x + 4 x - 20 x + 108 x + 162 x 15 14 13 12 11 10 9 + 89 x + 294 x + 580 x + 614 x + 378 x - 179 x - 255 x 8 7 6 5 4 3 2 + 870 x + 2011 x + 1897 x + 1044 x + 405 x + 140 x + 45 x + 10 x + 1) 5 4 3 2 20 x (x + 1) (x - x + 3 x + 2 x - 2 x - 1) 18 17 16 15 14 13 12 11 (x - 3 x + 4 x + 3 x + x - 26 x - 12 x + 71 x 10 9 8 7 6 5 4 - 118 x - 469 x - 481 x - 382 x - 554 x - 753 x - 688 x 3 2 - 430 x - 180 x - 45 x - 5). I have been calling these "strong divisibility sequences". What's the accepted term? The degree of these polynomials, 0,0,0,1,1,2,3,2,5,6,7,9,10,12,14,14,18,..., grows nonlinearly, in fact as eight interlaced quadratic progressions. This offers a good OEIS signature, but how to tabulate the actual, jagged, nontriangular coefficient array? --rwg
This is my fifth attempt to reply from laughably flaky Squirrelmail. (Each time with loss of draft and recipient list, and several minutes of unresponsiveness.) Apologies for any duplication.
AhA!! Well done Bill! At last someone is begining to see the light. But one doesn't need to be so esoteric as to go into elliptic divisibility (= Somos??) sequences. You can get all second order linear recurrences (which are divisibility sequences if you take a(0)=0) by tipping up Omar Khayyam's triangle (see A011973 in OEIS): (0) 1 1 1 1 1 2 1 3 1 1 4 3 1 5 6 1 1 6 10 4 1 7 15 10 1 1 8 21 20 5 1 9 28 35 15 1 1 10 36 56 35 6 1 11 45 84 70 21 1 .. .. .. .. .. .. ..
and then get a TWO-parameter family by multiplying the columns from right to left by a^0, a^1, a^2, ... and the left-falling diagonals downwards by b^0, b^1, b^2, ..., giving the sequence, which I'll call S(a,b), of polynomials
0, 1, a, a^2 + b, a^3 + 2ab, a^4 + 3a^2* b + b^2, a^5 + 4a^3*b + b^2, a^6 + 5a^4*b + 6a^2*b^2 + b^3, a^7 + 6a^5*b + 10a^3*b^2 + 4b^3, ...
which factor into algebraic & primitive parts, just like Bill's. I believe that S(2x,-1) are (one sort of) Chebyshev polynomials.
The lists that I gave earlier were for S(k,1) and S(k,-1). All the sequences for all a & b have almost all of the properties mentioned in my checklist, and people are beginning to add more.
Every second term of the sequence is S(a^2+2b,-b^2), every third term is S(a^3+3ab,b^3), etc.
Note that S(-a,b) = (-1)^(n+1) * S(a,b).
Everyone knows that S(1,1) is the Fibs, but it may not be so well known that S(4,1) is F(3n)/2, for example. Can anyone state the theorem behind item 13 in the checklist?
Just in case there's anyone rueing having deleted my earlier message, I'd be happy to resend to them personally. I don't want to bother the whole list, but I would like to enlist help from a few keen types to do the enormous amount of spring cleaning that needs to be done.
As another example, here is a list of the rather motley collection of present titles of S(k,-2):
...... k=-3 missing (-1)^(n+1) * S(3,-2) k=-2 A108520 see k=2 below. k=-1 A001607 a(n)=-a(n-1)-2a(n-2) k=0 missing 0,1,0,-2,0,4,0,-8,0,16,0,... k=1 A107920 Lucas & Lehmer numbers ... k=2 missing, though is (-1)^(n+1) * A108520, which is Expansion of 1/(1+2x+2x^2) k=3) A000225 Mersenne numbers, repeated as k=3) A168604 No. of ways of partitioning the multiset ... k=4 A007070 a(n)=4a(n-1)-2a(n-2) k=5 A107839 a(n)=5a(n-1)-2a(n-2),a(0)=1,a(1)=5 (though these initial values are ``wrong'') [cf. A005824 and A109165 which each interweave this seq with another. Also A159289, which has curious initial values.] k=6 A154244 a(n)=((3+sqrt7)^n-(3-sqrt7)^n)/2*sqrt7 k=7 missing from here on. ...........
Thankyou again for your patience, if you got this far. R.
On Sun, 13 Dec 2009, rwg@sdf.lonestar.org wrote:
NJAS> Here is the preface to a file:
This is a list, in OEIS numerical order, of divisibility sequences.
I confess to having looked up only a few of the 60-odd sequences, but it may come as news to some that there are elliptic divisibility sequences of polynomials. E.g. s[1]:s[2]:1,s[3]:-1,s[4]:x,s[n+4] :=(s[n+3]*s[n+1] + s[n+2]^2)/s[n] gives 1 1 2 1 3 - 1 4 x 5 x + 1 2 6 x - x - 1 3 7 - (x + x + 1) 8 - x (3 x + 2) 5 4 2 9 x - x + 3 x + 3 x + 1 5 4 3 2 10 - (x + 1) (x - x + 3 x + 2 x - 2 x - 1) 7 6 5 4 3 2 11 - (x - 2 x + 3 x + 9 x + 5 x + 3 x + 3 x + 1) 2 6 4 3 2 12 - x (x - x - 1) (x + 2 x + 5 x + 9 x + 9 x + 3) . . . where s[<prime>] is irreducible, etc. --rwg (Sorry about ascribing your words to Neil.) Note that these polynomials satisfy a stronger condition than gcd(s[a],s[b]) = s[gcd(a,b)], namely a|n -> s[a]|s[n]. E.g.,
2 5 4 2 18 (x - x - 1) (x - x + 3 x + 3 x + 1) 13 12 11 10 9 8 7 6 5 (x + x + 7 x + 19 x + 25 x + 78 x + 133 x + 108 x + 79 x 4 3 2 + 65 x + 24 x - 6 x - 6 x - 1) 22 21 20 19 18 17 16 19 - (x - 5 x + 15 x + 4 x - 20 x + 108 x + 162 x 15 14 13 12 11 10 9 + 89 x + 294 x + 580 x + 614 x + 378 x - 179 x - 255 x 8 7 6 5 4 3 2 + 870 x + 2011 x + 1897 x + 1044 x + 405 x + 140 x + 45 x + 10 x + 1) 5 4 3 2 20 x (x + 1) (x - x + 3 x + 2 x - 2 x - 1) 18 17 16 15 14 13 12 11 (x - 3 x + 4 x + 3 x + x - 26 x - 12 x + 71 x 10 9 8 7 6 5 4 - 118 x - 469 x - 481 x - 382 x - 554 x - 753 x - 688 x 3 2 - 430 x - 180 x - 45 x - 5). I have been calling these "strong divisibility sequences". What's the accepted term? The degree of these polynomials, 0,0,0,1,1,2,3,2,5,6,7,9,10,12,14,14,18,..., grows nonlinearly, in fact as eight interlaced quadratic progressions. This offers a good OEIS signature, but how to tabulate the actual, jagged, nontriangular coefficient array? --rwg
participants (5)
-
Cordwell, William R -
Eugene Salamin -
Richard Guy -
rwg@sdf.lonestar.org -
victor miller