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 11-20 of 43 results. Next

A106535 Numbers k such that the smallest x > 1 for which Fibonacci(x) == 0 mod k is x = k - 1.

Original entry on oeis.org

11, 19, 31, 59, 71, 79, 131, 179, 191, 239, 251, 271, 311, 359, 379, 419, 431, 439, 479, 491, 499, 571, 599, 631, 659, 719, 739, 751, 839, 971, 1019, 1039, 1051, 1091, 1171, 1259, 1319, 1399, 1439, 1451, 1459, 1499, 1531, 1559, 1571, 1619, 1759, 1811, 1831
Offset: 1

Views

Author

Peter K. Pearson (ppearson+att(AT)spamcop.net), May 06 2005

Keywords

Comments

This is a sister sequence to A000057, because this sequence, since {k : A001177(k) = k-1}, might be called a subdiagonal sequence of A001177, and {k : A001177(k) = k+1}, which might be called a superdiagonal sequence of A001177. Sequences A000057 and A106535 are disjoint. Is this sequence the set of all divisors of some family of sequences, like A000057 is? - Art DuPre, Jul 11 2012
Are all members of this sequence prime? Using A069106, any composite members must exceed 89151931. - Robert Israel, Oct 13 2015
From Jianing Song, Jul 02 2019: (Start)
Yes, all terms are primes. See a brief proof below.
Also, if p == 1 (mod 4) then b(p) divides (p-Legendre(p,5))/2. So terms in this sequence are congruent to 11 or 19 modulo 20.
Primes p such that ord(-(3+sqrt(5))/2,p) = p-1, where ord(z,p) is the smallest integer k > 0 such that (z^k-1)/p is an algebraic integer. (End)
Comments from Amiram Eldar, Jan 30 2022 (Start)
Sequence A003147, "Primes p with a Fibonacci primitive root", is defined in the paper: Daniel Shanks, Fibonacci primitive roots, Fibonacci Quarterly, Vol. 10, No. 2 (1972), pp. 163-168, and 181.
A second paper on this subject Daniel Shanks and Larry Taylor, An Observation of Fibonacci Primitive Roots, Fibonacci Quarterly, Vol. 11, No. 2 (1973), pp. 159-160,
deals with terms p == 3 (mod 4) of A003147, i.e., the intersection of A003147 and A002145 (or A004767).
It states that if g is a Fibonacci primitive root of a prime p such that p == 3 (mod 4) then g-1 and g-2 are also primitive roots of p.
The first 2000 terms of (A003147 intersect A002145) agree with the present sequence, although the definitions are quite different. Are these two sequences the same? (End)

Crossrefs

Similar sequences that give primes p such that A001177(p) = (p-1)/s: this sequence (s=1), A308795 (s=2), A308796 (s=3), A308797 (s=4), A308798 (s=5), A308799 (s=6), A308800 (s=7),A308801 (s=8), A308802 (s=9).

Programs

  • GAP
    Filtered([2..2000], n -> Fibonacci(n-1) mod n = 0 and Filtered( [2..n-2], x -> Fibonacci(x) mod n = 0 ) = [] );
    
  • Maple
    A106535 := proc(n)
            option remember;
            if n = 1 then
                    11;
            else
                    for a from procname(n-1)+1 do
                            if A001177(a) = a-1 then
                                    return a;
                            end if;
                    end do:
            end if;
    end proc: # R. J. Mathar, Jul 09 2012
    # Alternative:
    fmod:= proc(a,b) local A;
      uses LinearAlgebra[Modular];
      A:= Mod(b, <<1,1>|<1,0>>,integer[8]);
      MatrixPower(b,A,a)[1,2];
    end proc:
    filter:= proc(n)
      local cands;
      if fmod(n-1,n) <> 0 then return false fi;
      cands:= map(t -> (n-1)/t, numtheory:-factorset(n-1));
      andmap(c -> (fmod(c,n) > 0), cands);
    end proc:
    select(filter, [$2..10^4]); # Robert Israel, Oct 13 2015
  • Mathematica
    f[n_] := Block[{x = 2}, While[Mod[Fibonacci@ x, n] != 0, x++]; x];Select[Range@ 1860, f@ # == # - 1 &] (* Michael De Vlieger, Oct 13 2015 *)
  • PARI
    isok(n) = {x = 2; while(fibonacci(x) % n, x++); x == n-1;} \\ Michel Marcus, Oct 20 2015

Formula

{n: A001177(n) = n-1}. - R. J. Mathar, Jul 09 2012

Extensions

Corrected by T. D. Noe, Oct 25 2006

A001578 Smallest primitive prime factor of Fibonacci number F(n), or 1 if F(n) has no primitive prime factor.

Original entry on oeis.org

1, 1, 2, 3, 5, 1, 13, 7, 17, 11, 89, 1, 233, 29, 61, 47, 1597, 19, 37, 41, 421, 199, 28657, 23, 3001, 521, 53, 281, 514229, 31, 557, 2207, 19801, 3571, 141961, 107, 73, 9349, 135721, 2161, 2789, 211, 433494437, 43, 109441, 139, 2971215073, 1103, 97, 101
Offset: 1

Views

Author

Keywords

Comments

A prime factor of F(n) is called primitive if it does not divide F(r) for any r < n.
A Fibonacci number can have more than one primitive factor; the primitive factors of F(19) are 37 and 113.
From Robert Israel, Oct 13 2015: (Start)
Since gcd(F(n),F(k)) = F(gcd(n,k)), the non-primitive prime factors of F(n) are factors of F(k) for some proper divisors k of n.
Since prime p divides F(p-1) if p == 1 or 4 (mod 5), F(p+1) if p == 2 or 3 mod 5, F(p) if p = 5, we have a(n) >= n-1 if a(n) > 1.
a(n) = n-1 iff n=2 or n-1 is in A000057.
a(n) = n+1 iff n+1 is a prime in A106535. (End)

References

  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000045, A000057, A106535, A086597 (number of primitive prime factors in F(n)), A061488 (1's omitted), A262341 (largest primitive prime factor of F(n)).

Programs

  • Maple
    for n from 1 to 350 do
      f:= combinat:-fibonacci(n);
      if not isprime(n) then
        for k in map(t -> n/t, numtheory:-factorset(n)) do
           fk:= combinat:-fibonacci(k);
           g:= igcd(f,fk);
           while g > 1 do
             f:= f/g;
             g:= igcd(f,fk);
           od
        od
      fi;
      if f = 1 then A[n]:= 1; next fi;
      F:= map(t -> t[1],ifactors(f,easy)[2]);
      p:= select(type, F,integer);
      if nops(p) >= 1 then A[n]:= min(p); next fi;
      A[n]:= min(numtheory:-factorset(f));
    od:
    seq(A[i],i=1..350); # Robert Israel, Oct 13 2015
  • Mathematica
    prms={}; Table[f=First/@FactorInteger[Fibonacci[n]]; p=Complement[f, prms]; prms=Join[prms, p]; If[p=={}, 1, First[p]], {n, 50}]
  • PARI
    a(n) = {my(v = vector(n, k, fibonacci(k))); my(vf = vector(n, k, factor(v[k])[,1]~)); for (k=1, n-1, vf[n] = setminus(vf[n], vf[k]);); if (#vf[n], vecmin(vf[n]), 1);} \\ Michel Marcus, May 11 2021

Formula

a(n) = 1 if and only if n = 1, 2, 6, or 12, by Carmichael's theorem. - Jonathan Sondow, Dec 07 2017

Extensions

Edited by T. D. Noe, Apr 15 2004
Definition clarified at the suggestion of Joerg Arndt by Jonathan Sondow, Oct 13 2015

A253208 a(n) = 4^n + 3.

Original entry on oeis.org

4, 7, 19, 67, 259, 1027, 4099, 16387, 65539, 262147, 1048579, 4194307, 16777219, 67108867, 268435459, 1073741827, 4294967299, 17179869187, 68719476739, 274877906947, 1099511627779, 4398046511107, 17592186044419, 70368744177667, 281474976710659
Offset: 0

Views

Author

Vincenzo Librandi, Dec 29 2014

Keywords

Comments

Subsequence of A226807.

Crossrefs

Cf. Numbers of the form k^n+k-1: A000057 (k=2), A168607 (k=3), this sequence (k=4), A242329 (k=5), A253209 (k=6), A253210 (k=7), A253211 (k=8), A253212 (k=9), A253213 (k=10).

Programs

  • Magma
    [4^n+3: n in [0..30]];
    
  • Mathematica
    Table[4^n + 3, {n, 0, 30}] (* or *) CoefficientList[Series[(4 - 13 x) / ((1 - x) (1 - 4 x)), {x, 0, 40}], x]
  • PARI
    a(n)=4^n+3 \\ Charles R Greathouse IV, Oct 07 2015

Formula

G.f.: (4 - 13*x)/((1 - x)*(1 - 4*x)).
a(n) = 5*a(n-1) - 4*a(n-2) for n > 1.
From Elmo R. Oliveira, Nov 14 2023: (Start)
a(n) = 4*a(n-1) - 9 with a(0) = 4.
E.g.f.: exp(4*x) + 3*exp(x). (End)

A064414 Fix a > 0, b > 0, k > 0 and define G_1 = a, G_2 = b, G_k = G_(k-1) + G_(k-2); sequence gives numbers m such that there exists (a, b) where G_k is divisible by m.

Original entry on oeis.org

1, 2, 3, 4, 6, 7, 9, 14, 23, 27, 43, 49, 67, 81, 83, 86, 98, 103, 127, 134, 163, 167, 206, 223, 227, 243, 254, 283, 326, 343, 367, 383, 443, 446, 463, 467, 487, 503, 523, 529, 547, 566, 587, 607, 643, 647, 683, 686, 727, 729, 734, 787, 823, 827, 863, 883, 887
Offset: 1

Views

Author

Naohiro Nomoto, Oct 15 2001

Keywords

Comments

From Logan J. Kleinwaks, Oct 29 2017: (Start)
The squares of this sequence are the squares in A232656.
Conjecture: these are the numbers j such that j^2 = Sum_{d|j} phi(d)*A001177(d), where phi = Euler's totient function (A000010). See A232656. (End)

Examples

			If a = 1, b = 4, then G_k is (1, 4, 5, 9, 14, 23, ...) and no G_k is a multiple of 11. Therefore 11 is not in the sequence.
		

References

  • Teruo Nishiyama, Fibonacci numbers, Suuri-Kagaku, No. 285, March 1987, 67-69, (in Japanese).

Crossrefs

Prime terms are in A000057.

Programs

  • Mathematica
    g[a_, b_, k_] := Fibonacci[k-2]*a + Fibonacci[k-1]*b; ok[n_] := Catch[ Do[ test = Catch[ Do[ If[ Divisible[g[a, b, k], n], Throw[True]], {k, 1, 2*n}]]; If[test == Null, Throw[False]], {a, 1, Floor[Sqrt[n]]}, {b, 1, Floor[Sqrt[n]]}]] ; Reap[ Do[ If[ok[n] == Null, Print[n]; Sow[n]], {n, 1, 1000}]][[2, 1]] (* Jean-François Alcover, Jul 19 2012 *)

Extensions

More terms from David Wasserman, Jul 18 2002
Name edited by David A. Corneth, Oct 30 2017

A079346 Primes p such that F(p-(p/5)) is the first Fibonacci number that p divides.

Original entry on oeis.org

2, 3, 5, 7, 11, 19, 23, 31, 43, 59, 67, 71, 79, 83, 103, 127, 131, 163, 167, 179, 191, 223, 227, 239, 251, 271, 283, 311, 359, 367, 379, 383, 419, 431, 439, 443, 463, 467, 479, 487, 491, 499, 503, 523, 547, 571, 587, 599, 607, 631, 643, 647, 659, 683, 719, 727, 739, 751, 787, 823, 827
Offset: 1

Views

Author

Jon Perry, Jan 04 2003

Keywords

Comments

The n-th prime p is in this sequence iff A001602(n) = p-(5/p) (that is the maximum possible value of A001602(n)).

Examples

			7 belongs to this sequence since (7/5) = -1, F(8) = 21 and 7 does not divide F(1) to F(7).
		

Crossrefs

Union of A000057, A106535 and {5}.

Programs

  • PARI
    forprime (p=2,500, wss=p-kronecker(5,p); for(n=1, wss, if( fibonacci(n)%p==0, if( n==wss, print1(p","), break) ) ))

Extensions

Corrected and edited by Max Alekseyev, Nov 23 2007

A213892 Fixed points of a sequence h(n) defined by the minimum number of 3's in the relation n*[n,3,3,...,3,n] = [x,...,x] between simple continued fractions.

Original entry on oeis.org

2, 7, 19, 31, 47, 67, 71, 83, 151, 163, 167, 223, 227, 271, 307, 331, 359, 379, 431, 463, 479, 487, 499, 631, 643, 683, 691, 743, 787, 811, 839, 863, 947, 967, 1019, 1051, 1087, 1103, 1123, 1163, 1259, 1279, 1307, 1319, 1399, 1423, 1451, 1471
Offset: 1

Views

Author

Art DuPre, Jun 23 2012

Keywords

Comments

In a variant of A213891, multiply n by a number with simple continued fraction [n,3,3,...,3,n] and increase the number of 3's until the continued fraction of the product has the same first and last entry (called x in the NAME). Examples are
2 * [2, 3, 3, 2] = [4, 1, 1, 1, 1, 4],
3 * [3, 3, 3] = [9, 1, 9],
4 * [4, 3, 3, 3, 3, 3, 4] = [17, 4, 1, 2, 1, 4, 17],
5 * [5, 3, 3, 5] = [26, 1, 1, 26],
6 * [6, 3, 3, 3, 3, 3, 6] = [37, 1, 4, 2, 4, 1, 37],
7 * [7, 3, 3, 3, 3, 3, 3, 3, 7] = [51, 8, 2, 1, 2, 8, 51].
The number of 3's needed defines the sequence h(n) = 2, 1, 5, 2, 5, 7, 5, 9, 2, ... (n>=2).
The current sequence contains the fixed points of h, i.e., those n where h(n)=n.
We conjecture that this sequence contains prime numbers analogous to the sequence of prime numbers A000057, in the sense that, instead of referring to the Fibonacci sequences (sequences satisfying f(n) = f(n-1) + f(n-2) with arbitrary positive integer values for f(1) and f(2)) it refers to the sequences satisfying f(n) = 3*f(n-1) + f(n-2), A006190, A003688, A052924, etc. This would mean that a prime is in the sequence if and only if it divides some term in each of the sequences satisfying f(n) = 3*f(n-1) + f(n-2).
The above sequence h() is recorded as A262213. - M. F. Hasler, Sep 15 2015

Crossrefs

Programs

  • Mathematica
    f[m_, n_] := Block[{c, k = 1}, c[x_, y_] := ContinuedFraction[x FromContinuedFraction[Join[{x}, Table[m, {y}], {x}]]]; While[First@ c[n, k] != Last@ c[n, k], k++]; k]; Select[Range[2, 1000], f[3, #] == # &] (* Michael De Vlieger, Sep 16 2015 *)
  • PARI
    {a(n) = local(t, m=1); if( n<2, 0, while( 1,
       t = contfracpnqn( concat([n, vector(m,i,3), n]));
       t = contfrac(n*t[1,1]/t[2,1]);
       if(t[1]
    				

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

A213896 Fixed points of a sequence h(n) defined by the minimum number of 7's in the relation n*[n,7,7,...,7,n] = [x,...,x] between simple continued fractions.

Original entry on oeis.org

2, 3, 19, 31, 67, 79, 103, 127, 139, 151, 167, 179, 191, 263, 283, 359, 383, 443, 463, 479, 491, 503, 571, 631, 691, 787, 827, 883, 919, 1019, 1087, 1171, 1291, 1303, 1307, 1327, 1399, 1423, 1451, 1487
Offset: 1

Views

Author

Art DuPre, Jun 23 2012

Keywords

Comments

In a variant of A213891, multiply n by a number with simple continued fraction [n,7,7,..,7,n] and increase the number of 7's until the continued fraction of the product has the same first and last entry (called x in the NAME). Examples are
2 * [2, 7, 7, 2] = [4, 3, 1, 1, 3, 4],
3 * [3, 7, 7, 7, 3] = [9, 2, 2, 1, 1, 1, 2, 2, 9] ,
4 * [4, 7, 7, 7, 7, 7, 4] = [16, 1, 1, 3, 1, 1, 1, 6, 1, 1, 1, 3, 1, 1, 16],
5 * [5, 7, 7, 5] = [25, 1, 2, 2, 1, 25] ,
6 * [6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 6] = [36, 1, 5, 3, 1, 4, 10, 1, 2, 2, 4, 2, 2, 1, 10, 4, 1, 3, 5, 1, 36],
7 * [7, 7, 7] = [49, 1, 49] .
The number of 7's needed defines the sequence h(n) = 2, 3, 5, 2, 11, 1, 5, 11, 2,... (n>=2).
The current sequence contains the fixed points of h, i. e., those n where h(n)=n.
We conjecture that this sequence contains prime numbers analogous to the sequence of prime numbers A000057, in the sense that, instead of referring to the Fibonacci sequences (sequences satisfying f(n)=f(n-1)+f(n-2) with arbitrary positive integer values for f(1) and f(2)) it refers to the sequences satisfying f(n)=7*f(n-1)+f(n-2), A054413, A015453, etc. This would mean that a prime is in the sequence A213896 if and only if it divides some term in each of the sequences satisfying f(n)=7*f(n-1)+f(n-2).
The above sequence h() is recorded as A262217. - M. F. Hasler, Sep 15 2015

Crossrefs

Programs

  • PARI
    {a(n) = local(t, m=1); if( n<2, 0, while( 1,
       t = contfracpnqn( concat([n, vector(m,i,7), n]));
       t = contfrac(n*t[1,1]/t[2,1]);
       if(t[1]
    				

A213897 Fixed points of a sequence h(n) defined by the minimum number of 8's in the relation n*[n,8,8,...,8,n] = [x,...,x] between simple continued fractions.

Original entry on oeis.org

3, 7, 23, 31, 71, 107, 131, 139, 163, 199, 211, 227, 283, 347, 367, 379, 419, 431, 439, 487, 499, 503, 547, 571, 607, 619, 643, 691, 719, 751, 787, 811, 823, 827, 907, 911, 983, 991, 1031, 1051, 1091, 1151, 1163, 1231, 1303, 1319, 1367, 1399, 1423, 1439, 1459, 1499
Offset: 1

Views

Author

Art DuPre, Jun 24 2012

Keywords

Comments

In a variant of A213891, multiply n by a number with simple continued fraction [n,8,8,..,8,n] and increase the number of 8's until the continued fraction of the product has the same first and last entry (called x in the NAME). Examples are
2 * [2, 8, 2] = [4, 4, 4],
3 * [3, 8, 8, 8, 3] = [9, 2, 1, 2, 2, 2, 1, 2, 9],
4 * [4, 8, 4] = [16, 2, 16],
5 * [5, 8, 8, 5] = [25, 1, 1, 1, 1, 1, 1, 25],
6 * [6, 8, 8, 8, 6] = [36, 1, 2, 1, 4, 1, 2, 1, 36],
7 * [7, 8, 8, 8, 8, 8, 8, 8, 7] = [49, 1, 6, 4, 3, 2, 1, 2, 1, 2, 3, 4, 6, 1, 49].
The number of 8's needed defines the sequence h(n) = 1, 3, 1, 2, 3, 7, 1, 11, 5,.. (n>=2).
The current sequence contains the fixed points of h, i. e., those n where h(n)=n.
We conjecture that this sequence contains prime numbers analogous to the sequence of prime numbers A000057, in the sense that, instead of referring to the Fibonacci sequences (sequences satisfying f(n)=f(n-1)+f(n-2) with arbitrary positive integer values for f(1) and f(2)) it refers to the sequences satisfying f(n)=8*f(n-1)+f(n-2), A041025, A015454, etc. This would mean that a prime is in the sequence A213897 if and only if it divides some term in each of the sequences satisfying f(n)=8*f(n-1)+f(n-2).
The sequence h() is recorded as A262218. - M. F. Hasler, Sep 15 2015

Crossrefs

Programs

  • Mathematica
    f[m_, n_] := Block[{c, k = 1}, c[x_, y_] := ContinuedFraction[x FromContinuedFraction[Join[{x}, Table[m, {y}], {x}]]]; While[First@ c[n, k] != Last@ c[n, k], k++]; k]; Select[Range[2, 1000], f[8, #] == # &] (* Michael De Vlieger, Sep 16 2015 *)
  • PARI
    {a(n) = local(t, m=1); if( n<2, 0, while( 1,
       t = contfracpnqn( concat([n, vector(m,i,8), n]));
       t = contfrac(n*t[1,1]/t[2,1]);
       if(t[1]
    				

A213898 Fixed points of a sequence h(n) defined by the minimum number of 9's in the relation n*[n,9,9,...,9,n] = [x,...,x] between simple continued fractions.

Original entry on oeis.org

2, 11, 31, 43, 47, 67, 79, 103, 127, 199, 211, 223, 263, 307, 311, 383, 431, 439, 463, 467, 499, 523, 563, 571, 587, 691, 719, 751, 811, 839, 863, 883, 911, 967, 991, 1051, 1063, 1087, 1091, 1123, 1151, 1231, 1307, 1327, 1399, 1447, 1451, 1459, 1483, 1499
Offset: 1

Views

Author

Art DuPre, Jun 24 2012

Keywords

Comments

In a variant of A213891, multiply n by a number with simple continued fraction [n,9,9,..,9,n] and increase the number of 9's until the continued fraction of the product has the same first and last entry (called x in the NAME). Examples are
2 * [2, 9, 9, 2] = [4, 4, 1, 1, 4, 4],
3 * [3, 9, 3] = [9, 3, 9],
4 * [4, 9, 9, 9, 9, 9, 4] = [16, 2, 3, 1, 1, 1, 1, 8, 1, 1, 1, 1, 3, 2, 16] ,
5 * [5, 9, 9, 9, 9, 5] = [25, 1, 1, 4, 1, 1, 1, 1, 1, 1, 4, 1, 1, 25],
6 * [6, 9, 9, 9, 9, 9, 6] = [36, 1, 1, 1, 13, 6, 13, 1, 1, 1, 36],
7 * [7, 9, 9, 9, 9, 9, 7] = [49, 1, 3, 3, 6, 1, 6, 3, 3, 1, 49].
The number of 9's needed defines the sequence h(n) = 2, 1,5, 4, 5, 5, 5, 1, 14,... (n>=2).
The current sequence contains the fixed points of h, i. e., those n where h(n)=n.
We conjecture that this sequence contains prime numbers analogous to the sequence of prime numbers A000057, in the sense that, instead of referring to the Fibonacci sequence (sequences satisfying f(n)=f(n-1)+f(n-2) with arbitrary positive integer values for f(1) and f(2)) it refers to the sequences satisfying f(n)=9*f(n-1)+f(n-2) like A099371, A015455 etc. This would mean that a prime is in the sequence A213898 if and only if it divides some term in each of the sequences satisfying f(n)=9*f(n-1)+f(n-2).

Crossrefs

Programs

  • Mathematica
    f[m_, n_] := Block[{c, k = 1}, c[x_, y_] := ContinuedFraction[x FromContinuedFraction[Join[{x}, Table[m, {y}], {x}]]]; While[First@ c[n, k] != Last@ c[n, k], k++]; k]; Select[Range[2, 1000], f[9, #] == # &] (* Michael De Vlieger, Sep 16 2015 *)
  • PARI
    {a(n) = local(t, m=1); if( n<2, 0, while( 1,
       t = contfracpnqn( concat([n, vector(m,i,9), n]));
       t = contfrac(n*t[1,1]/t[2,1]);
       if(t[1]
    				
Previous Showing 11-20 of 43 results. Next