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

A337372 Primitively primeshift-abundant numbers: Numbers that are included in A246282 (k with A003961(k) > 2k), but none of whose proper divisors are.

Original entry on oeis.org

4, 6, 9, 10, 14, 15, 21, 35, 39, 49, 57, 69, 91, 125, 242, 275, 286, 325, 338, 363, 418, 425, 442, 475, 494, 506, 561, 575, 598, 646, 682, 715, 722, 725, 754, 775, 782, 806, 845, 847, 867, 874, 925, 957, 962, 1023, 1025, 1045, 1054, 1058, 1066, 1075, 1105, 1118, 1175, 1178, 1221, 1222, 1235, 1265, 1309, 1325, 1334, 1353
Offset: 1

Views

Author

Antti Karttunen, Aug 27 2020

Keywords

Comments

Numbers k whose only divisor in A246282 is k itself, i.e., A003961(k) > 2k, but for none of the proper divisors d|k, dA003961(d) > 2d.
Question: Do the odd terms in A326134 all occur here? Answer is yes, if the following conjecture holds: This is a subsequence of A263837, nonabundant numbers. In other words, we claim that any abundant number k (A005101) has A337345(k) > 1 and thus is a term of A341610. (The conjecture indeed holds. See the proof below).
From Antti Karttunen, Dec 06 2024: (Start)
Observation 1: The thirteen initial terms (4, 6, 9, ..., 69, 91) are only semiprimes in A246282, all other semiprimes being in A246281 (but none in A341610), and there seems to be only 678 terms m with A001222(m) = 3, from a(14) = 125 to the last one of them, a(2691) = 519963. There are more than 150000 terms m with A001222(m) = 4. In general, there should be only a finite number of terms m for any given k = A001222(m). Compare for example with A287728.
Observation 2: The intersection with A005101 (and thus also with A091191) is empty, which then implies the claims made in the sequences A378662, A378664, from which further follows that there are no 1's present in any of these sequences: A378658, A378736, A378740.
(End)
Proof of the latter observation by Jianing Song, Dec 11 2024: (Start)
Let's write p' for the next prime after the prime p. Also, write Q(n) = A003961(n)/sigma(n) which is multiplicative.
Proposition: For n > 1 not being a prime nor twice a prime, n has a factor p such that Q(n) > p'/p.
This implies that if n is abundant [including any primitively abundant n in A091191], then n has a factor p such that A003961(n/p)/(n/p) = (A003961(n)/n)/(p'/p) > sigma(n)/n [which is > 2 because n is abundant], so n/p is in A246282, meaning that n cannot be in this sequence.
Proof. We see that 1 <= Q(p) <= Q(p^2) <= ..., which implies that if n verifies the proposition, then every multiple of n also verifies it. Since n = p^2 > 4 and n = 8 verify the proposition, it suffices to consider the case where n = pq is the product of two distinct odd primes. Suppose WLOG that p < q, so q >= p', then using q/(q+1) >= p'/(p'+1) we have
Q(n) = p'q'/((p+1)(q+1)) >= p'^2*q'/(q(p+1)(p'+1)) > (p'^2-1)*q'/(q(p+1)(p'+1)) = (p'-1)/(p+1) * q'/q >= q'/q.
(End)

Examples

			14 = 2*7 is in the sequence as setting every prime to the next larger prime gives 3*11 = 33 > 28 = 2*14. Doing so for any proper divisor d of 14 gives a number < 2 * d. - _David A. Corneth_, Dec 07 2024
		

Crossrefs

Setwise difference A246282 \ A341610.
Positions of ones in A337345 and in A341609 (characteristic function).
Subsequence of A263837 and thus also of A341614.
Cf. also A005101, A091191, A326134.
Cf. also A337543.

Programs

  • Mathematica
    Block[{a = {}, b = {}}, Do[If[2 i < Times @@ Map[#1^#2 & @@ # &, FactorInteger[i] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] - Boole[i == 1], AppendTo[a, i]; If[IntersectingQ[Most@ Divisors[i], a], AppendTo[b, i]]], {i, 1400}]; Complement[a, b]] (* Michael De Vlieger, Feb 22 2021 *)
  • PARI
    A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    A252742(n) = (A003961(n) > (2*n));
    A337346(n) = sumdiv(n,d,(dA252742(d));
    isA337372(n) = ((1==A252742(n))&&(0==A337346(n)));
    
  • PARI
    is_A337372 = A341609;
    
  • PARI
    \\ See Corneth link

Formula

{k: 1==A337345(k)}.

A337345 Number of divisors d of n for which A003961(d) > 2*d, where A003961 is fully multiplicative with a(p) = nextprime(p).

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Aug 27 2020

Keywords

Comments

Number of terms of A246282 that divide n.
Number of divisors d of n for which A048673(d) > d.

Crossrefs

Inverse Möbius transform of A252742.
Cf. A003961, A048673, A246282, A337346, A337372 (positions of ones), A341609 (their characteristic function), A341610 (positions of terms > 1), A378658 [= a(A091191(n))], A378662, A378663.
Cf. also A080224, A337541, A341620.

Programs

  • PARI
    A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    A337345(n) = sumdiv(n,d,A003961(d)>(d+d));

Formula

a(n) = Sum_{d|n} A252742(d).
a(n) = A337346(n) + A252742(n).
From Antti Karttunen, Dec 10 2024: (Start)
a(n) = 1 <=> A341609(n) = 1.
a(n) = A378662(n) + A080224(n) = A378663(n) + A341620(n).
(End)

A378662 Number of divisors d of n such that sigma(d) <= 2*d < A003961(d), where A003961 is fully multiplicative with a(p) = nextprime(p).

Original entry on oeis.org

0, 0, 0, 1, 0, 1, 0, 2, 1, 1, 0, 2, 0, 1, 1, 3, 0, 2, 0, 2, 1, 0, 0, 3, 0, 0, 2, 3, 0, 3, 0, 4, 0, 0, 1, 3, 0, 0, 1, 3, 0, 3, 0, 2, 3, 0, 0, 4, 1, 2, 0, 2, 0, 3, 0, 4, 1, 0, 0, 4, 0, 0, 3, 5, 0, 1, 0, 2, 1, 3, 0, 4, 0, 0, 2, 2, 0, 2, 0, 4, 3, 0, 0, 5, 0, 0, 0, 3, 0, 5, 1, 2, 0, 0, 0, 5, 0, 3, 2, 3, 0, 1, 0, 3, 4
Offset: 1

Views

Author

Antti Karttunen, Dec 06 2024

Keywords

Comments

Number of terms of A341614 that divide n.
Claim: a(n) > 0 if and only if A003961(n) > 2*n [i.e., n is in A246282]. That a(n) must be zero when n is in A246281 is obvious, as is also that a(n) > 0 when n is a term of A341614 [as then A378664(n) = n], but that a(n) > 0 for all abundant numbers (A005101) is slightly less clear. So the claim boils down to this: All abundant numbers have at least one (by necessity a proper) divisor d|n such that it is in A341614, i.e., sigma(d) <= 2*d < A003961(d), i.e., that for abundant numbers n, A337345(n) is always strictly greater than A080224(n). Equivalently, of the all nonabundant divisors d of an abundant number, at least one is primeshift-abundant, i.e., A003961(d) > 2*d. This has been proved Dec 11 2024 by Jianing Song in A337372. The claim given in A378658 also follows from that proof.

Crossrefs

Inverse Möbius transform of A341612.
Cf. A246281 (positions of 0's), A246282 (of terms > 0).
Cf. also A337372, A378658.

Programs

  • Mathematica
    Table[Length@ Select[Divisors[n], DivisorSigma[1, #] <= 2 # < (Times @@ Map[Power @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Prime[PrimePi[p] + 1], e}] - Boole[# == 1]) &], {n, 105}] (* Michael De Vlieger, Dec 06 2024 *)
  • PARI
    A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    A341612(n) = ((sigma(n)<=(2*n))&&((2*n)<A003961(n)));
    A378662(n) = sumdiv(n,d,A341612(d));

Formula

a(n) = Sum_{d|n} A341612(d).
a(n) = A337345(n) - A080224(n).
a(n) = A080225(n) + A378663(n).

A378740 Distinct values of A378664(k) in the order of appearance, when k ranges over those primitively abundant numbers k for which A378664(k) is less than the largest proper divisor of k.

Original entry on oeis.org

6, 68, 170, 290, 646, 682, 754, 370, 410, 430, 470, 530, 4756, 6844, 8236, 30566, 10730, 11890, 12470, 43214, 45182, 46342, 47386, 15170, 15370, 54094, 17110, 17690, 62186, 62894, 19430, 75686, 79178, 39530, 89914, 41890, 43070, 95282, 97966, 46610, 103402, 47570, 106018, 107602, 107666, 48970, 109798, 111386, 51830
Offset: 1

Views

Author

Antti Karttunen, Dec 07 2024

Keywords

Comments

Term 1 does not occur in this sequence, thus all terms are in A341614, and also in A246282. See A378658, A378662, A378664, A378736 and A337372 for a proof.

Crossrefs

Cf. A246282, A337372, A378658, A378662, A378664, A378736, A378738, A378739, A378741 (corresponding numbers of A091191).

Programs

  • PARI
    \\ Uses program from A378738:
    memoA378740 = Map();
    k=0; n=0; while(k<200, n++; if(is_A378738(n), t=A378664(n); if(!mapisdefined(memoA378740, t), mapput(memoA378740, t, n); k++; print1(t, ", "); write("b378740.txt", k, " ", t))));

Formula

a(n) = A378664(A378741(n)), a(n)| A378741(n).
Showing 1-4 of 4 results.