A131992 -id:A131992 - cp's OEIS frontend

A006093 a(n) = prime(n) - 1.

1, 2, 4, 6, 10, 12, 16, 18, 22, 28, 30, 36, 40, 42, 46, 52, 58, 60, 66, 70, 72, 78, 82, 88, 96, 100, 102, 106, 108, 112, 126, 130, 136, 138, 148, 150, 156, 162, 166, 172, 178, 180, 190, 192, 196, 198, 210, 222, 226, 228, 232, 238, 240, 250, 256, 262, 268, 270

Offset: 1

N. J. A. Sloane

These are also the numbers that cannot be written as i*j + i + j (i,j >= 1). - Rainer Rosenthal, Jun 24 2001; Henry Bottomley, Jul 06 2002
The values of k for which Sum_{j=0..n} (-1)^j*binomial(k, j)*binomial(k-1-j, n-j)/(j+1) produces an integer for all n such that n < k. Setting k=10 yields [0, 1, 4, 11, 19, 23, 19, 11, 4, 1, 0] for n = [-1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9], so 10 is in the sequence. Setting k=3 yields [0, 1, 1/2, 1/2] for n = [-1, 0, 1, 2], so 3 is not in the sequence. - Dug Eichelberger (dug(AT)mit.edu), May 14 2001
n such that x^n + x^(n-1) + x^(n-2) + ... + x + 1 is irreducible. - Robert G. Wilson v, Jun 22 2002
Records for Euler totient function phi.
Together with 0, n such that (n+1) divides (n!+1). - Benoit Cloitre, Aug 20 2002; corrected by Charles R Greathouse IV, Apr 20 2010
n such that phi(n^2) = phi(n^2 + n). - Jon Perry, Feb 19 2004
Numbers having only the trivial perfect partition consisting of a(n) 1's. - Lekraj Beedassy, Jul 23 2006
Numbers n such that the sequence {binomial coefficient C(k,n), k >= n } contains exactly one prime. - Artur Jasinski, Dec 02 2007
Record values of A143201: a(n) = A143201(A001747(n+1)) for n > 1. - Reinhard Zumkeller, Aug 12 2008
From Reinhard Zumkeller, Jul 10 2009: (Start)
The first N terms can be generated by the following sieving process:
  start with {1, 2, 3, 4, ..., N - 1, N};
  for i := 1 until SQRT(N) do
  (if (i is not striked out) then
  (for j := 2 * i + 1 step i + 1 until N do
  (strike j from the list)));
  remaining numbers = {a(n): a(n) <= N}. (End)
a(n) = partial sums of A075526(n-1) = Sum_{1..n} A075526(n-1) = Sum_{1..n} (A008578(n+1) - A008578(n)) = Sum_{1..n} (A158611(n+2) - A158611(n+1)) for n >= 1. - Jaroslav Krizek, Aug 04 2009
A006093 U A072668 = A000027. - Juri-Stepan Gerasimov, Oct 22 2009
A171400(a(n)) = 1 for n <> 2: subsequence of A171401, except for a(2) = 2. - Reinhard Zumkeller, Dec 08 2009
Numerator of (1 - 1/prime(n)). - Juri-Stepan Gerasimov, Jun 05 2010
Numbers n such that A002322(n+1) = n. This statement is stronger than repeating the property of the entries in A002322, because it also says in reciprocity that this sequence here contains no numbers beyond the Carmichael numbers with that property. - Michel Lagneau, Dec 12 2010
a(n) = A192134(A095874(A000040(n))); subsequence of A192133. - Reinhard Zumkeller, Jun 26 2011
prime(a(n)) + prime(k) < prime(a(k) + k) for at least one k <= a(n): A212210(a(n),k) < 0. - Reinhard Zumkeller, May 05 2012
Except for the first term, numbers n such that the sum of first n natural numbers does not divide the product of first n natural numbers; that is, n*(n + 1)/2 does not divide n!. - Jayanta Basu, Apr 24 2013
BigOmega(a(n)) equals BigOmega(a(n)*(a(n) + 1)/2), where BigOmega = A001222. Rationale: BigOmega of the product on the right hand side factorizes as BigOmega(a/2) + Bigomega(a+1) = BigOmega(a/2) + 1 because a/2 and a + 1 are coprime, because BigOmega is additive, and because a + 1 is prime. Furthermore Bigomega(a/2) = Bigomega(a) - 1 because essentially all 'a' are even. - Irina Gerasimova, Jun 06 2013
Record values of A060681. - Omar E. Pol, Oct 26 2013
Deficiency of n-th prime. - Omar E. Pol, Jan 30 2014
Conjecture: All the sums Sum_{k=s..t} 1/a(k) with 1 <= s <= t are pairwise distinct. In general, for any integers d >= -1 and m > 0, if Sum_{k=i..j} 1/(prime(k)+d)^m = Sum_{k=s..t} 1/(prime(k)+d)^m with 0 < i <= j and 0 < s <= t then we must have (i,j) = (s,t), unless d = m = 1 and {(i,j),(s,t)} = {(4,4),(8,10)} or {(4,7),(5,10)}. (Note that 1/(prime(8)+1)+1/(prime(9)+1)+1/(prime(10)+1) = 1/(prime(4)+1) and Sum_{k=5..10} 1/(prime(k)+1) = 1/(prime(4)+1) + Sum_{k=5..7} 1/(prime(k)+1).) - Zhi-Wei Sun, Sep 09 2015
Numbers n such that (prime(i)^n + n) is divisible by (n+1), for all i >= 1, except when prime(i) = n+1. - Richard R. Forberg, Aug 11 2016
a(n) is the period of Fubini numbers (A000670) over the n-th prime. - Federico Provvedi, Nov 28 2020

Archimedeans Problems Drive, Eureka, 40 (1979), 28.
Harvey Dubner, Generalized Fermat primes, J. Recreational Math., 18 (1985): 279-280.
M. Gardner, The Colossal Book of Mathematics, pp. 31, W. W. Norton & Co., NY, 2001.
M. Gardner, Mathematical Circus, pp. 251-2, Alfred A. Knopf, NY, 1979.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe, Table of n, a(n) for n = 1..10000
Thomas F. Bloom, Unit fractions with shifted prime denominators, arXiv:2305.02689 [math.NT], 2023.
R. P. Boas & N. J. A. Sloane, Correspondence, 1974
Harvey Dubner, Generalized Fermat primes, J. Recreational Math. 18.4 (1985-1986), 279. (Annotated scanned copy)
Armel Mercier, Problem E 3065, American Mathematical Monthly, 1984, p. 649.
Armel Mercier, S. K. Rangarajan, J. C. Binz and Dan Marcus, Problem E 3065, American Mathematical Monthly, No. 4, 1987, pp. 378.
Poo-Sung Park, Additive uniqueness of PRIMES-1 for multiplicative functions, arXiv:1708.03037 [math.NT], 2017.
J. R. Rickard and J. J. Hitchcock, Problem Drive 4, Archimedeans Problems Drive, Eureka, 40 (1979), 28-29, 40. (Annotated scanned copy)
Index entries for sequences generated by sieves

a(n) = K(n, 1) and A034693(K(n, 1)) = 1 for all n. The subscript n refers to this sequence and K(n, 1) is the index in A034693. - Labos Elemer
Cf. A000040, A034694. Different from A075728.
Complement of A072668 (composite numbers minus 1), A072670(a(n))=0.
Essentially the same as A039915.
Cf. A084920, A006093, A050997, A008864, A060800, A131991, A131992, A131993.
Cf. A101301 (partial sums), A005867 (partial products).
Column 1 of the following arrays/triangles: A087738, A249741, A352707, A378979, A379010.
The last diagonal of A162619, and of A174996, the first diagonal in A131424.
Row lengths of irregular triangles A086145, A124223, A212157.

Filtered([1..280],IsPrime)-1; # Muniru A Asiru, Nov 25 2018

a006093 = (subtract 1) . a000040  -- Reinhard Zumkeller, Mar 06 2012

[NthPrime(n)-1: n in [1..100]]; // Vincenzo Librandi, Nov 17 2015

for n from 2 to 271 do if (n! mod n^2 = n*(n-1) and (n<>4) then print(n-1) fi od; # Gary Detlefs, Sep 10 2010
# alternative
A006093 := proc(n)
    ithprime(n)-1 ;
end proc:
seq(A006093(n),n=1..100) ; # R. J. Mathar, Feb 06 2019

Table[Prime[n] - 1, {n, 1, 30}] (* Vladimir Joseph Stephan Orlovsky, Apr 27 2008 *)
a[ n_] := If[ n < 1, 0, -1 + Prime @ n] (* Michael Somos, Jul 17 2011 *)
Prime[Range[60]] - 1 (* Alonso del Arte, Oct 26 2013 *)

isA006093(n) = isprime(n+1) \\ Michael B. Porter, Apr 09 2010

A006093(n) = prime(n)-1 \\ Michael B. Porter, Apr 09 2010

\\ Sieve as described in Rainer Rosenthal's comment.
m=270;s=vector(m);for(i=1,m,for(j=i,m,k=i*j+i+j;if(k<=m,s[k]=1)));for(k=1,m,if(s[k]==0,print1(k,", "))); \\ Hugo Pfoertner, May 14 2019

from sympy import prime
for n in range(1,100): print(prime(n)-1, end=', ') # Stefano Spezia, Nov 30 2018

a(n) = (p-1)! mod p where p is the n-th prime, by Wilson's theorem. - Jonathan Sondow, Jul 13 2010
a(n) = A000010(prime(n)) = A000010(A006005(n)). - Antti Karttunen, Dec 16 2012
a(n) = A005867(n+1)/A005867(n). - Eric Desbiaux, May 07 2013
a(n) = A000040(n) - 1. - Omar E. Pol, Oct 26 2013
a(n) = A033879(A000040(n)). - Omar E. Pol, Jan 30 2014

Correction for change of offset in A158611 and A008578 in Aug 2009 Jaroslav Krizek, Jan 27 2010
Obfuscating comments removed by Joerg Arndt, Mar 11 2010
Edited by Charles R Greathouse IV, Apr 20 2010

A008864 a(n) = prime(n) + 1.

Original entry on oeis.org

3, 4, 6, 8, 12, 14, 18, 20, 24, 30, 32, 38, 42, 44, 48, 54, 60, 62, 68, 72, 74, 80, 84, 90, 98, 102, 104, 108, 110, 114, 128, 132, 138, 140, 150, 152, 158, 164, 168, 174, 180, 182, 192, 194, 198, 200, 212, 224, 228, 230, 234, 240, 242, 252, 258, 264, 270, 272, 278, 282, 284

Offset: 1

N. J. A. Sloane, R. K. Guy

Sum of divisors of prime(n). - Labos Elemer, May 24 2001
For n > 1, there are a(n) more nonnegative Hurwitz quaternions than nonnegative Lipschitz quaternions, which are counted in A239396 and A239394, respectively. - T. D. Noe, Mar 31 2014
These are the numbers which are in A239708 or in A187813, but excluding the first 3 terms of A187813, i.e., a number m is a term if and only if m is a term > 2 of A187813, or m is the sum of two distinct powers of 2 such that m - 1 is prime. This means that a number m is a term if and only if m is a term > 2 such that there is no base b with a base-b digital sum of b, or b = 2 is the only base for which the base-b digital sum of m is b. a(6) is the only term such that a(n) = A187813(n); for n < 6, we have a(n) > A187813(n), and for n > 6, we have a(n) < A187813(n). - Hieronymus Fischer, Apr 10 2014
Does not contain any number of the format 1 + q + ... + q^e, q prime, e >= 2, i.e., no terms of A060800,  A131991, A131992, A131993 etc. [Proof: that requires 1 + p = 1 + q + ... + q^e, or p = q*(1 + ... + q^(e-1)). This is not solvable because the left hand side is prime, the right hand side composite.] - R. J. Mathar, Mar 15 2018
1/a(n) is the asymptotic density of numbers whose prime(n)-adic valuation is odd. - Amiram Eldar, Jan 23 2021

C. W. Trigg, Problem #1210, Series Formation, J. Rec. Math., 15 (1982), 221-222.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
R. P. Boas and N. J. A. Sloane, Correspondence, 1974
N. J. A. Sloane and Brady Haran, Eureka Sequences, Numberphile video (2021).

Cf. A000040, A060800, A131991, A131992, A131993, A141468.
Cf. A007953, A079696, A187813, A239703, A239708.
Column 1 of A341605, column 2 of A286623 and of A328464.
Partial sums of A125266.

a008864 = (+ 1) . a000040
-- Reinhard Zumkeller, Sep 04 2012, Oct 08 2012

[NthPrime(n)+1: n in [1..70]]; // Vincenzo Librandi, Jul 30 2016

A008864:=n->ithprime(n)+1; seq(A008864(n), n=1..50); # Wesley Ivan Hurt, Apr 11 2014

Prime[Range[70]]+1 (* Vladimir Joseph Stephan Orlovsky, Apr 27 2008 *)

forprime(p=2,1e3,print1(p+1", ")) \\ Charles R Greathouse IV, Jun 16 2011

A008864(n) = (1+prime(n)); \\ Antti Karttunen, Mar 14 2021

[nth_prime(n) +1 for n in (1..70)] # G. C. Greubel, May 20 2019

a(n) = prime(n) + 1 = A000040(n) + 1.
a(n) = A000005(A034785(n)) = A000203(A000040(n)). - Labos Elemer, May 24 2001
a(n) = A084920(n) / A006093(n). - Reinhard Zumkeller, Aug 06 2007
A239703(a(n)) <= 1. - Hieronymus Fischer, Apr 10 2014
From Ilya Gutkovskiy, Jul 30 2016: (Start)
a(n) ~ n*log(n).
Product_{n>=1} (1 + 2/(a(n)*(a(n) - 2))) = 5/2. (End)

A030514 a(n) = prime(n)^4.

Original entry on oeis.org

16, 81, 625, 2401, 14641, 28561, 83521, 130321, 279841, 707281, 923521, 1874161, 2825761, 3418801, 4879681, 7890481, 12117361, 13845841, 20151121, 25411681, 28398241, 38950081, 47458321, 62742241, 88529281, 104060401, 112550881, 131079601, 141158161

Offset: 1

Jeff Burch

Numbers with 5 divisors (1, p, p^2, p^3, p^4, where p is the n-th prime). - Alexandre Wajnberg, Jan 15 2006
Subsequence of A036967. - Reinhard Zumkeller, Feb 05 2008
The n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime. - Omar E. Pol, May 06 2008
The general product formula for even s is: product_{p = A000040} (p^s-1)/(p^s+1) = 2*Bernoulli(2s)/( binomial(2s, s)*Bernoulli^2(s)), where the infinite product is over all primes. Here, with s = 4, product_{n = 1, 2, ...} (a(n)-1)/(a(n)+1) = 6/7. In A030516, where s = 6, the product of the ratios is 691/715. For s = 8, the 8th row in A120458, the corresponding product of ratios is 7234/7293. - R. J. Mathar, Feb 01 2009
Except for the first three terms, all others are congruent to 1 mod 240. - Robert Israel, Aug 29 2014

R. J. Mathar, Table of n, a(n) for n = 1..1000
Eric Weisstein's World of Mathematics, Prime Power.
OEIS Wiki, Index entries for number of divisors
Index to sequences related to prime signature

Cf. A030078, A085964, A131991, A131992, A000005, A000040, A001248, A157290, A215267.

Cf. A258601.

a030514 = (^ 4) . a000040
a030514_list = map (^ 4) a000040_list
-- Reinhard Zumkeller, Jun 03 2015

[NthPrime(n)^4: n in [1..100] ]; // Vincenzo Librandi, Apr 22 2011

map(p -> p^4, select(isprime,[2,seq(2*i+1,i=1..100)])); # Robert Israel, Aug 29 2014

Array[Prime[#]^4 &, 5!] (* Vladimir Joseph Stephan Orlovsky, Sep 01 2008 *)

a(n)=prime(n)^4 \\ Charles R Greathouse IV, Mar 21 2013

[p**4 for p in prime_range(100)]
# Zerinvary Lajos, May 15 2007

a(n) = A000040(n)^(5-1) = A000040(n)^4, where 5 is the number of divisors of a(n). - Omar E. Pol, May 06 2008
A000005(a(n)) = 5. - Alexandre Wajnberg, Jan 15 2006
A056595(a(n)) = 2. - Reinhard Zumkeller, Aug 15 2011
Sum_{n>=1} 1/a(n) = P(4) = 0.0769931397... (A085964). - Amiram Eldar, Jul 27 2020
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(4)/zeta(8) = 105/Pi^4 (A157290).
Product_{n>=1} (1 - 1/a(n)) = 1/zeta(4) = 90/Pi^4 (A215267). (End)

Description corrected by Eric W. Weisstein

A050997 Fifth powers of primes.

Original entry on oeis.org

32, 243, 3125, 16807, 161051, 371293, 1419857, 2476099, 6436343, 20511149, 28629151, 69343957, 115856201, 147008443, 229345007, 418195493, 714924299, 844596301, 1350125107, 1804229351, 2073071593, 3077056399, 3939040643, 5584059449, 8587340257, 10510100501

Offset: 1

Eric W. Weisstein

Numbers k such that A062799(k) = 5.

Let r(n) = (a(n)+1)/(a(n)-1) if a(n) mod 4 = 3, (a(n)-1)/(a(n)+1) otherwise; then Product_{n>=1} r(n) = (31/33) * (244/242) * (3124/3126) * (16808/16806) * ... = 246016/259875. - Dimitris Valianatos, Mar 09 2020

T. D. Noe, Table of n, a(n) for n = 1..1000
Xavier Gourdon and Pascal Sebah, Some Constants from Number theory.
Eric Weisstein's World of Mathematics, MathWorld: Prime Power.
Index to sequences related to prime signature

Cf. A000040, A001248, A030078, A030514, A085965, A131992, A131993, A013663, A157291.

Cf. A258602.

a050997 = (^ 5) . a000040
a050997_list = map (^ 5) a000040_list
-- Reinhard Zumkeller, Jun 03 2015

[p^5: p in PrimesUpTo(300)]; // Vincenzo Librandi, Mar 27 2014

Array[Prime[ # ]^5 &, 30] (* Vladimir Joseph Stephan Orlovsky, May 01 2008 *)

PARI
```
vector(66,n,prime(n)^5)
```

A056595(a(n)) = 3. - Reinhard Zumkeller, Aug 15 2011
Sum_{n>=1} 1/a(n) = P(5) = 0.0357550174... (A085965). - Amiram Eldar, Jul 27 2020
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(5)/zeta(10) (A157291).
Product_{n>=1} (1 - 1/a(n)) = 1/zeta(5) = 1/A013663. (End)

A060800 a(n) = p^2 + p + 1 where p runs through the primes.

Original entry on oeis.org

7, 13, 31, 57, 133, 183, 307, 381, 553, 871, 993, 1407, 1723, 1893, 2257, 2863, 3541, 3783, 4557, 5113, 5403, 6321, 6973, 8011, 9507, 10303, 10713, 11557, 11991, 12883, 16257, 17293, 18907, 19461, 22351, 22953, 24807, 26733, 28057, 30103, 32221

Offset: 1

Jason Earls, Apr 27 2001

Terms are divisible by 3 iff p is of the form 6*m+1 (A002476). - Michel Marcus, Jan 15 2017

			a(3) = 31 because 5^2 + 5 + 1 = 31.

Harry J. Smith, Table of n, a(n) for n = 1..1000
R. J. Mathar, No common terms in the sequences sigma(p^i) and sigma(p^(i+1)) as p runs through the primes.

Cf. A001248, A131991, A131992, A131993, A253905.

Cf. A008864, A000203. - Zak Seidov, Feb 13 2016

[p^2+p+1: p in PrimesUpTo(200)]; // Vincenzo Librandi, Mar 20 2014

A060800:= n -> map (p -> p^(2)+p+1, ithprime(n)):
seq (A060800(n), n=1..41); # Jani Melik, Jan 25 2011

#^2 + # + 1&/@Prime[Range[200]] (* Vincenzo Librandi, Mar 20 2014 *)

{ n=0; forprime (p=2, prime(1000), write("b060800.txt", n++, " ", p^2 + p + 1); ) } \\ Harry J. Smith, Jul 13 2009

a(n) = A036690(n) + 1.
a(n) = 1 + A008864(n)*A000040(n) = (A030078(n) - 1)/A006093(n). - Reinhard Zumkeller, Aug 06 2007
a(n) = sigma(prime(n)^2) = A000203(A000040(n)^2). - Zak Seidov, Feb 13 2016
a(n) = A000203(A001248(n)). - Michel Marcus, Feb 15 2016
Product_{n>=1} (1 - 1/a(n)) = zeta(3)/zeta(2) (A253905). - Amiram Eldar, Nov 07 2022

More terms from Larry Reeves (larryr(AT)acm.org), May 03 2001

A131991 a(n) = 1 + prime(n) + prime(n)^2 + prime(n)^3.

Original entry on oeis.org

15, 40, 156, 400, 1464, 2380, 5220, 7240, 12720, 25260, 30784, 52060, 70644, 81400, 106080, 151740, 208920, 230764, 305320, 363024, 394420, 499360, 578760, 712980, 922180, 1040604, 1103440, 1236600, 1307020, 1455780, 2064640, 2265384

Offset: 1

Reinhard Zumkeller, Aug 06 2007

Number of points and lines for the prime(n)-Cremona-Richmond configuration. - Carlos Segovia Gonzalez, Jul 30 2020

			a(4)=400 because the 4th prime is 7, 7^3=343, 7^2=49, and 343+49+7+1=400.

Vincenzo Librandi, Table of n, a(n) for n = 1..1000
C. Segovia and M. Winklmeier, Calculating the dimension of the universal embedding of the symplectic dual polar space using languages, arXiv:1312.4315 [math.CO], 2013-2019.
C. Segovia and M. Winklmeier, Calculating the dimension of the universal embedding of the symplectic dual polar space using languages, The Elec. Jour. of Comb. 27(4) (2020), P4.39.

Cf. A030078, A008864, A131992, A131993.

Cf. A000040, A000203. - Zak Seidov, Feb 13 2016

[1+NthPrime(n)+NthPrime(n)^2+NthPrime(n)^3: n in [1..40]]; // Vincenzo Librandi, Dec 27 2010

A131991:= n -> map (p -> p^(3)+p^(2)+p+1, ithprime(n)):
seq (A131991(n), n=1..32); # Jani Melik, Jan 25 2011

#^3 + #^2 + # + 1 &/@Prime[Range[100]] (* Vincenzo Librandi, Mar 20 2014 *)

a(n) = 1 + A060800(n)*A000040(n).
a(n) = (A030514(n) - 1)/A006093(n).
a(n) = A000203(A000040(n)^3). - Zak Seidov, Feb 13 2016

A131993 1 + prime(n) + prime(n)^2 + prime(n)^3 + prime(n)^4 + prime(n)^5.

Original entry on oeis.org

63, 364, 3906, 19608, 177156, 402234, 1508598, 2613660, 6728904, 21243690, 29583456, 71270178, 118752606, 150508644, 234330768, 426237714, 727250580, 858672906, 1370581548, 1830004056, 2101864254, 3116505840, 3987077724

Offset: 1

Reinhard Zumkeller, Aug 06 2007

nonn

a(n) = 1 + A131992(n)*A000040(n).

Ivan Panchenko, Table of n, a(n) for n = 1..1000

Cf. A050997, A008864, A060800, A131991.

[1+(&+[NthPrime(n)^(k): k in [1..5]]): n in [1..100]]; // Berselli - Librandi, Apr 20 2011

Total[#^Range[0,5]]&/@Prime[Range[30]]  (* Harvey P. Dale, Apr 20 2011 *)

a(n) = (A030516(n) - 1)/A006093(n).

A190527 Primes of the form p^4 + p^3 + p^2 + p + 1, where p is prime.

Original entry on oeis.org

31, 2801, 30941, 88741, 292561, 732541, 3500201, 28792661, 39449441, 48037081, 262209281, 1394714501, 2666986681, 3276517921, 4802611441, 5908670381, 12936304421, 16656709681, 19408913261, 24903325661, 37226181521, 43713558101, 52753304641, 64141071121, 96427561501, 100648118041

Offset: 1

Bernard Schott, Dec 20 2012

nonn

These primes are generated by exactly A065509, cf. 2nd formula.
These numbers are repunit primes 11111_p, so they are Brazilian primes (A085104).
When p^4 + p^3 + p^2 + p + 1 = sigma(p^4) is prime, then it equals A193574(p^4), so that this sequence is a subsequence of A193574; by definition it is also a subsequence of A053699 and A131992. - Hartmut F. W. Hoft, May 05 2017

			a(3) = 30941 = 11111_13 = 13^4 + 13^3 + 13^2 + 13^1 + 1 is prime.

Vincenzo Librandi, Table of n, a(n) for n = 1..1100

Cf. A049409 (n^4 + ... + 1 is prime), A065509 (primes among these n), A193574.
Subsequence of A088548 (primes n^4 + ... + 1) and A085104 ("Brazilian" primes, of the form 1 + n + n^2 + ... + n^k).
Intersection of A000040 (primes) and A131992 (p^4 + ... + 1), subsequence of A053699 (n^4 + ... + 1).

[p: p in PrimesUpTo(600) | IsPrime(p) where p is p^4 +p^3+p^2+p+1]; // Vincenzo Librandi, May 06 2017

a190527[n_] := Select[Map[(Prime[#]^5-1)/(Prime[#]-1)&, Range[n]], PrimeQ]
a190527[100] (* data *) (* Hartmut F. W. Hoft, May 05 2017 *)
Select[#^4 + #^3 + #^2 + # + 1 &/@Prime[Range[100]], PrimeQ] (* Vincenzo Librandi, May 06 2017 *)

[q|p<-primes(100),ispseudoprime(q=(p^5-1)\(p-1))]
A190527_vec(N)=[(p^5-1)\(p-1)|p<-A065509_vec(N)] \\ M. F. Hasler, Mar 03 2020

a(n) = A193574(A065509(n)^4). - Hartmut F. W. Hoft, May 08 2017

a(n) = A053699(A065509(n)) = A000203(A065509(n)^4). - M. F. Hasler, Mar 03 2020

a(7) corrected and a(18)-a(26) added by Hartmut F. W. Hoft, May 05 2017

Edited by M. F. Hasler, Mar 06 2020

A319076 Square array T(n,k) read by antidiagonal upwards in which column k lists the partial sums of the powers of the k-th prime, n >= 0, k >= 1.

Original entry on oeis.org

1, 3, 1, 7, 4, 1, 15, 13, 6, 1, 31, 40, 31, 8, 1, 63, 121, 156, 57, 12, 1, 127, 364, 781, 400, 133, 14, 1, 255, 1093, 3906, 2801, 1464, 183, 18, 1, 511, 3280, 19531, 19608, 16105, 2380, 307, 20, 1, 1023, 9841, 97656, 137257, 177156, 30941, 5220, 381, 24, 1, 2047, 29524, 488281, 960800, 1948717

Offset: 0

Omar E. Pol, Sep 09 2018

T(n,k) is also the sum of the divisors of the n-th nonnegative power of the k-th prime, n >= 0, k >= 1.

			The corner of the square array is as follows:
         A126646 A003462 A003463  A023000    A016123    A091030     A091045
A000012        1,      1,      1,       1,         1,         1,          1, ...
A008864        3,      4,      6,       8,        12,        14,         18, ...
A060800        7,     13,     31,      57,       133,       183,        307, ...
A131991       15,     40,    156,     400,      1464,      2380,       5220, ...
A131992       31,    121,    781,    2801,     16105,     30941,      88741, ...
A131993       63,    364,   3906,   19608,    177156,    402234,    1508598, ...
.......      127,   1093,  19531,  137257,   1948717,   5229043,   25646167, ...
.......      255,   3280,  97656,  960800,  21435888,  67977560,  435984840, ...
.......      511,   9841, 488281, 6725601, 235794769, 883708281, 7411742281, ...
...

Rows 0-5: A000012, A008864, A060800, A131991, A131992, A131993.
Columns 1-15: A126646, A003462, A003463, A023000, A016123, A091030, A091045, A218722, A218726, A218732, A218734, A218740, A218744, A218746, A218750.
Main diagonal gives A319074.
Cf. A000079, A000244, A000351, A000420, A001020, A001022, A001026, A001029, A009967, A009973, A009975, A009981, A009985, A009987, A009991.
Cf. A000040, A000203, A006093, A319075.

T(n, k) = sigma(prime(k)^n); \\ Michel Marcus, Sep 13 2018

T(n,k) = A000203(A000040(k)^n).
T(n,k) = Sum_{j=0..n} A000040(k)^j.
T(n,k) = Sum_{j=0..n} A319075(j,k).
T(n,k) = (A000040(k)^(n+1) - 1)/(A000040(k) - 1).
T(n,k) = (A000040(k)^(n+1) - 1)/A006093(k).

cp's OEIS Frontend

A006093 a(n) = prime(n) - 1.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

GAP

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

Formula

Extensions

A008864 a(n) = prime(n) + 1.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Sage

Formula

A030514 a(n) = prime(n)^4.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Sage

Formula

Extensions

A050997 Fifth powers of primes.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

Formula

A060800 a(n) = p^2 + p + 1 where p runs through the primes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI