introduzione

La sequenza ECG inizia con 1 e 2, quindi la regola è che il termine successivo è il numero intero positivo più piccolo non già presente nella sequenza e il cui fattore comune con l'ultimo termine è maggiore di 1 (non sono coprimi).

I primi termini sono:

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

Si chiama ECG perché il grafico dei termini è abbastanza simile a un ECG.

È la sequenza A064413 nell'OEIS .

Sfida

Devi scrivere una funzione che accetta un intero n come input e genera quanti dei primi n termini della sequenza sono maggiori di n .

Poiché la regola della sequenza inizia con il terzo termine, l'intero di input deve essere maggiore o uguale a 3. Ad esempio, dato l'input 10l'output è 1perché il settimo termine è 12e nessuno degli altri primi dieci termini supera 10.

Casi test

3 -> 1

10 -> 1

100 -> 9

1000 -> 70

Regole

• Per numeri inferiori a 3, la funzione può generare 0 o un codice di errore.
• Nessun'altra regola particolare tranne: è il codice golf, più corto è, meglio è!

Gelatina , 20 19 18 byte

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

Questo è un programma completo.

Provalo online!

Come funziona

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.

Si noti che la sequenza generata è $[1, \mathbf{0}, 2, 4, 6, 3, 9, 12, 8, 10, 5, 15, \dots]$ . Poiché la chiamata del collegamento helper $n$ volte genera una sequenza di lunghezza $n + 1$ , lo $0$ viene praticamente ignorato.

Perl 6 , 66 63 59 58 byte

-4 byte grazie a Jo King

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

Provalo online!

Troppo lento su TIO per n = 1000.

JavaScript (ES6), 107 106 105 byte

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])

Provalo online!

Come?

$C$

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

$a$$[2,1]$$a[0]$

$k$$0$

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

a.indexOf(k) è uguale a:

• $-1$$k$$a$
• $0$$k$ è uguale all'ultimo valore (nel qual caso non è necessariamente coprimi con esso)
• $i\ge 1$

a.indexOf(k) + C(k, a[0])$0$$k$$a$$k$$-1+true=0$

Haskell, 89 82 byte

Modifica: -7 byte grazie a @ 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]

Provalo online!

Buccia , 16 byte

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

Provalo online!

Spiegazione

#>¹↑¡§ḟȯ←⌋→`-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

MATL , 29 byte

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

Provalo online!

Spiegazione:

	#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

APL (Dyalog Unicode) , ^{SBCS da} 39 byte

-2 byte grazie a ngn, -1 byte usando il controllo condizionale appropriato.

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

Provalo online!

JavaScript, 93 91 87 byte

Genera un errore di overflow per 0o 1, output 0per 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)

Provalo online

APL (NARS), caratteri 121, byte 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
∇

test in meno di un minuto qui in tempo di esecuzione:

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

Naturale non esiste alcun controllo per tipo e gamma ...

Japt, 23 21 byte

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

Provalo

@_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

Python 3 , 153 byte

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

Provalo online! (Avviso: la valutazione richiede circa 30 secondi)