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.

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