Warunki sekwencji EKG

13

Wprowadzenie

Sekwencja EKG zaczyna się od 1 i 2, a następnie reguła jest taka, że ​​następny warunek jest najmniejszą dodatnią liczbą całkowitą, która nie jest jeszcze w sekwencji i której wspólny czynnik z ostatnim terminem jest większy niż 1 (nie są koprime).

Pierwsze warunki to:

1, 2, 4, 6, 3, 9, 12, 8, 10, 5, 15, ...

Nazywa się EKG, ponieważ wykres jego warunków jest dość podobny do EKG.

Jest to sekwencja A064413 w OEIS .

Wyzwanie

Musisz napisać funkcję, która przyjmuje liczbę całkowitą n jako dane wejściowe i wyjściowe, o ile n pierwszych składników sekwencji jest większych niż n .

Ponieważ reguła sekwencji zaczyna się od trzeciego członu, wejściowa liczba całkowita musi być większa lub równa 3. Na przykład, dane wejściowe dają 10wynik, 1ponieważ siódmy człon jest 12i żaden z pozostałych dziesięciu pierwszych członów nie przekracza 10.

Przypadki testowe

3 -> 1

10 -> 1

100 -> 9

1000 -> 70

Zasady

  • W przypadku liczb całkowitych niższych niż 3 funkcja może wygenerować 0 lub kod błędu.
  • Żadnych innych szczególnych zasad z wyjątkiem: to jest golf golfowy, im krótszy, tym lepiej!
David
źródło
Czy możemy zastosować indeksowanie 0, przy 1czym jest to 0 wyraz sekwencji, a zatem, na przykład, 1510 wyraz, a nie 5?
Kudłaty
@Shaggy Myślę, że to sprawiedliwe, aby użyć tego jako metody matematycznej, ale tak naprawdę zmieni wynik przypadków testowych, a nawet samą zadaną funkcję. Dlatego uważam, że nie powinieneś tego robić. Przepraszam.
David
oeis.org/A064413/graph - OEIS może pisać wykresy? Schludny.
Magic Octopus Urn

Odpowiedzi:

7

Galaretka , 20 19 18 bajtów

S‘gṪ’ɗƇḟ¹Ṃṭ
1Ç¡>¹S

To jest pełny program.

Wypróbuj online!

Jak to działa

1Ç¡>¹S       Main link. Argument: n (integer)

1            Set the return value to 1.
 Ç¡          Call the helper link n times.
   >¹        Compare the elements of the result with n.
     S       Take the sum, counting elements larger than n.


S‘gṪ’ɗƇḟ¹Ṃṭ  Helper link. Argument: A (array or 1)

S            Take the sum of A.
 ‘           Increment; add 1.
     ɗƇ      Drei comb; keep only elements k of [1, ..., sum(A)+1] for which the
             three links to the left return a truthy value.
  g              Take the GCD of k and all elements of A.
   Ṫ             Tail; extract the last GCD.
    ’            Decrement the result, mapping 1 to 0.
       ḟ¹    Filterfalse; remove the elements that occur in A.
         Ṃ   Take the minimum.
          ṭ  Tack; append the minimum to A.

Zauważ, że wygenerowana sekwencja to . Ponieważ wywołanie linii pomocniczej razy generuje sekwencję o długości , jest praktycznie ignorowane.[1,0,2,4,6,3,9,12,8,10,5,15,]nn+10

Dennis
źródło
6

Perl 6 , 66 63 59 58 bajtów

-4 bajty dzięki Jo King

{sum (1,2,{+(1...all *gcd@_[*-1]>1,*∉@_)}...*)[^$_]X>$_}

Wypróbuj online!

Zbyt wolne dla TIO dla n = 1000.

nwellnhof
źródło
@JoKing Po tym, jak zdałem sobie sprawę, że first &f,1..*można to przepisać +(1...&f), Twoja sztuczka skrzyżowania pomogła.
nwellnhof
4

JavaScript (ES6), 107 106 105 bajtów

f=(n,a=[2,1],k=3)=>a[n-1]?0:a.indexOf(k)+(C=(a,b)=>b?C(b,a%b):a>1)(k,a[0])?f(n,a,k+1):(k>n)+f(n,[k,...a])

Wypróbuj online!

W jaki sposób?

Funkcja pomocnicza zwraca wartość true, jeśli dwie podane liczby całkowite nie są koprime:C

C = (a, b) => b ? C(b, a % b) : a > 1

Tablica jest inicjowana na i przechowuje wszystkie napotkane do tej pory wartości w odwrotnej kolejności. Dlatego ostatnią wartością jest zawsze .a[2,1]a[0]

Aby wiedzieć, czy kwalifikuje się jako następny termin sekwencji, sprawdzamy, czy następujące wyrażenie jest równe :k0

a.indexOf(k) + C(k, a[0])

a.indexOf(k) jest równy albo:

  • 1 jeśli nie znaleziono wka
  • 0 jeśli jest równe ostatniej wartości (w takim przypadku niekoniecznie jest to coprime z nią)k
  • niektóre inaczeji1

a.indexOf(k) + C(k, a[0])0kak1+true=0

Arnauld
źródło
4

Haskell, 89 82 bajtów

Edycja: -7 bajtów dzięki @ H.PWiz

f l=[n|n<-[3..],all(/=n)l,gcd(l!!0)n>1]!!0:l
g n=sum[1|x<-iterate f[2]!!(n-2),x>n]

Wypróbuj online!

nimi
źródło
82 bajty
H.PWiz
@ H.PWiz: ah, to sprytne. Dzięki!
nimi
4

Łuska , 16 bajtów

#>¹↑¡§ḟȯ←⌋→`-Nḣ2

Wypróbuj online!

Wyjaśnienie

#>¹↑¡§ḟȯ←⌋→`-Nḣ2  Implicit input, say n=10
              ḣ2  Range to 2: [1,2]
    ¡             Construct an infinite list, adding new elements using this function:
                   Argument is list of numbers found so far, say L=[1,2,4]
             N     Natural numbers: [1,2,3,4,5,6,7...
           `-      Remove elements of L: K=[3,5,6,7...
      ḟ            Find first element of K that satisfies this:
                    Argument is a number in K, say 6
     §    →         Last element of L: 4
         ⌋          GCD: 2
       ȯ←           Decrement: 1
                    Implicitly: is it nonzero? Yes, so 6 is good.
                  Result is the EKG sequence: [1,2,4,6,3,9,12...
   ↑              Take the first n elements: [1,2,4,6,3,9,12,8,10,5]
#                 Count the number of those
 >¹               that are larger than n: 1
Zgarb
źródło
3

MATL , 29 bajtów

qq:2:w"GE:yX-y0)yZdqg)1)h]G>z

Wypróbuj online!

Wyjaśnienie:

	#implicit input, n, say 10
qq:	#push 1:8
2:	#push [1 2]. Stack: {[1 .. 8], [1 2]}
w	#swap top two elements on stack
"	#begin for loop (do the following n-2 times):
 GE:	#push 1...20. Stack: {[1 2], [1..20]}
 y	#copy from below. Stack:{[1 2], [1..20], [1 2]}
 X-	#set difference. Stack: {[1 2], [3..20]}
 y0)	#copy last element from below. Stack:{[1 2], [3..20], 2}
 yZd	#copy from below and elementwise GCD. Stack:{[1 2], [3..20],[1,2,etc.]}
 qg)	#select those with gcd greater than 1. Stack:{[1 2], [4,6,etc.]}
 1)	#take first. Stack:{[1 2], 4}
 h	#horizontally concatenate. Stack:{[1 2 4]}
 ]	#end of for loop
G>z	#count those greater than input
	#implicit output of result
Giuseppe
źródło
proszę wyjaśnić, dlaczego podwajasz dane wejściowe (z GE:)?
David
2
a(n)2na(n)n2n=1000while
3

APL (Dyalog Unicode) , 39 bajtów SBCS

-2 bajty dzięki ngn, -1 bajt przy użyciu właściwego sprawdzania warunkowego.

{+/⍵<⍵⍴3{(1=⍺∨⊃⌽⍵)∨⍺∊⍵:⍵∇⍨⍺+1⋄⍵,⍺}⍣⍵⍳2}

Wypróbuj online!

jastrząb
źródło
przekazuje swój lewy argument do funkcji operandu, więc nie ma takiej potrzeby . także nie będzie wiązać się z rzeczą po prawej, ponieważ zaczyna się od funkcji ( ), więc nie ma takiej potrzeby .
ngn
2

JavaScript, 93 91 87 bajtów

Zgłasza błąd przepełnienia dla 0lub 1, wyniki 0dla 2.

n=>(g=x=>n-i?g[++x]|(h=(y,z)=>z?h(z,y%z):y)(x,c)<2?g(x):(g[c=x]=++i,x>n)+g(2):0)(i=c=2)

Wypróbuj online

Kudłaty
źródło
2

APL (NARS), znaki 121, bajty 242

∇r←a w;i;j;v
r←w⋄→0×⍳w≤2⋄i←2⋄r←⍳2⋄v←1,1,(2×w)⍴0
j←¯1+v⍳0
j+←1⋄→3×⍳1=j⊃v⋄→3×⍳∼1<j∨i⊃r⋄r←r,j⋄i+←1⋄v[j]←1⋄→2×⍳w>i
r←+/w<r
∇

przetestuj tutaj w czasie krótszym niż minuta:

  a¨3 10 100 1000 2000
1 1 9 70 128 

Naturalnie nie ma sprawdzania typu i zasięgu ...

RosLuP
źródło
1

Japt, 23 21 bajtów

@_jX ªAøZ}f}gA=ì)Aè>U

Spróbuj

@_jX ªAøZ}f}gA=ì)Aè>U
                          :Implicit input of integer U
             A            :10
               ì          :Digit array
              =           :Reassign to A
@          }g             :While the length of A < U+1, take the last element as X,
                          :pass it through the following function & push the result to A
 _       }f               :  Find the first integer Z >= 0 that returns falsey
  jX                      :    Is Z co-prime with X?
     ª                    :    OR
      AøZ                 :    Does A contain Z?
                )         :End loop
                 Aè>U     :Count the elements in A that are greater than U
Kudłaty
źródło
1

Python 3 , 153 bajty

import math
def f(n):
 l=[1];c=0
 for _ in range(n-1):l+=[min(*range(2,n*4),key=lambda x:n*8if x in l or math.gcd(x,l[-1])<2else x)];c+=l[-1]>n
 return c

Wypróbuj online! (Ostrzeżenie: ocena zajmuje około 30 sekund)

Czarna sowa Kai
źródło