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.

A347391 Numbers k such that sigma(k) is either their sibling in Doudna tree (A005940), or one of the sibling's descendants.

Original entry on oeis.org

3, 4, 5, 15, 20, 189, 945, 2125, 6375, 9261, 46305, 401625, 19679625
Offset: 1

Views

Author

Antti Karttunen, Aug 30 2021

Keywords

Comments

Numbers k > 1 such that nearest common ancestor of k and sigma(k) in Doudna tree is the parent of k, and sigma(k) is not a descendant of k.
Any hypothetical odd term x in A005820 (triperfect numbers) would also be a member of this sequence. This is illustrated in the following diagram which shows how the neighborhood of such x would look like in the Doudna tree (A005940). If m (the parent of x, x = A003961(m), m = A064989(x)) is even, then x is a multiple of 3, while if m is odd, then 3 does not divide x. Because the abundancy index decreases when traversing leftwards in the Doudna tree, m must be a term of A068403. Both x and m would also need to be squares, by necessity.
.
<--A003961-- m ---(*2)--->
.............../ \...............
/ \
/ \
x 2m
/ \ / \
etc.../ \.....2x sigma(x) = 3x..../ \.....4m
/ \ / \ / \
etc. etc. etc. \ / etc.
\ /
6x 9x = sigma(2x)
/ \ / \
etc. \ etc. etc.
\
12x = sigma(3x) if m odd.
.
From the diagram we also see that 2x would then need to be a term of A347392 (as well as that of A159907 and also in A074388, thus sqrt(x) should be a term of A097023), and furthermore, if x is not a multiple of 3 (i.e., when m is odd), then sigma(3*x) = 4*sigma(x) = 4*(3*x), thus 3*x = sigma(x) would be a term of A336702 (particularly, in A027687) and x would be a term of A323653.
Moreover, any odd square x in this sequence (for which sigma(x) would also be odd), would have an abundancy index of at least three (sigma(x)/x >= 3). See comments in A347383.
Note how 401625 = 6375 * 63 = 945 * 425, 46305 = 945 * 49, 9261 = 189 * 49, 6375 = 2125 * 3, 945 = 189 * 5 = 15 * 63 and 9261*2125 = 19679625. It seems that when the multiplicands are coprime, then they are both terms of this sequence, e.g. 2125 and 3, 189 and 5, 2125 and 9261.
From Antti Karttunen, Jul 10 2024: (Start)
Regarding the observation above, for two coprime odd numbers x, y, if both are included here because sigma(x) = 2^a * A064989(x) and sigma(y) = 2^b * A064989(y), then also their product x*y is included because in that case sigma(x*y) = 2^(a+b) * A064989(x*y).
Also, for two coprime odd numbers x, y, if both are included here because sigma(x) = A065119(i) * x and sigma(y) = A065119(j) * y, then also their product x*y is included because sigma(x*y) = A065119(k) * x*y, where A065119(k) = A065119(i)*A065119(j). The existence of such numbers (that would include odd triperfect and odd 6-perfect numbers, see A046061) is so far hypothetical, none is known.
It is not possible that the odd x is in this sequence if sigma(x) = k*A003961^e(x) and e = A061395(k)-2 >= 1.
Note that all odd terms < 2^33 here are some of the exponentially odd divisors of 19679625 (see A374199, also A374463 and A374464).
(End)
Question: from a(6) = 189 onward, are the rest of terms all in A347390?
Conjecture: sequence is finite.
If it exists, a(14) > 2^33.

Examples

			Sigma(3) = 4 is located as the sibling of 3 in the Doudna-tree (see the illustration in A005940), thus 3 is included in this sequence.
Sigma(4) = 7 is located as a grandchild of 3 (which is the sibling of 4) in the Doudna-tree, thus 4 is included in this sequence.
Sigma(5) = 6 is located as the sibling of 5 in the Doudna-tree, thus 5 is included in this sequence.
189 (= 3^3 * 7) is a term, as sigma(189) = 320, and 320 occurs as a descendant of 80 (which is the right sibling of 189) in the Doudna tree, as illustrated below:
.
             40
            /  \
   A003961 /    \ *2
          /      \
        189       80
        / \      / \
     etc   etc etc  160
                   / \
                 etc  320
                     / \
                   etc. etc.
.
945 (= 3^3 * 5 * 7) is a term, as sigma(945) = 1920, and 1920 occurs as a descendant of 240, which is the right sibling of 945 in the Doudna tree, as illustrated below:
            120
            /  \
   A003961 /    \ *2
          /      \
        945       240
        / \      / \
     etc   etc  etc  480
                   / \
                 etc  960
                     / \
                   etc. 1920
                        / \
                     etc. etc.
		

Crossrefs

Programs

  • PARI
    isA347391(n) = (1==A347381(n));
    
  • PARI
    A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
    A252463(n) = if(!(n%2),n/2,A064989(n));
    isA347391(n) = if(1==n,0,my(m=A252463(n), s=sigma(n)); while(s>m, if(s==n, return(0)); s = A252463(s)); (s==m));

A348738 Numbers k for which A326042(k) < k, where A326042(n) = A064989(sigma(A003961(n))).

Original entry on oeis.org

2, 3, 5, 6, 7, 8, 10, 11, 13, 14, 15, 17, 19, 20, 21, 22, 23, 24, 26, 27, 28, 29, 30, 31, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 45, 46, 47, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 80, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 97
Offset: 1

Views

Author

Antti Karttunen, Nov 02 2021

Keywords

Comments

Claim: If there is an odd term y of A336702 larger than one, and it is the least one of such terms, then it should satisfy condition that for all nontrivial unitary divisor pairs d and x/d of x = A064989(y) [with gcd(d,x/d) = 1, 1 < d < x], the other divisor should reside in this sequence, and the other divisor in A348739. Proof: Applying A064989 to the odd terms of A336702 gives the fixed points of A326042. Suppose there are other odd terms in A336702 in addition to its initial 1, and let y be the least of these odd terms > 1 and x = A064989(y). Because A326042 (from here on indicated with f) is multiplicative, it follows that if we take any two nontrivial unitary divisors a and b of x, with x = a*b, gcd(a,b) = 1, 1 < a,b < x, then f(a)*f(b) = f(x) = x. Because f(x)/x = 1, we must have f(a)/a * f(b)/b = 1, as also the ratio f(n)/n is multiplicative. But f(a)/a and f(b)/b cannot be equal to 1, because then a and b would also be fixed by f, which contradicts our assumption that x were the least such fixed point larger than one. Therefore f(a) < a and f(b) > b, or vice versa. See also the comments in A348930, A348933.
Moreover, all odd perfect numbers (a subsequence of A336702), if such numbers exist, should also satisfy the same condition, regardless of whether they are the least of such numbers or not, because having a non-deficient proper divisor will push the abundancy index (ratio sigma(n)/n) of any number over 2. That is, for any such pair of nontrivial unitary divisors d and x/d, both A003961(d) and A003961(x/d) should be deficient, i.e., neither one should be in A337386. See also the condition given in A347383.
Terms that occur also in A337386 are: 120, 240, 360, 420, 480, 504, 540, 600, 630, ...

Crossrefs

Positions of positive terms in A348736, positions of 1's in A348737 (characteristic function).
Almost complement of A348739.
Subsequences: A000040, A374464 (after its initial 1).
Cf. also A348930, A348933.

Programs

  • Mathematica
    f1[2, e_] := 1; f1[p_, e_] := NextPrime[p, -1]^e; s1[1] = 1; s1[n_] := Times @@ f1 @@@ FactorInteger[n]; f2[p_, e_] := NextPrime[p]^e; s2[1] = 1; s2[n_] := Times @@ f2 @@@ FactorInteger[n]; Select[Range[100], s1[DivisorSigma[1, s2[#]]] < # &] (* Amiram Eldar, Nov 04 2021 *)
  • PARI
    A003961(n) = my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); \\ From A003961
    A064989(n) = {my(f); f = factor(n); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f)};
    A326042(n) = A064989(sigma(A003961(n)));
    isA348738(n) = (A326042(n)
    				

A374463 Numbers k such that the odd parts of sigma(k) and A064989(k) are equal, where A064989 shifts the prime factorization one step towards lower primes, and sigma is the sum of divisors function.

Original entry on oeis.org

1, 3, 5, 15, 189, 945, 2125, 6375, 9261, 42550, 46305, 127650, 401625, 936100, 1191400, 2808300, 3574200, 8041950, 19679625, 22469750, 58378600, 58506250, 67409250, 175135800, 175518750, 176922900, 394055550, 494334500
Offset: 1

Views

Author

Antti Karttunen, Jul 11 2024

Keywords

Comments

Odd terms after 1 form a subsequence of A347391.
If x and y are included, and they are coprime (gcd(x,y) = 1), then x*y is also included.

Examples

			945 = 3^3 * 5 * 7 is included as A064989(945) = 2^3 * 3 * 5 = 120, and sigma(945) = 1920 = 2^7 * 3 * 5, with A000265(120) = A000265(1920) = 15.
		

Crossrefs

Programs

Formula

{k | A000265(sigma(k)) = A000265(A064989(k))}.
Showing 1-3 of 3 results.