Liczby ubogie w czynniki

20

Jeśli dodatnia liczba całkowita ma (ściśle) mniej czynników pierwszych (bez liczenia wielokrotności) niż jej następca i jej poprzednik, nazwiemy ją liczbą złą .N>2

Innymi słowy, i , w którym to liczba unikalnych głównych czynników .ω(N)<ω(N1)ω(N)<ω(N+1)ω(N)N

Zadanie

Możesz wybrać jeden z następujących formatów We / Wy:

  • Weź liczbę całkowitą i wyślij liczbę o niskiej wartości. W przypadku wybrania tego, może być indeksowane 0 lub 1.NNthN
  • Weź dodatnią liczbę całkowitą i wyślij pierwsze liczby N o niskiej liczbie czynników.NN
  • Wydrukuj sekwencję w nieskończoność.

Możesz przyjmować dane wejściowe i dostarczać dane wyjściowe za pomocą dowolnej standardowej metody , w dowolnym języku programowania , zwracając uwagę, że te luki są domyślnie zabronione. To jest kod golfowy, więc wygrywa najkrótsza przesłanka zgodna z zasadami.

Nie dołączę osobnych przypadków testowych, ponieważ metody konkurowania są różne, ale możesz odwołać się do pierwszych 100 terminów tej sekwencji, czyli OEIS A101934 :

11, 13, 19, 23, 25, 27, 29, 37, 41, 43, 47, 49, 53, 59, 61, 64, 67, 71, 73, 79, 81, 83, 89, 97, 101, 103, 107, 109, 113, 121, 125, 131, 137, 139, 149, 151, 155, 157, 163, 167, 169, 173, 179, 181, 191, 193, 197, 199, 211, 221, 223, 227, 229, 233, 239, 241, 243, 251, 259, 263, 265, 269, 271, 277, 281, 283, 289, 293, 307, 309, 311, 313, 317, 331, 337, 341, 343, 347, 349, 353, 359, 361, 365, 367, 371, 373, 379, 383, 389, 397, 401, 407, 409, 419, 421, 431, 433, 439, 441, 443

Na przykład występuje w tej sekwencji, ponieważ (5), (2 i 13) i (2 i 3), więc i .25ω(25)=1ω(26)=2ω(24)=2ω(25)<ω(24)ω(25)<ω(26)

Pan Xcoder
źródło
Czy mogę wyprowadzić wiodące n = przed każdą wartością?
Steadybox
@ Steadybox Sketchy, ale pozwolę na to: - /
Mr. Xcoder
Dodałem go jako alternatywną wersję.
Steadybox

Odpowiedzi:

7

Brachylog , 21 bajtów

⟨+₁≡-₁⟩{ḋdl}ᵐ⌋>~↰₂?ẉ⊥

Wypróbuj online!

Drukuje nieskończenie.

Wyjaśnienie

⟨+₁≡-₁⟩                  Fork: The output of the fork is [Input + 1, Input -1]
       {   }ᵐ            Map:
        ḋdl                (predicate 2) the output is the length of the prime decomposition
                           of the input with no duplicates
             ⌋           Take the minimum result of that map
              >          This minimum is bigger than…
               ~↰₂?      …the output of predicate 2 with Input as input
                  ?ẉ     Write Input followed by a new line if that's the case
                    ⊥    False: backtrack and try another value for Input
Fatalizować
źródło
5

Galaretka , 13 12 bajtów

Cr~ÆvÐṂN⁼Wø#

Wyświetla pierwsze n liczb o niskiej liczbie czynników.

Wypróbuj online!

Jak to działa

Cr~ÆvÐṂN⁼Wø#  Main link. No arguments.

          ø   Wrap the links to the left into a chain and begin a new chain.
           #  Read an integer n from STDIN and call the chain to the left with
              arguments k = 0, 1, 2, ... until n of them return a truthy value.
              Return those n values of k as an array.
C                 Complement; yield -k+1.
  ~               Bitwise NOT; yield -k-1.
 r                Range; yield [-k+1, -k, -k-1].
     ÐṂ           Yield those elements of [-k+1, -k, -k-1] for which the link to
                  the left returns the minimal value.
   Æv                 Count the number of unique prime factors.
                      Note that, for a negative argument, Æv counts -1 as well, and
                      0 is counted as a/the factor of 0. Negating the the arguments
                      eliminates the edge case 1 (no factors), which would be a
                      false positive otherwise.
                  For a factor-poor number, this yields [-k].
       N          Take the negatives of the resulting integers.
         W        Wrap; yield [k].
        ⁼         Test the results to both sides for equality.
Dennis
źródło
5

Python 2 , 123 119 bajtów

q=lambda n:sum(n%i<all(i%j for j in range(2,i))for i in range(2,n+1))
i=2
while 1:
 i+=1
 if q(i-1)>q(i)<q(i+1):print i

Wypróbuj online!

Pręt
źródło
@FryAmTheEggman dzięki! Nawet jeśli nie skorzystałem z twojej sugestii, zainspirowało mnie to do innego podejścia, które pozwoliło zaoszczędzić 4 bajty: D
Rod
Ładny! Byłem pewien, że istnieje sposób na uniknięcie dwóch brzydkich
lambd
4

MATL , 26 24 22 bajtów

`T@3:q+YFg!sdZSd0>?@QD

Drukuje sekwencję w nieskończoność.

Wypróbuj online!

Wyjaśnienie

`         % Do...while loop
  T       %   Push true. Will be used as loop condition
  @       %   Push (1-based) iteration index, k 
  3:q     %   Push [1 2 3] minus 1, that is, [0 1 2]
  +       %   Add, element-wise. Gives [k k+1 k+2]
  YF      %   Exponents of prime-factor decomposition. Gives a 3-row matrix
  g       %   Convert to logical: non-zero numbers become 1
  !s      %   Transpose, sum of each column. Gives a row vector of 3 elements, 
          %   which are the number of unique prime factors of k, k+1 and k+2 
  d       %   Consecutive differences. Gives a row vector of 2 elements
  ZS      %   Sign: replaces each number by -1, 0 or 1
  d       %   Consecutive difference. Gives a single number
  0>      %   Is it positive?
  ?       %   If so
    @Q    %     Push k+1
    D     %     Display
          %   End (implicit)
          % End (implicit). The stack contains true, which (is consumed and)
          % causes an infinite loop
Luis Mendo
źródło
3

Łuska , 22 bajty

f(ΠtSM<←ṙ1mȯLup§…←→)tN

Drukuje sekwencję w nieskończoność, wypróbuj ją online lub wyświetl pierwszą N !

Alternatywnie §oΛ>←t można użyć zamiast ΠtSM<←.

Wyjaśnienie

f(                  )tN  -- filter the tail of the naturals ([2,3…]) by:
  ΠtSM<←ṙ1m(Lup)§…←→     -- | takes a number as argument, example 11
                §…       -- | range of..
                  ←      -- | | the argument decremented (10)
                   →     -- | | to the argument incremented (12)
                         -- | : [10,11,12]
          m(   )         -- | map the following (example on 12) ..
              p          -- | | prime factors: [2,2,3]
             u           -- | | deduplicate: [2,3]
            L            -- | | length: 2
                         -- | : [2,1,2]
        ṙ1               -- | rotate by 1: [1,2,2]
    SM<                  -- | map the function (< X) over the list where X is ..
       ←                 -- | | the first element (1)
                         -- | : [0,1,1]
   t                     -- | tail: [1,1]
  Π                      -- | product: 1
                         -- : [11,13,19,23,25,27,29…
ბიმო
źródło
3

Pyth , 14 bajtów

.f!-.ml{Pb}tZh

Wypróbuj tutaj!

Początkowo była to sugestia dotycząca odpowiedzi Dopappa , ale powiedzieli mi, żebym opublikował ją osobno.

Jak to działa?

.f! -. ml {Pb} tZh | Pełny program Pobiera dane wejściowe ze STDIN, wysyła pierwsze N ​​do STDOUT.

.f | (var: Z) Wyprowadza pierwsze N ​​wartości, które spełniają predykat.
          } tZh | Tworzy listę [Z - 1, Z, Z + 1].
    .m | (var: b) Weź elementy o minimalnej wartości funkcji.
        Pb | Czynniki pierwsze b.
      l {| Deduplikuj, uzyskaj długość.
               | W przypadku liczby ubogiej w czynnik daje to wynik w postaci singletonu.
   - | Usuń z tego Z.
  ! | Logiczna negacja.
Pan Xcoder
źródło
3

Haskell, 105 86 bajtów

Dzięki @Wheat Wizard, @Bruce Forte i @Laikoni za uratowanie 19 bajtów.

[n|n<-[2..],d n<d(n-1),d n<d(n+1)] d x=[1|n<-[1..x],x`rem`n<1,all((>0).rem n)[2..n-1]]

Enderperson1010
źródło
podczas używania rem ==0i /=0można je zastąpić odpowiednio za pomocą <1i >0.
Pszenica Wizard
Nie ma potrzeby letdefiniowania djako funkcji pomocniczej jest w porządku (patrz przewodnik po regułach gry w golfa ). Również summoże być pominięty, porównanie działa tak samo na listach. 86 bajtów: Wypróbuj online!
Laikoni
2

Oktawa ,  87   83  79 bajtów

Podziękowania dla szarlatana @Cows za uratowanie bajtu i podziękowania @Luis Mendo za uratowanie trzech sześciu bajtów!

f=@(n)nnz(unique(factor(n)));n=9;while++n if[f(n-1) f(n+1)]>f(n)disp(n);end;end

Drukuje sekwencję w nieskończoność.

Wypróbuj online!

73 bajty z wiodącymi n =przed każdą wartością:

f=@(n)nnz(unique(factor(n)));n=9;while++n if[f(n-1) f(n+1)]>f(n)n end;end

Wypróbuj online!

Steadybox
źródło
Myślę, że funkcja fmoże stać się f=@(n)length(unique(factor(n)))o jeden bajt mniej.
Kritixi Lithos
2

05AB1E , 14 13 bajtów

Wysyła n-tą liczbę o niskiej liczbie czynników (1-indeksowana)

µ3LN+Íf€gÀć›P

Wypróbuj online!

Wyjaśnienie

µ                 # loop over N (1,2,3, ...) until counter equals input
 3L               # push the range [1,2,3]
   N+             # add current N
     Í            # subtract 2
      f           # get unique prime factors of each
       €g         # get length of each factor list
         À        # rotate left
          ć       # extract the head
           ›      # check if the remaining elements are strictly greater
            P     # product
                  # if 1, increment counter
                  # implicitly output final N
Emigna
źródło
1
Właśnie chciałem zaproponować przejście na µ , więc chyba wskażę moją alternatywę - N<N>Ÿmoże zastąpić 3LN+Í, jeśli to pomoże.
Pan Xcoder,
@ Mr.Xcoder: Podobnie, ®XŸN+działa również. Lub 0®X)N+w takim przypadku Ànie byłoby to konieczne. Niestety wszystkie kończą na tej samej liczbie bajtów.
Emigna
1

Pyth, 30 25 bajtów

#=hTI&>l{PhTKl{PT>l{PtTKT

To jest mój pierwszy prawdziwy golf Pyth, więc wszelkie komentarze są bardzo mile widziane.

Ogromne podziękowania dla Xcodera!

Wyjaśnienie

#                         | Loop until error
 =hT                      | Add one to T (initially 10)
    I&                    | If both...
      >l{PhTKl{PT         | The number of unique prime factors of T+1 is greater than that of T (store in K) 
                 >l{PtTK  | And the number of unique prime factors of T-1 is greater than K (from before)
                        T | Then implicitly print T

TIO .

Daniel
źródło
14 bajtów: .f!-.ml{Pb}tZh(drukuje pierwsze n) ( .fpobiera pierwsze n wartości, które spełniają warunek [1,2,3,...]i używa zmiennej Z, }tZhgeneruje zakres liczb całkowitych [Z - 1 ... Z + 1], .mzwraca listę elementów o minimalnej wartości funkcji (z b), l{Pbpobiera liczbę różnych dzielników, -odrzuca Zz listy, !stosuje logiczną negację)
Pan Xcoder
1
@ Mr.Xcoder, wow, nie sądzę, żebym to zrobił w najbliższym czasie! Czy nie powiedziałbyś, że zasługuje na własną odpowiedź?
Daniel
@Dopapp Ok, więc opublikowałem go osobno. +1 Witamy w Pyth golfa!
Pan Xcoder
25 bajtów przy użyciu twojego podejścia.
Pan Xcoder
Nie. his +1, tis -1, while Kjest zmienną, która jest przypisywana bez =. Na przykład K4przypisuje Kdo 4. Następnie możesz uzyskać do niego dostęp za pomocą K.
Pan Xcoder,
1

JavaScript (ES6), 94 bajty

Zwraca N-ty czynnik o niskiej wartości, indeksowany 0.

f=(i,n=9)=>(P=(n,i=k=1)=>++k>n?0:n%k?P(n,1):i+P(n/k--,0))(n)>P(++n)&P(n)<P(n+1)&&!i--?n:f(i,n)

Wypróbuj online!

W jaki sposób?

Najpierw definiujemy funkcję P (), która zwraca liczbę unikatowych czynników pierwszych danej liczby całkowitej.

P = (                   // P = recursive function taking:
  n,                    //   n = input number to test
  i =                   //   i = flag to increment the number of prime factors
  k = 1                 //   k = current divisor
) =>                    //
  ++k > n ?             // increment k; if k is greater than n:
    0                   //   stop recursion
  :                     // else:
    n % k ?             //   if k is not a divisor of n:
      P(n, 1)           //     do a recursive call with n unchanged and i = 1
    :                   //   else:
      i +               //     add i and
      P(n / k--, 0)     //     do a recursive call with n / k and i = 0; decrement k

Kod owijania brzmi teraz:

f = (                   // given:
  i,                    //   i = input index
  n = 9                 //   n = counter
) =>                    //
  P(n) > P(++n) &       // increment n; if P(n - 1) is greater than P(n)
  P(n) < P(n + 1) &&    // and P(n) is less than P(n + 1)
  !i-- ?                // and we have reached the correct index:
    n                   //   return n
  :                     // else:
    f(i, n)             //   go on with the next value
Arnauld
źródło
1

Japt , 29 27 26 bajtów

Nie do końca zadowolony z tego, ale przynajmniej jest lepszy niż moja pierwsza próba, która miała ponad 40 bajtów!

Zwraca Nliczbę th w sekwencji, indeksowaną 1.

È=Jõ_+X k â ÊÃé)e>Xo)«´U}a

Spróbuj


Wyjaśnienie

Domniemane wprowadzenie liczby całkowitej U.

È                       }a

Zwraca pierwszą liczbę całkowitą, Xktóra zwraca wartość true po przejściu przez następującą funkcję.

=Jõ             )

Przypisz tablicę [-1,0,1]do X.

 _+X      Ã

Przekaż każdy element tej tablicy przez funkcję, która najpierw dodaje bieżącą wartość X.

k â Ê

Uzyskaj długość ( Ê) niepowtarzalnych ( â) czynników pierwszych ( k) wyniku.

é

Obróć wynikowy układ o jeden w prawo.

e>Xo)

Pop ( o) ostatni element z Xi sprawdź, czy wszystkie pozostałe elementy są większe od niego.

«´U

Jeśli tak, zmniejsz Ui sprawdź, czy jest równe 0.

Kudłaty
źródło
1

Python 3 , 97 bajtów

n,g=9,lambda m,k=2,p=1:m//k and(m%k<p%k)+g(m,k+1,p*k*k)
while 1:n+=1;g(n-1)>g(n)<g(n+1)!=print(n)

Teoretycznie drukuje sekwencję w nieskończoność. W praktyce gostatecznie przekracza limit rekurencji.

Wypróbuj online!

Dennis
źródło
1

C (gcc) , 126 bajtów

j,t;o(n){for(j=t=0;j++<n;)if(n%j<1)for(t--;j>1>n%j;)n/=j;t=~t;}
m(J){for(J=3;;J++)if(o(J-1)>o(J)&o(J)<o(J+1))printf("%d,",J);}

Wypróbuj online!

Jonathan Frech
źródło
0

Czysty , 130 123 117 bajtów

import StdEnv
?n=[1\\q<-[p\\p<-[2..n]|and[gcd p i<2\\i<-[2..p-1]]]|n rem q<1]
f=[n\\n<-[3..]| ?n<min(?(n-1))(?(n+1))]

Odpowiada nieskończonej liczbie terminów sekwencji. Ponieważ wszystkie są zagnieżdżone, nie można bardzo dobrze korzystać z redukcji wykresów, dlatego jest dość powolny, nawet przy tak złym algorytmie.

Wypróbuj online!

Obrzydliwe
źródło
0

APL NARS, 124 bajty, 62 znaki

{(⍵>1E4)∨⍵≤0:¯1⋄x←9+⍳10×⍵⋄⍵↑(∊{v←⍴∘∪∘π¨⍵+2-⍳3⋄2=+/v>2⌷v}¨x)/x}

Powinien zwracać odpowiedź do 1E4, a następnie zwracać błąd -1; przypuszcza, że ​​9..10x instrument ma wystarczającą liczbę odpowiednich liczb; test:

  z←{(⍵>1E4)∨⍵≤0:¯1⋄x←9+⍳10×⍵⋄⍵↑(∊{v←⍴∘∪∘π¨⍵+2-⍳3⋄2=+/v>2⌷v}¨x)/x}  
  z 0
¯1
  z 1
11 
  z 20
11 13 19 23 25 27 29 37 41 43 47 49 53 59 61 64 67 71 73 79 
RosLuP
źródło
4
około 150 bajtów?
Kudłaty
@ Shaggy tak, to była jedna przybliżona wartość; ale + - jest odpowiedni dla
golfisty
Według moich obliczeń wynik tutaj to 147 bajtów, a nie 152 (liczba znaków nie ma znaczenia). Dalsza lektura: codegolf.meta.stackexchange.com/q/9428/58974
Kudłaty
@ Shaggy liczba 152 byłaby wielkością w bajtach jednego pliku zawierającego tylko ten kod (skopiuj przeszłość, zapisz ten kod w dokumencie notatnika (okna) i obserwuj „właściwość” tego pliku)
RosLuP
Nie liczymy tutaj bajtów.
Shaggy