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

A123317 Smallest prime power m such that n+m is a prime number.

Original entry on oeis.org

1, 1, 2, 1, 2, 1, 4, 3, 2, 1, 2, 1, 4, 3, 2, 1, 2, 1, 4, 3, 2, 1, 8, 5, 4, 3, 2, 1, 2, 1, 16, 5, 4, 3, 2, 1, 4, 3, 2, 1, 2, 1, 4, 3, 2, 1, 32, 5, 4, 3, 2, 1, 8, 5, 4, 3, 2, 1, 2, 1, 256, 5, 4, 3, 2, 1, 4, 3, 2, 1, 2, 1, 16, 5, 4, 3, 2, 1, 4, 3, 2, 1, 128, 5, 4, 3, 2, 1, 8, 7, 16, 5, 4, 3, 2, 1, 4, 3, 2, 1, 2, 1
Offset: 1

Views

Author

Reinhard Zumkeller, Sep 27 2006

Keywords

Examples

			n=23: 23+1=3*2^3, 23+2=5^2, 23+3=13*2, 23+2^2=3^3, 23+5=7*2^2, 23+7=5*3*2, but 23+8=31=A000040(11), therefore a(23)=8;
n=24: 24+1=5^2, 24+2=13*2, 24+3=3^3, 24+2^2=7*2^2, but 24+5=29=A000040(10), therefore a(24)=5;
the smallest occurring proper odd prime power is 9=3^2:
n=118: 118+1=17*7, 118+2=5*3*2^3, 118+3=11^2, 118+2^2=61*2, 118+5=41*3, 118+7=5^3, 118+2^3=7*2*3^2, but 118+3^2=127=A000040(31), therefore a(118)=9.
		

Crossrefs

Programs

  • Maple
    A123317 := proc(n)
        local m ;
        m :=1 ;
        if isprime(n+m) then
            return m ;
        end if;
        for m from 2 do
            if nops(numtheory[factorset](m)) = 1 then
                if isprime(n+m) then
                    return m;
                end if;
            end if;
        end do:
    end proc:
    seq(A123317(n),n=1..102) ; # R. J. Mathar, Aug 09 2019

Formula

A123318(n) = n + a(n);
a(A006093(n)) = 1; a(A040976(n)) = 2 for n>2.

A172183 a(n) is the smallest prime of the form p^q+n, where p and q are prime, or zero if no such prime exists.

Original entry on oeis.org

5, 11, 7, 13, 13, 31, 11, 17, 13, 19, 19, 37, 17, 23, 19, 41, 8209, 43, 23, 29, 29, 31, 31, 73, 29, 53, 31, 37, 37, 79, 0, 41, 37, 43, 43, 61, 41, 47, 43, 67, 73, 67, 47, 53, 53, 71, 79, 73, 53, 59, 59, 61, 61, 79, 59, 83, 61, 67, 67, 109, 0, 71, 67, 73, 73, 191, 71, 193, 73, 79
Offset: 1

Views

Author

Cheng Zhang (cz1(AT)rice.edu), Jan 28 2010, Mar 02 2010

Keywords

Comments

If n mod 6 = 1, both p and q must be 2, and a(n)=0 if n + 4 is not a prime. The values of a(n) for n=257,297,353,383,557 are either greater than 176 000 or 0. Several large entries: a(87) = 2^25633 + 87, a(717) = 2^3217 + 717, a(773) = 2^2539 + 773, a(927) = 2^1117 + 927.

Examples

			a(1)=5 because 5=2^2+1 is the smallest prime of the form p^q+1. a(2)=11 because 11=3^2+2. a(3)=7, because 7=2^2+3. a(17)=8209, because 8209=2^13+17. a(31)=0, because p^q+31 cannot be a prime.
		

Crossrefs

Programs

  • Mathematica
    For[l = {}; n = 1, n <= 70, n++, found = False; If[Mod[n, 2] == 0, For[rm = Infinity; i = 1, i < 100, i++, For[j = 1, j < 100, j++, p = Prime[i]; q = Prime[j]; r = p^q + n; If[r >= rm, Break[], If[PrimeQ[r], rm = r; found = True]]; ]; ], (* if n is odd, r=2^q+n *) If[Mod[n, 6] == 1, r = 4 + n; If[PrimeQ[r], found = True], For[j = 1, j < 1000, j++, q = Prime[j]; r = 2^q + n; If[PrimeQ[r], found = True; rm = r; Break[]]; ]; ]; ]; If[ ! found, rm = 0]; l = Append[l, rm]; ]; l
Showing 1-2 of 2 results.