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.

A335037 a(n) is the number of divisors of n that are themselves divisible by the product of their digits.

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 4, 3, 3, 2, 6, 1, 3, 4, 4, 1, 5, 1, 4, 3, 3, 1, 8, 2, 2, 3, 4, 1, 6, 1, 4, 3, 2, 3, 8, 1, 2, 2, 5, 1, 5, 1, 4, 5, 2, 1, 8, 2, 3, 2, 3, 1, 5, 3, 5, 2, 2, 1, 8, 1, 2, 4, 4, 2, 5, 1, 3, 2, 4, 1, 10, 1, 2, 4, 3, 3, 4, 1, 5, 3, 2, 1, 7, 2, 2, 2, 5
Offset: 1

Views

Author

Bernard Schott, Jun 03 2020

Keywords

Comments

Inspired by A332268.
A number that is divisible by the product of its digits is called Zuckerman number (A007602); e.g., 24 is a Zuckerman number because it is divisible by 2*4=8 (see links).
a(n) = 1 iff n = 1 or n is prime not repunit >= 13.
a(n) = 2 iff n is prime = 2, 3, 5, 7 or a prime repunit.
Numbers 1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 12, 15, 24, 111111111111111111111 (repunit with 19 times 1's) have all divisors Zuckerman numbers. The sequence of numbers with all Zuckerman divisors is infinite iff there are infinitely many repunit primes (see A004023).

Examples

			For n = 4, the divisors are 1, 2, 4 and they are all Zuckerman numbers, so a(4) = 3.
For n = 14, the divisors are 1, 2, 7 and 14. Only 1, 2 and 7 are Zuckerman numbers, so a(14) = 3.
		

Crossrefs

Similar with: A001227 (odd divisors), A087990 (palindromic divisors), A087991 (non-palindromic divisors), A242627 (divisors < 10), A332268 (Niven divisors).

Programs

  • Mathematica
    zuckQ[n_] := (prodig = Times @@ IntegerDigits[n]) > 0&& Divisible[n, prodig]; a[n_] := Count[Divisors[n], ?(zuckQ[#] &)]; Array[a, 100] (* _Amiram Eldar, Jun 03 2020 *)
  • PARI
    iszu(n) = my(p=vecprod(digits(n))); p && !(n % p);
    a(n) = sumdiv(n, d, iszu(d)); \\ Michel Marcus, Jun 03 2020

Formula

Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{n>=1} 1/A007602(n) = 3.26046... . - Amiram Eldar, Jan 01 2024

A337941 Numbers whose divisors are all Zuckerman numbers (A007602).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 11, 12, 15, 24, 1111111111111111111, 11111111111111111111111
Offset: 1

Views

Author

Bernard Schott, Oct 01 2020

Keywords

Comments

Inspired by A337741.
Zuckerman numbers are numbers that are divisible by the product of their digits (see link).
The next term is the repunit prime R_317 which is too large to include in the data.
Primes in this sequence are 2, 3, 5, 7 and all the repunit primes (see A004023).
This sequence is infinite if and only if there are infinitely many repunit primes.

Examples

			6 is a term since all the divisors of 6, i.e., 1, 2, 3 and 6, are Zuckerman numbers.
		

Crossrefs

Subsequence of A007602.
Similar sequences: A062687, A190217, A308851, A329419, A337741.
Cf. A004022 (subsequence of prime repunits).

Programs

  • Mathematica
    zuckQ[n_] := (prod = Times @@ IntegerDigits[n]) > 0 && Divisible[n, prod]; Select[Range[24], AllTrue[Divisors[#], zuckQ] &] (* Amiram Eldar, Oct 01 2020 *)
  • PARI
    isok(m) = {fordiv(m, d, my(p=vecprod(digits(d))); if (!p || (d % p), return (0))); return (1);} \\ Michel Marcus, Oct 05 2020

A340796 a(n) is the smallest number with exactly n divisors that are Brazilian.

Original entry on oeis.org

1, 7, 14, 24, 40, 48, 60, 84, 140, 144, 120, 168, 252, 700, 240, 336, 560, 360, 420, 672, 1120, 2304, 960, 720, 1008, 1080, 840, 2184, 1800, 1260, 2016, 5376, 8960, 2160, 1680, 2880, 4032, 3600, 7056, 19600, 3960, 2520, 3360, 6480, 9072, 9900, 6300, 11520, 16128
Offset: 0

Views

Author

Bernard Schott, Jan 21 2021

Keywords

Comments

Primes can be partitioned into Brazilian primes and non-Brazilian primes. If two distinct primes each larger than 11 are in the same category then the larger one has a multiplicity that is smaller than or equal to that of the smaller prime. - David A. Corneth, Jan 24 2021

Examples

			Of the eight divisors of 24, three are Brazilian numbers: 8, 12 and 24, and there is no smaller number with three Brazilian divisors, hence a(3) = 24.
		

Crossrefs

Similar with: A087997 (palindromes), A333456 (Niven), A335038 (Zuckerman).

Programs

  • Mathematica
    brazQ[n_] := Module[{b = 2, found = False}, While[b < n - 1 && Length[Union[IntegerDigits[n, b]]] > 1, b++]; b < n - 1]; d[n_] := DivisorSum[n, 1 &, brazQ[#] &]; m = 30; s = Table[0, {m}]; c = 0; n = 1; While[c < m, i = d[n]; If[i < m && s[[i + 1]] == 0, c++; s[[i + 1]] = n]; n++]; s (* Amiram Eldar, Jan 21 2021 *)
  • PARI
    isokb(n) = for(b=2, n-2, d=digits(n, b); if(vecmin(d)==vecmax(d), return(1))); \\ A125134
    isok(k, n) = sumdiv(k, d, isokb(d)) == n;
    a(n) = my(k=1); while (!isok(k, n), k++); k; \\ Michel Marcus, Jan 23 2021

Extensions

More terms from Amiram Eldar, Jan 21 2021

A357300 a(n) is the smallest number m with exactly n divisors whose first digit equals the first digit of m.

Original entry on oeis.org

1, 10, 100, 108, 120, 180, 1040, 1020, 1170, 1008, 1260, 1680, 10010, 10530, 10200, 10260, 10560, 10800, 11340, 10920, 12600, 10080, 15840, 18480, 15120, 102060, 104400, 101640, 100320, 102600, 100980, 117600, 114660, 107100, 174240, 113400, 105840, 100800, 120120, 143640
Offset: 1

Views

Author

Bernard Schott, Sep 23 2022

Keywords

Comments

a(m) <= a(551) = 18681062400 for m < 555. All terms with values up to 2*10^10 start with 1. Do there exist a(n) starting with any other digit? - Charles R Greathouse IV, Sep 25 2022

Examples

			Of the twelve divisors of 108, four have their first digit equals to the first digit of 108: 1, 12, 18 and 108, and there is no such smaller number, hence a(4) = 108.
		

Crossrefs

Cf. A335491 (with last digit), A206287, A355592, A357299.
Similar, but with: A333456 (Niven numbers), A335038 (Zuckerman numbers).

Programs

  • Mathematica
    f[n_] := IntegerDigits[n][[1]]; s[n_] := Module[{fn = f[n]}, DivisorSum[n, 1 &, f[#] == fn &]]; seq[len_, nmax_] := Module[{v = Table[0, {len}], c = 0, n = 1, i}, While[c < len && n < nmax, i = s[n]; If[i <= len && v[[i]] == 0, c++; v[[i]] = n]; n++]; v]; seq[40, 10^6] (* Amiram Eldar, Sep 23 2022 *)
  • PARI
    f(n) = my(fd=digits(n)[1]); sumdiv(n, d, digits(d)[1] == fd); \\ A357299
    a(n) = my(k=1); while (f(k)!=n, k++); k; \\ Michel Marcus, Sep 23 2022
    
  • PARI
    v=vector(1000); v[1]=r=1; forfactored(n=2, 10^11, t=a(n[1],n[2],r); if(t>r && v[t]==0, v[t]=n[1]; print(t" "n[1]" = "n[2]); while(v[r],r++); r--)) \\ Charles R Greathouse IV, Sep 25 2022

Extensions

More terms from Michel Marcus, Sep 23 2022

A335491 a(n) is the smallest number m with exactly n divisors whose last digit equals the last digit of m.

Original entry on oeis.org

1, 11, 40, 60, 160, 120, 640, 240, 360, 480, 8064, 600, 18144, 1920, 1440, 1200, 72576, 1800, 52416, 2400, 5760, 30720, 183456, 3600, 12960, 122880, 9000, 9600, 602784, 7200, 445536, 8400, 92160, 798336, 51840, 12600, 2159136, 576576, 368640, 16800, 2935296, 28800
Offset: 1

Views

Author

Bernard Schott, Jun 11 2020

Keywords

Comments

a(n) exists for any n >= 1. Indeed, the number 5*2^n (see A020714), n >= 1, has exactly n divisors (5*2^1, 5*2^2, ..., 5*2^n), with the last digit 0. - Marius A. Burtea, Jun 12 2020
It's always the case when a(n) ends in 0 then a(n) = 10 * A005179(n). Proof: Let v be the list of divisors of a(n) that end in 0. We then have |v| = n and lcm(v) = a(n) as a(n) is in v and all other terms in v divide a(n). We then have lcm(v)/10 = a(n)/10 where a(n)/10 has exactly n divisors. The least positive integer that has exactly n divisors is A005179(n). - Bernard Schott and David A. Corneth, Jun 12 2020
For some p_i^e_i||a(n) and p_j^e_j||a(n) where p_i and p_j are primes such that p_i < p_j, 10 | (p_j - p_i) and t^k||u denotes t^k|u but t^(k + 1) doesn't divide u i.e. gcd(t, u/t^k) = 1 denotes then e_i >= e_j. For example k * 11^2 * 31^3 where gcd(k, 11*31) = 1 can't be a term as the multiplicity of 11 is less than the multiplicity of 31. - David A. Corneth, Jun 12 2020

Examples

			Of the twelve divisors of 60, four have their last digit equals to the last digit of 60: 10, 20, 30 and 60, and there is no smaller number k with four divisors whose last digit equals the last digit of k, hence a(4) = 60.
		

Crossrefs

Similar with: A333456 (Niven numbers), A335038 (Zuckerman numbers).

Programs

  • Magma
    a:=[]; for n in [1..30] do k:=1; while #[d:d in Divisors(k)|k mod 10 eq d mod 10] ne n do k:=k+1; end while; Append(~a,k); end for; a; // Marius A. Burtea, Jun 12 2020
  • Mathematica
    d[n_] := DivisorSum[n, 1 &, Mod[# - n, 10] == 0 &]; mx = 20; c = 0; n = 1; s = Table[0, {mx}]; While[c < mx, i = d[n]; If[i <= mx && s[[i]] == 0, c++; s[[i]] = n]; n++]; s (* Amiram Eldar, Jun 12 2020 *)
  • PARI
    f(n) = my(u=n%10); sumdiv(n, d, (d%10) == u);
    a(n) = my(k=1); while(f(k) != n, k++); k; \\ Michel Marcus, Jun 12 2020
    

Extensions

Corrected and extended by Marius A. Burtea, Jun 12 2020

A340638 Integers whose number of divisors that are Zuckerman numbers sets a new record.

Original entry on oeis.org

1, 2, 4, 6, 12, 24, 72, 144, 360, 432, 1080, 2016, 2160, 6048, 8064, 15120, 24192, 48384, 88704, 120960, 241920, 266112, 532224, 1064448, 1862784, 2661120, 3725568, 5322240, 7451136, 10450944, 19160064, 20901888, 28740096, 38320128, 57480192, 99283968, 114960384
Offset: 1

Views

Author

Bernard Schott, Jan 14 2021

Keywords

Comments

A Zuckerman number is a number that is divisible by the product of its digits (A007602).
The terms in this sequence are not necessarily Zuckerman numbers. For example a(7) = 72 has product of digits = 14 and 72/14 = 36/7 = 5.142...
The first seven terms are the first seven terms of A087997, then A087997(8) = 66 while a(8) = 144.

Examples

			The 8 divisors of 24 are all Zuckerman numbers, and also, 24 is the smallest integer that has at least 8 divisors that are Zuckerman numbers, hence 24 is a term.
		

Crossrefs

Subsequence of A335038.
Similar for palindromes (A093036), repdigits (A340548), repunits (A340549), Niven numbers (A340637).

Programs

  • Mathematica
    zuckQ[n_] := (prod = Times @@ IntegerDigits[n]) > 0 && Divisible[n, prod]; s[n_] := DivisorSum[n, 1 &, zuckQ[#] &]; smax = 0; seq = {}; Do[s1 = s[n]; If[s1 > smax, smax = s1; AppendTo[seq, n]], {n, 1, 10^5}]; seq (* Amiram Eldar, Jan 14 2021 *)
  • PARI
    isokz(n) = iferr(!(n % vecprod(digits(n))), E, 0); \\ A007602
    lista(nn) = {my(m=0); for (n=1, nn, my(x = sumdiv(n, d, isokz(d));); if (x > m, m = x; print1(n, ", ")););} \\ Michel Marcus, Jan 15 2021

Extensions

More terms from David A. Corneth and Amiram Eldar, Jan 15 2021
Showing 1-6 of 6 results.