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

A369054 Number of representations of n as a sum (p*q + p*r + q*r) with three odd primes p <= q <= r.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 1
Offset: 0

Views

Author

Antti Karttunen, Jan 20 2024

Keywords

Comments

Number of solutions to n = x', where x' is the arithmetic derivative of x (A003415), and x is a product of three odd primes (not all necessarily distinct, A046316).
See the conjecture in A369055.

Examples

			a(27) = 1 as 27 can be expressed in exactly one way in the form (p*q + p*r + q*r), with p, q, r all being 3 in this case, as 27 = (3*3 + 3*3 + 3*3).
a(311) = 5 as 311 = (3*5 + 3*37 + 5*37) = (3*7 + 3*29 + 7*29) = (3*13 + 3*17 + 13*17) = (5*7 + 5*23 + 7*23) = (7*11 + 7*13 + 11*13). Expressed in the terms of arithmetic derivatives, of the A099302(311) = 8 antiderivatives of 311 [366, 430, 494, 555, 609, 663, 805, 1001], only the last five are products of three odd primes: 555 = 3*5*37, 609 = 3*7*29, 663 = 3*13*17, 805 = 5*7*23, 1001 = 7 * 11 * 13.
		

Crossrefs

Cf. A369055 [quadrisection, a(4n-1)], and its trisections A369460 [= a((12*n)-9)], A369461 [= a((12*n)-5)], A369462 [= a((12*n)-1)].
Cf. A369251 (positions of terms > 0), A369464 (positions of 0's).
Cf. A369063 (positions of records), A369064 (values of records).
Cf. A369241 [= a(2^n - 1)], A369242 [= a(n!-1)], A369245 [= a(A006862(n))], A369247 [= a(3*A057588(n))].

Programs

  • PARI
    \\ Use this for building up a list up to a certain n. We iterate over weakly increasing triplets of odd primes:
    A369054list(up_to) = { my(v = [3,3,3], ip = #v, d, u = vector(up_to)); while(1, d = ((v[1]*v[2]) + (v[1]*v[3]) + (v[2]*v[3])); if(d > up_to, ip--, ip = #v; u[d]++); if(!ip, return(u)); v[ip] = nextprime(1+v[ip]); for(i=1+ip,#v,v[i]=v[i-1])); };
    v369054 = A369054list(100001);
    A369054(n) = if(!n,n,v369054[n]);
    
  • PARI
    \\ Use this for computing the value of arbitrary n. We iterate over weakly increasing pairs of odd primes:
    A369054(n) = if(3!=(n%4),0, my(v = [3,3], ip = #v, r, c=0); while(1, r = (n-(v[1]*v[2])) / (v[1]+v[2]); if(r < v[2], ip--, ip = #v; if(1==denominator(r) && isprime(r),c++)); if(!ip, return(c)); v[ip] = nextprime(1+v[ip]); for(i=1+ip,#v,v[i]=v[i-1])));

Formula

a(n) = Sum_{i=1..A002620(n)} A369058(i)*[A003415(i)==n], where [ ] is the Iverson bracket.
For n >= 2, a(n) <= A099302(n).

A369245 Number of representations of the n-th Euclid number, A002110(n) + 1, as a sum (p*q + p*r + q*r) with three odd primes p <= q <= r. (Definition implies that p=3 and q > 3).

Original entry on oeis.org

0, 0, 0, 0, 1, 2, 2, 1, 0, 0, 1, 1, 0, 1, 2, 0, 1
Offset: 0

Views

Author

Antti Karttunen, Jan 22 2024

Keywords

Comments

Number of representations of the n-th Euclid number, A002110(n) + 1, as a sum of the form 3*(p+q) + p*q, where p and q are odd primes.
Question: Will there be an eventual growth spurt for this sequence? Even though all solutions must be multiples of 3 (but not of 9), because A006862(n) == 1 (mod 3), for n > 1, and the solutions belong to a set listed by A369461.
Similar sequence A369242 grows more vigorously because A033312(n) == -1 (mod 3) for n >= 3, thus allowing non-multiples of 3 as solutions. See comments in A369252.

Examples

			a(4) = 1 as there exists a natural number 399 = 3 * 7 * 19, whose arithmetic derivative (indicated with 399', see A003415) is computed as ((3*7) + (3*19) + (7*19)) = 211 = 1 + prime(4)# = A006862(4), and because 399 is the unique term in A046316 that satisfies the condition.
a(17) >= 1 because there exists (at least one) solution k = 4903038892893242229501 = 3 * 17 * 96138017507710631951 with A003415(k) = 1+A002110(17).
For other cases, see examples in A369246.
		

Crossrefs

Cf. also A116979, A369000, A369239 for similar counts, also A369241, A369242 and A369247.

Programs

  • PARI
    \\ Needs also program from A369054.
    A002110(n) = prod(i=1,n,prime(i));
    A369245(n) = A369054(A002110(n)+1);
    
  • PARI
    \\ Optimized version of above, employs the fact that solutions must all be multiples of 3. Outputs also terms for A369246.
    search_for_3k1_cases(n) = if(3!=(n%4),0, my(p = 5, q, c=0); while(1, q = (n-(3*p)) / (3+p); if(q < p, return(c), if(1==denominator(q) && isprime(q),c++; write("b369246_by_search_order_to.txt", n, " ", 3*p*q))); p = nextprime(1+p)));
    A002110(n) = prod(i=1,n,prime(i));
    A369245(n) = search_for_3k1_cases(A002110(n)+1);

Formula

a(n) = A369054(A006862(n)).

A369241 Number of representations of 2^n - 1 as a sum (p*q + p*r + q*r) with three odd primes p <= q <= r.

Original entry on oeis.org

0, 0, 0, 0, 0, 0, 0, 0, 2, 1, 1, 0, 2, 1, 3, 0, 2, 2, 6, 0, 13, 1, 13, 0, 15, 0, 24, 1, 49, 4, 47, 0, 156, 6, 129, 0, 441, 1, 616
Offset: 0

Views

Author

Antti Karttunen, Jan 21 2024

Keywords

Comments

Any solutions for odd cases must have p = 3, with q and r > 3, because A000225(2n-1) == 1 (mod 3), while on even n, 2^n - 1 is a multiple of 3. This explains why the odd bisection grows much more sluggishly than the even bisection.
Question 2: Is there an infinite number of 0's in this sequence? See also comments in A369055.

Crossrefs

Programs

  • PARI
    A369054(n) = if(3!=(n%4),0, my(v = [3,3], ip = #v, r, c=0); while(1, r = (n-(v[1]*v[2])) / (v[1]+v[2]); if(r < v[2], ip--, ip = #v; if(1==denominator(r) && isprime(r),c++)); if(!ip, return(c)); v[ip] = nextprime(1+v[ip]); for(i=1+ip,#v,v[i]=v[i-1])));
    search_for_3k1_cases(n) = if(3!=(n%4), 0, my(p = 5, q, c=0); while(1, q = (n-(3*p)) / (3+p); if(q < p, return(c), if(1==denominator(q) && isprime(q), c++; write("b369241_by_solutions_of_odd_bisection_to.txt", n, " ", 3*p*q))); p = nextprime(1+p)));
    A369241(n) = if(n%2, search_for_3k1_cases((2^n)-1), A369054((2^n)-1));

Formula

a(n) = A369054(A000225(n)).
For n >= 2, a(n) = A369055(2^(n-2)).
Showing 1-3 of 3 results.