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.

A373145 a(n) = gcd(A003415(n), A276085(n)), where A003415 is the arithmetic derivative and A276085 is the primorial base log-function.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 3, 2, 7, 1, 4, 1, 1, 8, 4, 1, 1, 1, 8, 2, 1, 1, 1, 2, 1, 3, 32, 1, 1, 1, 5, 2, 1, 12, 6, 1, 1, 8, 1, 1, 1, 1, 4, 1, 1, 1, 2, 2, 1, 4, 8, 1, 1, 8, 1, 2, 1, 1, 2, 1, 1, 17, 6, 6, 1, 1, 8, 2, 1, 1, 1, 1, 1, 1, 16, 6, 1, 1, 2, 4, 1, 1, 2, 2, 1, 8, 1, 1, 1, 20, 4, 2, 1, 12, 1, 1, 1, 1, 14, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Antti Karttunen, May 26 2024

Keywords

Crossrefs

Cf. A368998 (positions of even terms), A368999 (of odd terms), A373144 (of multiples of 3).
Cf. also A327858.

Programs

  • PARI
    A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
    A002110(n) = prod(i=1,n,prime(i));
    A276085(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*A002110(primepi(f[k, 1])-1)); };
    A373145(n) = gcd(A003415(n), A276085(n));

Formula

a(n) = gcd(A003415(n), A373146(n)) = gcd(A276085(n), A373146(n)).
For n > 1, a(n) = gcd(A276085(n), A373147(n)) = gcd(A003415(n), A373148(n)).

A369650 Numbers k such that A003415(k) = A276085(k), where A003415 is the arithmetic derivative, and A276085 is the primorial base log-function.

Original entry on oeis.org

1, 2, 10, 15, 28, 5005
Offset: 1

Views

Author

Antti Karttunen, Feb 05 2024

Keywords

Comments

Intersection with A048103 gives the fixed points (1, 2, 10, 15, 5005, ...) of A327859. Question: Does that set preclude nonsquarefree numbers? Certainly it does not contain any multiples of 9. See also comments in A328110.
If k == 2 (mod 4), then both A003415(k) and A276085(k) are odd, and the latter is of the form 4m+1 (if k has an odd number of prime factors), or of the form 4m+3 (if k has an even number of prime factors). Therefore, for k of the form 4m+2 to be included in this sequence, a necessary condition is that it must be either in the intersection of A026424 and A358772 (like, for example, 2 is) or in A369668 (the intersection of A028260 and A358774), like for example, 10 is.
If k is odd, then A276085(k) is even, and for A003415(k) to be even with k odd, then k has to be in A046337 (odd numbers with an even number of prime factors, counted with multiplicity). But A276085(A046337(n)) == 0 (mod 4) for all n, so also A003415(k) has to be a multiple of 4, so k has to be in A360110 (itself a subsequence of A369002), like for example k=15 and k=5005 are.
If it exists, a(7) > 2^19.

Examples

			As 5005 = 5*7*11*13, A003415(5*7*11*13) = (5*7*11) + (5*7*13) + (5*11*13) + (7*11*13) = 2556 = 2^2 * 3^2 * 71 = A276085(5005) = A002110(2) + A002110(3) + A002110(4) + A002110(5) [as 5, 7, 11 and 13 are prime(3) .. prime(6)], therefore 5005 is included in this sequence.
		

Crossrefs

Positions of 0's in A373146.
Intersection of A373487 and A373490.
Cf. also A351228.

Programs

  • PARI
    A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
    A276085(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*prod(i=1, primepi(f[k, 1]-1), prime(i))); };
    isA369650(n) = (A003415(n) == A276085(n));
Showing 1-2 of 2 results.