login
Numbers located at angle turns in a pentagonal spiral.
2

%I #33 Sep 08 2022 08:45:56

%S 1,2,3,4,5,7,9,11,12,14,17,20,23,24,27,31,35,39,40,44,49,54,59,60,65,

%T 71,77,83,84,90,97,104,111,112,119,127,135,143,144,152,161,170,179,

%U 180,189,199,209,219,220,230,241,252,263,264,275,287,299,311,312,324,337,350,363,364,377,391,405,419,420,434,449,464,479,480

%N Numbers located at angle turns in a pentagonal spiral.

%C The link illustrates with three figures:

%C Figure 1 contains the numbers located at angle turns in the pentagonal spiral;

%C Figure 2 contains the primes in the pentagonal spiral;

%C Figure 3 shows a variety of sequences that are associated with the numbers on the lines and diagonals in the pentagonal spiral. For example, the sequence A033537 given by the formula n(2n+5) generates {0, 7, 18, 33, 52, 75, ...} and the corresponding line in the spiral passes through {7, 18, 33, 52, 75, ...}.

%H Vincenzo Librandi, <a href="/A188551/b188551.txt">Table of n, a(n) for n = 1..5000</a>

%H Michel Lagneau, <a href="/A188551/a188551_3.pdf">Illustration of the numbers in the pentagonal spiral</a>

%H <a href="/index/Rec#order_11">Index entries for linear recurrences with constant coefficients</a>, signature (1,0,0,0,2,-2,0,0,0,-1,1).

%F From _R. J. Mathar_, Apr 12 2011: (Start)

%F a(n) = a(n-1) + 2*a(n-5) - 2*a(n-6) - a(n-10) + a(n-11).

%F G.f.: x*(1+x)*(1+x^2)*(x^2-x+1)*(x^3-x-1) / ((x^4+x^3+x^2+x+1)^2*(x-1)^3 ). (End)

%p with(numtheory):

%p T:=array(1..300): k:=1:

%p for n from 1 to 50 do:

%p x1:= 2*n^2 -1: T[k]:=x1:

%p x2:= (n+1)*(2*n-1): T[k+1]:=x2:

%p x3:=2*n^2+2*n-1: T[k+2]:=x3:

%p x4:= 2*n*(n+1): T[k+3]:=x4:

%p x5:=n*(2*n+3): T[k+4]:=x5:

%p k:=k+5:

%p od:

%p for p from 1 to 250 do:

%p z:= T[p]:

%p printf(`%d, `, z):

%p od:

%t CoefficientList[Series[(1 + x) (1 + x^2) (x^2 - x + 1) (x^3 - x - 1) / ((x^4 + x^3 + x^2 + x + 1)^2 (x - 1)^3), {x, 0, 80}], x] (* _Vincenzo Librandi_, Aug 18 2018 *)

%t LinearRecurrence[{1,0,0,0,2,-2,0,0,0,-1,1},{1,2,3,4,5,7,9,11,12,14,17},80] (* _Harvey P. Dale_, Jun 17 2021 *)

%o (Magma) I:=[1,2,3,4,5,7,9,11,12,14,17]; [n le 11 select I[n] else Self(n-1)+2*Self(n-5)-2*Self(n-6)-Self(n-10)+Self(n-11): n in [1..90]]; // _Vincenzo Librandi_, Aug 18 2018

%K nonn,easy

%O 1,2

%A _Michel Lagneau_, Apr 04 2011