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-10 of 11 results. Next

A378738 Primitively abundant numbers k for which A378665(k) > A378664(k).

Original entry on oeis.org

66, 102, 174, 186, 222, 246, 258, 282, 318, 354, 366, 402, 426, 438, 474, 498, 534, 582, 606, 618, 642, 654, 678, 748, 762, 786, 822, 834, 894, 906, 942, 978, 1002, 1038, 1074, 1086, 1146, 1158, 1182, 1194, 1266, 1338, 1362, 1374, 1398, 1434, 1446, 1506, 1542, 1578, 1614, 1626, 1662, 1686, 1698, 1758, 1842, 1866
Offset: 1

Views

Author

Antti Karttunen, Dec 07 2024

Keywords

Comments

Subsequence of A378737: 1496 is its first term that does not occur here.
Equal to primitively abundant numbers k such that A032742(k) > A378664(k), because for primitively abundant numbers the greatest non-abundant divisor is the largest proper divisor, A378665(k) = A032742(k).
Question: What is the asymptotic density of these numbers among A091191? Does it tend to 1?
Conjecture: A001222(a(n)) = 3 <=> 3|a(n).

Examples

			Examples given in A378737 for 66, 748, 1866, and 1870 all work also here, because those four numbers are all in A091191.
		

Crossrefs

Intersection of A091191 and A378737.
Cf. A001222, A032742, A294930, A337372, A341612, A341614, A378664, A378665, A378735, A378736, A378739 [= A378664(a(n))], A378741 (subsequence), A378742 (subsequence after its initial term).

Programs

  • 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)));
    A378664(n) = { fordiv(n,d,if(A341612(n/d), return(n/d))); (1); };
    A032742(n) = if(1==n,n,n/vecmin(factor(n)[,1]));
    is_A091191(n) = if(sigma(n)<=2*n, 0, fordiv(n,d,if(d2*d, return(0))); (1));
    is_A378738(n) = (is_A091191(n) && (A378664(n)!=A032742(n)));

Formula

{k such that A294930(k) = 1 and A032742(k) > A378664(k)}.

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).

A378737 Abundant numbers k for which A378665(k) > A378664(k).

Original entry on oeis.org

66, 102, 174, 186, 222, 246, 258, 282, 318, 354, 366, 402, 426, 438, 474, 498, 534, 582, 606, 618, 642, 654, 678, 748, 762, 786, 822, 834, 894, 906, 942, 978, 1002, 1038, 1074, 1086, 1146, 1158, 1182, 1194, 1266, 1338, 1362, 1374, 1398, 1434, 1446, 1496, 1506, 1542, 1578, 1614, 1626, 1662, 1686, 1698, 1758, 1842
Offset: 1

Views

Author

Antti Karttunen, Dec 06 2024

Keywords

Comments

For most of these numbers k, A378664(k) = 6. Note that A003961(6) = A003961(2*3) = 3*5 = 15 > 2*6, while sigma(6) = 12, making 6 non-abundant. Note that there seems to be only a finite amount (namely 13, see A337372) of such "stopper semiprimes" that would prevent of A378736 obtaining value 1.
Other possible values that A378664 obtains on these numbers are for example 68, 136, 170, 256, 290, 370, 410, 430, 470, 530, 646, 682, 754, 1276, 1292, 1364, 1508, 1628, 1804, 1892, 2068, 2332, 4756, 6844, 10846, 15334, etc. See A378740, which contains some of these.

Examples

			66 is a term as A378665(66) = 33, but A378664(66) = 6.
748 is a term as A378665(748) = 374, but A378664(748) = 68.
1866 is a term as A378665(1866) = 933, but A378664(1866) = 6.
1870 is a term as A378665(1870) = 935, but A378664(1870) = 170.
		

Crossrefs

Subsequence of A005101.
Cf. A378738 (subsequence).

Programs

Formula

{A005101(i) for such indices i where A378735(i) > A378736(i)}.

A378741 Terms k of A378738 for which A378664(k) obtains novel values.

Original entry on oeis.org

66, 748, 1870, 3190, 3230, 3410, 3770, 4070, 4510, 4730, 5170, 5830, 52316, 75284, 90596, 152830, 182410, 202130, 211990, 216070, 225910, 231710, 236930, 257890, 261290, 270470, 290870, 300730, 310930, 314470, 330310, 378430, 395890, 434830, 449570, 460790, 473770, 476410, 489830, 512710, 517010, 523270, 530090
Offset: 1

Views

Author

Antti Karttunen, Dec 07 2024

Keywords

Comments

a(n) is a multiple of A378740(n), by definition.

Crossrefs

Cf. A378740 [= A378664(a(n))].
Subsequence of A378738, which is a subsequence of A091191.

Programs

A378739 A378664(k) computed for those primitively abundant numbers k for which A378664(k) is less than the largest proper divisor of k.

Original entry on oeis.org

6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 68, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 170, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 290, 646, 6, 6, 6, 6, 682, 6
Offset: 1

Views

Author

Antti Karttunen, Dec 07 2024

Keywords

Comments

Of the initial 20000 terms, 19803 are 6's.

Crossrefs

Cf. A091191, A378664, A378738, A378740 (distinct terms in the order of appearance), A378741.

Formula

a(n) = A378664(A378738(n)).

A378742 Primitively abundant numbers k for which A378664(k) = 6, where A378664 is the greatest divisor d of n such that sigma(d) <= 2*d < A003961(d), or 1 if no such divisor exists.

Original entry on oeis.org

12, 66, 102, 174, 186, 222, 246, 258, 282, 318, 354, 366, 402, 426, 438, 474, 498, 534, 582, 606, 618, 642, 654, 678, 762, 786, 822, 834, 894, 906, 942, 978, 1002, 1038, 1074, 1086, 1146, 1158, 1182, 1194, 1266, 1338, 1362, 1374, 1398, 1434, 1446, 1506, 1542, 1578, 1614, 1626, 1662, 1686, 1698, 1758, 1842, 1866
Offset: 1

Views

Author

Antti Karttunen, Dec 07 2024

Keywords

Comments

Apparently all the terms are of the form 6*p, where p is any prime except one of the 3, 5, 7, 13, 19, 23.

Crossrefs

Cf. A378664.
After the initial term, a subsequence of A378738.
Subsequence of A008588 and of A091191.

Programs

  • 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)));
    A378664(n) = { fordiv(n,d,if(A341612(n/d), return(n/d))); (1); };
    is_A091191(n) = if(sigma(n)<=2*n, 0, fordiv(n,d,if(d2*d, return(0))); (1));
    is_A378742(n) = (is_A091191(n) && (A378664(n)==6));

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)}.

A341614 Numbers k such that sigma(k) <= 2k < A003961(k).

Original entry on oeis.org

4, 6, 8, 9, 10, 14, 15, 16, 21, 27, 28, 32, 35, 39, 44, 45, 49, 50, 52, 57, 63, 64, 68, 69, 75, 76, 81, 91, 92, 98, 99, 105, 110, 116, 117, 124, 125, 128, 130, 135, 136, 147, 148, 152, 153, 154, 164, 165, 170, 171, 172, 175, 182, 184, 188, 189, 190, 195, 207, 212, 225, 230, 231, 232, 236, 238, 242, 243, 244, 245, 248, 250, 255, 256
Offset: 1

Views

Author

Antti Karttunen, Feb 22 2021

Keywords

Crossrefs

Intersection of A263837 and A246282 (nonabundant numbers in A246282).
Union of A000396 and A341615.
Union of A337372 and A341611 (see also A341610).
Cf. A341612 (characteristic function), A326134 (a subsequence).
Cf. also A378662, A378664.

Programs

  • Mathematica
    Select[Range[256], DivisorSigma[1, #] <= 2 # <
    Times @@ Map[#1^#2 & @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Prime[PrimePi@ p + 1], e}] &] (* Michael De Vlieger, Feb 22 2021 *)

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).

A378736 Greatest divisor d of the n-th abundant number such that sigma(d) <= 2*d < A003961(d).

Original entry on oeis.org

6, 9, 10, 8, 15, 9, 10, 21, 16, 27, 28, 15, 6, 35, 9, 39, 16, 28, 44, 45, 32, 50, 6, 52, 27, 28, 57, 15, 63, 44, 69, 35, 16, 75, 52, 32, 81, 28, 6, 44, 45, 6, 64, 98, 99, 50, 68, 52, 105, 27, 110, 6, 32, 76, 117, 16, 6, 63, 6, 130, 44, 135, 136, 92, 35, 6, 32, 147, 75, 152, 153, 154, 52, 6, 64, 81, 165, 28, 170, 171
Offset: 1

Views

Author

Antti Karttunen, Dec 06 2024

Keywords

Comments

There are no 1's in this sequence. See A378662, A378664 and A337372 for a proof.

Crossrefs

Programs

  • PARI
    A003961(n) = { my(f = factor(n)); for (i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); };
    A294935(n) = (sigma(n)<=(2*n));
    A341612(n) = ((sigma(n)<=(2*n))&&((2*n)<A003961(n)));
    A378664(n) = { fordiv(n,d,if(A341612(n/d), return(n/d))); (1); };
    k=0; n=0; while(k<20000, n++; if(!A294935(n), k++; print1(A378664(n),", ")));

Formula

a(n) = A378664(A005101(n)).
a(n) <= A378735(n).

Extensions

Unnecessary escape-clause removed from the definition by Antti Karttunen, Dec 12 2024
Showing 1-10 of 11 results. Next