รับจำนวนธรรมชาติ , กลับ -th นายกคิวบา $n$ $n$

Cubimes Primes

Cuban prime เป็นจำนวนเฉพาะของแบบฟอร์ม

พี = \frac{x^{3} - Y^{3}}{x - Y}

$p = \frac{x^3-y^3}{x-y}$

โดยที่และหรือ $y>0$ $x = 1+y$ $x = 2+y$

รายละเอียด

คุณอาจใช้การจัดทำดัชนีตาม 0 หรือ 1 สิ่งที่เหมาะกับคุณที่สุด
คุณสามารถกลับนายก -th ให้ดัชนีหรือแรกช่วงเวลาในการสั่งซื้อที่เพิ่มขึ้นหรือมิฉะนั้นคุณสามารถกลับรายการอนันต์ / เครื่องกำเนิดไฟฟ้าที่ผลิตเฉพาะในการสั่งซื้อที่เพิ่มขึ้น $n$ $n$ $n$

กรณีทดสอบ

คำศัพท์สองสามคำแรกมีดังต่อไปนี้:

(#1-13)   7, 13, 19, 37, 61, 109, 127, 193, 271, 331, 397, 433, 547,
(#14-24) 631, 769, 919, 1201, 1453, 1657, 1801, 1951, 2029, 2269, 2437,
(#25-34) 2791, 3169, 3469, 3571, 3889, 4219, 4447, 4801, 5167, 5419,
(#35-43) 6211, 7057, 7351, 8269, 9241, 10093, 10267, 11719, 12097,
(#44-52) 12289, 13267, 13669, 13873, 16651, 18253, 19441, 19927, 20173

คำศัพท์เพิ่มเติมสามารถพบได้ใน OEIS: พวกเขาจะแบ่งออกเป็นสองลำดับขึ้นอยู่กับ wheterหรือ : A002407และA002648 $x = 1+y$ $x = 2+y$

— flawr
แหล่งที่มา

2

เราสามารถคืนค่า n จำนวนแรกที่ไม่ได้จัดเรียงได้หรือไม่?

— J42161217

@ J42161217 ไม่สิ่งที่ควรจะอยู่ในลำดับที่เพิ่มขึ้น

— ข้อบกพร่อง

23

JavaScript (V8) , 54 ไบต์

โปรแกรมเต็มรูปแบบที่พิมพ์ช่วงเวลาคิวบาตลอดไป

for(x=0;;){for(k=N=~(3/4*++x*x);N%++k;);~k||print(-N)}

ลองออนไลน์!

^{หมายเหตุ: หากคุณไม่มีกระดาษที่ไม่มีที่สิ้นสุดในเครื่องพิมพ์ของคุณอย่าพยายามเรียกใช้สิ่งนี้ในคอนโซลเบราว์เซอร์ของคุณซึ่งprint()อาจมีความหมายแตกต่างกัน}

JavaScript (ES6), 63 61 60 59 ไบต์

ส่งคืนไพรม์ Cuban -th ที่สร้างดัชนี 1 ตำแหน่ง $n$

f=(n,x)=>(p=k=>N%++k?p(k):n-=!~k)(N=~(3/4*x*x))?f(n,-~x):-N

ลองออนไลน์!

อย่างไร?

สิ่งนี้มีพื้นฐานมาจากข้อเท็จจริงที่ว่าจำนวนเฉพาะของคิวบาเป็นช่วงเวลาของแบบฟอร์ม:

{พี}_{n} = ⌊ \frac{3 n^{2}}{4} ⌋ + 1, n \geq 3

$p_n=\left\lfloor\frac{3n^2}{4}\right\rfloor+1,\;n\ge3$

สูตรข้างต้นสามารถเขียนเป็น:

{พี}_{n} = {\begin{cases} \frac{3 n^{2} + 1}{4} ถ้า n แปลก \\ \frac{3 n^{2} + 4}{4} ถ้า n เป็นคู่ \end{cases}

$p_n=\begin{cases} \dfrac{3n^2+1}{4}\;\text{ if }n\text{ is odd}\\ \dfrac{3n^2+4}{4}\;\text{ if }n\text{ is even} \end{cases}$

หรือสำหรับใด ๆ : $y>0$

{พี}_{2 Y + 1} = \frac{3 (2 Y + 1)^{2} + 1}{4} = 3 Y^{2} + 3 Y + 1

$p_{2y+1}=\dfrac{3(2y+1)^2+1}{4}=3y^2+3y+1$

p_{2 y + 2} = \frac{3 (2 y + 2)^{2} + 4}{4} = 3 Y^{2} + 6 Y + 4

$p_{2y+2}=\dfrac{3(2y+2)^2+4}{4}=3y^2+6y+4$

ซึ่งคือสำหรับและตามลำดับ $\dfrac{x^3-y^3}{x-y}$ $x=y+1$ $x=y+2$

— Arnauld
แหล่งที่มา

7

05AB1E , 16 12 9 ไบต์

สร้างรายการที่ไม่มีที่สิ้นสุด
ที่บันทึกไว้ 4 ไบต์กับเควิน Cruijssen ของ พอร์ตของสูตร Arnaulds
บันทึกอีก 3 ไบต์ด้วยGrimy

∞n3*4÷>ʒp

ลองออนไลน์!

คำอธิบาย

∞          # on the list of infinite positive integers
 n3*4÷>    # calculate (3*N^2)//4+1 for each
       ʒp  # and filter to only keep primes

— Emigna
แหล่งที่มา

คุณได้พิมพ์ผิดในคำอธิบายของคุณ: " ใส่สำเนาของN^2+3ในกอง " 3*N^2ควรจะเป็น นอกจากนี้เหตุผลที่)แทน¯? เพราะพิมพ์ง่ายกว่า? และด้วยเหตุผลบางอย่างฉันรู้สึกว่าNnN‚3*¬sO‚สามารถมีขนาดสั้นลง 1 ไบต์ แต่ฉันไม่เห็น Nn3*DN3*+‚เล็กน้อยทางเลือกเท่ากับไบต์คือ แต่ฉันอาจเห็นสิ่งที่ไม่มี .. ;) คำตอบที่ดีดังนั้น +1 จากฉัน

— Kevin Cruijssen

1

จริง ๆ แล้วฉันพยายามที่จะพอร์ตคำตอบของฉันเพื่อ 05AB1E แต่ล้มเหลวอย่างน่าสังเวช : D

— Arnauld

1

ที่จริงแล้วการสร้างรายการไม่มีที่สิ้นสุดสะดวกกว่า: 9 ไบต์ด้วยwithn3 * 4 ÷> ʒp

— Grimmy

1

ตกลงฉันไม่คุ้นเคยกับสเป็คที่ขัดแย้งกับตัวเอง :-)

— WGroleau

6

@WGroleau ฉันคิดว่าคุณไม่เคยพัฒนาซอฟต์แวร์อย่างมืออาชีพมาก่อน ฉันกังวลมากขึ้นเมื่อฉันได้รับสเป็คที่ไม่ขัดแย้งกับตัวเอง

— MikeTheLiar

7

R , 75 73 ไบต์

n=scan()
while(F<n)F=F+any(!(((T<-T+1)*1:4-1)/3)^.5%%1)*all(T%%(3:T-1))
T

ลองออนไลน์!

-2 ไบต์โดยสังเกตว่าฉันสามารถลบวงเล็บถ้าฉันใช้*แทน&(ลำดับความสำคัญที่แตกต่างกัน)

แสดงผลnคิวบาที่สำคัญ (ดัชนี 1)

มันใช้ความจริง (ที่ให้ไว้ใน OEIS) ว่าช่วงเวลาคิวบามีรูปแบบ $p=1+3n^2$ หรือ $4p=1+3n^2$ สำหรับบาง $n$ , เช่น $n=\sqrt{\frac{a\cdot p-1}{3}}$ เป็นจำนวนเต็มสำหรับหรือ4 $a=1$ $a=4$

เคล็ดลับคือว่าไม่มีนายกสามารถของรูปแบบที่ $2p=1+3n^2$ หรือ $3p=1+3n^2$ (*) เพื่อให้เราสามารถประหยัด 2 ไบต์โดยการตรวจสอบสูตรสำหรับ ( ) แทนที่จะ ( ) $a\in\{1, 2, 3, 4\}$ 1:4 $a\in\{1, 4\}$ c(1,4)

โค้ดเวอร์ชันที่ไม่อัปเดตเล็กน้อย:

# F and T are implicitly initialized at 0 and 1
# F is number of Cuban primes found so far
# T is number currently being tested for being a Cuban prime
n = scan()                       # input
while(F<n){
  T = T+1                        # increment T 
  F = F +                        # increment F if
    (!all(((T*1:4-1)/3)^.5 %% 1) # there is an integer of the form sqrt(((T*a)-1)/3)
     & all(T%%(3:T-1)))          # and T is prime (not divisible by any number between 2 and T-1)
  }
T                                # output T

(*) ไม่สำคัญอาจเป็นในรูปแบบ $3p=1+3n^2$ อื่น $1=3(p-n^2)$ จะหารด้วย3 $3$

ไม่มีที่สำคัญอื่น ๆ กว่า $p=2$ (ซึ่งไม่ได้เป็นนายกคิวบา) สามารถของรูปแบบที่ $2p=1+3n^2$ : $n$ จะต้องแปลกคือ $n=2k+1$ 1การขยายให้ $2p=4+12k(k+1)$ ดังนั้น $p=2+6k(k+1)$ และ $p$ จะเท่ากัน

— Robin Ryder
แหล่งที่มา

สิ่งที่เกี่ยวกับการหลีกเลี่ยงการวนซ้ำโดยใช้ขอบเขตบนบนคิวบาที่สำคัญ?

— ซีอาน

@ ซีอานฉันคิดเกี่ยวกับสิ่งนั้น แต่ไม่สามารถเกิดขึ้นได้ด้วยข้อ จำกัด ดังกล่าว คุณมีหรือไม่

— Robin Ryder

5

ภาษา Wolfram (Mathematica) , 66 65 56 bytes

(f=1+⌊3#/4#⌋&;For[n=i=0,i<#,PrimeQ@f@++n&&i++];f@n)&

ลองออนไลน์!

J42161217 -1โดยใช้⌊ ⌋แทนFloor[ ]
attinat
- -1โดยใช้⌊3#/4#⌋แทน⌊3#^2/4⌋
- -8สำหรับFor[n=i=0,i<#,PrimeQ@f@++n&&i++]แทนn=2;i=#;While[i>0,i-=Boole@PrimeQ@f@++n]

— speedstyle
แหล่งที่มา

1

65 ไบต์ ยินดีต้อนรับสู่ ppcg คำตอบแรกที่ดี! +1

— J42161217

ขอบคุณ! (ผู้แฝงตัวมานาน) ฉันไม่สามารถแยกคำตอบที่มีอยู่ของคุณได้ดังนั้นฉันจึงเขียนของตัวเองและมันก็สั้นลงเล็กน้อย ฉันอาจทำ Python ด้วย

— speedstyle

2

56 bytes

— attinat

@attinat ฉันคิดว่าสูตรของ Arnauld ใช้ได้กับ n> 2 เท่านั้นฉันจึงไม่ได้ขึ้นต้นด้วย 0 - แม้ว่าในตัวอย่างของคุณจะใช้ได้กับ n ทั้งหมด (เพราะเริ่ม 1 1 4 7 13 ... ดังนั้นช่วงเวลานั้นคือ 13 13 .. )

— speedstyle

3

Java 8, 94 88 86 84 ไบต์

v->{for(int i=3,n,x;;System.out.print(x<1?++n+" ":""))for(x=n=i*i++*3/4;~n%x--<0;);}

-6 ไบต์โดยใช้Java prime-checker ของ@SaraJดังนั้นอย่าลืม upvote เธอ!
-2 ไบต์ขอบคุณที่@ OlivierGrégoire ตั้งแต่การตรวจสอบจำนวนครั้งแรกที่เราเป็น7เราสามารถวางต่อท้ายจากซาร่าลากตัวตรวจสอบซึ่งก็คือการบอกเลิกห่วงสำหรับ%n -2 ไบต์ขอบคุณ@ OlivierGrégoireโดยตอบคำตอบของ@Arnauldn=1

ส่งออกที่คั่นด้วยช่องว่างอย่างไม่มีกำหนด

ลองออนไลน์

คำอธิบาย (ของรุ่น 86 ไบต์เก่า): สิ่งที่ต้องทำ: อัปเดตคำอธิบาย

$p_n=\left\lfloor\frac{3n^2}{4}\right\rfloor+1,\;n\ge3$

v->{                     // Method with empty unused parameter and no return-type
  for(int i=3,           //  Loop-integer, starting at 3
          n,x            //  Temp integers
      ;                  //  Loop indefinitely:
      ;                  //    After every iteration:
       System.out.print( //     Print:
        n==x?            //      If `n` equals `x`, which means `n` is a prime:
         n+" "           //       Print `n` with a space delimiter
        :                //      Else:
         ""))            //       Print nothing
    for(n=i*i++*3/4+1,   //   Set `n` to `(3*i^2)//4+1
                         //   (and increase `i` by 1 afterwards with `i++`)
        x=1;             //   Set `x` to 1
        n%++x            //   Loop as long as `n` modulo `x+1`
                         //   (after we've first increased `x` by 1 with `++x`)
             >0;);}      //   is not 0 yet
                         //   (if `n` is equal to `x`, it means it's a prime)

— Kevin Cruijssen
แหล่งที่มา

ฉันไม่คิดว่ามันจะเป็นไปได้ แต่วิธีอื่นในการค้นหา Cuban primes ใช้สูตรนี้: v->{for(int n=7,i=3,p,x,d,r=0;;i+=++r%2*3,n+=i,System.out.print(x>1?x+" ":""))for(x=n,d=1;++d<n;x=x%d<1?0:n);}อาจมีบางคนสามารถใช้สิ่งนี้เพื่อเล่นกอล์ฟ? ฉันทำไม่ได้

— Olivier Grégoire

1

@ OlivierGrégoireคุณสามารถตีกอล์ฟของคุณได้มากขึ้นโดยการลบสิ่งที่ไม่ได้ใช้,pและเปลี่ยนi+=++r%2*3,n+=iไปn+=i+=++r%2*3แต่หลังจากนั้นฉันจะจบที่ 106 ไบต์ โดยใช้ Java 11 String#repeatกับนายก regex 105 v->{for(int n=7,i=3,r=0;;n+=i+=++r%2*3)if(!"x".repeat(n).matches(".?|(..+?)\\1+"))System.out.println(n);}ไบต์:

— Kevin Cruijssen

ใช่ฉันเดาว่ามันไม่สามารถเล่นกอล์ฟได้มากนักแม้จะมีข้อผิดพลาด (ตอนนี้ชัดเจน) ของฉัน ขอบคุณที่ให้การเดินทาง;)

— Olivier Grégoire

@ OlivierGrégoireอาจจะดีสำหรับคุณเช่นกัน แต่ดูเหมือนว่าจะมีการตรวจสอบวนซ้ำที่สั้นกว่าใน Java ดูการแก้ไขและคำตอบเฉพาะของ SaraJ

— Kevin Cruijssen

ฉันอาจจะผิด แต่%nไม่จำเป็นต้องใช้ล่าสุดมันคืออะไร?

— Olivier Grégoire

2

ภาษา Wolfram (Mathematica) , 83 ไบต์

(t=1;While[Length[l=Select[Join@@Array[{(v=3#^2+1)+3#,v}&,t++],PrimeQ]]<#];Sort@l)&

ลองออนไลน์!

— J42161217
แหล่งที่มา

2

เยลลี่ 12 ไบต์

²×3:4‘
ÇẒ$#Ç

ลองออนไลน์!

ขึ้นอยู่กับ@ วิธี ใช้ n ใน stdin และคืนค่าช่วงเวลาคิวบาจำนวนมาก

— นิคเคนเนดี
แหล่งที่มา

1

ภาษา Wolfram (Mathematica) , 83 ไบต์

วิธีการแก้ปัญหานี้จะส่งผลให้นายกรัฐมนตรีคิวบาที่ n ได้รับประโยชน์เพิ่มเติมจากความรวดเร็วและจดจำผลลัพธ์ก่อนหน้านี้ทั้งหมดในสัญลักษณ์ f

(d:=1+3y(c=1+y)+3b c;e:=If[PrimeQ@d,n++;f@n=d];For[n=y=b=0,n<#,e;b=1-b;e,y++];f@#)&

ลองออนไลน์!

— Kelly Lowder
แหล่งที่มา

1

ช่องว่าง 180 ไบต์

[S S S T    S N
_Push_2][S N
S _Duplicate][N
S S N
_Create_Label_OUTER_LOOP][S N
N
_Discard_top_stack][S S S T N
_Push_1][T  S S S _Add][S N
S _Duplicate][S N
S _Duplicate][T S S N
_Multiply][S S S T  T   N
_Push_3][T  S S N
_Multiply][S S S T  S S N
_Push_4][T  S T S _Integer_divide][S S S T  N
_Push_1][T  S S S _Add][S S S T N
_Push_1][S N
S _Duplicate_1][N
S S S N
_Create_Label_INNER_LOOP][S N
N
_Discard_top_stack][S S S T N
_Push_1][T  S S S _Add][S N
S _Duplicate][S N
S _Duplicate][S T   S S T   T   N
_Copy_0-based_3rd][T    S S T   _Subtract][N
T   S T N
_Jump_to_Label_PRINT_if_0][S T  S S T   S N
_Copy_0-based_2nd][S N
T   _Swap_top_two][T    S T T   _Modulo][S N
S _Duplicate][N
T   S S S N
_Jump_to_Label_FALSE_if_0][N
S N
S N
_Jump_to_Label_INNER_LOOP][N
S S T   N
_Create_Label_PRINT][T  N
S T _Print_as_integer][S S S T  S T S N
_Push_10_(newline)][T   N
S S _Print_as_character][S N
S _Duplicate][N
S S S S N
_Create_Label_FALSE][S N
N
_Discard_top_stack][S N
N
_Discard_top_stack][N
S N
N
_Jump_to_Label_OUTER_LOOP]

เพิ่มตัวอักษรS(ช่องว่าง), T(แท็บ) และN(บรรทัดใหม่) เป็นการเน้นเท่านั้น
[..._some_action]เพิ่มเป็นคำอธิบายเท่านั้น

เอาท์พุทคั่นด้วย newline ไปเรื่อย

ลองใช้ออนไลน์ (ด้วยพื้นที่ว่างเปล่าแท็บและบรรทัดใหม่เท่านั้น)

คำอธิบายในรหัสเทียม:

$p_n=\left\lfloor\frac{3n^2}{4}\right\rfloor+1,\;n\ge3$

Integer i = 2
Start OUTER_LOOP:
  i = i + 1
  Integer n = i*i*3//4+1
  Integer x = 1
  Start INNER_LOOP:
    x = x + 1
    If(x == n):
      Call function PRINT
    If(n % x == 0):
      Go to next iteration of OUTER_LOOP
    Go to next iteration of INNER_LOOP

function PRINT:
  Print integer n
  Print character '\n'
  Go to next iteration of OUTER_LOOP

— Kevin Cruijssen
แหล่งที่มา

1

Python 3 , 110 108 102 ไบต์

วิธีที่คล้ายกับคำตอบ Mathematicaของฉัน(เช่นisPrime(1+⌊¾n²⌋) else n++) โดยใช้ตัวตรวจสอบนายกคนนี้ golfedและส่งกลับตัวสร้างแบบไม่ จำกัด ที่ไม่ระบุชื่อ

from itertools import*
(x for x in map(lambda n:1+3*n**2//4,count(2)) if all(x%j for j in range(2,x)))

ลองออนไลน์!

mypetlion -2เนื่องจากเครื่องกำเนิดไฟฟ้านิรนามที่ไม่ระบุชื่อได้รับอนุญาตมากกว่าชื่อที่กำหนด
-6โดยเริ่มต้นcountที่ 2 ₊₁ดังนั้นand x>1ในตัวตรวจสอบหลักที่ฉันยืมมานั้นไม่จำเป็น_-7

— speedstyle
แหล่งที่มา

คำตอบที่เข้าสู่ตัวแปรมักไม่ถือว่าเป็นรูปแบบที่ถูกต้องของ "เอาท์พุท" คุณสามารถปรับปรุงคำตอบของคุณใหม่เพื่อให้ผลลัพธ์เป็นเอาต์พุตไปยัง stdout หรือส่งคืนโดยฟังก์ชันได้หรือไม่?

— mypetlion

1

g=ตั้งแต่ฟังก์ชั่นที่ไม่ระบุชื่อที่ได้รับอนุญาตและความท้าทายชัดเจนให้กำเนิดอนันต์ผมได้ลบ ฉันได้รวมเฉพาะในสถานที่แรกเพราะมันได้รับอนุญาตให้ภาพอย่างรวดเร็วใน TIO print(next(g) for i in range(52))กับ

— speedstyle

1

Japt , 14 13 ไบต์

ที่ดัดแปลงมาจากสูตรของ Arnauld 1 การจัดทำดัชนี

@µXj}f@Ò(X²*¾

ลองมัน

บันทึก 1 ไบต์ด้วย EmbodimentOfIgnorance

— มีขนดก
แหล่งที่มา

13 ไบต์ ไม่ได้ทดสอบอย่างละเอียดว่า

— ศูนย์รวมแห่งความไม่รู้

ขอบคุณ @EmbodimentofIgnorance ฉันลองแล้ว แต่มันใช้งานไม่ได้ (ผลัดกันออกมาผมลืม

— Shaggy

1

แร็กเก็ต , 124 ไบต์

(require math)(define(f n[i 3])(let([t(+(exact-floor(* 3/4 i i))1)][k(+ 1 i)])(if(prime? t)(if(= 0 n)t(f(- n 1)k))(f n k))))

ลองออนไลน์!

ส่งคืนไพรม์คิวบ์ n-th ที่มีค่าดัชนี 0

ใช้สูตรคำตอบ JavaScriptของ@ Arnauld

— Galen Ivanov
แหล่งที่มา

1

Python 3 , 83 ไบต์

พิมพ์ช่วงเวลาคิวบาตลอดไป

P=k=1
while 1:P*=k*k;x=k;k+=1;P%k>0==((x/3)**.5%1)*((x/3+.25)**.5%1-.5)and print(k)

ลองออนไลน์!

$x = 1+y$ $x=2+y$

พี = \frac{(1 + Y)^{3} - Y^{3}}{(1 + Y) - Y} = 1 + 3 Y + 3 Y^{2} \Leftrightarrow Y = - \frac{1}{2} \pm \sqrt{\frac{1}{4} + \frac{พี - 1}{3}}

$p=\frac{(1+y)^3-y^3}{(1+y)-y} = 1 + 3y +3y^2 \Leftrightarrow y = -\frac{1}{2}\pm\sqrt{\frac{1}{4}+\frac{p-1}{3}}$

พี = \frac{(2 + Y)^{3} - Y^{3}}{(1 + Y) - Y} = 4 + 6 Y + 3 Y^{2} \Leftrightarrow Y = - 1 \pm \sqrt{\frac{พี - 1}{3}}

$p=\frac{(2+y)^3-y^3}{(1+y)-y} = 4 + 6y +3y^2 \Leftrightarrow y = -1 \pm\sqrt{\frac{p-1}{3}}$

y

$y$

\pm

$\pm$

- 1

$-1$

— OVS
แหล่งที่มา

1

Perl 6 , 33 31 ไบต์

-2 ไบต์ขอบคุณ Grimy

{grep &is-prime,1+|¾*$++²xx*}

ลองออนไลน์!

บล็อกโค้ดแบบไม่ระบุชื่อที่ส่งคืนรายการอนันต์ช่วงเวลาของคิวบา สิ่งนี้ใช้สูตรของ Arnauldเพื่อสร้างคิวบิกคิวบาที่เป็นไปได้จากนั้น&is-primeกรอง

คำอธิบาย:

{                           }  # Anonymous code block
 grep &is-prime,               # Filter the primes from
                         xx*   # The infinite list
                   ¾*          # Of three quarters
                     $++²      # Of an increasing number squared
                1+|            # Add one by ORing with 1

— ล้อเล่น
แหล่งที่มา

1

1+0+|สามารถเป็นเพียง1+|

— Grimmy

0

Pari / GP , 51 ไบต์

โดยใช้สูตรของ Arnauld

n->a=0;for(i=1,n,until(isprime(p=3*a^2\4+1),a++));p

ลองออนไลน์!

— alephalpha
แหล่งที่มา

0

APL (NARS), 98 ตัวอักษร, 196 ไบต์

r←h w;y;c;v
r←c←y←0⋄→4
→3×⍳∼0πv←1+3×y×1+y+←1⋄r←v⋄→0×⍳w≤c+←1
→2×⍳∼0πv+←3×y+1⋄c+←1⋄r←v
→2×⍳w>c

เยื้อง:

r←h w;y;c;v
r←c←y←0⋄→4
    →3×⍳∼0πv←1+3×y×1+y+←1⋄r←v⋄→0×⍳w≤c+←1
    →2×⍳∼0πv+←3×y+1⋄c+←1⋄r←v
    →2×⍳w>c

ทดสอบ:

  h ¨1..20
7 13 19 37 61 109 127 193 271 331 397 433 547 631 769 919 1201 1453 1657 1801 
  h 1000
25789873
  h 10000
4765143511

มันขึ้นอยู่กับ: ถ้า y ใน N, Cuban Prime ที่เป็นไปได้หนึ่งอันคือ

S1=1+3y(y+1)

Cuban Prime ที่เป็นไปได้ต่อไปจะเป็น

S2=3(y+1)+S1

— RosLuP
แหล่งที่มา