A329011 -id:A329011 - cp's OEIS frontend

A329012 a(n) = p(1,n), where p(x,n) is the strong divisibility sequence of polynomials based on sqrt(5) as in A327322.

1, 7, 52, 406, 16496, 27664, 1663936, 2081968, 18513664, 833245952, 16665967616, 13888655872, 1666655481856, 8333310963712, 55555495903232, 104166621927424, 16666663803355136, 9259258622967808, 1666666620853682176, 4166666620853682176, 55555555311219638272

Offset: 1

Clark Kimberling, Nov 23 2019

a(n) is a strong divisibility sequence; i.e., gcd(a(h),a(k)) = a(gcd(h,k)). Conjecture: there is no upper bound for the number of consecutive equal digits among numbers in this sequence, as suggested, for example, by 34 straight 1's in a(96) and 38 straight 6's in a(97).

			See Example in A327322.

Cf. A327320, A327321, A329011, A329013.

c[poly_] := If[Head[poly] === Times, Times @@ DeleteCases[(#1 (Boole[MemberQ[#1, x] || MemberQ[#1, y] || MemberQ[#1, z]] &) /@Variables /@ #1 &)[List @@ poly], 0], poly];
r = Sqrt[5]; f[x_, n_] := c[Factor[Expand[(r x + r)^n - (r x - 1/r)^n]]];
Flatten[Table[CoefficientList[f[x, n], x], {n, 1, 12}]];  (* A327322 *)
Table[f[x, n] /. x -> 0, {n, 1, 30}]   (* A329011 *)
Table[f[x, n] /. x -> 1, {n, 1, 30}]   (* A329012 *)
Table[f[x, n] /. x -> 2, {n, 1, 30}]   (* A329013 *)
(* Peter J. C. Moses, Nov 01 2019 *)

A329013 a(n) = p(2,n), where p(x,n) is the strong divisibility sequence of polynomials based on sqrt(5) as in A327322.

Original entry on oeis.org

1, 12, 147, 1836, 116721, 301644, 27679401, 52496748, 704739609, 47763633852, 1436395799961, 1798109838252, 323942200421841, 2430837436077972, 24315999958264707, 68401618078375404, 16418241358998948801, 13682794309260216588, 3694504558135555477881

Offset: 1

Clark Kimberling, Nov 23 2019

a(n) is a strong divisibility sequence; i.e., gcd(a(h),a(k)) = a(gcd(h,k)).

			See Example in A327322.

Cf. A327320, A327321, A329011, A329012.

c[poly_] := If[Head[poly] === Times, Times @@ DeleteCases[(#1 (Boole[MemberQ[#1, x] || MemberQ[#1, y] || MemberQ[#1, z]] &) /@Variables /@ #1 &)[List @@ poly], 0], poly];
r = Sqrt[5]; f[x_, n_] := c[Factor[Expand[(r x + r)^n - (r x - 1/r)^n]]];
Flatten[Table[CoefficientList[f[x, n], x], {n, 1, 12}]];  (* A327322 *)
Table[f[x, n] /. x -> 0, {n, 1, 30}]   (* A329011 *)
Table[f[x, n] /. x -> 1, {n, 1, 30}]   (* A329012 *)
Table[f[x, n] /. x -> 2, {n, 1, 30}]   (* A329013 *)
(* Peter J. C. Moses, Nov 01 2019 *)

cp's OEIS Frontend

A329012 a(n) = p(1,n), where p(x,n) is the strong divisibility sequence of polynomials based on sqrt(5) as in A327322.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica