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 80 results. Next

A244052 Highly regular numbers a(n) defined as positions of records in A010846: a(1) = 1 and a(n) is the least number k > a(n-1) such that A010846(k) > A010846(a(n-1)).

Original entry on oeis.org

1, 2, 4, 6, 10, 12, 18, 24, 30, 42, 60, 84, 90, 120, 150, 180, 210, 330, 390, 420, 630, 840, 1050, 1260, 1470, 1680, 1890, 2100, 2310, 2730, 3570, 3990, 4620, 5460, 6930, 8190, 9240, 10920, 11550, 13650, 13860
Offset: 1

Views

Author

Michael De Vlieger, Jun 18 2014

Keywords

Comments

Analogous to highly divisible numbers (A002182).

Examples

			a(2) = 2 because already for k = 2 A010846(2) = 2 > A010846(1) = 1.
a(3) = 4 because for k = 3  A010846(3) = 2 = A010846(2), but
for k = 4 A010846(4) = 3 > A010846(2) = 2.
a(4) = 6 because for k = 5 A010846(5) = 2 < A010846(4) = 3, but
A010846(6) = 5 > A010846(4) = 3.
		

Crossrefs

Programs

  • Mathematica
    Function[w, Map[Position[w, #][[1, 1]] &, Union@ Rest@ FoldList[Max, 0, w]]]@ Table[Count[Range@ n, k_ /; PowerMod[n, Floor@ Log2@ n, k] == 0], {n, 10^3}] (* simplest, or *)
    f[n_] := If[n == 1, 1, Length@ Function[w, ToExpression@ StringJoin["Module[{n = ", ToString@ n, ", k = 0}, Flatten@ Table[k++, ",  Most@ Flatten@ Map[{#, ", "} &, #], "]]"] &@ MapIndexed[Function[p, StringJoin["{", ToString@ Last@ p, ", 0, Log[", ToString@ First@ p, ", n/(", ToString@ InputForm[Times @@ Map[Power @@ # &, Take[w, First@ #2 - 1]]], ")]}"]]@ w[[First@ #2]] &, w]]@Map[{#, ToExpression["p" <> ToString@ PrimePi@ #]} &, FactorInteger[n][[All, 1]]]]; Function[w, Map[Position[w, #][[1, 1]] &, Union@ Rest@ FoldList[Max, 0, w]]]@ Array[f, 14000] (* Michael De Vlieger, Mar 08 2017, more efficient *)

Extensions

Edited, giving new name and example. - Wolfdieter Lang, Jun 29 2014

A372192 a(n) = A010846(n) - A373738(n).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 1, 0, 1, 1, 1, 0, 1, 0, 1, 0, 1, 0, 6, 0, 0, 1, 2, 1, 1, 0, 2, 1, 1, 0, 6, 0, 1, 1, 1, 0, 1, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 7, 0, 1, 1, 0, 1, 7, 0, 2, 1, 7, 0, 1, 0, 2, 0, 2, 1, 7, 0, 1, 0, 1, 0, 8, 1, 1, 2
Offset: 1

Views

Author

Michael De Vlieger, Aug 06 2024

Keywords

Crossrefs

Programs

  • Mathematica
    {0}~Join~Table[Count[Range[n], _?(Divisible[n, Times @@ FactorInteger[#][[All, 1]]] &)] - Floor[(1/PrimeNu[n]!)*Times @@ Map[Log[#, n] + 1 &, FactorInteger[n][[All, 1]]]], {n, 2, 120}]

A372322 a(n) = A010846(A372111(n)).

Original entry on oeis.org

1, 2, 5, 6, 5, 11, 18, 8, 16, 22, 5, 28, 13, 33, 23, 38, 11, 26, 12, 9, 58, 28, 80, 5, 30, 55, 19, 27, 19, 56, 37, 21, 27, 87, 44, 44, 48, 38, 18, 58, 42, 5, 110, 26, 112, 140, 38, 45, 32, 144, 102, 59, 5, 139, 225, 39, 44, 22, 180, 86, 114, 34, 23, 133, 41, 115
Offset: 1

Views

Author

Michael De Vlieger, May 05 2024

Keywords

Comments

Let r(x) = A010846(x), the number of m <= x such that rad(m) | x, where rad = A007947.
Let row k of A162306 contain { m : rad(m) | k, m <= k }. Thus r(k) is the length of row k of A162306.
a(n) is the length of row A372111(n) of A162306.
Analogous to A371909, which instead regards A109890 and A109735.

Examples

			Let s(x) = A372111(x) and let r(x) = A010846(x).
a(1) = 1 since r(s(1)) = r(1) = 1.
a(2) = 2 since r(s(2)) = r(3) = 2. For prime p, r(p) = card({1, p}) = 2.
a(3) = 5 since r(s(3)) = r(6) = 5. r(6) = card({1, 2, 3, 4, 6}) = 5.
a(4) = 6 since r(s(4)) = r(10) = 6. r(10) = card({1, 2, 4, 5, 8, 10}) = 6.
a(5) = 5 since r(s(5)) = r(15) = 5. r(15) = card({1, 3, 5, 9, 15}) = 5.
a(6) = 11 since r(s(6)) = r(24) = 11. r(24) = card({1, 2, 3, 4, 6, 8, 9, 12, 16, 18, 24}) = 11, etc.
		

Crossrefs

Programs

  • Mathematica
    nn = 68; c[_] := False;
    rad[x_] := rad[x] = Times @@ FactorInteger[x][[All, 1]];
    f[x_] := Select[Range[x], Divisible[x, rad[#]] &];
    Array[Set[{a[#], c[#]}, {#, True}] &, 2]; s = a[1] + a[2];
    {1}~Join~Reap[Do[r = f[s]; k = SelectFirst[r, ! c[#] &];
      Sow[Length[r]]; c[k] = True;
      s += k, {i, 3, nn}] ][[-1, 1]]

A384000 Smallest number k with n distinct prime factors such that A010846(k) = A024718(n) (a tight lower bound), or -1 if such k does not exist.

Original entry on oeis.org

1, 2, 6, 1001, 268801, 3433936673, 2603508937756211
Offset: 0

Views

Author

Michael De Vlieger, May 19 2025

Keywords

Comments

These numbers k have the smallest A010846(k) for a number with n distinct prime factors.
a(7) <= 1398483454696343742813089 = 1049 * 2819 * 3319 * 3433 * 3457 * 3463 * 3467.
a(8) <= 32829974457045619959776094471833047127947.

Examples

			Table of a(n), n = 0..6, showing prime decomposition and cardinality of row a(n) of A162306, c(n) = A010846(a(n)) = A024718(n).
n               a(n)   c(n)    prime factors of a(n)        a(n)
----------------------------------------------------------------------
0                  1     1     -
1                  2     2     2                            A000040(1)
2                  6     5     2,   3                       A138109(1)
3               1001    15     7,  11,  13                  A383177(1)
4             268801    50    13,  23,  29,  31             A383178(2)
5         3433936673   176    41,  83,  97, 101, 103        A383179(209)
6   2603508937756211   638   163, 373, 439, 457, 461, 463
Tables of terms m in r(a(n)) = row a(n) of A162306, writing instead only exponents i of prime power factors p^i | m for  each p | a(n), written in order of the prime base:
For n = 2, i.e., squarefree semiprime k in A138109 (that achieves the lower bound), we have the following ordered exponent combinations in a rank-2 table:
  00  10  20
  01  11
Thus row 6 of A162306 has the following elements:
   1   2   4
   3   6
For n = 3, i.e., sphenic k in A383177 (that achieves the lower bound), we have the following ordered exponent combinations in a rank-3 table:
  000 100 200 300     001 101 201     002
  010 110 210         011 111
  020 120
Thus row 1001 of A162306 has the following elements:
    1   7  49 343      13   91 637    169
   11  77 539         141 1001
  121 857
		

Crossrefs

A301892 a(n) = A010846(A002182(n)).

Original entry on oeis.org

1, 2, 3, 5, 8, 11, 14, 15, 26, 36, 44, 49, 58, 76, 131, 156, 174, 206, 266, 308, 339, 388, 428, 460, 766, 550, 568, 979, 1124, 1238, 1411, 1548, 1659, 1754, 1983, 2048, 2160, 3689, 4211, 4617, 5245, 5731, 6135, 6482, 7308, 7539, 7949, 8477, 9198, 9681, 10306
Offset: 1

Views

Author

Michael De Vlieger, Mar 28 2018

Keywords

Comments

We define an "n-regular" number as 1 <= m <= n such that m | n^e with integer e >= 0. The divisor d is a special case of regular number m such that d | n^e with e = 0 or e = 1. Regular numbers m can exceed n; we are concerned only with regulars m <= n herein.
Since highly composite numbers represent those numbers that set records in the divisor counting function A000005, and since the divisor is a special case of regular number, this sequence applies the "regular counting function" A010846 to terms in A002182.
Only 13 HCNs less than 36 * 10^6 are also "highly regular", i.e., appear in A244052. The largest HCN that is also highly regular is 27720, the 25th HCN and the 47th highly regular number.
Only 2 and 6 set records for the ratio A010846(n)/A000005(n).
Conjectures:
Let "tier" t consist of all terms A002110(t) <= m < A002110(t + 1) in A244052, where all such m in tier t have A001221(m) = t. The intersection of A002182 and A244052 is finite, consisting of 13 terms: {1, 2, 4, 6, 12, 24, 60, 120, 180, 840, 1260, 1680, 27720}. All of these terms are also in A060735 and not in A288813, as the latter are squarefree and have "gaps" among prime divisors. This intersection has the following number of terms in the "tiers" 0 through 5 of A244052: {1, 2, 3, 3, 3, 1}. If we look at A060735 as a number triangle T(n,k) = k * A002110(n) with 1 <= k < prime(n + 1), the terms are:
{0, 1},
{{1,1}, {1,2}},
{{2,1}, {2,2}, {2,4}},
{{3,2}, {3,4}, {3,6}},
{{4,4}, {4,6}, {4,8}},
{5,12}.

Examples

			A002182(4) = 6. There are five numbers 1 <= m <= 6 such that m divides an integer power of 6: {1, 2, 3, 4, 6}. Thus, a(4) = 5.
		

Crossrefs

Programs

  • Mathematica
    With[{s = Array[DivisorSigma[0, #] &, 10^6]}, Map[With[{n = FirstPosition[s, #][[1]]}, Count[Range@ n, _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)]] &, Union@ FoldList[Max, s]]]
  • PARI
    a010846(n) = sum(k=1, n, if(gcd(n, k)-1, 0, moebius(k)*(n\k))) \\ after Benoit Cloitre in A010846
    r=0; for(x=1, oo, if(numdiv(x) > r, print1(a010846(x), ", "); r=numdiv(x))) \\ Felix Fröhlich, Mar 30 2018

A380483 a(n) = least k such that A010846(6*prime(k)) = n, or -1 if no solution exists.

Original entry on oeis.org

3, 4, -1, -1, 5, 6, -1, 7, 8, 9, -1, 10, -1, 12, 13, 14, 16, 17, 19, 21, 23, 24, 25, 29, 31, 32, 35, 38, 40, 44, 48, 54, 55, 62, 67, 69, 73, 77, 84, 93, 98, 106, 119, 124, 130, 136, 151, 165, 173, 184, 191, 211, 219, 232, 243, 270, 296, 310, 328, 343, 378, 399, 422
Offset: 18

Views

Author

Michael De Vlieger, Jul 08 2025

Keywords

Comments

Conjecture: n = 30 is the largest number for which a(n) = -1.

Examples

			Let f = A010846.
For k >= 3 and m = 6*prime(k), f(k) = [k=3] + f(6) + Sum_{j=0..floor(log_3 m)} floor(log_2 m/3^j), with Iverson brackets. (This, since for k = 3, prime(k)^2 < 2*3*prime(k) < prime(k)^3, but for k > 3, 2*3*prime(k) < prime(k)^2.)
a(18) = 3 since 2*3*prime(3) = 30 is the smallest sphenic number k, and f(30) = A363061(3) = 18. This is to say that row 30 of A162306 = {m : rad(m) | 30} = {1,2,3,4,5,6,8,9,10,12,15,16,18,20,24,25,27,30} has cardinality 18.
a(19) = 4 since 2*3*prime(4) = 42 and f(42) = 19. This is to say that row 30 of A162306 = {m : rad(m) | 42} = {1,2,3,4,6,7,8,9,12,14,16,18,21,24,27,28,32,36,42} has cardinality 19.
a(m) = -1, m = 20..21, since for k >= 3, there exists no solution to f(6*prime(k)) = m.
a(22) = 5 since 2*3*prime(5) = 66 and f(66) = 22.
a(23) = 6 since 2*3*prime(6) = 78 and f(78) = 23.
a(24) = = -1 since for k >= 3, there exists no solution to f(6*prime(k)) = 24.
a(25) = 7 since 2*3*prime(7) = 102 and f(102) = 25.
a(26) = 8 since 2*3*prime(8) = 114 and f(114) = 26.
a(27) = 9 since 2*3*prime(9) = 138 and f(138) = 27.
a(28) = = -1 since for k >= 3, there exists no solution to f(6*prime(k)) = 28.
a(29) = 10 since 2*3*prime(10) = 174 and f(174) = 29.
a(30) = = -1 since for k >= 3, there exists no solution to f(6*prime(k)) = 30.
Note: 2*3*prime(11) = 186; f(186) = f(174) = 29.
a(31) = 12 since 2*3*prime(12) = 222 and f(222) = 31, etc.
		

Crossrefs

Programs

  • Mathematica
    k = 3; Table[While[Set[t, Boole[k == 3] + 5 + Sum[Floor[1 + Log2[#/3^j]], {j, 0, Floor[Log[3, #]]}] &[6*Prime[k] ] ] < n, k++]; If[t == n, k, -1], {n, 18, 30}]

A383180 Irregular table T(n,k) = A010846(A019565(2^n + k)).

Original entry on oeis.org

1, 2, 2, 5, 2, 6, 5, 18, 2, 6, 5, 19, 5, 20, 16, 68, 2, 7, 6, 22, 5, 21, 18, 77, 5, 22, 17, 79, 16, 74, 60, 283, 2, 7, 6, 23, 5, 23, 18, 80, 5, 22, 18, 82, 16, 78, 62, 295, 5, 24, 19, 87, 16, 82, 64, 315, 15, 80, 62, 316, 55, 290, 226, 1161
Offset: 0

Views

Author

Michael De Vlieger, May 09 2025

Keywords

Examples

			Triangle begins:
  0: 1;
  1: 2;
  2: 2, 5;
  3: 2, 6, 5, 18;
  4: 2, 6, 5, 19, 5, 20, 16, 68;
  5: 2, 7, 6, 22, 5, 21, 18, 77, 5, 22, 17, 79, 16, 74, 60, 283;
   ...
		

Crossrefs

Programs

  • Mathematica
    (* Load the "theta" program at the Mathematica link in A369609, then: *)
    f[x_] := Times @@ Prime@ Position[Reverse@ IntegerDigits[x, 2], 1][[All, 1]]; Table[theta[f[2^n + k] ], {n, 0, 7}, {k, 0, 2^n - 1}]

Formula

T(0,0) = 1.
T(n,0) = 2.
T(n,2^(n-1)-1) = A363061(n).

A384004 a(n) = smallest k such that A010846(k) = n.

Original entry on oeis.org

1, 2, 4, 8, 6, 10, 22, 12, 44, 18, 24, 50, 98, 36, 48, 54, 224, 30, 42, 70, 108, 66, 78, 162, 102, 60, 138, 84, 174, 260, 132, 90, 126, 228, 354, 120, 234, 168, 350, 306, 150, 516, 408, 180, 252, 552, 696, 294, 240, 336, 612, 378, 270, 1416, 300, 702, 1332, 360
Offset: 1

Views

Author

Michael De Vlieger, Jun 10 2025

Keywords

Comments

For n > 2, a(n) is composite, since A010846(p) = 2 for prime p.
For n <= 3, a(n) = 2^n; for n > 3, a(n) < 2^n, and a(n) is in A024619.
Smallest k with omega(k) = i is A002110(i).
Conjecture: there are only 8 powerful terms (i.e., in A001694) in the sequence.

Examples

			Table of n, a(n) for n=1..10, showing row a(n) of A162306, replacing lpf(a(n)) with p, and A119288(a(n)) with q. Note: A010846(n) is the length of row n of A162306.
 n  a(n)  row n of A162306
----------------------------------------------------------
 1:   1   {1}
 2:   2   {1, p}
 3:   4   {1, p, p^2}
 4:   8   {1, p, p^2, p^3}
 5:   6   {1, p, q, p^2, p*q}
 6:  10   {1, p, p^2, q, p^3, p*q}
 7:  22   {1, p, p^2, p^3, q, p^4, p*q}
 8:  12   {1, p, q, p^2, p*q, p^3, q^2, p^2*q}
 9:  44   {1, p, p^2, p^3, q, p^4, p*q, p^5, p^2*q}
10:  18   {1, p, q, p^2, p*q, p^3, q^2, p^2*q, p^4, p*q^2}
		

Crossrefs

Programs

  • Mathematica
    (* First, load the theta program from the algorithms linked in A369609, then: *)
    nn = 2310; t[_] := 0; u = 1; Do[(If[t[#] == 0, t[#] = n]; If[# == u, While[t[u] != 0, u++]]) &[theta[n]], {n, nn}]; Array[t, u - 1]

A384960 a(n) = smallest sphenic number k such that A010846(k) = n.

Original entry on oeis.org

1001, 105, 231, 30, 42, 70, 110, 66, 78, 170, 102, 114, 138, 370, 174, 826, 222, 246, 258, 318, 354, 402, 438, 498, 534, 582, 654, 762, 786, 894, 978, 1038, 1158, 1338, 1506, 1542, 1758, 1986, 2082, 2202, 2334, 2598, 2922, 3126, 3462, 3918, 4098, 4398, 4614, 5262
Offset: 15

Views

Author

Michael De Vlieger, Jul 06 2025

Keywords

Comments

a(1) = A384000(3) = 1001; A010846(1001) = A024718(3) = 15; 1001 is the smallest number k with 3 distinct prime factors that has the smallest possible number of terms in row k of A162306, i.e., m <= k such that rad(m) | k.
For n > 30, 6 | a(n).

Examples

			Table of a(n) indicating prime factors and S, where S = {ceiling(log_p a(n))} for all primes p that divide a(n), in order of the magnitude of p.
                                Prime power factor
                                    1111223344455
 n  m=a(n) pi(facs(m))    S     23571379391713739
-------------------------------------------------
15   1001   4.5.6       4.3.3   ...111
16    105   2.3.4       5.3.3   .111
17    231   2.4.5       5.3.3   .1.11
18     30   1.2.3       5.4.3   111
19     42   1.2.4       6.4.2   11.1
20     70   1.3.4       7.3.3   1.11
21    110   1.3.5       7.3.2   1.1.1
22     66   1.2.5       7.4.2   11..1
23     78   1.2.6       7.4.2   11...1
24    170   1.3.7       8.4.2   1.1...1
25    102   1.2.7       7.5.2   11....1
26    114   1.2.8       7.5.2   11.....1
27    138   1.2.9       8.5.2   11......1
28    370   1.3.12      9.4.2   1.1........1
29    174   1.2.10      8.5.2   11.......1
30    826   1.4.17     10.4.2   1..1............1
31    222   1.2.12      8.5.2   11.........1
32    246   1.2.13      8.6.2   11..........1
33    258   1.2.14      9.6.2   11...........1
34    318   1.2.16      9.6.2   11.............1
		

Crossrefs

Programs

  • Mathematica
    (* See Mathematica code link for function definitions for kappaomega and theta *)
    s =  kappaomega[3, 6000]; t = Map[theta, s];
    Map[s[[FirstPosition[t, #][[1]] ]] &, Union[t]]

A275581 Numbers n such that A010846(n) >= n/2.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 10, 12, 18, 30
Offset: 1

Views

Author

Michael De Vlieger, Dec 25 2016

Keywords

Comments

Let 1 <= r <= n be a number "regular to" n, that is, a product of prime divisors p that also divide n.
This sequence includes numbers n such that there are at least as many regulars r among the range 1 <= m <= n as m that are nonregular to n.
Divisors d are a special case of regular r such that d also divides n in addition to all prime divisors p of d also dividing n.
With the exception of r = 1, all regular numbers r are part of the cototient of n.
For prime n, A010846(n) = A000005(n). Composite n > 4 can have nondivisor r. The only prime numbers in a(n) are {2,3}.
Regular r divides some power of n, but reliably divides n^r. When r is a divisor d of n, r | n.
Because the divisors of n are a subset of the regulars of n, it is interesting to compare A010846(n) with A000005(n) = tau(n):
A010846(n) = n for n = (1,2). Tau(n) = n for n = (1,2).
n/2 < A010846(n) < n for n = (3,4,6,10,12,18,30). n/2 < tau(n) < n for n = (3,4,6).
A010846(n) = n/2 for n = 8. Tau(n) = n/2 for n = (8,12).
Sequence A020490 includes numbers n such that tau(n) >= phi(n). The number 24 is in A020490 but not in a(n). This is because A010846(24) = 11 (cf. A162306(24) = {1,2,3,4,6,8,9,12,16,18,24}).
Numbers n such that tau(n) >= n/2: {1,2,3,4,6,8,12}.

Crossrefs

Programs

  • Mathematica
    Select[Range[10^3], 2 Count[Range@ #, k_ /; PowerMod[#, k, k] == 0] >= # &] (* or *)
    Select[Range[10^3], Function[n, 2 Total[MoebiusMu[#] Floor[n/#] &@ Select[Range@ n, CoprimeQ[#, n] &]] >= n]] (* or *)
    Select[Range[10^3], Function[n, 2 (1 + Count[Range@ n, m_ /; SubsetQ[FactorInteger[n][[All, 1]], FactorInteger[m][[All, 1]]]]) >= n]]
Showing 1-10 of 80 results. Next