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.

User: Richard R. Forberg

Richard R. Forberg's wiki page.

Richard R. Forberg has authored 69 sequences. Here are the ten most recent ones:

A376256 Numbers which are the minimum of a cycle in the map x -> phi(sigma(x)).

Original entry on oeis.org

1, 2, 4, 8, 12, 16, 48, 72, 128, 240, 432, 576, 720, 1200, 1728, 1800, 6912, 10368, 15552, 27648, 32768, 41472, 67392, 142560, 184320, 272160, 326592, 712800, 1140480, 1190400, 1658880, 3345408, 3571200, 5702400, 6220800, 10222080, 14859936, 29719872, 40255488, 50319360, 113218560, 118879488
Offset: 1

Author

Richard R. Forberg, Sep 16 2024

Keywords

Comments

No further terms < 254731536.
The status of 254731536 is unknown, but conjectured not a term.
Additional terms include 2142720000, 5033164800, 150493593600, 3852635996160.
See further cycles in the linked document here below which contains 422 cycles. It includes the 80 cycles complied by Jud McCrainie in a linked document at A095955.

Crossrefs

Union of A001229, A373435, A373453, A373454, etc.

Programs

  • PARI
    \\ Naive program, assumes eventual termination (ok upto 254731535).
    isok(n)={my(M=Map(),p=n); while(!mapisdefined(M,p) && p>=n, mapput(M,p,1); p=eulerphi(sigma(p))); p==n} \\ Andrew Howroyd, Sep 19 2024

A375920 Decimal expansion of Sum_{n>=2} (zeta(n)^n - 1).

Original entry on oeis.org

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

Author

Richard R. Forberg, Sep 02 2024

Keywords

Comments

It is interesting to note that this sum is very close in value to 1/3 of Product_{n>=2} zeta(n)^n, A375887, where that factor's first 30 digits are: 0.333319653211135001436063576617.

Examples

			3.255474995780369262094368665069015138075282643803397585341859272265720258...
		

Crossrefs

Cf. A375887 (Product_{n>=2} zeta(n)^n), A021002 (Product_{n>2} zeta(n)), A093720 (Sum_{n>=2} zeta(n)/n!), A013661 (zeta(2)).

Programs

  • Maple
    evalf(Sum(Zeta(n)^n - 1, n = 2 .. infinity), 120); # Vaclav Kotesovec, Sep 02 2024
  • Mathematica
    RealDigits[N[Sum[Zeta[n]^n - 1, {n, 2, 1000}], 150]][[1]]

A375887 Decimal expansion of Product_{n>=2} zeta(n)^n.

Original entry on oeis.org

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

Author

Richard R. Forberg, Sep 01 2024

Keywords

Comments

It is interesting to note that this product is very close in value to 3 * Sum_{n>=2} (zeta(n)^n-1), A375920, where that factor's first 30 digits are: 3.00012312615292744064909403341.

Examples

			9.766825821453289699230682695640792162028987950967280928488833051400227...
		

Crossrefs

Cf. A375920,(Sum_{n>=2} (zeta(n)^n-1)), A021002 (Product_{n>=2} zeta(n)), A093720 (Sum_{n >= 2} zeta(n)/n!), A013661 (zeta(2)).

Programs

  • Maple
    evalf(Product(Zeta(n)^n, n = 2 .. infinity), 150); # Vaclav Kotesovec, Sep 02 2024
  • Mathematica
    RealDigits[N[Product[Zeta[n]^n, {n, 2, 500}], 150]][[1]]
  • PARI
    prodinf(k = 2, zeta(k)^k) \\ Amiram Eldar, Sep 02 2024

A341183 Numbers k such that (30*k)^4 + 1 is prime.

Original entry on oeis.org

3, 6, 7, 18, 20, 23, 27, 31, 39, 40, 55, 56, 57, 68, 75, 80, 88, 90, 101, 103, 113, 115, 116, 124, 126, 129, 132, 136, 139, 148, 153, 163, 165, 173, 194, 203, 211, 215, 218, 232, 260, 277, 284, 285, 288, 290, 293, 302, 309, 336, 344, 347, 357, 358, 362, 368, 377, 390, 394, 397, 401
Offset: 1

Author

Richard R. Forberg, Feb 06 2021

Keywords

Comments

This is the third sequence in a family defined by this same condition, when generalized as (30*k)^(2^m).
A111175 and A138220 are the first and second sequences in this family, with m = 0 and m = 1, respectively.

Crossrefs

Programs

  • Mathematica
    (* This program gives some results for each member of the family for m up to 10. *)
    Do[resultk = {};
    Do[
      r = (30 k)^(2^m) + 1;
      If[PrimeQ[r], AppendTo[resultk, k]];
      , {k, 1, 1000}];
    If[Length[resultk] > 0, Print[{m, 2^m, resultk}]], {m, 0, 10}]
  • PARI
    isok(k) = isprime((30*k)^4 + 1); \\ Michel Marcus, Feb 13 2021

A274610 Values of c such that p^2-c and p^2+c are both positive primes, for the special case when p^2-c = 3, or c = p^2-3, where p is a prime.

Original entry on oeis.org

1, 22, 118, 358, 1678, 3478, 5038, 6238, 10198, 17158, 19318, 22798, 32758, 36478, 52438, 58078, 72358, 109558, 143638, 177238, 212518, 229438, 248998, 259078, 326038, 410878, 434278, 502678, 737878, 863038, 982078, 1062958, 1100398, 1274638, 1442398, 1515358, 1559998, 1635838, 1907158, 2042038
Offset: 1

Author

Richard R. Forberg, Jun 30 2016

Keywords

Comments

All values of a(n), except {1, 22}, equal 118 mod 120.
Corresponding values of p = sqrt(c+3)begin {2,5,11,19,41,59,71,79,101...}.
This relates to a comment at A047222 regarding c values for the general case of p^2-c and p^2+c both being positive primes.

Programs

  • Mathematica
    result = {}; Do[c = Prime[i]^2 - 3; If[PrimeQ[Prime[i]^2 + c], AppendTo[result, {Prime[i], c}]], {i, 1, 1000}]; result
  • PARI
    lista(nn) = forprime(p=2, nn, if (isprime(2*p^2-3), print1(p^2-3, ", "))); \\ Michel Marcus, Jul 05 2016

A274609 Primes p such that both 2p-1 and 2p^2-2p+1 are prime.

Original entry on oeis.org

2, 3, 31, 331, 1171, 2011, 2281, 3181, 4621, 4861, 6151, 6211, 6961, 7951, 8521, 9151, 11251, 12211, 13411, 15661, 17491, 18121, 19141, 20641, 22531, 23071, 23581, 24631, 25411, 26041, 26161, 26431, 26641, 27091, 27271, 27361, 27691, 28201, 28621, 29221, 31891, 33151, 34261, 35491, 36451
Offset: 1

Author

Richard R. Forberg, Jun 30 2016

Keywords

Comments

All values of a(n), except {2,3}, are equal to 1 mod 30.
These are also primes p such that both p^2+c and p^2-c are positive primes, for some c, when c is a square, since that requires c = (p-1)^2. Corresponding c values begin {1, 4, 900, 108900, ...}. This relates to a comment at A047222.

Examples

			31^2 - 30^2 = 61 and 31^2 + 30^2 = 1861 are both prime.
		

Crossrefs

Cf. A047222.

Programs

  • Mathematica
    result = {}; Do[If[PrimeQ[2*Prime[i] - 1] && PrimeQ[2*Prime[i]^2 - 2*Prime[i] + 1], AppendTo[result, Prime[i]]], {i, 1, 10000}]; result
    Select[Prime[Range[4000]],AllTrue[{2#-1,2#^2-2#+1},PrimeQ]&] (* Harvey P. Dale, Dec 26 2022 *)
  • PARI
    is(n)=isprime(2*n-1) && isprime(2*n^2-2*n+1) && isprime(n) \\ Charles R Greathouse IV, Jul 15 2016

A270424 Numbers m such that m^2 is the sum of the squares of two or more consecutive primes.

Original entry on oeis.org

586, 6088, 8174, 11585, 11707, 270106, 288818, 375661, 724909, 732910, 937423, 1141509, 1326970, 1619934, 1776809, 1930140, 2239367, 2489647, 3063687, 3649371, 3790381, 3941615, 4193988, 4821615, 4887146, 5572173, 6047246, 6192322, 8088524, 9158347
Offset: 1

Author

Richard R. Forberg, Mar 30 2016

Keywords

Comments

m^2 = Sum_{i=k..j} prime(i)^2 is a square, for some k,j, j > k.
The 30 numbers given above are the only m values for all possible summations where the resulting m^2 < 10^14 (m <10^7). This requires searching from k values up to ~482,000, but with decreasing j-k ranges for efficiency.
Values of k that yield results begin: 13, 37, 101, 183, 235, 588, 805, 891, 1066, ... but do not correspond fully to the order of the m values shown.
Number of sequential summands (i.e., j-k+1) vary widely, with the smallest being 28 and largest being 10360, for those m values listed above.
Also note j-k+1 mod 8 = {0,1,4}, as expected, since prime(i)^2 mod 24 = 1, for i > 2.

Examples

			586 is in the sequence because 586^2 = 343396 = Sum_{i=13..40} prime(i)^2.
		

Crossrefs

Programs

  • Mathematica
    lim = 20000^2; L={}; P=Prime[Range[2 + PrimePi@ Sqrt[lim/2]]]^2; i = 1; While[ P[[i]] + P[[i+1]] <= lim, s = P[[i]]; j = i+1; While[(s += P[[j++]]) <= lim,If[IntegerQ@ Sqrt@ s, AppendTo[L, Sqrt@ s]]]; i++]; Union@L (* Giovanni Resta, Apr 13 2016 *)
    result = {}; k = 3; While[k <= 481167, resultk = {}; sump = 0;
    count = 0; i = k; While[sump < 10^14, sump += Prime[i]^2;
      If[Mod[i - k + 1, 8] == 1 || Mod[i - k + 1, 8] == 0 ||
        Mod[i - k + 1, 8] == 4, If[i != k && IntegerQ[Sqrt[sump]], count++;
        AppendTo[resultk, {k, i - k + 1, sump}]]]; i++];
    If[count > 0, AppendTo[result, resultk]]; k++]; result (* Only for k>2, so as to use index values to reduce repeated checking Sqrt - Richard R. Forberg, Apr 14 2016 *)

A270361 Odd primes p for which there exists an odd prime q < p such that p*q - 1 is a square.

Original entry on oeis.org

13, 29, 53, 61, 73, 89, 97, 109, 137, 149, 157, 173, 233, 241, 277, 317, 349, 353, 373, 389, 397, 409, 433, 461, 521, 541, 569, 593, 617, 641, 653, 661, 673, 701, 709, 733, 757, 769, 773, 821, 829, 853, 877, 881, 929, 937
Offset: 1

Author

Richard R. Forberg, Mar 15 2016

Keywords

Comments

Conjecture: For any odd prime p there is at most one odd prime q, with q < p, for which p*q-1 is square.
(Note: If p were not restricted to being prime, there could be multiple primes q which make p*q-1 a square, with increasing multiplicities for larger p. The upper limit on the multiplicity of solutions, for prime and nonprime p and q values, is given by A006278. Note, however that those limits allow for solutions where q=2, and therefore two solutions when p and q are prime.)
a(n) is a subsequence of the Pythagorean Primes (A002144), which are of the form 4k+1.
The density of a(n) values among the primes declines with increasing n. For example, a(n) is about 22% of the first 1000 primes, and drops to about 15% of "incremental" primes around prime(10000). The density continues to fall among even larger primes. Twice those percentages apply as a portion of A002144.
All values of q also belong to A002144. It appears the set of q values "intends to" fully comprise A002144. This is notable because p values comprise an increasingly sparse subsequence within A002144, and each p value has just one q value.
The ability to fully comprise A002144 with q values is further challenged by the fact that for any given q value (i.e., any term of A002144) multiple values of p > q can be found such that p*q-1 is square. Thus q values are "promiscuous", and apparently without bounds on the number of p values they can serve.
Contrast this with primes p and q such that p*q+1 is square. The result are the Twin Primes (A001359 and A006512), arranged in a simple one-to-one correspondence, with p = q+2.

Examples

			13 is in this sequence because 13*5 - 1 = 64, which is square, with 5 < 13.
		

Crossrefs

Cf. A002144.

Programs

  • Mathematica
    result = {}; Do[p = Prime[i]; Do[q = Prime[j]; r = p*q - 1;
      If[Mod[r, 8] == 1 || Mod[r, 8] == 0 || Mod[r, 8] == 4,
       If[IntegerQ[Sqrt[r]], AppendTo[result, p]]], {j, 2, i - 1}],
    {i, 3, 1000}]; result
  • PARI
    lista(nn) = {forprime(p=3, nn, ok = 0; forprime (q=3, p-1, if (issquare(p*q-1), ok = 1; break);); if (ok, print1(p, ", ")););} \\ Michel Marcus, Apr 06 2016
  • Python
    from gmpy2 import is_prime, is_square
    for p in range(3, 10 ** 4, 2):
        flag = 0
        if is_prime(p):
            for q in range(3, p, 2):
                if is_square(p * q - 1) and is_prime(q):
                    flag = 1
                    break
        if flag:
            print(p, end=", ")
    # Soumil Mandal, Apr 07 2016
    

A268593 The largest of three distinct primes, a > b > c > 0, such that each of the following four numbers are also prime: a+b+c, a+b-c, a-b+c, a-b-c, regardless of sign.

Original entry on oeis.org

13, 17, 19, 23, 37, 43, 47, 53, 67, 73, 83, 97, 107, 113, 127, 137, 139, 157, 167, 181, 193, 197, 199, 223, 227, 229, 233, 251, 263, 269, 277, 283, 317, 337, 347, 349, 373, 379, 383, 397, 409, 431, 433, 439, 443, 461
Offset: 1

Author

Richard R. Forberg, Feb 07 2016

Keywords

Comments

A set of seven primes defined by the rules above may be called a "Prime Septet". For a given "a" value, from a(n) above, there may be multiple such septets.
These septets are found by searching all combinations of three distinct primes for a given candidate "a" value.
Once found, the following three properties of these septets are observed without exception for "a" values up to prime(800), resulting in these conjectures:
1. a-b-c = +-3.
2. b > a/2.
3. 5 <= c < a/2.
4. If a prime p >= 3 belongs to one Prime Septet, then it belongs to an infinite number of Prime Septets.
From the first conjecture it follows that no values for "a", "b" or "c" belong to A172256.
The first three conjectures can be used to accelerate the search for prime septets. Further observations below are from accelerated searches.
There are a total of 43 "missing" values in a(n), compared to the full odd prime complement of A172256, starting with 3, 5, 7, 11, 29, 31, 41, 71, ...., 10677.
There are a total of 25 missing "b" values compared to the full odd prime complement of A172256. The last missing value is 3067.
It appears the values of "c" comprise the full odd prime complement of A172256.
The three sets of values for a+b+c, a+b-c and a-b+c, by contrast, split their membership between A172256 and its odd prime complement.
Positive vs. negative values for a-b-c occur with approximately equal frequency, with positive values at 48.05% of the total for all 1106 septets with a <= prime(1000); reaching 48.30% of the total for all 65821 septets with a <= prime(10000); and reaching the 48.97% of the total for all 195359 septets with a <= prime(20000). No other members of these septets are negative.
In the first 1000 primes, 502 primes play the role of an "a" value, with an average of 4 Prime Septets subordinate to that "a" value.
In the 1000 primes from prime(9001) to prime(10000), 376 primes play the role of an "a" value, with an average of 26 Prime Septets subordinate to each "a" value, with none having only 1,2,3,4 or 5 septets.
It appears the highest "a" value having only one subordinate septet is prime(2127)=18587, where {a,b,c} = {18587, 18427, 157}. This was checked for candidate "a" values up to prime(20000). By prime(10^6) it is typical for each "a" value to support 500 to 1500 Prime Septets.
It might seem that all odd primes should belong to at least one Prime Septet. This is not the case.
The primes not included in any Prime Septet may be called "Lonely Primes". The first few such odd primes are: 151, 179, 239, 293, 313. They are a subset of A172256.
Within the first 1000 odd primes there are 259 primes that are "Lonely". Within the first 10000 odd primes there are 3978 such primes. In the 1000 primes from prime(9001) to prime(10000), 455 of them are "Lonely".
It appears likely that "Lonely Primes" are a majority among all primes.

Examples

			The first Prime Septet is based on {a,b,c} = {13,11,5}, with a+b+c =29,
a+b-c=19, a-b+c=7, a-b-c =-3. The first "a" value associated with more than one septet is 47, with {a,b,c} values of {47,37,13} and {47,43,7}.
		

Crossrefs

Programs

  • PARI
    listp(nn) = {forprime(p=5, nn, ok = 0; forprime(q=3, p-1, forprime(r=2, q-1, if (isprime(abs(p+q+r)) && isprime(abs(p+q-r)) && isprime(abs(p-q+r)) && isprime(abs(p-q-r)), ok = 1; break);); if (ok, break);); if (ok, print1(p, ", ")););} \\ Michel Marcus, Mar 13 2016

A248716 Numerator of (1/e)*Sum_{k>=0} (1/k!)*(Sum_{j=0..k} j^n).

Original entry on oeis.org

2, 3, 17, 27, 293, 791, 10583, 25685, 448303, 251411, 4503535, 6331107, 4436875097, 3335427631, 19619696071, 75379875277, 3019260651391, 16773385986619, 3047463007411973, 2732480436961811, 398377271835431771, 173581842021095897, 1716900426430701553, 35001773773285490879, 6684326532123939298051
Offset: 0

Author

Richard R. Forberg, Dec 28 2014

Keywords

Comments

Denominators are given by A130190, which in that entry are associated with the denominators of the z-sequence for a certain Sheffer matrix (triangle), but also apply here. See also the formula given there for the denominator of a certain sum involving the Stirling2 numbers.
Note that a(0) = 2 uses 0^0 := 1. For n >= 1 use triangle A079618 and the formula (1/e)*Sum_{k>=0} ((k+1)^n)/k! = Bell(n+1) = A000110(n+1). - Wolfdieter Lang, Feb 03 2015
If instead of Sum_{j=0..k} j^n one uses the sum with falling factorials, namely F(k, n) := Sum_{j=0..k} A008279(j, n) = A008279(k+1, n+1)/(n+1) the result for the rationals R(n) = (1/e)*Sum_{k>=0} (1/k!)*F(k, n) becomes very simple, namely R(n) = (n+2)/(n+1), n >= 0. - Wolfdieter Lang, Feb 03 2015

Examples

			Terms up to n = 10, with denominators, are 2/1, 3/2, 17/6, 27/4, 293/15, 791/12, 10583/42, 25685/24, 448303/90, 251411/10, 4503535/33, ... .
From _Wolfdieter Lang_, Feb 03 2015: (Start)
With triangle A079618, A064538 and the Bell numbers A000110 the rationals r(n) are:
n=4: (1/30)*(-1*1 + 0*2 + 10*5 + 15*15 + 6*52) = 293/15.
n=9: (1/20)*(0*1 + (-3)*2 +  0*5 + 10*15 + 0*52 + (-14)*203 + 0*877 + 15*4140 + 10*21147 + 2*115975) = 251411/10.
(End)
		

Crossrefs

Programs

  • Mathematica
    Numerator[Table[(1/Exp[1])*Sum[Sum[j^n/k!, {j, 0, k}], {k, 0, Infinity}],
    {n, 0, 100}]]

Formula

a(n) = numerator(r(n)) with the rationals r(n) = (1/e)*Sum_{k>=0} (1/k!)*(Sum_{j=0..k} j^n), n >= 0, where 0^0 := 1.
a(n) = numerator(r(n)), with r(n) = (1/A064538(n))*Sum_{k=0..n} T(n+1,k+1)*Bell(k+1), with T(n,k) = A079618(n,k) and Bell(n) = A000110(n), for n >= 1. a(0) = 2 using 0^0 := 1. See comments above. - Wolfdieter Lang, Feb 03 2015

Extensions

Edited. Comment and formula rewritten. Cross references added. - Wolfdieter Lang, Feb 03 2015