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-6 of 6 results.

A175734 Largest n-digit number with 3 divisors.

Original entry on oeis.org

9, 49, 961, 9409, 97969, 994009, 9840769, 99460729, 999002449, 9998200081, 99996985729, 999966000289, 9999995824729, 99999820000081, 999997874844049, 9999997800000121, 99999977853408361, 999999874000003969, 9999999828172082689, 99999999340000001089
Offset: 1

Views

Author

Jaroslav Krizek, Aug 24 2010

Keywords

Comments

Because every number with exactly three divisors must be the square of a prime, the last digit of each term in this sequence is either 1 or 9. - Harvey P. Dale, Aug 18 2011 [with thanks to Ant King and Christopher Hunt Gribble]

Crossrefs

Programs

Formula

a(n) = A132153(n)^2.

Extensions

a(14)-a(20) from Harvey P. Dale, Aug 18 2011

A131581 The next prime greater than the square root of 10^n.

Original entry on oeis.org

2, 5, 11, 37, 101, 317, 1009, 3163, 10007, 31627, 100003, 316241, 1000003, 3162283, 10000019, 31622777, 100000007, 316227767, 1000000007, 3162277669, 10000000019, 31622776621, 100000000003, 316227766069, 1000000000039
Offset: 0

Views

Author

Robert G. Wilson v, Jan 12 2007

Keywords

Comments

The difference between a(n) and floor(sqrt(10^n)): 1, 2, 1, 6, 1, 1, 9, 1, 7, 5, 3, 14, 3, 6, 19, 1, 7, 1, 7, 9, ....
Values for which the difference between a(n) and floor(sqrt(10^n)) equals one: 0, 2, 4, 5, 7, 15, 17, 25, 145, 1079, ..., (1350). Only even terms are 0, 2 & 4; all the rest must be odd.

Crossrefs

Programs

  • Mathematica
    NextPrim[n_] := Block[{k = n + 1}, While[ !PrimeQ[k], k++ ]; k]; Table[ NextPrim@ Floor@ Sqrt[10^n], {n, 0, 25}]
    Table[NextPrime[Sqrt[10^n]],{n,0,30}] (* Harvey P. Dale, Aug 15 2017 *)

Formula

a(n) = sqrt(A175733(n+1)). [From Jaroslav Krizek, Aug 24 2010]

A136582 Sqrt(10)-primes: primes obtained by concatenating the first digits in the decimal expansion of sqrt(10).

Original entry on oeis.org

3, 31, 3162277, 316227766016837933, 316227766016837933199889354443271
Offset: 1

Views

Author

Lekraj Beedassy, Jan 09 2008

Keywords

Comments

n such that floor(sqrt(10^(2*n-1))) is prime (1, 2, 7, 18, 33, ...) are given in A136583.
This sequence is the list of prime terms in A017934.
The next term has 206 digits. - Harvey P. Dale, Dec 06 2023

Crossrefs

Programs

  • Magma
    // by Jason Kimberley, Aug 2011
    for n in [1..499 by 2] do
      f := Isqrt(10^n);
      if IsPrime(f) then
        printf "%o,", f;
      end if;
    end for;
  • Mathematica
    Module[{nn=50,sq10},sq10=RealDigits[Sqrt[10],10,nn][[1]];Select[FromDigits/@Table[Take[sq10,n],{n,nn}],PrimeQ]] (* Harvey P. Dale, Dec 06 2023 *)

Formula

a(n) = A017934(2*A136583(n)-1).

A175733 a(n) is the smallest n-digit number with 3 divisors.

Original entry on oeis.org

4, 25, 121, 1369, 10201, 100489, 1018081, 10004569, 100140049, 1000267129, 10000600009, 100008370081, 1000006000009, 10000033772089, 100000380000361, 1000000025191729, 10000001400000049, 100000000621806289, 1000000014000000049, 10000000055856073561
Offset: 1

Views

Author

Jaroslav Krizek, Aug 24 2010

Keywords

Crossrefs

Programs

  • Mathematica
    Table[Prime[1 + PrimePi[Sqrt[10^n]]]^2, {n, 0, 25}] (* T. D. Noe, Aug 19 2011 *)
  • Python
    from math import sqrt
    from sympy import nextprime
    def a(n): return nextprime(int(sqrt(10**n+1)))**2
    print([a(n) for n in range(20)]) # Michael S. Branicky, Apr 25 2021

Formula

a(n) = A131581(n-1)^2.

A136583 n such that floor(sqrt(10^(2*n-1))) is (probably) prime.

Original entry on oeis.org

1, 2, 7, 18, 33, 206, 468, 1061, 6831, 40377
Offset: 1

Views

Author

Lekraj Beedassy, Jan 09 2008

Keywords

Comments

Number of digits of sqrt(10)-primes (A136582).
The n such that A017934(2*n-1) is (probably) prime.

Crossrefs

Programs

  • Magma
    for n in [1..10^6] do if IsPrime(Isqrt(10^(2*n-1))) then printf "%o, ", n; end if; end for; // Jason Kimberley, Sep 03 2011
  • Mathematica
    rd = RealDigits[Sqrt[10], 10, 10^5][[1]]; Do[ If[ PrimeQ@ FromDigits@ Take[rd, n], Print@n], {n, 10^5}] (* Robert G. Wilson v, Jan 20 2008 *)

Extensions

a(6) - a(8) from Robert G. Wilson v, Jan 20 2008
Probable terms a(9) and a(10) from Jason Kimberley, Aug 19 and Sep 03 2011

A133225 Largest prime <= 2^((n+1)/2).

Original entry on oeis.org

2, 2, 3, 5, 7, 11, 13, 19, 31, 43, 61, 89, 127, 181, 251, 359, 509, 719, 1021, 1447, 2039, 2887, 4093, 5791, 8191, 11579, 16381, 23167, 32749, 46337, 65521, 92681, 131071, 185363, 262139, 370723, 524287, 741431, 1048573, 1482907, 2097143, 2965819
Offset: 1

Views

Author

Anthony C Robin, Jan 03 2008

Keywords

Comments

If one is trying to decide whether an n+1 digit binary number is prime, this is the largest prime for which one needs to test divisibility. For example a six digit number like 110101 must be below 64, so only primes up to 7 are needed to test divisibility. Compare with sequence A132153.

Crossrefs

Cf. A132153.

Programs

Formula

a(n) = A007917[A017910(n+1)]. - R. J. Mathar

Extensions

Showing 1-6 of 6 results.