cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 21-30 of 32 results. Next

A309593 Numbers k with 4 zeros in a fundamental period of A006190 mod k.

Original entry on oeis.org

5, 10, 13, 25, 26, 29, 37, 41, 50, 58, 65, 73, 74, 82, 89, 97, 109, 125, 130, 137, 145, 146, 149, 157, 169, 178, 181, 185, 193, 194, 197, 205, 218, 229, 233, 241, 250, 269, 274, 281, 290, 293, 298, 314, 317, 325, 338, 349, 353, 362, 365, 370, 373, 377
Offset: 1

Views

Author

Jianing Song, Aug 10 2019

Keywords

Comments

Numbers k such that A322906(k) = 4.
Also numbers k such that A214027(k) is odd.

Crossrefs

Cf. A322907.
Let {x(n)} be a sequence defined by x(0) = 0, x(1) = 1, x(n+2) = m*x(n+1) + x(n). Let w(k) be the number of zeros in a fundamental period of {x(n)} modulo k.
| m=1 | m=2 | m=3
-----------------------------+----------+---------+----------
The sequence {x(n)} | A000045 | A000129 | A006190
The sequence {w(k)} | A001176 | A214027 | A322906
Primes p such that w(p) = 1 | A112860* | A309580 | A309586
Primes p such that w(p) = 2 | A053027 | A309581 | A309587
Primes p such that w(p) = 4 | A053028 | A261580 | A309588
Numbers k such that w(k) = 1 | A053031 | A309583 | A309591
Numbers k such that w(k) = 2 | A053030 | A309584 | A309592
Numbers k such that w(k) = 4 | A053029 | A309585 | this seq
* and also A053032 U {2}

Programs

  • PARI
    for(k=1, 400, if(A322906(k)==4, print1(k, ", ")))

A086598 Number of distinct prime factors in Lucas(n).

Original entry on oeis.org

1, 0, 1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 3, 1, 2, 3, 1, 1, 3, 1, 2, 3, 3, 2, 3, 3, 2, 3, 2, 2, 4, 1, 2, 3, 3, 4, 4, 1, 2, 4, 3, 1, 5, 2, 4, 6, 3, 1, 4, 2, 4, 4, 3, 1, 4, 4, 2, 4, 3, 3, 6, 1, 2, 6, 2, 5, 5, 2, 2, 5, 4, 1, 4, 2, 3, 7, 2, 4, 4, 1, 2, 5, 4, 2, 6, 4, 2, 5, 3, 2, 6, 3, 3, 4, 4, 5, 4, 2, 4, 7, 4, 3, 6, 3, 4, 9
Offset: 0

Views

Author

T. D. Noe, Jul 24 2003

Keywords

Comments

Interestingly, the Lucas numbers separate the primes into three disjoint sets: (A053028) primes that do not divide any Lucas number, (A053027) primes that divide Lucas numbers of even index and (A053032) primes that divide Lucas numbers of odd index.

Crossrefs

Cf. A000204 (Lucas numbers), A086599 (number of prime factors, counting multiplicity), A086600 (number of primitive prime factors).

Programs

  • Magma
    [#PrimeDivisors(Lucas(n)): n in [1..100]]; // Vincenzo Librandi, Jul 26 2017
  • Mathematica
    Lucas[n_] := Fibonacci[n+1] + Fibonacci[n-1]; Table[Length[FactorInteger[Lucas[n]]], {n, 150}]
  • PARI
    a(n)=omega(fibonacci(n-1)+fibonacci(n+1)) \\ Charles R Greathouse IV, Sep 14 2015
    

Formula

a(n) = Sum{d|n and n/d odd} A086600(d) + 1 if 6|n, a Mobius-like transform

Extensions

a(0)=1 prepended by Max Alekseyev, Jun 15 2025

A122487 2 together with odd primes p that divide Fibonacci[(p+1)/2].

Original entry on oeis.org

2, 13, 17, 37, 53, 73, 97, 113, 137, 157, 173, 193, 197, 233, 257, 277, 293, 313, 317, 337, 353, 373, 397, 433, 457, 557, 577, 593, 613, 617, 653, 673, 677, 733, 757, 773, 797, 853, 857, 877, 937, 953, 977, 997, 1013, 1033, 1093, 1097, 1117, 1153, 1193, 1213
Offset: 1

Views

Author

Alexander Adamchuk, Sep 16 2006

Keywords

Comments

Primes of the form 2x^2+2xy+13y^2. Discriminant = -100. - T. D. Noe, May 02 2008
Primes of the form a^2 + b^2 such that a^2 == b^2 (mod 5). - Thomas Ordowski, May 18 2015

Crossrefs

Programs

  • Mathematica
    Select[Prime[Range[1000]],IntegerQ[Fibonacci[(#1+1)/2]/#1]&]
  • PARI
    is(n)=my(k=n%20); (k==13||k==17||k==2) && isprime(n) \\ Charles R Greathouse IV, May 18 2015

Formula

Except for 2, the primes are congruent to {13, 17} (mod 20). - T. D. Noe, May 02 2008
2 together with all primes p == {13, 17} (mod 20). - Thomas Ordowski, May 18 2015

Extensions

Definition changed by T. D. Noe, May 02 2008

A106299 Primes that do not divide any term of the Lucas 3-step sequence A001644.

Original entry on oeis.org

2, 103, 199, 211, 421, 757, 883, 907, 991, 1021, 1123, 1237, 1543, 1567, 1621, 1699, 1753, 1873, 2113, 2539, 2731, 2797, 2803, 3391, 3433, 3463, 3499, 3613, 3631, 3793, 3853, 3919, 4093, 4591, 4723, 4933, 4951, 4987, 5107, 5179, 5527, 5791, 5839, 6073
Offset: 1

Views

Author

T. D. Noe, May 02 2005

Keywords

Comments

If a prime p divides a term a(k) of this sequence, then k must be less than the period of the sequence mod p. Hence these primes are found by computing A001644(k) mod p for increasing k and stopping when either A001644(k) mod p = 0 or the end of the period is reached. Interestingly, for all of these primes except 211, the period of the sequence A001644(k) mod p is (p-1)/d, where d is a small integer. The only other exceptional primes less than 1000000 are 23977 and 47093.

Crossrefs

Cf. A053028 (primes not dividing any Lucas number), A106300 (primes not dividing any Lucas 4-step number), A106301 (primes not dividing any Lucas 5-step number).

Programs

  • Mathematica
    n=3; lst={}; Table[p=Prime[i]; a=Join[Table[ -1, {n-1}], {n}]; a=Mod[a, p]; a0=a; While[s=Mod[Plus@@a, p]; a=RotateLeft[a]; a[[n]]=s; !(a==a0 || s==0)]; If[s>0, AppendTo[lst, p]], {i, 1000}]; lst

A112312 Least index k such that the n-th prime divides the k-th tribonacci number.

Original entry on oeis.org

4, 8, 15, 6, 9, 7, 29, 19, 30, 78, 15, 20, 36, 83, 30, 34, 65, 69, 101, 133, 32, 19, 271, 110, 20, 187, 14, 185, 106, 173, 587, 80, 12, 35, 11, 224, 72, 38, 42, 315, 101, 26, 73, 172, 383, 27, 84, 362, 35, 250, 37, 29, 507, 305, 55, 38, 178, 332, 62, 537, 778, 459, 31, 124
Offset: 1

Views

Author

Jonathan Vos Post, Nov 29 2005

Keywords

Comments

The tribonacci numbers are indexed so that trib(0) = trib(1) = 0, trib(2) = 1, for n>2: trib(n) = trib(n-1) + trib(n-2) + trib(n-3). See A112618 for another version.
Brenner proves that every prime divides some tribonacci number T(n). For the similar 3-step Lucas sequence A001644, there are primes (A106299) that do not divide any term.

Examples

			a(1) = 4 because prime(1) = 2 and tribonacci( 4) = 2.
a(2) = 8 because prime(2) = 3 and tribonacci( 8) = 24 = 3 * 2^3.
a(3) = 15 because prime(3) = 5 and tribonacci(15) = 1705 = 5 *(11 * 31).
a(4) = 6 because prime(4) = 7 and tribonacci( 6) = 7.
a(5) = 9 because prime(5) = 11 and tribonacci( 9) = 44 = 11 * 4.
a(6) = 7 because prime(6) = 13 and tribonacci( 7) = 13.
a(7) = 29 because prime(7) = 17 and tribonacci(29) = 8646064 = 17 *(2^4 * 7 * 19 * 239).
		

Crossrefs

Cf. also A112618 = this sequence minus 1.

Programs

  • Mathematica
    a[0] = a[1] = 0; a[2] = 1; a[n_] := a[n] = a[n - 1] + a[n - 2] + a[n - 3]; f[n_] := Module[{k = 2, p = Prime[n]}, While[Mod[a[k], p] != 0, k++ ]; k]; Array[f, 64] (* Robert G. Wilson v *)

Formula

a(n) = minimum k such that prime(n) | A000073(k) and A000073(k) >= prime(n). a(n) = minimum k such that A000040(n) | A000073(k) and A000073(k) >= A000040(n).

Extensions

Corrected and extended by Robert G. Wilson v, Dec 01 2005

A194363 Lucas entry points: smallest m >= 0 such that the n-th prime divides Lucas(m), or -1 if there is no such m.

Original entry on oeis.org

0, 2, -1, 4, 5, -1, -1, 9, 12, 7, 15, -1, 10, 22, 8, -1, 29, -1, 34, 35, -1, 39, 42, -1, -1, 25, 52, 18, -1, -1, 64, 65, -1, 23, -1, 25, -1, 82, 84, -1, 89, 45, 95, -1, -1, 11, 21, 112, 114, 57, -1, 119, 60, 125, -1, 44, -1, 135, -1, 14, 142, -1, 22, 155, -1
Offset: 1

Views

Author

T. D. Noe, Oct 09 2011

Keywords

Comments

The -1 terms are for the primes in A053028. Note that 2 divides the zeroth Lucas number. In the plots, the uppermost line consists of the odd primes in A000057. Note that when a(n) > 0, then a(n) = A001602(n)/2.

Crossrefs

Cf. A000204 (Lucas numbers), A001602 (Fibonacci entry points), A223486 (Lucas entry points), A000040 (prime numbers).

Programs

  • Mathematica
    lim = 100; luc = LucasL[Range[0, Prime[lim]]]; Table[s = Select[Range[p], Mod[luc[[#]], p] == 0 &, 1]; If[s == {}, -1, s[[1]] - 1], {p, Prime[Range[lim]]}]

Formula

a(n) = A223486(A000040(n)). - Jon Maiga, Jul 01 2021

A223486 Lucas entry points: a(n) = least k such that n divides Lucas number L_k (=A000032(k), for k >= 0), or -1 if there is no such k.

Original entry on oeis.org

0, 0, 2, 3, -1, 6, 4, -1, 6, -1, 5, -1, -1, 12, -1, -1, -1, 6, 9, -1, -1, 15, 12, -1, -1, -1, 18, -1, 7, -1, 15, -1, -1, -1, -1, -1, -1, 9, -1, -1, 10, -1, 22, 15, -1, 12, 8, -1, 28, -1, -1, -1, -1, 18, -1, -1, -1, 21, 29, -1, -1, 15, -1, -1, -1, -1, 34, -1
Offset: 1

Views

Author

Casey Mongoven, Mar 20 2013

Keywords

Comments

If one takes L_k, for k >= 1, that is A000204, then a(1) = 1 and a(2) = 3 followed by the given numbers. This fits then with A106291(n) = A253808(n)*a(n), n >= 1 (where in A253808 a negative entry at position n indicates, as in the present sequence, that the Lucas numbers are not divisible by n. For odd primes not dividing any Lucas numbers see A053028. No power 2^m, m >= 3 divides any Lucas number, see, e.g., Vajda, p. 81). - Wolfdieter Lang, Jan 20 2015

Examples

			a(9) = 6 because L_6 = 18 is the first number in the Lucas sequence (A000032) that 9 divides.
		

References

  • A. Brousseau, Fibonacci and Related Number Theoretic Tables. Fibonacci Association, San Jose, CA, 1972, p. 25.
  • S. Vajda, Fibonacci and Lucas numbers and the Golden Section, Ellis Horwood Ltd., Chichester, 1989.

Crossrefs

Cf. A000032, A000204, A001177, A194363, A053028 (primes not dividing any Lucas numbers), A106291, A253808.

Programs

  • Mathematica
    test[n_] := Module[{a, b, t, cnt = 1}, {a, b} = {2, 1}; While[cnt++; t = b; b = Mod[a + b, n]; a = t; ! (b == 0 || {a, b} == {2, 1})]; If[b == 0, cnt, -1]]; Join[{0, 0}, Table[test[i], {i, Range[3, 100]}]] (* T. D. Noe, Mar 22 2013 *)

Extensions

Edited. Added "k >= 0" in the name and added cross references. - Wolfdieter Lang, Jan 20 2015

A106300 Primes that do not divide any term of the Lucas 4-step sequence A073817.

Original entry on oeis.org

2789, 3847, 4451, 4751, 5431, 6203, 8317, 9533, 9629, 9907, 10093, 11839, 13903, 13907, 14207, 15823, 16319, 16759, 19543, 20939, 21379, 21859, 25303, 26683, 29483, 30871, 31267, 31699, 32003, 32771, 33967, 34963, 36229, 37061, 39983
Offset: 1

Views

Author

T. D. Noe, May 02 2005

Keywords

Comments

If a prime p divides a term a(k) of this sequence, then k must be less than the period of the sequence mod p. Hence these primes are found by computing A073817(k) mod p for increasing k and stopping when either A073817(k) mod p = 0 or the end of the period is reached. Interestingly, for all of these primes, the period of the sequence A073817(k) mod p appears to be (p-1)/d, where d is a small integer.

Crossrefs

Cf. A053028 (primes not dividing any Lucas number), A106299 (primes not dividing any Lucas 3-step number), A106301 (primes not dividing any Lucas 5-step number).

Programs

  • Mathematica
    n=4; lst={}; Table[p=Prime[i]; a=Join[Table[ -1, {n-1}], {n}]; a=Mod[a, p]; a0=a; While[s=Mod[Plus@@a, p]; a=RotateLeft[a]; a[[n]]=s; !(a==a0 || s==0)]; If[s>0, AppendTo[lst, p]], {i, 10000}]; lst

A106301 Primes that do not divide any term of the Lucas 5-step sequence A074048.

Original entry on oeis.org

2, 691, 3163, 4259, 5419, 6637, 6733, 14923, 25111, 27947, 29339, 34123, 34421, 34757, 42859, 55207, 57529, 59693, 61643, 68897, 70249, 75991, 82763, 83177, 85607, 86441, 87103, 93169, 93283, 98573, 106121, 106433, 114847, 129589, 132313
Offset: 1

Views

Author

T. D. Noe, May 02 2005

Keywords

Comments

If a prime p divides a term a(k) of this sequence, then k must be less than the period of the sequence mod p. Hence these primes are found by computing A074048(k) mod p for increasing k and stopping when either A074048(k) mod p = 0 or the end of the period is reached. Interestingly, for all of these primes, the period of the sequence A074048(k) mod p appears to be (p-1)/d, where d is a small integer.

Crossrefs

Cf. A053028 (primes not dividing any Lucas number), A106299 (primes not dividing any Lucas 3-step number), A106300 (primes not dividing any Lucas 4-step number).

Programs

  • Mathematica
    n=5; lst={}; Table[p=Prime[i]; a=Join[Table[ -1, {n-1}], {n}]; a=Mod[a, p]; a0=a; While[s=Mod[Plus@@a, p]; a=RotateLeft[a]; a[[n]]=s; !(a==a0 || s==0)]; If[s>0, AppendTo[lst, p]], {i, 10000}]; lst

A266587 Smallest index of a Lucas number (A000032) that is divisible by prime(n), if it exists, or 0 if it does not exist (for n > 1).

Original entry on oeis.org

0, 2, 0, 4, 5, 0, 0, 9, 12, 7, 15, 0, 10, 22, 8, 0, 29, 0, 34, 35, 0, 39, 42, 0, 0, 25, 52, 18, 0, 0, 64, 65, 0, 23, 0, 25, 0, 82, 84, 0, 89, 45, 95, 0, 0, 11, 21, 112, 114, 57, 0, 119, 60, 125, 0, 44, 0, 135, 0, 14, 142, 0, 22, 155, 0, 0, 55, 0, 58, 87, 0, 179, 184, 0, 189, 192, 0, 0, 50, 102, 209, 0, 215, 0, 219, 222, 112, 0, 23, 232, 234, 239, 244, 245
Offset: 1

Views

Author

Richard R. Forberg, Jan 01 2016

Keywords

Comments

These a(n) values may be called Lucas "entry points".
a(1)=0, corresponding to prime(1)= 2, must be viewed as a special case here, because 0, in this case only, is the correct smallest index.
For all other zeros in a(n) (e.g., for a(3), a(6), a(7) corresponding to prime(3)=5, prime(6)=13, prime(7)=17), there are no Lucas numbers divisible by these primes. The full set of primes not divisible into any Lucas number are given by A053028.
For a(n) > 0, a(n) is always equal to either (prime(n)-1)/k or (prime(n)+1)/k, for some k >= 1. This is similar to the Fibonacci entry points given by A001602.
For each value of a(n) > 0 there is an infinite, periodic subsequence within the Lucas numbers divisible by prime(n), given by Lucas(a(n) + 2*i*a(n)), for all i >= 0. Again this is analogous to the A001602 for Fibonacci. The periodicity of the subsequence is 2*a(n), twice the entry point vs. equal to the entry point for Fibonacci.
Conjecture: Infinite Lucas subsequences divisible by the powers of odd primes, p(n), for which a(n) > 0, are given by:
Lucas(a(n) + a(n)*(p(n)-1)*(Sum_{j=1..m-1} p(n)^(j-1)) + 2*i*a(n)*p(n)^(m-1)) is divisible by p(n)^m, where p(n) > 2, i >= 0, m > 1.
Note: the formula above also works for m=1 if the "Sum" is assumed to be zero when the upper summation index limit is less than initial summation index. See second Mathematica example below which works in this way for all m >= 1, to demonstrate the rule for p = p(n) with corresponding a = a(n).
The divisibility of Lucas numbers by powers of 2 is limited to 2^1 and 2^2, as follows: Lucas(3*i) is divisible by 2 and Lucas(3+6i) is divisible by 4, for all i >= 0.

Examples

			For prime(10) = 29, we get a(10) = 7, because Lucas(7)= 29 is the first Lucas number divisible by 29. Also note 7 = (29-1)/4.
For prime(11) = 31, we get a(11) = 15, because Lucas(15) = 1364 is the first Lucas number divisible by 31. Also note 15 = (31-1)/2.
		

Crossrefs

Essentially the same as A194363.

Programs

  • Mathematica
    result={}; Do[iresult=0; Do[If[Divisible[LucasL[i], Prime[k]], iresult=i; Break[]], {i, 1, 2000}]; AppendTo[result, iresult], {k, 2, 200}]; result
    p = 23; a = 12; m = 4; Table[Divisible[LucasL[a + a*(p - 1)*Sum[p^(j - 1), {j, 1, m - 1}] + 2a*i*p^(m - 1)], p^m], {i, 1, 100}]
Previous Showing 21-30 of 32 results. Next