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.

User: Eric Chen

Eric Chen's wiki page.

Eric Chen has authored 88 sequences. Here are the ten most recent ones:

A347818 Smallest n-digit brilliant number.

Original entry on oeis.org

4, 10, 121, 1003, 10201, 100013, 1018081, 10000043, 100140049, 1000000081, 10000600009, 100000000147, 1000006000009, 10000000000073, 100000380000361, 1000000000000003, 10000001400000049, 100000000000000831, 1000000014000000049, 10000000000000000049, 100000000380000000361
Offset: 1

Author

Eric Chen, Sep 15 2021

Keywords

Comments

A brilliant number is a semiprime (products of two primes, A001358) whose two prime factors have the same number of decimal digits. For an n-digit brilliant number, the two prime factors must each have ceiling(n/2) decimal digits.
Since all brilliant numbers are semiprimes, a(n) >= A098449(n), also, a(n) = A098449(n) for n = 1, 2, 4, 16, 78, ..., are there infinitely many n such that a(n) = A098449(n)?

Examples

			a(6) =    100013 =   103 * 971.
a(7) =   1018081 =  1009 * 1009.
a(8) =  10000043 =  2089 * 4787.
a(9) = 100140049 = 10007 * 10007.
		

Programs

  • Mathematica
    Join[{4,10},Table[Module[{k=1},While[PrimeOmega[10^n+k]!=2||Length[ Union[ IntegerLength/@ FactorInteger[ 10^n+k][[;;,1]]]]!=1,k+=2];10^n+k],{n,2,20}]] (* Harvey P. Dale, Jan 09 2024 *)
  • PARI
    isA078972(n)=my(f=factor(n)); (#f[, 1]==1 && f[1, 2]==2) || (#f[, 1]==2 && f[1, 2]==1 && f[2, 2]==1 && #Str(f[1, 1])==#Str(f[2, 1]))
    A084476(n)=for(k=0,10^n,if(isA078972(10^(2*n-1)+k),return(k)))
    a(n)=if(n%2,nextprime(10^((n-1)/2))^2,10^(n-1)+A084476(n/2)) \\ after Charles R Greathouse IV in A078972

Formula

a(n) = 10^(n-1) + A083289(n).
a(2*n) = 10^(2*n-1) + A084476(n).
a(2*n+1) = A003617(n+1)^2.
a(n) >= A098449(n).

A347821 Smallest prime p such that n*p+1 is a perfect power, or 0 if no such p exists.

Original entry on oeis.org

3, 13, 5, 2, 3, 2801, 5, 3, 7, 50544702849929377, 13, 2, 2, 241, 13, 3, 19, 19, 17, 463, 3, 11, 89, 2, 23, 757, 29, 732541, 31, 917087137, 29, 7, 3
Offset: 1

Author

Eric Chen, Sep 25 2021

Keywords

Comments

For every n, all sufficiently large primes p such that n*p+1 is a perfect power are of the form ((n+1)^q-1)/n with q prime.
a(34) = (35^313-1)/34 is too large to include; it has 482 decimal digits.
a(35) - a(37) = {37, 61, 1483}.
a(38) = (39^349-1)/38 is too large to include; it has 554 decimal digits.
a(39) - a(100) = {5, 2, 43, 3500201, 5, 71, 43, 3851, 178481, 11, 47, 3221, 5, 178250690949465223, 2971, 127, 53, 3, 7, 3541, 61, 2, 59, 2, 61, 17, 3, 751410597400064602523400427092397, 21700501, 4831, 7, 19, 73, 5, 7, 5701, 73, 6007, 79, 39449441, 6481, 19, 79, 48037081, 6218272796370530483675222621221, 2, 3, 438668366137, 89, 5, 23, 331, 89, 654022685443, 11, 1001523179, 97, 3, 792806586866086631668831, 9901, 97, 10303}.
If n*p+1 = m^k, then n*p = m^k-1 = (m-1)*(m^(k-1) + m^(k-2) + ... + m + 1). If p >= n, then m^k = n*p+1 >= n^2+1 > n^2, and we have these three cases: Case 1: m-1 > n, then p can't be prime. Case 2: m-1 = n, this is A084738. Case 3: m-1 < n. If gcd(n, m-1) != m-1, then because m^(k-1) + m^(k-2) + ... + m + 1 > n, p can't be prime. This implies m-1 | n. The three cases means that we only need to check p < n and numbers m such that m-1 | n.
The first numbers n such that a(n) = 0 are {124, 215, 224, 242, ...}. a(268) is unknown; it is the smallest prime of the form (269^q - 1)/268 with prime q if such a prime exists (in which case it must be greater than (269^63659-1)/268), otherwise 0.

Crossrefs

Programs

  • PARI
    a(n)=forprime(p=2,2^32,if(ispower(n*p+1),return(p)))
    
  • PARI
    b(n)=forprime(p=2,2^16,if(ispseudoprime(q=((n+1)^p-1)/n),return(q)))
    a(n)=forprime(p=2,2^30,if(ispower(n*p+1),return(p)));b(n) \\ this program might be incorrect beyond a(300)

Formula

a(n) <= A084738(n+1) if A084738(n+1) > 0.

A347819 Minimal elements for the base-10 representations of the primes greater than 10.

Original entry on oeis.org

11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 227, 251, 257, 277, 281, 349, 409, 449, 499, 521, 557, 577, 587, 727, 757, 787, 821, 827, 857, 877, 881, 887, 991, 2087, 2221, 5051, 5081, 5501, 5581, 5801, 5851, 6469, 6949, 8501
Offset: 1

Author

Eric Chen, Sep 16 2021

Keywords

Comments

Sequence is finite with 77 terms, the largest being 5*10^30 + 27 (which can be written 5(0_28)27, where 0_28 means the string of 28 0's). See text file for proof (this file also has proofs for bases 2, 3, 4, 5, 6, 8, 12).
Minimal elements for the base b representations of the primes > b for other bases b: (see the text file for 9 <= b <= 16) (all written in base b)
b=2: {11}
b=3: {12, 21, 111}
b=4: {11, 13, 23, 31, 221}
b=5: {12, 21, 23, 32, 34, 43, 104, 111, 131, 133, 313, 401, 414, 3101, 10103, 14444, 30301, 33001, 33331, 44441, 300031, 10^95 + 13}
b=6: {11, 15, 21, 25, 31, 35, 45, 51, 4401, 4441, 40041}
b=7: {14, 16, 23, 25, 32, 41, 43, 52, 56, 61, 65, 113, 115, 131, 133, 155, 212, 221, 304, 313, 335, 344, 346, 364, 445, 515, 533, 535, 544, 551, 553, 1022, 1051, 1112, 1202, 1211, 1222, 2111, 3031, 3055, 3334, 3503, 3505, 3545, 4504, 4555, 5011, 5455, 5545, 5554, 6034, 6634, 11111, 11201, 30011, 30101, 31001, 31111, 33001, 33311, 35555, 40054, 100121, 150001, 300053, 351101, 531101, 1100021, 33333301, 5100000001, 33333333333333331} (conjectured, not proven)
b=8: {13, 15, 21, 23, 27, 35, 37, 45, 51, 53, 57, 65, 73, 75, 107, 111, 117, 141, 147, 161, 177, 225, 255, 301, 343, 361, 401, 407, 417, 431, 433, 463, 467, 471, 631, 643, 661, 667, 701, 711, 717, 747, 767, 3331, 3411, 4043, 4443, 4611, 5205, 6007, 6101, 6441, 6477, 6707, 6777, 7461, 7641, 47777, 60171, 60411, 60741, 444641, 500025, 505525, 3344441, 4444477, 5500525, 5550525, 55555025, 444444441, 744444441, 77774444441, 7777777777771, 555555555555525, (10^220-1)/9*40 + 7}.
Equivalently: primes > 10 such that no proper substring (i.e., deleting any positive number of digits) is again a prime > 10. - M. F. Hasler, May 03 2022

Examples

			277 is in this sequence because none of 2, 7, 27, 77 is a prime > 10.
857 is in this sequence because none of 8, 5, 7, 85, 87, 57 is a prime > 10.
991 is in this sequence because none of 9, 1, 99, 91 is a prime > 10.
149 is not in this sequence because 19 is subsequence of 149 and 19 is a prime > 10.
389 is not in this sequence because 89 is subsequence of 389 and 89 is a prime > 10.
439 is not in this sequence because 43 is subsequence of 439 and 43 is a prime > 10.
		

Crossrefs

Cf. A071062 (primes > 10 are not required).
Minimal sets for other sets: A071070 (for composites), A071071 (powers of 2), A071072 (multiples of 4), A071073 (multiples of 3), A111055 (primes of the form 4*n+1), A111056 (primes of the form 4*n+3), A114835 (palindromic primes), A130448 (minimal set of squares).

Programs

  • PARI
    a(n, k, b)=v=[]; for(r=1, length(digits(n, b)), if(r+length(digits(k, 2))-length(digits(n, b))>0 && digits(k, 2)[r+length(digits(k, 2))-length(digits(n, b))]==1, v=concat(v, digits(n, b)[r]))); fromdigits(v, b)
    iss(n, b)=for(k=1, 2^length(digits(n, b))-2, if(ispseudoprime(a(n, k, b)) && a(n, k, b)>b, return(0))); 1
    is(n, b=10)=isprime(n) && n>b && iss(n, b) \\ Test whether n is a minimal element for the base b representations of the primes > b. Default value b = 10 for this sequence.
    select( {is_A347819(n,b=10)=for(L=2, #n=digits(n,b), forvec(d=vector(L, i, [1,#n]), n[d[1]]&& isprime(fromdigits(vecextract(n,d),b))&& return(L==#n), 2))}, [1..8888]) \\ Better select among primes([1,N]). - M. F. Hasler, May 03 2022

Extensions

Edited by M. F. Hasler, May 03 2022

A347773 Square array read by antidiagonals downwards: T(n,k) is the smallest positive integer whose n-th power is the sum of k n-th powers of positive integers, or 0 if no such number exists.

Original entry on oeis.org

1, 2, 1, 3, 5, 1, 4, 3, 0, 1, 5, 2, 6, 0, 1, 6, 4, 7, 422481, 0, 1, 7, 3, 4, 353
Offset: 1

Author

Eric Chen, Sep 15 2021

Keywords

Comments

a(26) = T(5,3) is conjectured to be 0, but this has not been proved.
By Fermat's last theorem, T(n,2) = 0 for n > 2.
Euler's sum of powers conjecture is that T(n,k) = 0 for n > k > 1, but this conjecture is not true: T(4,3) = 422481, T(5,4) = 144, there are no known counterexamples for n >= 6.
There are no known 0s for k > 2.
Conjecture: If T(n,k) = 0, then T(r,k) = T(n,s) = T(r,s) = 0 for all r >= n, 2 <= s <= k.

Examples

			Table begins:
  n\k |  1   2       3    4   5   6     7     8
  ----+----------------------------------------
   1  |  1   2       3    4   5   6     7     8
   2  |  1   5       3    2   4   3     4     4
   3  |  1   0       6    7   4   3     5     2
   4  |  1   0  422481  353   5   3     9    13
   5  |  1   0       ?  144  72  12    23    14
   6  |  1   0       ?    ?   ?   ?  1141   251
   7  |  1   0       ?    ?   ?   ?   568   102
   8  |  1   0       ?    ?   ?   ?     ?  1409
T(2,5) = 4 because 4^2 = 1^2 + 1^2 + 1^2 + 2^2 + 3^2 and there is no smaller square that is the sum of 5 positive squares.
T(4,3) = 422481 because 422481^4 = 95800^4 + 217519^4 + 414560^4 and there is no smaller 4th power that is the sum of 3 positive 4th powers.
T(7,7) = 568 because 568^7 = 127^7 + 258^7 + 266^7 + 413^7 + 430^7 + 439^7 + 525^7 and there is no smaller 7th power that is the sum of 7 positive 7th powers.
		

Crossrefs

Cf. A007666 (main diagonal), A264764 (subdiagonal for k = n-1).
Cf. A175610 and A003828 (both for a(19) = T(4,3) = 422481).
Cf. A003294 and A039664 (both for a(25) = T(4,4) = 353).
Cf. A134341 (for a(33) = T(5,4) = 144).
Cf. A063922 and A063923 (both for a(41) = T(5,5) = 72).
Cf. A130012, A130022 (these two sequences are not rows of this table, since they require DISTINCT n-th powers, but this table does not have that requirement).

Programs

  • PARI
    /* return 0 instead of 1 for n=1, and oo loop when T(n, k)=0 */ A347773(p, n, s, m)={ /* Check whether s can be written as sum of n positive p-th powers not larger than m^p. If so, return the base a of the largest term a^p. */ s>n*m^p && return; n==1&&return(ispower(s, p, &n)*n); /* if s and m are not given, s>=n and m are arbitrary. */ !s&&for(m=round(sqrtn(n, p)), 9e9, A347773(p, n, m^p, m-1)&&return(m)); for(a=ceil(sqrtn(s\n, p)), min(sqrtn(max(0, s-n+1), p), m), A347773(p, n-1, s-a^p, a)&&return(a)); } /* after M. F. Hasler in A007666 */ /* Just enter "A347773(n, k)" to get T(n, k) */

Formula

T(n,1) = 1.
T(1,k) = k.
T(n,2) = 0 for n > 2.
T(n,n) = A007666(n).
T(n,n-1) = A264764(n).
T(3,k) <= A130012(k).
T(4,k) <= A130022(k).

A347772 Square array read by antidiagonals downwards: T(n,k) is the smallest prime p not dividing n such that (p-1) / ord_p(n) = k (n>=2, k>=1), or 0 if no such p exists.

Original entry on oeis.org

3, 7, 2, 43, 11, 0, 113, 67, 3, 2, 251, 13, 0, 11, 11, 31, 41, 17, 13, 19, 2, 1163, 61, 0, 101, 7, 3, 3, 73, 883, 31, 0, 5, 73, 17, 2, 397, 313, 0, 199, 31, 29, 13, 5, 7, 151, 271, 73, 827, 139, 1031, 113, 0, 3, 2, 331, 431, 0, 569, 463, 19, 251, 13, 103, 7, 5, 1753, 5743, 151, 487, 97, 43
Offset: 2

Author

Eric Chen, Sep 14 2021

Keywords

Examples

			Table begins:
  n\k |  1   2    3    4     5    6     7    8    9    10    11    12
  ----+--------------------------------------------------------------
   2  |  3   7   43  113   251   31  1163   73  397   151   331  1753
   3  |  2  11   67   13    41   61   883  313  271   431  5743   193
   4  |  0   3    0   17     0   31     0   73    0   151     0   433
   5  |  2  11   13  101     0  199   827  569  487    31  1453   181
   6  | 11  19    7    5    31  139   463   97   37   101   353   241
   7  |  2   3   73   29  1031   19    43  113  883   311   353  1453
   8  |  3  17   13  113   251    7  1163   89  109   431  1013   577
   9  |  2   5    0   13     0   67     0  313    0    41     0    61
  10  |  7   3  103   53    11   79   211   41   73   281   353    37
  11  |  2   7  193    5   191   19   379  449  199  1301  2531  1549
  12  |  5  23   19   37   271   13    29  193  487    11    89   373
  ...
		

Crossrefs

Row 2: A101208.
Row 3: A101209.
Row 10: A054471.
Column 1: A056619.

Programs

  • PARI
    a(m, n)=forprime(p=2, 2^40, if(gcd(m, p)==1 && znorder(Mod(m, p))==(p-1)/n, return(p)))
    is(m, n)=p=core(m); if(p>1 && p%4==1 && n%p==0 && n%2, return(1)); 0
    A347772(m, n)=if(is(m, n) || (issquare(m) && n%2 && (m%2==0 || n>1)), 0, a(m, n))

Formula

T(n,k) = 0 if n is square, k is odd, n > 1.
T(n,k) = 0 if n is even square, k is odd.
T(n,k) = 0 if (let n' be the squarefree part (A007913) of n) n' == 1 (mod 4), n' > 1, k is divisible by n', k is odd.
T(27,k) = 0 for k == 4 or 8 (mod 12).
T(n,k) == 1 mod k if nonzero.

A347771 Unitary nontotient numbers: values not in range of unitary totient function uphi(n).

Original entry on oeis.org

5, 9, 11, 13, 17, 19, 21, 23, 25, 27, 29, 33, 34, 35, 37, 38, 39, 41, 43, 45, 47, 49, 50, 51, 53, 55, 57, 59, 61, 65, 67, 68, 69, 71, 73, 74, 75, 76, 77, 79, 81, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 98, 99, 101, 103, 105, 107, 109, 110, 111, 113, 114, 115, 117, 118, 119, 121, 122, 123, 125, 129, 131, 133, 134, 135
Offset: 1

Author

Eric Chen, Sep 13 2021

Keywords

Comments

Numbers not appearing in A047994.
Indices of -1 in A135347.
Unitary version of A007617.
This sequence to A047994 is A007617 to A000010.
This sequence to A135347 is A007617 to A049283 (for the case that no such numbers exist, A135347 uses -1 and A049283 uses 0).
All odd numbers not of the form 2^k-1 (i.e. not in A000225) are in this sequence, since uphi(n) = A047994(n) is an even number unless n is a power of 2 (A000079), in this case uphi(n) = n-1.
The intersection of this sequence and A049225 is empty, since for squarefree numbers, all divisors are unitary divisors, note that the intersection of this sequence and A002202 is not empty, the number 110 is in both sequences.

Programs

  • Mathematica
    Select[Range[135], Length[invUPhi[#]] == 0 &] (* Amiram Eldar, Apr 01 2023, using the function invUPhi from A361966 *)
  • PARI
    A047994(n)=my(f=factor(n)~); prod(i=1, #f, f[1, i]^f[2, i]-1)
    is(n)=for(k=1,n^2,if(A047994(k)==n,return(0)));1 \\ after A047994

Formula

A361967(a(n)) = 0. - Amiram Eldar, Apr 01 2023

A347770 Conjectured list of numbers which are perfect, amicable, or sociable.

Original entry on oeis.org

6, 28, 220, 284, 496, 1184, 1210, 2620, 2924, 5020, 5564, 6232, 6368, 8128, 10744, 10856, 12285, 12496, 14264, 14288, 14316, 14536, 14595, 15472, 17296, 17716, 18416, 19116, 19916, 22744, 22976, 31704, 45946, 47616, 48976, 63020, 66928, 66992, 67095, 69615, 71145, 76084, 79750
Offset: 1

Author

Eric Chen, Sep 13 2021

Keywords

Comments

By definition, this is the union of A000396, A259180, and A122726. However, at present A122726 is not known to be complete. There is no proof that 564 (for example) is missing from this sequence. - N. J. A. Sloane, Sep 17 2021
Numbers m for which there exists k>=1 such that s^k(m) = m, where s is A001065.
Conjecture: There are no aliquot cycles containing even numbers and odd numbers simultaneously, i.e., every aliquot cycle either has only even numbers or has only odd numbers.

Examples

			Known aliquot cycles (sorted by smallest member):
{6}
{28}
{220, 284}
{496}
{1184, 1210}
{2620, 2924}
{5020, 5564}
{6232, 6368}
{8128}
{10744, 10856}
{12285, 14595}
{12496, 14288, 15472, 14536, 14264}
{14316, 19116, 31704, 47616, 83328, 177792, 295488, 629072, 589786, 294896, 358336, 418904, 366556, 274924, 275444, 243760, 376736, 381028, 285778, 152990, 122410, 97946, 48976, 45946, 22976, 22744, 19916, 17716}
{17296, 18416}
...
		

Extensions

Edited with new definition (pointing out that the list is only conjectured to be complete) by N. J. A. Sloane, Sep 17 2021

A347769 a(0) = 0; a(1) = 1; for n > 1, a(n) = A001065(a(n-1)) = sigma(a(n-1)) - a(n-1) (the sum of aliquot parts of a(n-1)) if this is not yet in the sequence; otherwise a(n) is the smallest number missing from the sequence.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 16, 15, 13, 14, 17, 18, 21, 19, 20, 22, 23, 24, 36, 55, 25, 26, 27, 28, 29, 30, 42, 54, 66, 78, 90, 144, 259, 45, 33, 31, 32, 34, 35, 37, 38, 39, 40, 50, 43, 41, 44, 46, 47, 48, 76, 64, 63, 49, 51, 52, 53, 56, 57, 58, 59, 60, 108, 172
Offset: 0

Author

Eric Chen, Sep 13 2021

Keywords

Comments

This sequence is a permutation of the nonnegative integers iff Catalan's aliquot sequence conjecture (also called Catalan-Dickson conjecture) is true.
a(563) = 276 is the smallest number whose aliquot sequence has not yet been fully determined.
As long as the aliquot sequence of 276 is not known to be finite or eventually periodic, a(563+k) = A008892(k).

Examples

			a(0) = 0, a(1) = 1;
since A001065(a(1)) = 0 has already appeared in this sequence, a(2) = 2;
since A001065(a(2)) = 1 has already appeared in this sequence, a(3) = 3;
...
a(11) = 11;
since A001065(a(11)) = 1 has already appeared in this sequence, a(12) = 12;
since A001065(a(12)) = 16 has not yet appeared in this sequence, a(13) = A001065(a(12)) = 16;
since A001065(a(13)) = 15 has not yet appeared in this sequence, a(14) = A001065(a(13)) = 15;
since A001065(a(14)) = 9 has already appeared in this sequence, a(15) = 13;
...
		

Crossrefs

Cf. A032451.
Cf. A001065 (sum of aliquot parts).
Cf. A003023, A044050, A098007, A098008: ("length" of aliquot sequences, four versions).
Cf. A007906.
Cf. A115060 (maximum term of aliquot sequences).
Cf. A115350 (termination of the aliquot sequences).
Cf. A098009, A098010 (records of "length" of aliquot sequences).
Cf. A290141, A290142 (records of maximum term of aliquot sequences).
Aliquot sequences starting at various numbers: A000004 (0), A000007 (1), A033322 (2), A010722 (6), A143090 (12), A143645 (24), A010867 (28), A008885 (30), A143721 (38), A008886 (42), A143722 (48), A143723 (52), A008887 (60), A143733 (62), A143737 (68), A143741 (72), A143754 (75), A143755 (80), A143756 (81), A143757 (82), A143758 (84), A143759 (86), A143767 (87), A143846 (88), A143847 (96), A143919 (100), A008888 (138), A008889 (150), A008890 (168), A008891 (180), A203777 (220), A008892 (276), A014360 (552), A014361 (564), A074907 (570), A014362 (660), A269542 (702), A045477 (840), A014363 (966), A014364 (1074), A014365 (1134), A074906 (1521), A143930 (3630), A072891 (12496), A072890 (14316), A171103 (46758), A072892 (1264460).

Programs

  • PARI
    A347769_list(N)=print1(0, ", "); if(N>0, print1(1, ", ")); v=[0, 1]; b=1; for(n=2, N, if(setsearch(Set(v), sigma(b)-b), k=1; while(k<=n, if(!setsearch(Set(v), k), b=k; k=n+1, k++)), b=sigma(b)-b); print1(b, ", "); v=concat(v, b))

A305548 a(n) = 27*n.

Original entry on oeis.org

0, 27, 54, 81, 108, 135, 162, 189, 216, 243, 270, 297, 324, 351, 378, 405, 432, 459, 486, 513, 540, 567, 594, 621, 648, 675, 702, 729, 756, 783, 810, 837, 864, 891, 918, 945, 972, 999, 1026, 1053, 1080, 1107, 1134, 1161, 1188, 1215, 1242, 1269, 1296, 1323, 1350, 1377, 1404, 1431, 1458, 1485, 1512
Offset: 0

Author

Eric Chen, Jun 05 2018

Keywords

Crossrefs

For a(n) = k*n: A001489 (k=-1), A000004 (k=0), A001477 (k=1), A005843 (k=2), A008585 (k=3), A008591 (k=9), A008607 (k=25), A252994 (k=26), this sequence (k=27), A135628 (k=28), A195819 (k=29), A249674 (k=30), A135631 (k=31), A174312 (k=32), A044102 (k=36), A085959 (k=37), A169823 (k=60), A152691 (k=64).

Programs

  • Mathematica
    Range[0,2000,27]
  • PARI
    a(n)=27*n

Formula

a(n) = 27*n.
a(n) = A008585(A008591(n)) = A008591(A008585(n)).
G.f.: 27*x/(x-1)^2.
From Elmo R. Oliveira, Apr 10 2025: (Start)
E.g.f.: 27*x*exp(x).
a(n) = 2*a(n-1) - a(n-2). (End)

A305534 Index of the smallest prime in the n-Fibonacci sequence, or the Lucas U(n,-1) sequence.

Original entry on oeis.org

3, 2, 2, 3, 2, 3, 2, 5, 29, 3, 2, 5, 2, 3, 23, 3, 2, 7, 2, 3, 29, 19, 2, 3, 83, 3, 53, 19, 2, 5, 2, 5, 5, 5479, 71, 3, 2, 17, 11, 3, 2, 37, 2, 31, 5, 11, 2, 5
Offset: 1

Author

Eric Chen, Jun 04 2018

Keywords

Comments

Smallest k such that the k-th Fibonacci polynomial evaluated at x=n is prime. (The first few Fibonacci polynomials are 1, x, x^2 + 1, x^3 + 2*x, x^4 + 3*x^2 + 1, x^5 + 4*x^3 + 3*x, ...)
All terms are primes, since if a divides b, then the a-th term of the n-Fibonacci sequence also divides the b-th term of the n-Fibonacci sequence.
Corresponding primes are 2, 2, 3, 17, 5, 37, 7, 4289, 726120289954448054047428229, 101, 11, 21169, 13, 197, 82088569942721142820383601, 257, 17, 34539049, 19, 401, ...
a(n) = 2 if and only if n is prime.
a(n) = 3 if and only if n^2 + 1 is prime (A005574), except n=2 (since 2 is the only prime p such that p^2 + 1 is also prime).
a(34) > 1024, does a(n) exist for all n >= 1? (However, 17 is the only prime in the first 1024 terms of the 4-Fibonacci sequence, and it seems that 17 is the only prime in the 4-Fibonacci sequence.)
a(35)..a(48) = 71, 3, 2, 17, 11, 3, 2, 37, 2, 31, 5, 11, 2, 5, a(50)..a(54) = 11, 11, 23, 2, 3, a(56) = 3, a(58)..a(75) = 5, 2, 47, 2, 5, 311, 13, 233, 3, 2, 5, 11, 5, 2, 7, 2, 3, 5. Unknown terms a(34), a(49), a(55), a(57), exceed 1024, if they exist.
a(49) > 20000, if it exists. - Giovanni Resta, Jun 06 2018

Crossrefs

Cf. A001605, A096650, A209493, which are the indices of the primes in the n-Fibonacci sequence for n = 1, 2, 3.
Cf. A005478, A086383, A201001, which are the primes in the n-Fibonacci sequence for n = 1, 2, 3.
Cf. A000045, A000129, A006190, A001076, A052918, A005668, A054413, A041025, A099371, A041041, A049666, A041061 (the n-Fibonacci sequence for n = 1 to 12).
Cf. A302990 (for n-step Fibonacci sequence instead of n-Fibonacci sequence).

Programs

  • PARI
    b(n,k)=([n,1;1,0]^k)[1,2]
    a(n)=for(k=1,2^12,if(ispseudoprime(b(n,k)),return(k)))

Extensions

a(34)-a(48) from Giovanni Resta, Jun 06 2018