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.

Showing 1-10 of 21 results. Next

A064768 Sum of primes dividing the repunit numbers with a prime number of digits [A031974] (with repetition).

Original entry on oeis.org

11, 40, 312, 4888, 534888, 265371785, 5365294080, 1111111111111111111, 11111111111111111111111, 77843964391, 57336415063797550469, 2212394296770453025145, 201763709900322803748658482662
Offset: 1

Views

Author

Jason Earls, Oct 18 2001

Keywords

Crossrefs

A068326 a(n) = Product_{k=1..n} A031974(k).

Original entry on oeis.org

1, 11, 1221, 13566531, 15073921825941, 167488020286558453130451, 186097800318379671476024826838541061, 2067753336870885217944964707274937180441685906828771, 2297503707634316908597988192875387398496877151223641424395368232574581
Offset: 0

Views

Author

Jani Melik, Feb 27 2002

Keywords

Crossrefs

Cf. A031974.

Programs

  • PARI
    a(n) = prod(k=1, n, (10^prime(k) - 1)/9); \\ Michel Marcus, Jan 04 2021

Extensions

More terms from Jason Earls, May 19 2002

A090707 Primes whose decimal representation is a valid number in base 4 and interpreted as such is again a prime.

Original entry on oeis.org

2, 3, 11, 13, 23, 31, 101, 103, 113, 131, 211, 223, 233, 311, 331, 1013, 1021, 1033, 1103, 1201, 1213, 1223, 1231, 1301, 2003, 2111, 2113, 2131, 2203, 2213, 2311, 2333, 3001, 3011, 3203, 3221, 3301, 3323, 10111, 10211, 10303, 10313, 10321, 10331
Offset: 1

Views

Author

Cino Hilliard, Jan 18 2004

Keywords

Examples

			13 is prime in decimal and also when considered as a number in base 4: 13 [base 4] = 7 [base 10] which is also prime.
		

Crossrefs

Programs

  • Mathematica
    Select[ FromDigits@# & /@ IntegerDigits[ Prime@ Range@ 270, 4], PrimeQ] (* Robert G. Wilson v, Jan 05 2014 *)
    FromDigits[#,10]&/@Select[Tuples[{0,1,2,3},5],AllTrue[{FromDigits[#,4],FromDigits[#,10]},PrimeQ]&] (* Harvey P. Dale, Jul 30 2021 *)
  • PARI
    forprime(p=2,1e4, if(isprime(t=fromdigits(digits(p,4))), print1(t", "))) \\ Charles R Greathouse IV, Apr 22 2015

Extensions

Name, example and offset corrected by M. F. Hasler, Jan 03 2014

A089971 Primes whose decimal representation also represents a prime in base 2.

Original entry on oeis.org

11, 101, 10111, 101111, 1011001, 1100101, 10010101, 10011101, 10100011, 10101101, 10110011, 11000111, 11100101, 100111001, 101001011, 101101111, 101111011, 101111111, 110111011, 111001001, 1000001011, 1001001011, 1001110111, 1010000011, 1010000111, 1010001101
Offset: 1

Views

Author

Cino Hilliard, Jan 18 2004

Keywords

Comments

See A065720 for the primes given by these terms considered as numbers written in base 2, i.e., the sequence with the definition "working in the opposite sense". - M. F. Hasler, Jan 05 2014
A subsequence of A020449. - M. F. Hasler, Jan 11 2014

Examples

			a(1)=11 is a prime and its decimal representation is also a valid base-2 representation (because all digits are < 2), and 11_2 = 3_10 is again a prime.
		

Crossrefs

Cf. A031974, A089981, A090707, A090708, A090709, A090710, A235394, A235395, A000040 and references therein.

Programs

  • Mathematica
    Select[ FromDigits@# & /@ IntegerDigits[ Prime@ Range@ 270, 2], PrimeQ] (* Robert G. Wilson v, Jan 05 2014 *)
  • PARI
    is_A089971(p)=vecmax(d=digits(p))<2&&isprime(vector(#d, i, 2^(#d-i))*d~)&&isprime(p) \\ "d" is implicitly declared local. Putting isprime(p) to the end improves performance when the function is applied to primes only or to very large numbers. - M. F. Hasler, Jan 05 2014
    
  • PARI
    fixBase(n, oldBase, newBase)=my(d=digits(n, oldBase), t=newBase-1); for(i=1, #d, if(d[i]>t, for(j=i, #d, d[j]=t); break)); fromdigits(d, newBase)
    list(lim)=my(v=List(), t); forprime(p=2, fixBase(lim\1, 10, 2), if(isprime(t=fromdigits(digits(p, 2), 10)), listput(v, t))); Vec(v) \\ Charles R Greathouse IV, Nov 07 2016
    
  • Python
    from sympy import isprime, primerange
    def aupto(limit):
        alst = []
        for p in primerange(2, limit+1):
            t = int(bin(p)[2:])
            if isprime(t): alst.append(t)
        return alst
    print(aupto(2**11)) # Michael S. Branicky, Aug 19 2021

Extensions

Definition and example reworded, offset corrected, and cross-references added by M. F. Hasler, Jan 05 2014

A089981 Primes whose decimal representation also represents a prime in base 3.

Original entry on oeis.org

2, 2111, 2221, 10211, 12011, 12211, 20201, 21011, 21101, 21211, 22111, 101021, 101111, 102101, 102121, 110221, 111121, 111211, 120011, 120121, 121001, 121021, 122011, 201101, 202001, 202021, 210011, 210101, 1000211, 1010201, 1012201
Offset: 1

Views

Author

Cino Hilliard, Jan 18 2004

Keywords

Comments

See A065721 for the primes given by these terms considered as numbers written in base 3, i.e., the sequence with the definition "working in the opposite sense". - M. F. Hasler, Jan 05 2014

Examples

			2111 is a prime and its decimal representation is also a valid base-3 representation (because all digits are < 3), and 2111[3] = 67[10] is again a prime. Therefore 2111 is in the sequence.
		

Crossrefs

Cf. A031974, A089971, A090707, A090708, A090709, A090710, A235394, A235395, A000040 and further references therein.

Programs

  • Mathematica
    Select[ FromDigits@# & /@ IntegerDigits[ Prime@ Range@ 270, 3], PrimeQ] (* Robert G. Wilson v, Jan 05 2014 *)
    FromDigits/@Select[Tuples[{0,1,2},7],AllTrue[{FromDigits[#],FromDigits[ #,3]},PrimeQ]&] (* Harvey P. Dale, Aug 15 2022 *)
  • PARI
    is_A089981(p)=vecmax(d=digits(p))<3&&isprime(vector(#d,i,3^(#d-i))*d~)&&isprime(p) \\ "d" is implicitly declared local. Putting isprime(p) to the end improves performance when the function is applied to primes only, as below, or to very large numbers. - M. F. Hasler, Jan 05 2014
    
  • PARI
    forprime(p=2,1e6,is_A089981(p)&&print1(p",")) \\ - M. F. Hasler, Jan 05 2014
    
  • PARI
    fixBase(n, oldBase, newBase)=my(d=digits(n, oldBase), t=newBase-1); for(i=1, #d, if(d[i]>t, for(j=i, #d, d[j]=t); break)); fromdigits(d, newBase)
    list(lim)=my(v=List(), t); forprime(p=2, fixBase(lim\1, 10, 3), if(isprime(t=fromdigits(digits(p, 3), 10)), listput(v, t))); Vec(v) \\ Charles R Greathouse IV, Nov 07 2016

Extensions

Definition and example reworded, offset corrected and cross-references added by M. F. Hasler, Jan 05 2014

A235394 Primes whose decimal representation is a valid number in base 8 and interpreted as such is again a prime.

Original entry on oeis.org

2, 3, 5, 7, 13, 23, 37, 53, 73, 103, 107, 131, 211, 227, 263, 277, 307, 337, 373, 401, 431, 433, 463, 467, 521, 541, 547, 557, 577, 631, 643, 661, 673, 701, 1013, 1063, 1151, 1153, 1201, 1223, 1327, 1423, 1451, 1453, 1531, 1567, 1613, 1627, 1663, 1721, 2011, 2017
Offset: 1

Views

Author

Robert G. Wilson v, Jan 09 2014

Keywords

Examples

			a(5) = 13_10 = prime(5), 13_8 = 3 + 1*8 = 11_10 = prime(4).
a(8) = 53_10 = prime(16), 53_8 = 3 + 5*8 = 43_10 = prime(14). - _Marius A. Burtea_, Jun 30 2019
		

Crossrefs

Programs

  • Magma
    [n:n in PrimesUpTo(2100)| Max(Intseq(n,10)) le 7 and IsPrime(Seqint(Intseq(Seqint(Intseq(n),8))))]; // Marius A. Burtea, Jun 30 2019
  • Mathematica
    Select[FromDigits@# & /@ IntegerDigits[ Prime@ Range@ 270, 8], PrimeQ]
  • PARI
    fixBase(n,oldBase,newBase)=my(d=digits(n,oldBase),t=newBase-1); for(i=1,#d, if(d[i]>t, for(j=i,#d, d[j]=t); break)); fromdigits(d,newBase)
    list(lim)=my(v=List(),t); forprime(p=2,fixBase(lim\1,10,8), if(isprime(t=fromdigits(digits(p,8),10)), listput(v,t))); Vec(v) \\ Charles R Greathouse IV, Nov 07 2016
    

A235395 Primes whose decimal representation is a valid number in base 9 and interpreted as such is again a prime.

Original entry on oeis.org

2, 3, 5, 7, 41, 47, 67, 131, 151, 241, 331, 337, 461, 557, 601, 641, 661, 751, 757, 827, 887, 1031, 1181, 1217, 1231, 1321, 1327, 1367, 1471, 1637, 1877, 2027, 2081, 2111, 2131, 2207, 2281, 2287, 2351, 2357, 2647, 2731, 2861, 3037, 3121, 3181, 3187, 3307, 3347
Offset: 1

Views

Author

Robert G. Wilson v, Jan 09 2014

Keywords

Crossrefs

Programs

  • Mathematica
    Select[FromDigits@# & /@ IntegerDigits[ Prime@ Range@ 270, 9], PrimeQ]
  • PARI
    fixBase(n, oldBase, newBase)=my(d=digits(n, oldBase), t=newBase-1); for(i=1, #d, if(d[i]>t, for(j=i, #d, d[j]=t); break)); fromdigits(d, newBase)
    list(lim)=my(v=List(), t); forprime(p=2, fixBase(lim\1, 10, 9), if(isprime(t=fromdigits(digits(p, 9), 10)), listput(v, t))); Vec(v) \\ Charles R Greathouse IV, Nov 07 2016

A090709 Primes whose decimal representation is a valid number in base 6 and interpreted as such is again a prime.

Original entry on oeis.org

2, 3, 5, 11, 31, 101, 151, 211, 241, 251, 331, 421, 431, 521, 1021, 1151, 1231, 1321, 2011, 2131, 2311, 2351, 2441, 2531, 2551, 3041, 3221, 3251, 3301, 3541, 4021, 4111, 4201, 4421, 4441, 4451, 5011, 5021, 5101, 5231, 5441, 5531, 10331, 11131, 11311
Offset: 1

Views

Author

Cino Hilliard, Jan 18 2004

Keywords

Examples

			31 is prime in decimal and a valid number in base 6: 31_6 = 19, a prime.
		

Crossrefs

Programs

  • Magma
    [n:n in PrimesUpTo(12000)| Max(Intseq(n,10)) le 5 and IsPrime(Seqint(Intseq(Seqint(Intseq(n),6))))]; // Marius A. Burtea, Jun 30 2019
  • Mathematica
    Select[ FromDigits@# & /@ IntegerDigits[ Prime@ Range@ 270, 6], PrimeQ] (* Robert G. Wilson v, Jan 05 2014 *)
    vn6pQ[n_]:=Module[{idn=IntegerDigits[n]},Max[idn]<6&&PrimeQ[ FromDigits[ idn,6]]]; Select[Prime[Range[1500]],vn6pQ] (* Harvey P. Dale, Jul 10 2015 *)
  • PARI
    fixBase(n, oldBase, newBase)=my(d=digits(n, oldBase), t=newBase-1); for(i=1, #d, if(d[i]>t, for(j=i, #d, d[j]=t); break)); fromdigits(d, newBase)
    list(lim)=my(v=List(), t); forprime(p=2, fixBase(lim\1, 10, 6), if(isprime(t=fromdigits(digits(p, 6), 10)), listput(v, t))); Vec(v) \\ Charles R Greathouse IV, Nov 07 2016
    
  • Python
    from gmpy2 import digits, is_prime
    A090709_list = [n for n in (int(digits(d,6)) for d in range(10**6) if is_prime(d)) if is_prime(n)] # Chai Wah Wu, Apr 09 2016
    

Extensions

Following suggestions by V.J. Pohjola and Donovan Johnson, name, example and offset corrected by M. F. Hasler, Jan 03 2014

A090710 Primes with digits less than 7 whose decimal representation is also a prime when interpreted in base 7.

Original entry on oeis.org

2, 3, 5, 23, 41, 43, 61, 113, 131, 241, 313, 401, 421, 443, 461, 463, 661, 1013, 1033, 1051, 1123, 1231, 1301, 1433, 1453, 1543, 1613, 2111, 2131, 2153, 2203, 2333, 2441, 2531, 2551, 3121, 3163, 3251, 3323, 3433, 3541, 4001, 4111, 4153, 4201, 4241, 4421
Offset: 1

Views

Author

Cino Hilliard, Jan 18 2004

Keywords

Comments

Note that the definition of, e.g., A090714 works "the other way round". - M. F. Hasler, Jan 03 2014

Examples

			23 is a prime and a valid number in base 7, and 23 [base 7] = 17 [base 10] is again a prime.
		

Crossrefs

Programs

  • Mathematica
    Select[ FromDigits@# & /@ IntegerDigits[ Prime@ Range@ 270, 7], PrimeQ] (* Robert G. Wilson v, Jan 05 2014 *)
    FromDigits/@Select[Tuples[{0,1,2,3,4,5,6},4],AllTrue[ {FromDigits[ #], FromDigits[ #,7]}, PrimeQ]&] (* The program uses the AllTrue function from Mathematica version 10 *) (* Harvey P. Dale, Dec 29 2015 *)
  • PARI
    is_A090710(p,b=7)=vecmax(d=digits(p))M. F. Hasler, Jan 03 2014
    
  • PARI
    fixBase(n, oldBase, newBase)=my(d=digits(n, oldBase), t=newBase-1); for(i=1, #d, if(d[i]>t, for(j=i, #d, d[j]=t); break)); fromdigits(d, newBase)
    list(lim)=my(v=List(), t); forprime(p=2, fixBase(lim\1, 10, 7), if(isprime(t=fromdigits(digits(p, 7), 10)), listput(v, t))); Vec(v) \\ Charles R Greathouse IV, Nov 07 2016

Extensions

Name, example and offset corrected by M. F. Hasler, Jan 03 2014

A090708 Primes whose decimal representation is a valid number in base 5 and interpreted as such is again a prime.

Original entry on oeis.org

2, 3, 23, 43, 131, 241, 313, 401, 1123, 1231, 1321, 2111, 2113, 2221, 2311, 3323, 4003, 4241, 4423, 10103, 10301, 10433, 11243, 11423, 12011, 12413, 13331, 14323, 14411, 20113, 20201, 20443, 21011, 21143, 21341, 21433, 22111, 22133, 22441, 23431, 24113, 24421, 24443, 30211, 31223
Offset: 1

Views

Author

Cino Hilliard, Jan 18 2004

Keywords

Examples

			23 is prime when read as base-10 number and also when read as base-5 number, 23 [base 5] = 13 [base 10].
		

Crossrefs

Programs

  • Magma
    [n:n in PrimesUpTo(32000)| Max(Intseq(n,10)) le 4 and IsPrime(Seqint(Intseq(Seqint(Intseq(n),5))))]; // Marius A. Burtea, Jun 30 2019
  • Mathematica
    Select[ FromDigits@# & /@ IntegerDigits[ Prime@ Range@ 270, 5], PrimeQ] (* Robert G. Wilson v, Jan 05 2014 *)
  • PARI
    fixBase(n, oldBase, newBase)=my(d=digits(n, oldBase), t=newBase-1); for(i=1, #d, if(d[i]>t, for(j=i, #d, d[j]=t); break)); fromdigits(d, newBase)
    list(lim)=my(v=List(), t); forprime(p=2, fixBase(lim\1, 10, 5), if(isprime(t=fromdigits(digits(p, 5), 10)), listput(v, t))); Vec(v) \\ Charles R Greathouse IV, Nov 07 2016
    

Extensions

Name, example and offset corrected by M. F. Hasler, Jan 03 2014
More terms from Alejandro J. Becerra Jr., Aug 13 2018
Showing 1-10 of 21 results. Next