Cyfrowe sekwencje produktów

22

Oto interesująca sekwencja odkryta przez Paula Loomisa, matematyka z Bloomsburg University. Z jego strony w tej sekwencji:

Zdefiniuj za
f(n) = f(n-1) + (the product of the nonzero digits of f(n-1))
f(0) = xpomocą xdowolnej dodatniej liczby całkowitej zapisanej w bazie 10.

Tak więc, zaczynając od f(0)=1, otrzymujesz następującą sekwencję
1, 2, 4, 8, 16, 22, 26, 38, 62, 74, 102, 104, ...

Do tej pory tak standardowe. Interesująca właściwość wchodzi w grę, gdy weźmiesz inną liczbę całkowitą jako punkt początkowy, ostatecznie sekwencja zbiegnie się w punkt wzdłuż powyższej x=1sekwencji. Na przykład, zaczynając od x=3wydajności
3, 6, 12, 14, 18, 26, 38, 62, 74, 102, ...

Oto kilka kolejnych sekwencji, z których każda pokazana jest tylko do osiągnięcia 102:

5, 10, 11, 12, 14, 18, 26, 38, 62, 74, 102, ...
7, 14, 18, 26, 38, 62, 74, 102, ...
9, 18, 26, 38, 62, 74, 102, ...
13, 16, 22, 26, 38, 62, 74, 102, ...
15, 20, 22, 26, 38, 62, 74, 102, ...
17, 24, 32, 38, 62, 74, 102, ...
19, 28, 44, 60, 66, 102, ...

Przypuszczał i udowodnił empirycznie x=1,000,000, że ta właściwość (tj. Że wszystkie liczby wejściowe są zbieżne w tej samej sekwencji) jest prawdziwa.

Wyzwanie

Biorąc pod uwagę dodatnią liczbę całkowitą wejściową, wypisz liczbę 0 < x < 1,000,000, w której f(x)sekwencja zbiega się z f(1)sekwencją. Na przykład x=5byłoby to 26, ponieważ jest to pierwsza wspólna liczba dla obu sekwencji.

 x output
 1 1
 5 26
19 102
63 150056

Zasady

  • Jeśli dotyczy, możesz założyć, że wejście / wyjście będzie pasować do rodzimego typu Integer w twoim języku.
  • Dane wejściowe i wyjściowe można podać dowolną dogodną metodą .
  • Dopuszczalny jest pełny program lub funkcja. Jeśli funkcja, możesz zwrócić dane wyjściowe zamiast je drukować.
  • Standardowe luki są zabronione.
  • To jest więc obowiązują wszystkie zwykłe zasady gry w golfa, a wygrywa najkrótszy kod (w bajtach).
AdmBorkBork
źródło

Odpowiedzi:

5

JavaScript (ES6), 81 67 bajtów

Zapisano 1 bajt dzięki @ l4m2

f=(n,x=1)=>x<n?f(x,n):x>n?f(+[...n+''].reduce((p,i)=>p*i||p)+n,x):n

Wypróbuj online!

Skomentował

f = (n,                   // n = current value for the 1st sequence, initialized to input
        x = 1) =>         // x = current value for the 2nd sequence, initialized to 1
  x < n ?                 // if x is less than n:
    f(x, n)               //   swap the sequences by doing a recursive call to f(x, n)
  :                       // else:
    x > n ?               //   if x is greater than n:
      f(                  //     do a recursive call with the next term of the 1st sequence:
        +[...n + '']      //       coerce n to a string and split it
        .reduce((p, i) => //       for each digit i in n:
          p * i || p      //         multiply p by i, or let p unchanged if i is zero
        ) + n,            //       end of reduce(); add n to the result
        x                 //       let x unchanged
      )                   //     end of recursive call
    :                     //   else:
      n                   //     return n
Arnauld
źródło
`` `f = (n, x = 1) => x <n? f (x, n): x> n? f (+ [... n + '']. zmniejsz ((p, i) = > p * i || p) + n, x): n ``
l4m2
4

Galaretka , 18 14 bajtów

ḊḢDo1P+Ʋ;µQƑ¿Ḣ

Dane wejściowe to tablica singletonów.

Wypróbuj online!

Jak to działa

ḊḢDo1P+Ʋ;µQƑ¿Ḣ  Main link. Argument: [n]

            ¿   While...
          QƑ      all elements of the return value are unique...
         µ          execute the chain to the left.
Ḋ                     Dequeue; remove the first item.
 Ḣ                    Head; extract the first item.
                      This yields the second item of the return value if it has
                      at least two elements, 0 otherwise.
       Ʋ              Combine the links to the left into a chain.
  D                     Take the decimal digits of the second item.
   o1                   Perform logical OR with 1, replacing 0's with 1's.
     P                  Take the product.
      +                 Add the product with the second item.
        ;             Prepend the result to the previous return value.
             Ḣ  Head; extract the first item.
Dennis
źródło
2

Python 2 , 111 95 93 bajtów

Używając rozpakowywania replace(*'01')jak w @Rod odpowiedz
-18 bajtów dzięki @Lynn

l=[1,input()]
while cmp(*l):l[0]+=eval('*'.join(`l[0]`.replace(*'01')));l.sort()
print l[0]  

Wypróbuj online!

Dead Possum
źródło
1
Ach, i stan pętli może być while cmp(*l)również!
Lynn
@Lynn Tak! Jeszcze raz dziękuję
Dead Possum
2

Haskell , 65 bajtów

(!1)
a!b|a>b=b!a|a<b=(a+product[read[d]|d<-show a,d>'0'])!b|1<2=a

Wypróbuj online!

Lynn
źródło
2

Python 2 , 78 bajtów

f=lambda a,b=1:a*(a==b)or f(*sorted([a+eval('*'.join(`a`.replace(*'01'))),b]))

Wypróbuj online!

Lynn
źródło
Pracowałem nad rozwiązaniem z lambda, ale utknąłem z shortcircuting na kilka minut, dobra robota!
Dead Possum
2

Łuska , 13 bajtów

→UΞm¡S+ȯΠf±dΘ

Pobiera dane wejściowe jako listę singletonów.

Wypróbuj online!

Wyjaśnienie

                 Implicit input, e.g 5
            Θ    Prepend a zero to get  [0,5]
   m             Map the following over [0,5]
    ¡              Iteratatively apply the following function, collecting the return values in a list
           d         Convert to a list of digits
         f±          keep only the truthy ones
       ȯΠ            then take the product
     S+              add that to the original number
                After this map, we have [[0,1,2,4,8,16,22,26,38,62...],[5,10,11,12,14,18,26,38,62,74...]]
  Ξ             Merge the sorted lists:  [0,1,2,4,5,8,10,11,12,14,16,18,22,26,26,38,38,62,62,74...]
 U              Take the longest unique prefix: [0,1,2,4,5,8,10,11,12,14,16,18,22,26]
→               Get the last element and implicitely output: 26
H.PWiz
źródło
1

Python 3 , 126 125 bajtów

m=[1]
n=[int(input())]
while not{*m}&{*n}:
 for l in m,n:l+=l[-1]+eval('*'.join(str(l[-1]).replace(*'01'))),
print({*m}&{*n})

Wypróbuj online!

Weź dane jako ciąg

Pręt
źródło
0

J , 50 bajtów

milcząca definicja funkcji stylu

[:{.@(e.~#])/[:(+[:*/@(*#])(#~10)&#:)^:(<453)"0,&1

jeśli argument (powiedzmy 63) został wklejony do wyrażenia REPL, może to być 45 np

{.(e.~#])/(+[:*/@(*#])(#~10)&#:)^:(<453)"0]1,63
  • ,&1 dołącz 1, aby wygenerować sekwencję wyszukiwania oraz sekwencję argumentów
  • ^:(<453)"0 iteruje każdy do momentu osiągnięcia 1 miliona w sekwencji 1
  • + [: */@(*#]) (#~10)&#: widelec dodaje do haka, który tworzy iloczyn cyfr
  • (e.~ # ])/ używa powtórz element, jeśli istnieje, aby uzyskać przecięcie list
  • {. zwraca tylko pierwszą wspólną wartość

Wypróbuj online!

jayprich
źródło
0

R , 110 86 bajtów

o=c(1,1:9);o=o%o%o%o%o;o=c(o%o%o)
x=c(1,n);while((x=sort(x))<x[2])x[1]=(x+o[x+1])[1]
x

TIO

poprzednia wersja 110:

f=function(x){if((x[1]=x[1]+(c((y=(y=c(1,1:9))%o%y%o%y)%o%y))[x[1]+1])==x[2]){x[1]}else{f(sort(x))}}
f(c(1,n))

TIO

jayprich
źródło