OFFSET
1,2
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
EXAMPLE
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.
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 *)
PROG
(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
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Bernard Schott, Sep 23 2022
EXTENSIONS
More terms from Michel Marcus, Sep 23 2022
STATUS
approved