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

A378664 Greatest divisor d of n such that sigma(d) <= 2*d < A003961(d), or 1 if no such divisor exists, where A003961 is fully multiplicative with a(p) = nextprime(p).

Original entry on oeis.org

1, 1, 1, 4, 1, 6, 1, 8, 9, 10, 1, 6, 1, 14, 15, 16, 1, 9, 1, 10, 21, 1, 1, 8, 1, 1, 27, 28, 1, 15, 1, 32, 1, 1, 35, 9, 1, 1, 39, 10, 1, 21, 1, 44, 45, 1, 1, 16, 49, 50, 1, 52, 1, 27, 1, 28, 57, 1, 1, 15, 1, 1, 63, 64, 1, 6, 1, 68, 69, 35, 1, 9, 1, 1, 75, 76, 1, 39, 1, 16, 81, 1, 1, 28, 1, 1, 1, 44, 1, 45, 91, 92
Offset: 1

Views

Author

Antti Karttunen, Dec 06 2024

Keywords

Comments

Largest term of {1} U A341614 that divides n.

Crossrefs

Positions of fixed points (where a(n)=n) is given by {1} U A341614.
Cf. A246281 (positions of 1's), A246282 (of terms > 0), A005101 (of terms that are neither 1 nor fixed points).

Programs

  • Mathematica
    Table[If[Length[#] == 0, 1, Max[#]] &@ Select[Divisors[n], DivisorSigma[1, #] <= 2 # < (Times @@ Map[Power @@ # &, FactorInteger[#] /. {p_, e_} /; e > 0 :> {Prime[PrimePi[p] + 1], e}] - Boole[# == 1]) &], {n, 92}] (* 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)));
    A378664(n) = { fordiv(n,d,if(A341612(n/d), return(n/d))); (1); };

Formula

a(n) <= A378665(n).

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

A378735 Greatest non-abundant divisor of the n-th abundant number.

Original entry on oeis.org

6, 9, 10, 8, 15, 9, 10, 21, 16, 27, 28, 15, 33, 35, 9, 39, 16, 28, 44, 45, 32, 50, 51, 52, 27, 28, 57, 15, 63, 44, 69, 35, 16, 75, 52, 32, 81, 28, 87, 44, 45, 93, 64, 98, 99, 50, 68, 52, 105, 27, 110, 111, 32, 76, 117, 16, 123, 63, 129, 130, 44, 135, 136, 92, 35, 141, 32, 147, 75, 152, 153, 154, 52, 159, 64, 81, 165
Offset: 1

Views

Author

Antti Karttunen, Dec 06 2024

Keywords

Crossrefs

Cf. also A378736 and the points of differences A378737.

Programs

  • PARI
    A294935(n) = (sigma(n)<=(2*n));
    A378665(n) = { fordiv(n,d,if(A294935(n/d), return(n/d))); (1); };
    k=0; n=0; while(k<105, n++; if(!A294935(n), k++; print1(A378665(n),", ")));

Formula

a(n) = A378665(A005101(n)).
a(n) >= A378736(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)}.
Showing 1-6 of 6 results.