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-2 of 2 results.

A327750 Numbers without zero digits such that after adding the product of its digits to it, a number with the same product of digits is obtained.

Original entry on oeis.org

28, 128, 214, 239, 266, 318, 326, 364, 494, 497, 563, 598, 613, 637, 695, 819, 1128, 1214, 1239, 1266, 1318, 1326, 1364, 1494, 1497, 1563, 1598, 1613, 1637, 1695, 1819, 2114, 2139, 2168, 2285, 2313, 2356, 2369, 2419, 2594, 2639, 2791, 3118, 3126, 3148, 3213, 3235
Offset: 1

Views

Author

Bernard Schott, Sep 24 2019

Keywords

Comments

The idea of this sequence comes from a problem in the annual Moscow Mathematical Olympiad (MMO) in 2003: Level A, problem 2. The problem only asks to find a ten-digit number that has the property of the name.
When an integer k belongs to this sequence, the integer 111..11//k obtained by concatenation // of 111..11 and k is also a term; hence, there are primitive terms as 28, 214, 239, 266, 318, 326, ... (A340908).
A subset of it is formed by the numbers 239, 326, 364, 497, 563, 598, 613, 637, 695, 819, 1239, 1326, 1364, 1497, 1563, 1598, 1613, 1637, 1695, 1819, 2139, 2313, 2356, 2369, ... for which the number obtained after adding the product of the digits has exactly the same digits (they are obtained by permuting the digits of the initial number). So, 239 + 2*3*9 = 239 + 54 = 293, 326 + 3*2*6 = 326 + 36 = 362, 3235 + 3*2*3*5 = 3235 + 90 = 3325, 23286 + 2*3*2*8*6 = 23286 + 576 = 23862. - Marius A. Burtea, Sep 24 2019
This subset is A247888. - Bernard Schott, Jul 22 2020

Examples

			28 + 2*8 = 44 and 2*8 = 4*4 hence 28 is a term.
326 + 3*2*6 = 362 and 3*2*6 = 3*6*2 hence 326 is another term.
		

Crossrefs

Subsequences: A247888, A340907, A340908 (primitives).

Programs

  • Magma
    [k:k in [1..3500]| not 0 in Intseq(k) and &*Intseq(k) eq &*(Intseq(k+&*Intseq(k)))]; // Marius A. Burtea, Sep 24 2019
    
  • Mathematica
    pd[n_] := Times @@ IntegerDigits[n]; aQ[n_] := (p = pd[n]) > 0 && pd[n + p] == p; Select[Range[5000], aQ] (* Amiram Eldar, Sep 24 2019 *)
  • PARI
    isok(n) = my(d = digits(n), p); vecmin(d) && (p=vecprod(d)) && (vecprod(digits(n+p)) == p); \\ Michel Marcus, Sep 24 2019
    
  • Python
    def test(n):
        m, p = n, 1
        while m > 0:
            m, p = m//10, p*(m%10)
        if p == 0:
            return 0
        m, q = n+p, 1
        while m > 0:
            m, q = m//10, q*(m%10)
        return p == q
    n, a = 0, 0
    while n < 100:
        a = a+1
        if test(a):
            n = n+1
        print(n,a) # A.H.M. Smeets, Sep 25 2019

Extensions

More terms from Amiram Eldar, Sep 24 2019

A340907 Numbers m without zero digits such that pod(q) = pod(k) = pod(m) where q = k + pod(k) and k = m + pod(m) where pod is the product of digits, A007954.

Original entry on oeis.org

262713, 267338, 283628, 342713, 351678, 432713, 451676, 516469, 516657, 516675, 622713, 634838, 651674, 716655, 728364, 851673, 857297, 916465, 1262713, 1267338, 1283628, 1342713, 1351678, 1432713, 1451676, 1516469, 1516657, 1516675, 1622713, 1634838, 1651674
Offset: 1

Views

Author

Bernard Schott, Jan 26 2021

Keywords

Comments

The idea of this sequence comes from a remark of Amiram Eldar in Discussion section of A327750 (m + pod(m) = k with pod(k) = pod(m)) in September 2019.
Question: is it possible to get a longer string of integers with this rule?
From David A. Corneth, Feb 01 2021: (Start)
The product of digits of a(n) is a multiple of 6. In terms below 10^10 however all products of digits of a(n) are a multiple of 36. Is that product a multiple of 36 for all a(n)?
The least term k such that k + 6 is here is k = 56516718. Are there consecutive terms that differ by less than 6? (End)

Examples

			262713 + pod(262713) = 262713 + 504 = 263217, whose product of digits is also 504, and 263217 + 504 = 263721 whose product of digits is again 504; hence, m=262713, k=263217, q=263721 and pod(m)=pod(k)=pod(q)=504, so 262713 is a term.
		

References

  • Roman Fedorov, Alexei Belov, Alexander Kovaldzhi, and Ivan Yashchenko, Moscow-Mathematical Olympiads, 2000-2005, Level A, Problem 2, 2003; MSRI, 2011, p. 15 and 97.

Crossrefs

Subsequence of A327750.

Programs

  • Mathematica
    pod[n_] := Times @@ IntegerDigits[n]; seqQ[n_] := Module[{p = pod[n], k, q}, k = n + p; q = k + pod[k]; p > 0 && Equal @@ {p, pod[k], pod[q]}]; Select[Range[2*10^6], seqQ] (* Amiram Eldar, Jan 26 2021 *)
  • PARI
    isok(m) = my(pm=vecprod(digits(m)), k=m+pm, pk=vecprod(digits(k)), q=k+pk, pq=vecprod(digits(q))); pm && (pm==pk) && (pk==pq); \\ Michel Marcus, Jan 26 2021

Extensions

Terms from Amiram Eldar, Jan 26 2021
Showing 1-2 of 2 results.