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

A244053 Let m = A244052(n) = n-th highly regular number; a(n) = number of numbers r <= m, all of whose prime divisors p also divide m.

Original entry on oeis.org

1, 2, 3, 5, 6, 8, 10, 11, 18, 19, 26, 28, 32, 36, 41, 44, 68, 77, 80, 96, 115, 131, 145, 156, 166, 174, 183, 192, 283, 295, 313, 322, 382, 395, 452, 463, 505, 519, 551, 567, 593, 629, 660, 691, 717, 743, 766, 1161, 1224, 1253, 1257, 1285, 1306, 1526
Offset: 1

Views

Author

Michael De Vlieger, Jun 18 2014

Keywords

Comments

Analogous to A002183.
Records transform of A010846. -Michael De Vlieger, Mar 08 2017

Examples

			a(5) = 6 since 6 is the fifth record value of A010846. The first record value is 1 set at position 1; the second is 2 set at position 2, the third is 3 set at position 4, the fourth is 5 set at position 6. Sequence A244052 records the positions of these record values.
		

Crossrefs

Programs

  • Mathematica
    a010846[n_] := Block[{pf, a}, a[x_] := First /@ FactorInteger@ x; pf = a@ n; If[n == 1, 1, 1 + Count[Range@ n, ?(SubsetQ[pf, a@ #] &)]]]; f[n] := Block[{t = {}, max = 0, x}, Do[If[(x = a010846@ i) > max, max = x; AppendTo[t, a010846[i]]], {i, n}]; t]; f@ 1000 (* Michael De Vlieger, Feb 10 2015 *)
    Union@ Rest@ FoldList[Max, 0, Array[Count[Range@ #, k_ /; PowerMod[#, Floor@ Log2@ #, k] == 0] &, 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]]]]; Union@ Rest@ FoldList[Max, 0, Array[f, 10^4]] (* Michael De Vlieger, Mar 08 2017, more efficient *)

A243103 Product of numbers m with 2 <= m <= n whose prime divisors all divide n.

Original entry on oeis.org

1, 2, 3, 8, 5, 144, 7, 64, 27, 3200, 11, 124416, 13, 6272, 2025, 1024, 17, 35831808, 19, 1024000, 3969, 247808, 23, 859963392, 125, 346112, 729, 2809856, 29, 261213880320000000, 31, 32768, 264627, 18939904, 30625, 26748301344768, 37, 23658496, 369603, 32768000000, 41
Offset: 1

Views

Author

Michael De Vlieger, Aug 19 2014

Keywords

Comments

This sequence is the product of n-regular numbers.
A number m is said to be "regular" to n or "n-regular" if all the prime factors p of m also divide n.
The divisor is a special case of a regular m such that m also divides n in addition to all of its prime factors p | n.
Analogous to A007955 (Product of divisors of n).
If n is 1 or prime, a(n) = n.
If n is a prime power, a(n) = A007955(n).
Note: b-file ends at n = 4619, because a(4620) has more than 1000 decimal digits.
Product of the numbers 1 <= k <= n such that (floor(n^k/k) - floor((n^k - 1)/k)) = 1. - Michael De Vlieger, May 26 2016

Examples

			a(12) = 124416 since 1 * 2 * 3 * 4 * 6 * 8 * 9 * 12 = 124416. These numbers are products of prime factors that are the distinct prime divisors of 12 = {2, 3}.
From _David A. Corneth_, Feb 09 2015: (Start)
Let p# be the product of primes up to p, A002110. Then
a(13#) ~= 8.3069582 * 10 ^ 4133
a(17#) ~= 1.3953000 * 10 ^ 22689
a(19#) ~= 3.8258936 * 10 ^ 117373
a(23#) ~= 6.7960327 * 10 ^ 594048
a(29#) ~= 1.3276817 * 10 ^ 2983168
a(31#) ~= 2.8152792 * 10 ^ 14493041
a(37#) ~= 1.9753840 * 10 ^ 69927040
Up to n = 11# already in the table.
(End)
		

Crossrefs

Cf. A162306 (irregular triangle of regular numbers of n), A010846 (number of regular numbers of n), A244974 (sum of regular numbers of n), A007955, A244052 (record transform of regular numbers of n).

Programs

  • Maple
    A:= proc(n) local F, S, s, j, p;
      F:= numtheory:-factorset(n);
      S:= {1};
      for p in F do
        S:= {seq(seq(s*p^j, j=0..floor(log[p](n/s))), s=S)}
      od;
      convert(S,`*`)
    end proc:
    seq(A(n), n=1..100); # Robert Israel, Feb 09 2015
  • Mathematica
    regularQ[m_Integer, n_Integer] := Module[{omega = First /@ FactorInteger @ m }, If[Length[Select[omega, Divisible[n, #] &]] == Length[omega], True, False]]; a20140819[n_Integer] := Times @@ Flatten[Position[Thread[regularQ[Range[1, n], n]], True]]; a20140819 /@ Range[41]
    regulars[n_] := Block[{f, a}, f[x_] := First /@ FactorInteger@ x; a = f[n];{1}~Join~Select[Range@ n, SubsetQ[a, f@ #] &]]; Array[Times @@ regulars@ # &, 12] (* Michael De Vlieger, Feb 09 2015 *)
    Table[Times @@ Select[Range@ n, (Floor[n^#/#] - Floor[(n^# - 1)/#]) == 1 &], {n, 41}] (* Michael De Vlieger, May 26 2016 *)
  • PARI
    lista(nn) = {vf = vector(nn, n, Set(factor(n)[,1])); vector(nn, n, prod(i=1, n, if (setintersect(vf[i], vf[n]) == vf[i], i, 1)));} \\ Michel Marcus, Aug 23 2014
    
  • PARI
    for(n=1, 100, print1(prod(k=1, n, k^(floor(n^k/k) - floor((n^k - 1)/k))),", ")) \\ Indranil Ghosh, Mar 22 2017
    
  • Python
    from sympy import primefactors
    def A243103(n):
        y, pf = 1, set(primefactors(n))
        for m in range(2,n+1):
            if set(primefactors(m)) <= pf:
                y *= m
        return y # Chai Wah Wu, Aug 28 2014
    
  • Scheme
    ;; A naive implementation, code for A123275bi given under A123275:
    (define (A243103 n) (let loop ((k n) (m 1)) (cond ((= 1 k) m) ((= 1 (A123275bi n k)) (loop (- k 1) (* m k))) (else (loop (- k 1) m)))))
    ;; Antti Karttunen, Mar 22 2017

Formula

a(n) = product of terms of n-th row of irregular triangle A162306(n,k).
a(n) = Product_{k=1..n} k^( floor(n^k/k)-floor((n^k -1)/k) ). - Anthony Browne, Jul 06 2016
From Antti Karttunen, Mar 22 2017: (Start)
a(n) = Product_{k=2..n, A123275(n,k)=1} k.
For n >= 1, A046523(a(n)) = A283990(n).
(End)

A299990 a(n) = A243822(n) - A000005(n).

Original entry on oeis.org

-1, -2, -2, -3, -2, -3, -2, -4, -3, -2, -2, -4, -2, -2, -3, -5, -2, -2, -2, -4, -3, -1, -2, -5, -3, -1, -4, -4, -2, 2, -2, -6, -2, 0, -3, -4, -2, 0, -2, -5, -2, 3, -2, -3, -4, 0, -2, -5, -3, 0, -2, -3, -2, 0, -3, -5, -2, 0, -2, 2, -2, 0, -4, -7, -3, 6, -2, -2
Offset: 1

Views

Author

Michael De Vlieger, Feb 25 2018

Keywords

Comments

Since A010846(n) = A000005(n) + A243822(n), this sequence examines the balance of the two components among "regular" numbers.
Value of a(n) is generally less frequently negative as n increases.
a(1) = -1.
For primes p, a(p) = -2 since 1 | p and the cototient is restricted to the divisor p.
For perfect prime powers p^e, a(p^e) = -(e + 1), since all m < p^e in the cototient of p^e that do not have a prime factor q coprime to p^e are powers p^k with 1 < p^k <= p^e; all such p^k divide p^e.
Generally for n with A001221(n) = 1, a(n) = -1 * A000005(n), since the cototient is restricted to divisors, and in the case of p^e > 4, divisors and numbers in A272619(p^e) not counted by A010846(p^e).
For m >= 3, a(A002110(m)) is positive.
For m >= 9, a(A244052(m)) is positive.

Examples

			a(6) = -3 since 6 has 4 divisors, and 4 | 6^2; A243822(6) = 1 and A000005(6) = 4; 1 - 4 = -3. Alternatively, A010846(6) = 5; 5 - 2*4 = -3.
a(30) = 2 since 30 has 8 divisors and the numbers {4, 8, 9, 12, 16, 18, 20, 24, 25, 27} divide 30^e with e > 1; A243822(30) = 10 and A000005(30) = 8; 10 - 8 = 2. Alternatively, A010846(30) = 18; 18 - 2*8 = 2.
Some values of a(n) and related sequences:
   n  a(n) A010846(n) A243822(n) A000005(n) A272618(n)
  ----------------------------------------------------
   1   -1          1          0          1  0
   2   -2          2          0          2  0
   3   -2          2          0          2  0
   4   -3          3          0          3  0
   5   -2          2          0          2  0
   6   -3          5          1          4  {4}
   7   -2          2          0          2  0
   8   -4          4          0          4  0
   9   -3          3          0          3  0
  10   -2          6          2          4  {4,8}
  11   -2          2          0          2  0
  12   -4          8          2          6  {8,9}
  ...
  30    2         18         10          8  {4,8,9,12,16,18,20,24,25,27}
  ...
  34    0          8          4          4  {4,8,16,32}
  ...
		

Crossrefs

Programs

  • Mathematica
    Table[Count[Range[n], _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)] - 2 DivisorSigma[0, n], {n, 68}]

Formula

a(n) = A010846(n) - 2*A000005(n).

A299991 Numbers n for which A243822(n) > A000005(n).

Original entry on oeis.org

30, 42, 60, 66, 70, 74, 78, 82, 84, 86, 90, 94, 98, 102, 106, 110, 114, 118, 120, 122, 126, 130, 132, 134, 138, 140, 142, 146, 150, 154, 156, 158, 162, 165, 166, 168, 170, 174, 178, 180, 182, 186, 190, 194, 195, 198, 202, 204, 206, 210, 214, 218, 220, 222, 226
Offset: 1

Views

Author

Michael De Vlieger, Feb 25 2018

Keywords

Comments

Composite numbers m have nondivisors k in the cototient such that k | n^e with e > 1. These k appear in row n of A272618 and are enumerated by A243822(n). These nondivisors k are a kind of "regular" number along with divisors d of n; both are listed in row n of A162306 and are together enumerated by A010846(n).
This sequence lists numbers that have more nondivisors k in the cototient of n than divisors d.
This sequence contains all n for which A299990(n) is positive.
The smallest odd term is 165.
For m >= 3, A002110(m) is in a(n).
For m >= 9, A244052(m) is in a(n).

Examples

			30 is the first term since it is the smallest number for which A243822(n) > A000005(n), alternatively, for which A010846(n) > 2*A000005(n).
		

Crossrefs

Programs

  • Mathematica
    Select[Range@ 226, Function[n, Count[Range[n], _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)] > 2 DivisorSigma[0, n]]]

A288813 Irregular triangle read by rows: T(m, k) is the list of squarefree numbers A002110(m) < t < 2*A002110(m) such that A001221(t) = m.

Original entry on oeis.org

3, 10, 42, 330, 390, 2730, 3570, 3990, 4290, 39270, 43890, 46410, 51870, 53130, 570570, 690690, 746130, 870870, 881790, 903210, 930930, 1009470, 11741730, 13123110, 14804790, 15825810, 16546530, 17160990, 17687670, 18888870, 281291010, 300690390, 340510170
Offset: 1

Views

Author

Michael De Vlieger, Jun 24 2017

Keywords

Comments

a(n) = terms t of row m of A288784 such that A002110(m) < t < 2*A002110(m).
The only odd term is 3; the only other term not ending in 10, 30, 70, or 90 in decimal is 42.
All terms t in row m have A001221(t) = m and at least one prime q coprime to t such that q < A006530(t).
Consider "tier" m and primorial p_m# = A002110(m), let "distension" i = pi(A006530(T(m, k))) - m and let "depth" j = m - pi(A053669(T(m, k))) + 1. Distension is the difference in the index of gpf(T(m, k)) and pi(m), while depth is the difference between the index of the least prime totative of T(m, k) and pi(m) + 1. We can calculate the maximum distension i given m and j via i_max = A020900(m - j + 1) - m - j + 1. This enables us to use permutations of 0 and 1 values in the notation A054841 and produce a(n) with some efficiency.
The most efficient method of generating a(n) is via f(x) = A287352(x), i.e., subtracting 1 from all values in row x of A287352. We use a pointer variable to direct increment on f(p_m#) = a constant array of m 1's, until we have exhausted producing terms p_m# < t < 2*p_m#. This enables the generation of T(m, k) for 1 <= m <= 100.

Examples

			Triangle begins:
n     a(n)
1:       3
2:      10
3:      42
4:     330     390
5:    2730    3570    3990    4290
6:   39270   43890   46410   51870   53130
7:  570570  690690  746130  870870  881790  903210  930930  1009470
       ...
		

Crossrefs

Programs

  • Mathematica
    Table[Function[P, Select[Range[P + 1, 2 P - 1], And[SquareFreeQ@ #, PrimeOmega@ # == n] &]]@ Product[Prime@ i, {i, n}], {n, 7}] // Flatten (* Michael De Vlieger, Jun 24 2017 *)
    f[n_] := Block[{P = Product[Prime@ i, {i, n}], lim, k = 1, c, w = ConstantArray[1, n]}, lim = 2 P; Sort@ Reap[Do[w = If[k == 1, MapAt[# + 1 &, w, -k], Join[Drop[MapAt[# + 1 &, w, -k], -k + 1], ConstantArray[1, k - 1]]]; c = Times @@ Map[If[# == 0, 1, Prime@ #] &, Accumulate@ w]; If[c < lim, Sow[c]; k = 1, If[k == n, Break[], k++]], {i, Infinity}] ][[-1, 1]] ]; Array[f, 9] // Flatten (* Michael De Vlieger, Jun 28 2017, faster *)
  • PARI
    primo(n) = prod(i=1, n, prime(i));
    row(n) = my(vrow = []); for (j=primo(n)+1, 2*primo(n)-1, if (issquarefree(j) && (omega(j)==n), vrow = concat(vrow, j))); vrow;
    tabf(nn) = for (n=1, nn, print(row(n))); \\ Michel Marcus, Jun 29 2017

A300156 Indices of records in A299990.

Original entry on oeis.org

1, 30, 42, 66, 78, 90, 102, 114, 138, 150, 210, 330, 390, 510, 570, 630, 870, 990, 1050, 1470, 1890, 2100, 2310, 2730, 3570, 3990, 4620, 5460, 6510, 6930, 8190, 9240, 10710, 11550, 13650, 16170, 19110, 20790, 23100, 24570, 25410, 30030, 39270, 43890, 46410
Offset: 1

Views

Author

Michael De Vlieger, Feb 26 2018

Keywords

Comments

A010846(n) = A000005(n) + A243822(n).
Successive terms in this sequence represent increasing differences A243822(n) - A000005(n).
A000079 = records in -1 * A299990, since A243822(p^e)=0 for e>=0, n = 2^k sets records in A000005(n). The corresponding records are in A000027.

Examples

			The first term is 1: A299990(1) = -1. The first term that exceeds -1 in A299990 corresponds to n = 30; A299990(30) = 2, thus 30 is the second term.
		

Crossrefs

Programs

  • Mathematica
    With[{s = Table[Count[Range[n], _?(PowerMod[n, Floor@ Log2@ n, #] == 0 &)] - 2 DivisorSigma[0, n], {n, 10^3}]}, Map[FirstPosition[s, #][[1]] &, Union@ FoldList[Max, s]]]

A292867 Indices of records in A243823.

Original entry on oeis.org

1, 8, 14, 16, 20, 22, 26, 28, 32, 38, 40, 44, 46, 48, 50, 52, 54, 56, 58, 62, 64, 68, 72, 78, 80, 86, 88, 92, 94, 96, 100, 108, 114, 122, 124, 126, 130, 132, 138, 144, 156, 160, 162, 174, 186, 192, 204, 216, 222, 228, 234, 240, 246, 252, 258, 264, 270, 276, 282
Offset: 1

Views

Author

Michael De Vlieger, Oct 02 2017

Keywords

Comments

Except for A292867(1) = 1, all terms are even.
Some conjectures:
1. The only prime powers p^e in this sequence are {8, 16, 32, 64}.
2. Squarefree terms m appear throughout. (There are 261 squarefree values among the first 1261 terms.)
3. Terms that set records for omega(m) are 1, followed by 2^e, with 3 <= e <= 6, then 2^e * 3 with 6 <= e <= 8, then 2^7 * A002110(k) with k >= 1.
4. Primorials A002110(n) for n >= 6 appear in this sequence. The first primorials in m are terms 6 through 8 of A002110 (i.e., 30030, 510510, 9699690) at n = 419, 774, 1258, respectively.
5. Outside of a(n) with 2 <= n <= 21 and n = {29, 30}, all terms of A244052 are also in this sequence. This observation applies to the smallest 104 terms of A244052.
6. For very large n, all terms are also in A244052. For small n, few terms of A244052 appear and are separated by many other numbers. Since numbers m in A244052 are products of k primes, many of which are the smallest primes, phi is minimized and A010846(m) becomes infinitesimal in comparison to m. Therefore A243823(m) is tantamount to the cototient of m. The size of n required to observe this agreement between this sequence and A244052 is unknown.

Examples

			8 is in the sequence since it is the first number n such that A243823(n) > 0. 14 appears immediately after 8 since A243823(14) = 3, and 3 is greater than the values that precede it.
Table of indices a(n) of records b(n) in A243823 = n - phi(n) - A010846(n) + 1:
.
   n  a(n) b(n) phi(a(n)) A010846(a(n))
  -------------------------------------
   1    1    0         1         1
   2    8    1         4         4
   3   14    3         6         6
   4   16    4         8         5
   5   20    5         8         8
   6   22    6        10         7
   7   26    8        12         7
   8   28    9        12         8
   9   32   11        16         6
  10   38   13        18         8
  11   40   14        16        11
  12   44   16        20         9
  13   46   17        22         8
  14   48   18        16        15
  15   50   19        20        12
  16   52   20        24         9
  17   54   21        18        16
  18   56   22        24        11
  19   58   23        28         8
  20   62   25        30         8
  ...
		

Crossrefs

Programs

A293555 Indices of records in A243822.

Original entry on oeis.org

1, 6, 10, 18, 30, 42, 60, 78, 84, 90, 126, 150, 210, 330, 390, 420, 630, 840, 990, 1050, 1470, 1890, 2100, 2310, 2730, 3570, 3990, 4620, 5460, 6930, 8190, 9240, 10920, 11550, 13650, 13860, 16170, 19110, 20790, 23100, 25410, 30030, 39270, 43890, 46410, 51870
Offset: 1

Views

Author

Michael De Vlieger, Oct 22 2017

Keywords

Comments

From Michael De Vlieger, Nov 17 2017: (Start)
Terms in a(n) appear in A244052 except {78, 126, 990, 19110, 6276270, ...}.
Primorials A002110(t) seem to divide this sequence into "tiers" thus: all terms A002110(t) <= m < A002110(t + 1), wherein A001221(m) = t as seen in A244052.
Terms in A244052 appear in a(n) except {2, 4, 12, 24, 120, 180, 1260, 1680, 18480, 27720, 360360, ...}. These numbers seem to have significantly more divisors than terms that are slightly greater or lesser in a(n).
Conjecture: all terms of a(n) with n > 92 also appear in A244052, and all terms in A244052 greater than a(92) = 6276270 appear in a(n).
(End)

Examples

			From _Michael De Vlieger_, Nov 17 2017: (Start)
Consider A243822(n), a function that counts numbers k < n such that k | n^e with e >= 2. The numbers k themselves appear in A272618(n).
a(1) = 1 since the number 1 has 0 such k. Primes p also have 0 such k, since 1 | p and all other numbers k < p are coprime to p. Prime powers p^e have 0 such k since any number k | n^e divides n^1.
a(2) = 6 since it is the smallest number to have 1 such k (i.e., 4 | 6^2). The numbers 7, 8, and 9 are prime powers having 0 such k.
a(3) = 10 since it has 2 such k (i.e., 4 | 10^2, 8 | 10^3), etc.
(End)
		

Crossrefs

Programs

A288784 Irregular triangle read by rows: T(n,m) is the list of numbers k*A002110(n) <= k*t < (k + 1)*A002110(n) such that A001222(k*t) = n, with 1 <= k < prime(n + 1).

Original entry on oeis.org

1, 2, 3, 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, 4290, 4620, 5460, 6930, 8190, 9240, 10920, 11550, 13650, 13860, 16170, 18480, 20790, 23100, 25410, 27720
Offset: 0

Views

Author

Michael De Vlieger, Jun 15 2017

Keywords

Comments

A060735 and A002110 are subsets.
This sequence is a necessary but insufficient condition for A244052. Terms that are in A060735 and A002110 are also in A244052. The first terms of this sequence that are not in A244052 are {3, 4290, 881790, 903210, 1009470, 17160990, 363993630, 380570190, 406816410, 434444010, ...}.
Primorial p_n# = A002110(n) is the smallest squarefree number with n prime factors. Consider the list of squarefree numbers t with n prime factors greater than and including A002110(n) but less than 2*A002110(n). Extend the list to include products k*t of this list with 1 <= k < prime(n+1) such that k*t < (k+1)*p_n#. This list contains squarefree numbers k*t with n distinct primes and presumes that the number (k+1)*p_n# serves as a "limit" beyond which k*t > (k+1)p_n# are not in the sequence.

Examples

			Triangle begins:
    n   T(n,m)
    0:   1;
    1:   2,  3,  4;
    2:   6, 10, 12, 18, 24;
    3:  30, 42, 60, 84, 90, 120, 150, 180;
       ...
		

Crossrefs

Programs

  • Mathematica
    Table[Function[P, Function[s, Flatten@ Map[Function[k, Select[k s, # < (k + 1) P &]], Range[1, Prime[n + 1] - 1]]]@ Select[Range[P, 2 P - 1], And[SquareFreeQ@ #, PrimeOmega@ # == n] &]]@ Product[Prime@ i, {i, n}], {n, 0, 5}] (* Michael De Vlieger, Jun 15 2017 *)

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
Showing 1-10 of 13 results. Next