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: }] == {} &]] (* From _Olivier Gérard

}] == {} &]] (* From _Olivier Gérard's wiki page.

}] == {} &]] (* From _Olivier Gérard has authored 3 sequences.

A202286 Smallest prime which is equal to a multiple of its reversal +- a prime smaller than itself in exactly n ways.

Original entry on oeis.org

31, 41, 71, 61, 6421, 8501, 8116001, 845534401
Offset: 1

Author

Claudio Meller, Dec 15 2011, a(3) from Olivier Gérard

Keywords

Comments

The sequence is finite, more specifically there cannot be more than 19 terms. Proof: To have p = k*R(p) -+ q, with q < p, we must have 0 < k = (p +- q) / R(p) < 2p / (p/10) = 20, since the prime p cannot end in 0 and therefore R(p) > p/10. (R(p) and p have the same number of digits.) Thus, for a given prime p, there cannot be more than 19 solutions (k=1..19) to p - k*R(p) = +- q, and therefore no a(n) beyond n=19. - M. F. Hasler, Mar 13 2012, improved following remarks from Hans Havermann, Mar 14 2012
A refined analysis shows that the maximal number is less. On one hand, if R(p) is odd, then only even k can yield a prime. Therefore the maximal number of solutions can only be obtained for p starting with an even digit <= 8, and therefore p/R(p) < 9, thus k < 2p/R(p) < 18, k <= 17. Moreover, R(p) is not a multiple of 3 (since p isn't), therefore 1/3 of the k-values lead to q == 0 (mod 3) and are excluded, which leaves at most 11 possibilities. Other k-values lead to q == 0 (mod 5), unless p starts with '5'. In any case there cannot be more than 9 solutions. - M. F. Hasler, Mar 14 2012
a(9) > 8*10^15, if it exists. - Giovanni Resta, Oct 29 2018

Examples

			a(4)=61 because 61 = 16 x 2 + 29, 61 = 16 x 3 + 13, 61 = 16 x 4 - 3, and 61 = 16 x 5 - 19.
		

Crossrefs

Programs

  • PARI
    A202286(n)={ forprime(p=1,default(primelimit), my(r=A004086(p)); 2*p > n*r & sum(k=1,(2*p-1)\r, isprime(abs(p-k*r)))==n & return(p))}  \\ M. F. Hasler, Mar 14 2012

Extensions

a(7)-a(8) from Hans Havermann, Mar 12 2012

A123161 Numbers of words on successive versions of the Google lookup page.

Original entry on oeis.org

13, 33, 53, 61, 37, 28
Offset: 1

Author

N. J. A. Sloane, based on a suggestion from Olivier Gérard, Jul 09 2008

Keywords

A116987 A variant of A014843 (see Mathematica code).

Original entry on oeis.org

2, 2, 2, 3, 1, 7, 2, 4, 1, 9, 6, 9, 14, 8, 5, 0, 14, 11, 19, 3, 20, 7, 18, 25, 5, 2, 24, 12, 5, 27, 9, 18, 0, 8, 0, 10, 26, 12, 3, 30, 21, 19, 24, 37, 27, 39, 17, 14, 36, 35, 24, 3, 32, 17, 42, 1, 47, 56, 44, 0, 11, 49, 21, 2, 55, 55, 17, 5, 61, 69, 55, 3, 14, 14, 59, 38, 22, 62
Offset: 3

Author

N. J. A. Sloane, based on correspondence from Olivier Gérard, Mar 15 2006

Keywords

Comments

Apparently the present definition of A014843 does not give the terms shown there. The present sequence is a failed attempt to reproduce them.

Programs

  • Mathematica
    Table[ Mod[Plus@@Table[Mod[Plus@@IntegerDigits[n,k],k],{k,2,n-1}],n],{n,3,100}]