Semut utama 🐜


50

"Semut utama" adalah hewan keras kepala yang menavigasi bilangan bulat dan membaginya sampai hanya tersisa bilangan prima!


Awalnya, kami memiliki array tak terbatas A yang berisi semua bilangan bulat> = 2: [2,3,4,5,6,.. ]

Membiarkan pmenjadi posisi semut pada array. Mulanya,p = 0 (array diindeks 0)

Setiap belokan, semut akan bergerak sebagai berikut:

  • jika A[p]prima, semut bergerak ke posisi berikutnya:p ← p+1
  • lain, jika A[p]adalah bilangan komposit, biarkan qmenjadi pembagi yang lebih kecil> 1. Kami membagi A[p]oleh q, dan kami tambahkan quntuk A[p-1]. Semut bergerak ke posisi sebelumnya:p ← p-1

Inilah langkah pertama untuk semut:

 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  ... 
          ^

Program Anda harus menampilkan posisi semut setelah nbergerak. (Anda bisa berasumsi n <= 10000)

Kasus uji:

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

Sunting. Anda juga dapat menggunakan daftar yang diindeks 1, dapat diterima untuk menampilkan hasil 1, 7, 10, 275, 513 untuk test case di atas.

Ini adalah kode-golf, jadi kode dengan kode terpendek dalam byte akan menang.


32
Jujur saya pikir ada semut di layar saya ketika saya melihat ini di Pertanyaan Jaringan Panas.
Kodos Johnson

14
Saya bertanya-tanya apakah urutannya didefinisikan dengan baik untuk sembarang besar n(atau apakah case komposit bisa mendorong semut ke kiri inisial 2).
Martin Ender

1
@ SupupChafouin sehingga output untuk kasus uji bisa: 1,7,10,275,513jika 1-pengindeksan dinyatakan? Atau apakah mereka masih harus mencocokkan output Anda.
Tom Carpenter

12
@ MartinEnder Pertanyaan terbuka lainnya adalah apakah prime> 7 akhirnya dapat ditinggalkan untuk selamanya.
Arnauld

2
@Arnauld Out sejauh memindahkan n = 1.000.000.000 (di mana p = 17156661), hubungan antara n dan p sangat dekat dengan p = n / (ln (n) * ln (ln (n))).
Penguino

Jawaban:


11

Alice , 45 byte

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

Cobalah online!

Sebagian besar implementasi langsung.

nWaktu Looping di Alice biasanya dilakukan dengan menekan n-1kali alamat kembali , kemudian kembali pada akhir setiap iterasi dengan k. Terakhir kali melalui loop, kinstruksi tidak punya tempat untuk kembali, dan eksekusi berlanjut.

Program ini menggunakan kinstruksi yang sama untuk berhenti lebih awal ketika nomornya prima. Akibatnya, iterasi terakhir akan selalu menggerakkan semut ke kiri. Untuk mengimbangi bug ini, kami melakukan n+1iterasi pada larik 1-diindeks, yang memberikan hasil persis seperti yang kita inginkan (dan memberikan kasing n=0gratis).


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

Cobalah online!

Ah, langka for- elselingkaran! The elseklausul hanya mengeksekusi jika fortubuh tidak keluar melalui break. Dalam kasus kami, ini berarti kami memeriksa semua qdan tidak menemukan satupun dari mereka untuk dibagi p.


7

Oktaf , 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

Cobalah online!

Kode ini akan menampilkan posisi dalam pengindeksan 1, jadi output untuk kasus uji adalah:

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

Versi ini menggunakan beberapa optimisasi Oktaf sehingga tidak kompatibel dengan MATLAB. Kode di bawah ini adalah versi yang kompatibel dengan 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

Menggunakan pengindeksan 1 seperti pada versi Oktaf. Saya sudah mengujinya terhadap semua kasus uji di MATLAB. Sebagai contoh, output pada 100000 adalah 3675 (satu-pengindeksan).

Versi komentar dari kode di atas:

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

Yang menarik, ini adalah posisi semut vs jumlah iterasi untuk 10000 nilai pertama dari n.

Posisi Semut

Tampaknya Semut mungkin akan cenderung tak hingga, tetapi siapa tahu, penampilan bisa menipu.


  • MATLAB: Disimpan 6 byte dengan foralih - alih whiledan menghapus tanda kurung dari if- Terima kasih @Giuseppe
  • MATLAB: Hemat 2 byte - Terima kasih @Sanchises
  • Oktaf: Hemat 10 byte dengan menggunakan Oktaf \=dan +=operasi - Terima kasih @Giuseppe
  • Oktaf: Hemat 2 byte dengan i++dan i--- Terima kasih @LuisMendo
  • Oktaf: Hemat 7 byte - Terima kasih @Sanchises

Untuk membuatnya bekerja di TIO, saya pikir Anda perlu endmencocokkan tanda tangan fungsi
Giuseppe

@Giuseppe Ah, baiklah. Dalam MATLAB, trailing endbersifat opsional.
Tom Carpenter

Anda dapat membuat fungsi anonim dengan menggunakan @ (n) di awal alih-alih menggunakan fungsi i = a (n)
Michthan

@Michthan tidak bisa melakukan itu di MATLAB. Saya tidak berpikir itu mungkin dalam Oktaf karena memiliki loop?
Tom Carpenter

1
Trailing endjuga opsional dalam Oktaf. Ini hanya diperlukan karena Anda memiliki kode setelah fungsi
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

Demo

NB: Anda mungkin harus menambah ukuran tumpukan standar mesin Anda agar bisa lulus semua test case.

Cobalah online!


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

Cobalah online! Contoh penggunaan: ([0]#[2..]!!) 10hasil 6(diindeks 0).

Fungsi #beroperasi pada dua daftar, bagian depan terbalik array [p-1, p-2, ..., 1]dan sisa array yang tak terbatas [p, p+1, p+2, ...]. Ini membangun daftar posisi yang tak terbatas, dari mana nposisi th dikembalikan diberi input n.

Pola ini ((a:b)#(p:q))mengikat pnilai posisi semut saat ini dan ake nilai posisi semut sebelumnya. badalah awalan dari array dari posisi 1 ke p-2dan qsisanya tak terbatas mulai dari posisi p+1.

Kami membangun daftar panggilan rekursif dengan cara berikut: Kami melihat setiap pembagi ddari p(yang lebih besar dari satu dan lebih kecil dari p) dalam urutan menaik dan menambahkan b#(a+d:div p d:q)untuk masing-masing, yaitu nilai saat pini dibagi oleh ddan bergerak semut satu langkah ke kiri tempat dditambahkan a. Kemudian kita tambahkan (p:a:b)#qke akhir daftar ini, yang menunjukkan semut bergerak satu langkah ke kanan.

Kami kemudian mengambil panggilan rekursif pertama dari daftar dan menambahkan posisi saat ini, yang bertepatan dengan panjang daftar awalan b. Karena pembagi berada dalam urutan menaik, memilih yang pertama dari daftar panggilan rekursif memastikan kami menggunakan yang terkecil. Selain itu, karena (p:a:b)#qditambahkan ke akhir daftar, itu hanya dipilih jika tidak ada pembagi dan pdengan demikian prima.

Suntingan:
-2 byte dengan mengalihkan daftar fungsi dari turun ke urutan naik.
-12 byte berkat ide Zgarb untuk mengindeks ke dalam daftar tanpa batas alih-alih menangani penghitung, dan dengan beralih ke pengindeksan 0.


2
96 byte dengan membangun daftar dan pengindeksan yang tak terbatas, alih-alih membawa-bawa penghitung.
Zgarb

1
@ Zgarb Terima kasih banyak! Itu bahkan hanya 94 byte ketika beralih ke pengindeksan-0.
Laikoni

5

TI-BASIC, 108 103 102 98 byte

Input dan output disimpan di Ans. Output diindeks 1.

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

Anda dapat melepas byte fPart(∟A(P)/F:dengan fPart(F¹∟A(P:. Hal yang sama di baris berikutnya.
Scott Milner

@ScottMilner Itu tidak selalu berhasil. not(fPart(7⁻¹7adalah 0 tetapi not(fPart(7/71.
kamoroso94

5

MATL , 41 byte

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

Output berbasis 1. Waktu program habis untuk test case terakhir dalam juru bahasa online.

Cobalah online!

Penjelasan

Program menerapkan prosedur sebagaimana dijelaskan dalam tantangan. Untuk melakukannya, ia menggunakan clipboard manual dan otomatis MATL yang luar biasa berat.

Pembagi terkecil diperoleh sebagai entri pertama dalam dekomposisi faktor prima.

"Kesenjangan" update dilakukan oleh Timpa yang sesuai masuknya berbagai A . The "add" update dilakukan oleh unsur-bijaksana menambah Sebuah array yang berisi angka nol kecuali pada posisi yang diinginkan.

: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

3

Pyth - 44 byte

Implementasi prosedur yang mudah.

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

Test Suite .


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

Cukup jelas (tidak begitu golf-ish). Jika Anda tidak menghitung f(n)=bagian, itu adalah 82 byte. Anda juga dapat memulai dengan n->(85 byte).

Ini adalah bahasa 1-diindeks.


Sunting: Modifikasi 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++))akan mencetak ilustrasi jalan semut (diberi terminal yang cukup lebar). Misalnya illustrate(150,25)akan memberikan 150 langkah pertama pada 25 kolom, seperti ini:

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


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


Cobalah online!

Martin Ender menyimpan 15 byte


Anda memiliki ruang nyasar di depan Divisors, Anda dapat menggunakan notasi infiks untuk Do, dan Anda bisa kembali tsebagai gantinya t-1(hasil berbasis 1).
Martin Ender

2

Python 3 , 158 149 133 byte

Ini adalah implementasi prosedural langsung dengan satu atau dua kebiasaan untuk memastikan kode bekerja untuk semua kasus uji. Saya menggunakan [*range(2,n+9)]untuk memastikan bahwa A cukup besar (kecuali untuk n<3, n+9lebih dari cukup). The elseklausul membagi tua A[p]oleh d, decrements p, dan kemudian menambahkan dke baru A[p], yang pasti praktek coding yang buruk. Kalau tidak, cukup mudah. Selamat datang saran bermain golf!

Sunting: -9 byte tanpa sympyterima kasih kepada Halvard Hummel. -14 byte dari Felipe Nardi Batista, -6 byte dari beberapa isyarat dari jawaban Jonathan Frech's Python 2

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)

Cobalah online!



148 byte dengan menjadikannya program penuh
Felipe Nardi Batista

if d-m:A[p]...dan else:p+=1untuk menyimpan byte
Felipe Nardi Batista

143 byte dengan menghapus elsepernyataan
Felipe Nardi Batista

setelah menghapus elsepernyataan tersebut, tidak ada perbedaan dalam byte ke versi fungsi
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;

Jalankan sebagai pipa dengan -Ratau coba online .

Output kosong untuk input 0; masukkan +setelah echountuk literal0

atau gunakan versi 1-diindeks ini (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;

2

R , 123 byte

Implementasi yang mudah. Ini disediakan sebagai fungsi, yang mengambil jumlah gerakan sebagai input dan mengembalikan posisi p.

Itu loop di atas urutan dan bergerak penunjuk maju dan mundur sesuai dengan aturan. Outputnya berbasis 0.

Catatan: untuk menemukan faktor prima terkecil dari bilangan x, ia menghitung modulus x relatif terhadap semua bilangan bulat dari 0 hingga x. Kemudian ekstrak angka-angka dengan modulus sama dengan 0, yang selalu [0,1, ..., x]. Jika angka ketiga tersebut bukan x, maka itu adalah faktor prima terkecil dari 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}

Cobalah online!


2

C (gcc), 152 148 byte

Diperkecil

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;}

Dibentuk dengan beberapa komentar

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;
}

Fungsi utama untuk pengujian

#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;
}

Untuk menunjukkan setiap langkah

  1. Deklarasikan tampilan () di dalam 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");
      }
      ...
  2. Tampilan panggilan ()

      A[i] = i + 2;
    display();
  3. Tampilan panggilan ()

      }
      display();
    }

Anda dapat memotong beberapa byte dengan mendeklarasikan A sebagai array dan menginisialisasi kontrol loop Anda sebelum loop jika memungkinkan, kan?
Pasang kembali Monica

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

Aduh, mengedit "keadaan" tidak ideal di Clojure. Anda harus meningkatkan eksponen untuk input yang lebih besar.


Mengapa Anda menggunakan pencocokan pola di loop? Anda harus dapat kehilangan beberapa byte tanpa itu.
clismique

Juga, Anda mungkin bisa mengubah firsthal itu menjadi somepernyataan.
clismique

Tanpa pencocokan pola, saya harus mengulangi recurdua kali, satu untuk setiap if-letcabang. Juga (dec i)akan diduplikasi. somebutuh predikat, saya bisa menggunakan +karena kita berurusan dengan angka tetapi ini adalah satu karakter lebih lama dari 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;}

Penjelasan:

Coba di sini.

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

1

Clojure, 198 193 191 byte

Ini harus sangat golf ...

#(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 : Disimpan 5 byte dengan mengubah (first(filter ...))ke(some ...)

Golf 2 : Disimpan 2 byte dengan mengubah (zero? ...)ke(= ... 0)


Pemakaian:

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

Kode tidak dikunci:

(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))))))))
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.