50

La "formica primaria" è un animale ostinato che naviga gli interi e li divide fino a quando rimangono solo i numeri primi!

Inizialmente, abbiamo un array infinito A contenente tutti gli interi> = 2: [2,3,4,5,6,.. ]

Lascia che psia la posizione della formica sull'array. Inizialmente, p = 0(l'array è indicizzato 0)

Ad ogni turno, la formica si sposterà come segue:

se A[p]è primo, la formica si sposta nella posizione successiva:p ← p+1
altrimenti, se A[p]è un numero composto, qsia il suo divisore più piccolo> 1. Dividiamo A[p]per qe aggiungiamo qa A[p-1]. La formica si sposta nella posizione precedente:p ← p-1

Ecco le prime mosse per la formica:

 2  3  4  5  6  7  8  9  ... 
 ^
 2  3  4  5  6  7  8  9  ... 
    ^
 2  3  4  5  6  7  8  9  ... 
       ^
 2  5  2  5  6  7  8  9  ... 
    ^
 2  5  2  5  6  7  8  9  ... 
       ^
 2  5  2  5  6  7  8  9  ... 
          ^
 2  5  2  5  6  7  8  9  ... 
             ^
 2  5  2  7  3  7  8  9  ... 
          ^

Il tuo programma dovrebbe produrre la posizione della formica dopo le nmosse. (puoi presumere n <= 10000)

Casi test:

0 => 0
10 => 6
47 => 9
4734 => 274
10000 => 512

Modificare. puoi anche utilizzare elenchi con 1 indice, è accettabile visualizzare i risultati 1, 7, 10, 275, 513 per il caso di test sopra riportato.

Questo è code-golf, quindi vince il codice con il codice più breve in byte.

code-golf sequence primes

— Arnaud
fonte

32

Onestamente ho pensato che ci fosse una formica sul mio schermo quando l'ho visto nelle Domande sulla rete calda.

— Kodos Johnson,

14

Mi chiedo se la sequenza sia ben definita per arbitrariamente grande n(o se il caso composito potrebbe mai spingere la formica a sinistra dell'iniziale 2).

— Martin Ender,

1

@SuperChafouin quindi i risultati di per i casi di test possono essere: 1,7,10,275,513se indicato 1-indexing? O avrebbero ancora bisogno di abbinare le tue uscite.

— Tom Carpenter,

12

@MartinEnder Un'altra domanda aperta è se un primo> 7 alla fine può essere lasciato indietro per sempre.

— Arnauld,

2

@Arnauld Out per quanto riguarda lo spostamento n = 1.000.000.000 (dove p = 17156661), la relazione tra n e p è molto vicina a p = n / (ln (n) * ln (ln (n))).

— Penguino,

11

Alice , 45 byte

/o
\i@/.&wqh!]k.&[&w;;?c]dt.n$k&;[.?~:![?+!kq

Provalo online!

Implementazione per lo più semplice.

I ntempi di ciclo in Alice vengono generalmente eseguiti spingendo i tempi di un indirizzo di ritorno n-1, quindi ritornando alla fine di ogni iterazione con k. L'ultima volta attraverso il ciclo, l' kistruzione non ha nessun posto in cui tornare e l'esecuzione procede.

Questo programma usa le stesse kistruzioni per fermarsi presto quando il numero è primo. Di conseguenza, l'iterazione finale sposta sempre la formica a sinistra. Per compensare questo errore, eseguiamo n+1iterazioni su un array a 1 indice, che fornisce esattamente il risultato desiderato (e fornisce il caso n=0gratuitamente).

— Nitrodon
fonte

7

Python 2 , 120 byte

p=0
A=range(2,input()+2)
for _ in A:
 for q in range(2,A[p]):
	if A[p]%q<1:A[p]/=q;p-=1;A[p]+=q;break
 else:p+=1
print p

Provalo online!

Ah, la rara for- elseciclo! La elseclausola viene eseguita solo se il forcorpo non viene chiuso tramite break. Nel nostro caso, ciò significa che abbiamo controllato tutti i messaggi di posta qelettronica e non ne abbiamo trovato nessuno da dividere p.

— Lynn
fonte

7

Ottava , 109 103 101 94 byte

function i=a(n)i=1;for l=z=2:n+1
if nnz(q=factor(z(i)))>1
z(i--)/=p=q(1);z(i--)+=p;end
i++;end

Provalo online!

Questo codice genererà la posizione in 1-indicizzazione, quindi gli output per i casi di test sono:

0 => 1
10 => 7
47 => 10
4734 => 275
10000 => 513

Questa versione utilizza alcune ottimizzazioni di Ottava, quindi non è compatibile con MATLAB. Il codice seguente è una versione compatibile MATLAB.

MATLAB, 130 123 118 117 byte

function i=a(n)
z=2:n+1;i=1;for l=z
q=factor(z(i));if nnz(q)>1
z(i)=z(i)/q(1);i=i-1;z(i)=z(i)+q(1);else
i=i+1;end
end

Utilizza l'indicizzazione 1 come con la versione Octave. L'ho testato su tutti i casi di test in MATLAB. Ad esempio, l'output a 100000 è 3675 (indicizzazione singola).

Una versione commentata del codice sopra:

function i=a(n)
    z=2:n+1;                %Create our field of numbers
    i=1;                    %Start of at index of 1 (MATLAB uses 1-indexing)
    for l=1:n               %For the required number of iterations
        q=factor(z(i));     %Calculate the prime factors of the current element
        if nnz(q)>1         %If there are more than one, then not prime
            z(i)=z(i)/q(1); %So divide current by the minimum
            i=i-1;          %Move back a step
            z(i)=z(i)+q(1); %And add on the minimum to the previous.
        else
            i=i+1;          %Otherwise we move to the next step
        end
    end

Per motivi di interesse, si tratta delle posizioni delle formiche rispetto al numero di iterazioni per i primi 10000 valori di n.

Sembra probabile che la formica tenderà probabilmente all'infinito, ma chissà, l'aspetto può essere ingannevole.

MATLAB: salvato 6 byte con forinvece di whilee rimuovendo parentesi da if- Grazie @Giuseppe
MATLAB: Salva 2 byte - Grazie @Sanchises
Octave: salva 10 byte usando Octave \=e +=operazioni - Grazie @Giuseppe
Octave: salva 2 byte con i++e i--- Grazie @LuisMendo
Octave: salva 7 byte - Grazie @Sanchises

— Tom Carpenter
fonte

Per farlo funzionare su TIO, penso che tu abbia bisogno enddi abbinare la firma della funzione

— Giuseppe,

@Giuseppe Ah, ok. In MATLAB il trailing endè facoltativo.

— Tom Carpenter,

puoi fare una funzione anonima usando @ (n) all'inizio invece di usare la funzione i = a (n)

— Michthan

@Michthan non può farlo in MATLAB. Non penso che sia possibile in Octave in quanto ha dei loop?

— Tom Carpenter,

1

Il trailing endè facoltativo anche in Octave. Qui è necessario solo perché hai il codice dopo la funzione

— Luis Mendo,

6

JavaScript (ES6), 91 byte

f=(n,a=[p=0])=>n--?f(n,a,(P=n=>++x<n?n%x?P(n):a[a[p]/=x,--p]+=x:p++)(a[p]=a[p]||p+2,x=1)):p

dimostrazione

NB: Potrebbe essere necessario aumentare le dimensioni dello stack predefinite del motore per far passare tutti i casi di test.

Provalo online!

— Arnauld
fonte

6

Haskell , 108 106 94 byte

([0]#[2..]!!)
(a:b)#(p:q)=length b:([b#(a+d:div p d:q)|d<-[2..p-1],mod p d<1]++[(p:a:b)#q])!!0

Provalo online! Esempio di utilizzo: ([0]#[2..]!!) 10rese 6(indicizzate 0).

La funzione #opera su due elenchi, il fronte inverso dell'array [p-1, p-2, ..., 1]e il resto infinito dell'array [p, p+1, p+2, ...]. Costruisce un elenco infinito di posizioni, da cui nviene restituita la th posizione data un input n.

Il modello si ((a:b)#(p:q))lega pal valore della posizione corrente della formica e aal valore della posizione precedente. bè il prefisso dell'array dalla posizione 1 a p-2e qil resto infinito a partire dalla posizione p+1.

Costruiamo un elenco di chiamate ricorsive nel modo seguente: Osserviamo ciascun divisore ddi p(che è più grande di uno e più piccolo di p) in ordine crescente e aggiungiamo b#(a+d:div p d:q)per ciascuno di essi, cioè il valore corrente pviene diviso per de la formica si sposta un passo a sinistra dove dviene aggiunto a. Quindi aggiungiamo (p:a:b)#qalla fine di questo elenco, che indica che la formica si muove di un passo verso destra.

Prendiamo quindi la prima di quelle chiamate ricorsive dall'elenco e anteponiamo la posizione corrente, che coincide con la lunghezza dell'elenco prefissi b. Poiché i divisori sono in ordine crescente, selezionando il primo dall'elenco delle chiamate ricorsive si assicura che utilizziamo il più piccolo. Inoltre, poiché (p:a:b)#qviene aggiunto alla fine dell'elenco, viene selezionato solo se non ci sono divisori ed pè quindi primo.

Modifica:
-2 byte cambiando l'elenco di funzioni dall'ordine decrescente a quello crescente.
-12 byte grazie all'idea di Zgarb di indicizzare in un elenco infinito invece di gestire un contatore e passando a 0-indicizzazione.

— Laikoni
fonte

2

96 byte creando un elenco infinito e indicizzando, invece di portare in giro il contatore.

— Zgarb,

1

@Zgarb Grazie mille! Sono anche solo 94 byte quando si passa all'indicizzazione 0.

— Laikoni,

5

TI-BASIC, 108 103 102 98 byte

Ingresso e uscita sono memorizzati in Ans. L'output è 1 indicizzato.

Ans→N
seq(X,X,2,9³→A
1→P
For(I,1,N
1→X:∟A(P→V
For(F,2,√(V
If X and not(fPart(V/F:Then
DelVar XV/F→∟A(P
P-1→P
F+∟A(P→∟A(P
End
End
P+X→P
End

— kamoroso94
fonte

Puoi togliere un byte fPart(∟A(P)/F:con fPart(F¹∟A(P:. Stessa cosa nella riga successiva.

— Scott Milner,

@ScottMilner Non sempre funziona. not(fPart(7⁻¹7è 0 ma not(fPart(7/7è 1.

— kamoroso94

5

MATL , 41 byte

:Q1y"XHyw)Zp?5MQ}1MtYf1)/H(8MyfHq=*+9M]]&

L'output è basato su 1. Il programma scade per l'ultimo caso di test nell'interprete online.

Provalo online!

Spiegazione

Il programma applica la procedura descritta nella sfida. Per fare ciò, fa un uso insolitamente pesante degli appunti manuali e automatici di MATL.

Il divisore più piccolo si ottiene come prima voce nella decomposizione del fattore primo.

L'aggiornamento "dividere" avviene sovrascrivendo il testo corrispondente della matrice A . L'aggiornamento "aggiungi" viene eseguito aggiungendo ad A un elemento un array che contiene zeri tranne nella posizione desiderata.

:Q        % Implicitly input n. Push array [2 3 ... n+1]. This is the initial array A. 
          % It contains all required positions. Some values will be overwritten
1         % Push 1. This is the initial value for p
y         % Duplicate from below
"         % For each loop. This executes the following n times.
          %   STACK (contents whosn bottom to top): A, p
  XH      %   Copy p into clipboard H
  y       %   Duplicate from below. STACK: A, p, A
  w       %   Swap. STACK: A, A, p
  )       %   Reference indexing. STACK: A, A[p]
  Zp      %   Isprime. STACK: A, false/true
  ?       %   If true (that is, if A[p] is prime)
    5M    %     Push p from automatic clipboard. STACK: A, p
    Q     %     Add 1. STACK: A, p+1
  }       %   Else (that is, if A[p] is not prime)
    1M    %     Push A[p] from automatic clipboard. STACK: A, A[p]
    t     %     Duplicate. STACK: A, A[p], A[p]
    Yf    %     Prime factors, with repetitions. STACK: A, A[p], prime factors of A[p]
    1)    %     Get first element, d. STACK: A, A[p], d
    /     %     Divide. STACK: A, A[p]/d
    H     %     Push p from clipboard H. STACK: A, A[p]/d, p
    (     %     Assignment indexing: write value. STACK: A with A[p] updated
    8M    %     Push d from automatic clipboard.
    y     %     Duplicate from below. STACK: A with A[p] updated, d, A with A[p] updated
    f     %     Find: push indices of nonzero entries.
          %     STACK: A with A[p] updated, d, [1 2 ... n]
    Hq    %     Push p from clipboard H, subtract 1.
          %     STACK: A with A[p] updated, d, [1 2 ... n], p-1
    =     %     Test for equality, element-wise.
          %     STACK: A with A[p] updated, d, [0 ... 0 1 0 ... 0]
    *     %     Multiply, element-wise. STACK: A with A[p] updated, [0 ... 0 d 0 ... 0]
    +     %     Add, element-wise. STACK: A with A[p-1] and A[p] updated
    9M    %     Push p-1 from automatic clipboard.
          %     STACK: A with A[p-1] and A[p] updated, p-1
  ]       %   End if. The stack contains the updated array and index
]         % End for each. Process the next iteration
&         % Specify that the following implicit display function should display only
          % the top of the stack. Implicitly display

— Luis Mendo
fonte

3

Pyth - 44 byte

Implementazione procedurale semplice.

K}2hhQVQJ@KZIP_J=hZ).?=GhPJ XKZ/JG=tZ XZKG;Z

Test Suite .

— Maltysen
fonte

3

PARI / GP, 87 byte

f(n)=A=[2..9^5];p=1;for(i=1,n,q=factor(A[p])[1,1];if(A[p]-q,A[p]/=q;p--;A[p]+=q,p++));p

Abbastanza autoesplicativo (non così golf-ish). Se non si conteggia la f(n)=parte, si tratta di 82 byte. Puoi anche iniziare con n->(85 byte).

È una lingua con 1 indice.

Modifica: la modifica illustrate(n,m)=A=[2..m+1];p=1;for(i=1,n,for(j=1,m,printf("%5s",if(j==p,Str(">",A[j],"<"),Str(A[j]," "))));print();q=factor(A[p])[1,1];if(A[p]!=q,A[p]/=q;p--;A[p]+=q,p++))stamperà un'illustrazione della camminata della formica (dato un terminale abbastanza largo). Ad esempio illustrate(150,25)fornirà i primi 150 passaggi su 25 colonne, in questo modo:

  > 2 <3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2> 3 <4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 3> 4 <5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2> 5 <2 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5> 2 <5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 2> 5 <6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 2 5> 6 <7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 2> 7 <3 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 2 7> 3 <7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 2 7 3> 7 <8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 2 7 3 7> 8 <9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 2 7 3> 9 <4 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 2 7> 6 <3 4 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 2> 9 <3 3 4 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5> 5 <3 3 3 4 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5> 3 <3 3 4 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3> 3 <3 4 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3> 3 <4 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 3> 4 <9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3> 5 <2 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5> 2 <9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 2> 9 <10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5> 5 <3 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5> 3 <10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5 3> 10 <11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5> 5 <5 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5 5> 5 <11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5 5 5> 11 <12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5 5 5 11> 12 <13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5 5 5> 13 <6 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5 5 5 13> 6 <13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5 5 5> 15 <3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5 5> 8 <5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5> 7 <4 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5 7> 4 <5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5 5> 9 <2 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 5> 8 <3 2 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3> 7 <4 3 2 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3 7> 4 <3 2 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3 3> 9 <2 3 2 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 3> 6 <3 2 3 2 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5> 5 <3 3 2 3 2 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5> 3 <3 2 3 2 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3> 3 <2 3 2 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3> 2 <3 2 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2> 3 <2 5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3> 2 <5 3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2> 5 <3 13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 5> 3 <13 14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 5 3> 13 <14 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 5 3 13> 14 <15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 5 3> 15 <7 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 5> 6 <5 7 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2> 7 <3 5 7 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7> 3 <5 7 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3> 5 <7 15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3 5> 7 <15 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3 5 7> 15 <16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3 5> 10 <5 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3> 7 <5 5 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3 7> 5 <5 16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3 7 5> 5 <16 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3 7 5 5> 16 <17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3 7 5> 7 <8 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3 7 5 7> 8 <17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3 7 5> 9 <4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3 7> 8 <3 4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7 3> 9 <4 3 4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2 7> 6 <3 4 3 4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 2> 9 <3 3 4 3 4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3> 5 <3 3 3 4 3 4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5> 3 <3 3 4 3 4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3> 3 <3 4 3 4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3> 3 <4 3 4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 3> 4 <3 4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3> 5 <2 3 4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5> 2 <3 4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2> 3 <4 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 3> 4 <17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2> 5 <2 17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 5> 2 <17 18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 5 2> 17 <18 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 5 2 17> 18 <19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 5 2> 19 <9 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 5 2 19> 9 <19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 5 2> 22 <3 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 5> 4 <11 3 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2> 7 <2 11 3 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7> 2 <11 3 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2> 11 <3 19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2 11> 3 <19 20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2 11 3> 19 <20 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2 11 3 19> 20 <21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2 11 3> 21 <10 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2 11> 6 <7 10 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2> 13 <3 7 10 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2 13> 3 <7 10 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2 13 3> 7 <10 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2 13 3 7> 10 <21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2 13 3> 9 <5 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2 13> 6 <3 5 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 2> 15 <3 3 5 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7> 5 <5 3 3 5 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5> 5 <3 3 5 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5> 3 <3 5 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3> 3 <5 21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 3> 5 <21 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 3 5> 21 <22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 3> 8 <7 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3> 5 <4 7 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 5> 4 <7 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3> 7 <2 7 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7> 2 <7 22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 2> 7 <22 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 2 7> 22 <23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 2> 9 <11 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7> 5 <3 11 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 5> 3 <11 23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 5 3> 11 <23 24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 5 3 11> 23 <24 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 5 3 11 23> 24 <25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 5 3 11> 25 <12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 5 3> 16 <5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 5> 5 <8 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 5 5> 8 <5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 5> 7 <4 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 5 7> 4 <5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7 5> 9 <2 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3 7> 8 <3 2 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5 3> 9 <4 3 2 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 5> 6 <3 4 3 2 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5> 7 <3 3 4 3 2 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7> 3 <3 4 3 2 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3> 3 <4 3 2 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 3> 4 <3 2 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3> 5 <2 3 2 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5> 2 <3 2 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2> 3 <2 5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3> 2 <5 12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3 2> 5 <12 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3 2 5> 12 <25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3 2> 7 <6 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3 2 7> 6 <25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3 2> 9 <3 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3> 5 <3 3 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3 5> 3 <3 25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3 5 3> 3 <25 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3 5 3 3> 25 <26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3 5 3> 8 <5 26
   2 5 5 5 3 3 2 3 5 3 3 5 2 7 5 7 3 5 2 3 5> 5 <4 5 26

— Jeppe Stig Nielsen
fonte

2

Python 2 , 142 127 byte

Risparmiato quindici byte grazie a Sherlock9 .

T=range(2,input()+2);p=0
for _ in T:
 m=T[p];d=min(k for k in range(2,m+1)if m%k<1);p+=1
 if d<m:T[p-1]/=d;p-=2;T[p]+=d
print p

Provalo online!

— Jonathan Frech
fonte

127 byte se non si utilizzaP(n)

— Sherlock9

126 byte

— Felipe Nardi Batista,

@FelipeNardiBatista Sfortunatamente, il tuo suggerimento non sembra funzionare per i casi di testn<=4

— Sherlock9

2

Mathematica, 118 103 byte

(s=Range[2,5^6];t=1;Do[If[PrimeQ@s[[t]],t++,s[[t]]/=(k=#2&@@ Divisors@s[[t]]);s[[t-1]]+=k;t--],#];t-1)&

Provalo online!

Martin Ender ha salvato 15 byte

— J42161217
fonte

Avresti avuto uno spazio randagio di fronte Divisors, puoi usare la notazione infix per Doe puoi semplicemente tornare al tposto di t-1(risultato basato su 1).

— Martin Ender,

2

Python 3 , 158 149 133 byte

Questa è un'implementazione procedurale semplice con una o due stranezze per assicurarsi che il codice funzioni per tutti i casi di test. Uso [*range(2,n+9)]per assicurare che A sia abbastanza grande (tranne che n<3, n+9è più che sufficiente). La elseclausola divide il vecchio A[p]per d, i decrementi pe quindi si aggiunge dal nuovo A[p], che è sicuramente una cattiva pratica di codifica. Altrimenti, abbastanza semplice. Suggerimenti di golf benvenuti!

Modifica: -9 byte senza sympygrazie a Halvard Hummel. -14 byte da Felipe Nardi Batista, -6 byte da alcuni spunti della risposta Python 2 di Jonathan Frech

p,_,*A=range(int(input())+2)
for _ in A:
 m=A[p];d=min(k for k in range(2,m+1)if m%k<1);p+=1
 if d<m:A[p-1]//=d;p-=2;A[p]+=d
print(p)

Provalo online!

— Sherlock9
fonte

145 byte

— Halvard Hummel,

148 byte rendendolo un programma completo

— Felipe Nardi Batista,

if d-m:A[p]...e else:p+=1per salvare un byte

— Felipe Nardi Batista,

143 byte rimuovendo la elsedichiarazione

— Felipe Nardi Batista,

dopo aver rimosso l' elseistruzione, non vi è alcuna differenza in byte rispetto alla versione della funzione

— Felipe Nardi Batista,

2

PHP, 102 + 1 byte

for($a=range(2,$argn);$argn--;$d<$a[+$p]?$a[$p--]/=$d+!$a[$p]+=$d:$p++)for($d=1;$a[+$p]%++$d;);echo$p;

Esegui come pipe -Ro provalo online .

Uscita vuota per input 0; inserire +dopo echoper un letterale0

o usa questa versione con 1 indice (103 + 1 byte):

for($a=range($p=1,$argn);$argn--;$d<$a[$p]?$a[$p--]/=$d+!$a[$p]+=$d:$p++)for($d=1;$a[$p]%++$d;);echo$p;

— Titus
fonte

2

R , 123 byte

Un'implementazione semplice. Viene fornito come una funzione, che accetta il numero di mosse come input e restituisce la posizione p.

Passa sopra la sequenza e sposta il puntatore avanti e indietro secondo le regole. L'output è basato su 0.

Nota: per trovare il fattore primo più piccolo di un numero x, calcola il modulo di x rispetto a tutti i numeri interi da 0 a x. Quindi estrae i numeri con modulo uguale a 0, che sono sempre [0,1, ..., x]. Se il terzo numero non è x, allora è il fattore primo più piccolo di x.

p=function(l){w=0:l;v=w+1;j=1;for(i in w){y=v[j];x=w[!y%%w][3]
if(x%in%c(NA,y))j=j+1
else{v[j]=y/x;j=j-1;v[j]=v[j]+x}}
j-2}

Provalo online!

— NOFP
fonte

2

C (gcc), 152 148 byte

minified

int f(int n){int*A=malloc(++n*4),p=0,i,q;for(i=0;i<n;i++)A[i]=i+2;for(i=1;i<n;i++){for(q=2;A[p]%q;q++);if(A[p++]>q){A[--p]/=q;A[--p]+=q;}}return p;}

Formulato con alcuni commenti

int f(int n) {
  int *A = malloc(++n * 4), p = 0, i, q;
  // Initialize array A
  for (i = 0; i < n; i++)
    A[i] = i + 2;
  // Do n step (remember n was incremented)
  for (i = 1; i < n; i++) {
    // Find smallest divisor
    for (q = 2; A[p] % q; q++)
      ;
    if (A[p++] > q) {
      A[--p] /= q;
      A[--p] += q;
    }
  }
  return p;
}

Funzione principale per i test

#include <stdlib.h>
#include <stdio.h>
int main(int argc, char **argv) {
  if (argc != 2)
    return 2;
  int n = atoi(argv[1]);
  int p = f(n);
  printf("%d => %d\n", n, p);
  return 0;
}

Per mostrare ogni passaggio

Dichiara display () all'interno di f ()

int f(int n) {
  int *A = malloc(++n * 4), p = 0, i, q;
  void display(void) {
    for (int i=0; i < p; i++) {
      printf(" %d", A[i]);
    }
    printf(" \033[1;31m%d\033[m", A[p]);
    if (p+1 < n)
      printf(" %d", A[p+1]);
    printf("\n");
  }
  ...

Visualizzazione chiamate ()
```
  A[i] = i + 2;
display();
```
Visualizzazione chiamate ()
```
  }
  display();
}
```

— user285259
fonte

Puoi radere alcuni byte dichiarando A come un array e inizializzando i tuoi controlli di loop prima dei loop, ove possibile, giusto?

— Ripristina Monica il

1

Clojure, 185 byte

#(loop[[n p][(vec(range 2 1e3))0]i %](if(= i 0)p(recur(if-let[q(first(for[i(range 2(n p)):when(=(mod(n p)i)0)]i))][(assoc n p(/(n p)q)(dec p)(+(n(dec p))q))(dec p)][n(inc p)])(dec i))))

Inoltre, modificare uno "stato" non è l'ideale in Clojure. Dovrai aumentare l'esponente per input più grandi.

— NikoNyrh
fonte

Perché hai usato il pattern matching nel loop? Dovresti riuscire a perdere qualche byte senza quello.

— clismique,

Inoltre, potresti essere in grado di cambiare la firstcosa in una somedichiarazione.

— clismique,

Senza pattern matching ho dovuto ripetere recurdue volte, una per ciascun if-letramo. Inoltre (dec i)sarebbe duplicato. someha bisogno di un predicato, potrei usare +come abbiamo a che fare con i numeri, ma questo è un personaggio più lungo di first. CMIIW

— NikoNyrh,

1

Java 8, 138 135 byte

n->{int a[]=new int[++n],s=0,p=0,j=0;for(;j<n;a[j++]=j+1);for(;++s<n;p++)for(j=1;++j<a[p];)if(a[p]%j<1){a[p--]/=j;a[p--]+=j;}return p;}

Spiegazione:

Provalo qui.

n->{                     // Method with integer as both parameter and return-type
  int a[]=new int[++n],  //  Integer-array with a length of `n+1`
      s=0,               //  Steps-counter (starting at 0)
      p=0,               //  Current position (starting at 0)
      j=0;               //  Index integer (starting at 0)
  for(;j<n;              //  Loop (1) from 0 to the input (inclusive due to `++n` above)
    a[j++]=j+1           //   And fill the array with 2 through `n+2`
  );                     //  End of loop (1)
  for(;++s<n;            //  Loop (2) `n` amount of steps:
      p++)               //    And after every iteration: increase position `p` by 1
    for(j=1;             //   Reset `j` to 1
        ++j<a[p];)       //   Inner loop (3) from 2 to `a[p]` (the current item)
      if(a[p]%j<1){      //    If the current item is divisible by `j`:
        a[p--]/=j;       //     Divide the current item by `j`
        a[p--]+=j;}      //     And increase the previous item by `j`
                         //     And set position `p` two steps back (with both `p--`)
                         //   End of inner loop (3) (implicit / single-line body)
                         //  End of loop (2) (implicit / single-line body)
  return p;              //  Return the resulting position `p`
}                        // End of method

— Kevin Cruijssen
fonte

1

Clojure, 198 193 191 byte

Questo deve essere fortemente golfato ...

#(loop[i(vec(range 2(+ % 9)))c 0 p 0](if(= % c)p(let[d(dec p)u(i p)f(some(fn[n](if(=(mod u n)0)n))(range 2(inc u)))e(= u f)](recur(if e i(assoc i d(+(i d)f)p(/ u f)))(inc c)(if e(inc p)d)))))

Golf 1 : salvato 5 byte cambiando (first(filter ...))in(some ...)

Golf 2 : salvato 2 byte cambiando (zero? ...)in(= ... 0)

Uso:

(#(...) 10000) => 512

Codice non golfato:

(defn prime-ant [n]
  (loop [counter 0
         pos 0
         items (vec (range 2 (+ n 9)))]
    (if (= n counter) pos
      (let [cur-item (nth items pos)
            prime-factor
            (some #(if (zero? (mod cur-item %)) %)
              (range 2 (inc cur-item)))
            equals? (= cur-item prime-factor)]
        (recur
          (inc counter)
          (if equals? (inc pos) (dec pos))
          (if equals? items
            (assoc items
              (dec pos) (+ (items (dec pos)) prime-factor)
              pos (/ cur-item prime-factor))))))))

— clismique
fonte

La formica principale 🐜

Alice , 45 byte

Python 2 , 120 byte

Ottava , 109 103 101 94 byte

MATLAB, 130 123 118 117 byte

JavaScript (ES6), 91 byte

dimostrazione

Haskell , 108 106 94 byte

TI-BASIC, 108 103 102 98 byte

MATL , 41 byte

Spiegazione

Pyth - 44 byte

PARI / GP, 87 byte

Python 2 , 142 127 byte

Mathematica, 118 103 byte

Python 3 , 158 149 133 byte

PHP, 102 + 1 byte

R , 123 byte

C (gcc), 152 148 byte

minified

Formulato con alcuni commenti

Funzione principale per i test

Per mostrare ogni passaggio

Clojure, 185 byte

Java 8, 138 135 byte

Clojure, 198 193 191 byte

Uso:

Codice non golfato: