[seqfan] Re: A269254

Bob Selcoe rselcoe at entouchonline.net
Mon Oct 23 03:50:18 CEST 2017


Hi Don, Brad and Seqfans,

Don - interesting!  As Brad showed, a(n) = -1 when n = m^2 - 2, m>= 2.  I 
was wondering if there might be other terms as well.

This might help.  m^2 -, m >= 2 == j^2 + 4j + 2, j >= 1.  So for these 
terms:

b(0) = 1;
b(1) = n+1 = j^2 + 4j + 3;
b(2) = n(n+1) - 1 = (j^2 + 4j +2)*(j^2 + 4j +3) - 1 = j^4 + 8j^3 + 21j^2 + 
20j +5

The coefficients of b(0), b(1) and b(2) are the 2nd, 6th and 10th rows of 
the Fibonacci coefficient triangle A011973; but note the degrees of the 
polynomials are in reverse order.  Generally, b(k) is obtained by the 
(4k+2)-th row of the triangle (in reverse order) for all k.  So:

b(3) = j^6 + 12j^5 + 55j^4 + 120j^3 + 126j^2 + 56j + 7;
b(4) = j^8 + 16j^7 + 105j^6 ... + 9;
etc.

These polynomials can be described as the product (c+d)*(c-d), where c(k) = 
(j+2)*c(k-1) - d(k-1) and d(k) = c(k-1). for example,  j=2; n=14:

                                    c         d

b(0)  =          1            1         0
b(1)  =        15           4          1
b(2)  =      209          15        4
b(3) =     2911         56       15
b(4) =   40545       209      56
b(5) = 546719       780    209
etc.

This satisfies that condition that all b(k) are non-prime for all j (unless 
I missed something in Brad's proof, this seems to be an alternate way of 
showing that a(n) = -1 for n = m^2 - 2).

So if a(n) = -1 when n != m^2 - 2, perhaps it can be shown that other 
polynomials correspond to other reverse-order rows in A011973 which generate 
only non-primes?

Regards,
Bob Selcoe


--------------------------------------------------
From: "Don Reble" <djr at nk.ca>
Sent: Sunday, October 22, 2017 4:15 PM
To: "Sequence Fanatics Discussion list" <seqfan at list.seqfan.eu>
Subject: [seqfan] Re: A269254

>> could a(110) be -1?
>
>    It could!
>    Trisect sequence-110 to get:
>
> s(0) -1 = -1 * 1
> s(3) 12209 = 29 * 421
> s(6) 16246150031 = 3191 * 5091241
> s(9) 21618264461738561 = 350981 * 61593831181
> s(12) 28766775971285404815839 = 38604719 * 745162164534481
> s(15) 38279085781688731361830743569 = 4246168109 * 9014971804944317941
> s(18) 50936831077090977385276030140145391
>      = 467039887271 * 109063128151054193913721
> s(21) 67780113009314371791483566295225436698401
>      = 51370141431701 * 1319445715356481833023876701
> s(24) 90192962978053418280696346184682355821321113279
>      = 5650248517599839 * 15962654155319589064868666412961
>
>    So far, this subsequence is the product of a Lucas sequence
>    [t(n+2) = 110 t(n+1) - t(n)] and an integer sequence.
>
>
> s(1) 1 = 1 * 1
> s(4) 1342879 = 139 * 9661
> s(7) 1786928798929 = 15289 * 116876761
> s(10) 2377812544869509551 = 1681651 * 1413975042901
> s(13) 3164083819079723345430241 = 184966321 * 17106269952137521
> s(16) 4210351415532437651518789281919 = 20344613659 * 
> 206951652466984684141
> s(19) 5602588318103384725927427610425725489
>      = 2237722536169 * 2503701074439310756598281
> s(22) 7455196197246420601834317660713681347165711
>      = 246129134364931 * 30289775391615129066341317381
> s(25) 9920405923784291913924768096855946930622570930881
>      = 27071967057606241 * 366445700184058757005286501075041
>
>    ... product of a Lucas sequence [t(n+2) = 110 t(n+1) - t(n)]
>    and an integer sequence.
>
>
> s(2) 111 = 111 * 1
> s(5) 147704481 = 111 * 1330671
> s(8) 196545921732159 = 111 * 1770683979569
> s(11) 261537761671184312049 = 111 * 2356196051091750559
> s(14) 348020453322798282592510671 = 111 * 3135319399304489032364961
> s(17) 463100376622786452935704990267521
>      = 111 * 4172075465070148224645990903311
> s(20) 616233778160295228874631761116689658399
>      = 111 * 5551655659101758818690376226276483409
> s(23) 820003801584096951829983459112209722751529809
>      = 111 * 7387421635892765331801652784794682186950719
> s(26) 1091154458653294057113443794307969480012661481283631
>      = 111 * 9830220348227874388409403552324049369483436768321
>
>    ... multiples of 111. (But there's another factorization.)
>
> -- 
> Don Reble  djr at nk.ca
>
>
> --
> Seqfan Mailing list - http://list.seqfan.eu/
> 


More information about the SeqFan mailing list