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

A023514 a(n) = sum of exponents in prime-power factorization of prime(n) + 1.

Original entry on oeis.org

1, 2, 2, 3, 3, 2, 3, 3, 4, 3, 5, 2, 3, 3, 5, 4, 4, 2, 3, 5, 2, 5, 4, 4, 3, 3, 4, 5, 3, 3, 7, 4, 3, 4, 4, 4, 2, 3, 5, 3, 5, 3, 7, 2, 4, 5, 3, 6, 4, 3, 4, 6, 3, 5, 3, 5, 5, 5, 2, 3, 3, 4, 4, 5, 2, 3, 3, 3, 4, 4, 3, 6, 5, 3, 4, 8, 4, 2, 3, 3, 5, 2, 7, 3, 5, 4, 5, 2, 4, 5, 5, 7, 4, 4, 5, 6, 4, 4, 3
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Maple
    with(numtheory):
    a:= n-> add(i[2], i=ifactors(ithprime(n)+1)[2]):
    seq(a(n), n=1..100);
  • Mathematica
    Array[Plus@@Last/@FactorInteger[Prime[ # ]+1]&,6! ] (* Vladimir Joseph Stephan Orlovsky, Feb 28 2010 *)
  • PARI
    a(n) = bigomega(prime(n)+1); \\ Michel Marcus, Jan 04 2016

Formula

a(n) = A001222(A000040(n)+1). - Zak Seidov, Jan 04 2016

A210936 Sum of prime factors of prime(n)-1 (counted with multiplicity).

Original entry on oeis.org

0, 2, 4, 5, 7, 7, 8, 8, 13, 11, 10, 10, 11, 12, 25, 17, 31, 12, 16, 14, 12, 18, 43, 17, 13, 14, 22, 55, 13, 15, 15, 20, 23, 28, 41, 15, 20, 14, 85, 47, 91, 15, 26, 15, 18, 19, 17, 42, 115, 26, 35, 26, 16, 17, 16, 133, 71, 16, 30, 18, 52, 77, 25, 38, 22, 83
Offset: 1

Views

Author

Paolo P. Lava, Mar 30 2012

Keywords

Comments

From an idea of Michael B. Porter.

Examples

			prime(10) = 29, and 29-1 = 28 = 2*2*7, so a(10) = 2+2+7 = 11.
		

Crossrefs

Programs

  • Maple
    with(numtheory);
    P:=proc(i)
    local a,k,n;
    for n from 1 to i do
      a:=ifactors(ithprime(n)-1)[2]; print(add(a[k][1]*a[k][2],k=1..nops(a)));
    od; end:
    # alternative
    A210936 := proc(n)
            local p,pplus,f ;
            p := ithprime(n) ;
            pplus := ifactors(p-1)[2] ;
            add(op(1,f)*op(2,f),f=pplus) ;
    end proc:
    seq(A210936(n),n=1..300) ; # R. J. Mathar, May 25 2022

Formula

a(n) = A001414(A006093(n)). - Michel Marcus, Oct 05 2013

A127305 Primes p such that p + (sum of the prime factors of p-1) + (sum of prime factors of p+1) is prime.

Original entry on oeis.org

2, 13, 47, 71, 101, 151, 193, 239, 241, 293, 331, 337, 359, 383, 397, 401, 421, 463, 487, 557, 577, 709, 773, 797, 821, 929, 1019, 1031, 1069, 1093, 1103, 1181, 1217, 1249, 1327, 1367, 1423, 1499, 1571, 1759, 1787, 1789, 1831, 1871, 1877, 1913, 1933, 2053
Offset: 1

Views

Author

J. M. Bergot, Mar 28 2007

Keywords

Comments

The primes are taken "with multiplicity".

Examples

			2 is a term, since 2 + 0 + 3 = 5 is a prime.
13 is a term since 13 + (2+2+3) + (2+7) = 29 is prime, i.e. the prime factors are added with multiplicity.
151 is prime, 150 = 2*3*5*5, 152 = 2*2*2*19. 151 + 2+3+5+5 + 2+2+2+19 = 191 is prime, hence 151 is a term.
		

Crossrefs

Programs

  • Magma
    [ p: p in PrimesInInterval(3, 2100) | IsPrime(&+[ &+[ k[1]*k[2]: k in Factorization(n)]: n in [p-1..p+1] ] ) ]; /* Klaus Brockhaus, Apr 06 2007 */
  • Maple
    spf:= proc(n) local t; add(t[1]*t[2], t = ifactors(n)[2]) end proc:
    filter:= proc(p) isprime(p) and isprime(p+spf(p-1)+spf(p+1)) end proc:
    select(filter, [$2..10000]); # Robert Israel, May 25 2022
  • Mathematica
    pspfQ[n_] := PrimeQ[n + Total[Flatten[Table[#[[1]], {#[[2]]}] & /@ Flatten[FactorInteger[n + {1,-1}], 1] ] ] ]; {2}~Join~Select[Prime[Range[400]], pspfQ] (* Harvey P. Dale, Jan 08 2015, corrected by Michael De Vlieger, May 25 2022 *)

Extensions

Edited and extended by Klaus Brockhaus, Apr 06 2007
Edited by N. J. A. Sloane, May 25 2022

A343416 a(n) = A001414(n) + A000203(n) + A001414(A000203(n)) + A000203(A001414(n)).

Original entry on oeis.org

1, 11, 15, 25, 22, 30, 29, 41, 44, 41, 42, 54, 50, 55, 56, 85, 61, 78, 68, 76, 70, 73, 80, 94, 90, 93, 73, 92, 99, 112, 105, 104, 97, 104, 99, 139, 134, 125, 116, 126, 137, 149, 146, 137, 119, 140, 154, 182, 117, 167, 146, 149, 172, 157, 131, 161, 151, 166, 191, 224, 218, 190, 150, 294, 155, 205
Offset: 1

Views

Author

J. M. Bergot and Robert Israel, Apr 14 2021

Keywords

Comments

If n = prime(k) is prime, a(n) = 3*prime(k)+1+A210934(k).

Examples

			For n = 3, A001414(3) = 3, A000203(3) = A000203(A001414(3)) = 4 and A001414(A000203(3)) = 4, so a(3) = 3+4+4+4 = 15.
		

Crossrefs

Programs

  • Maple
    spf:= proc(n) local t; add(t[1]*t[2],t=ifactors(n)[2]) end proc:
    f:= proc(n) local a,b;
      a:= spf(n);
      b:= numtheory:-sigma(n);
      a+b+spf(b)+numtheory:-sigma(a)
    end proc:
    map(f, [$1..100]);
Showing 1-4 of 4 results.