Quadrat special primes: Difference between revisions
Added Wren
(Added Wren) |
|||
Line 349:
Found 49 of the smallest primes < 16,000 such that the difference of successive terma are the smallest quadrat numbers
done...
</pre>
=={{header|Wren}}==
{{libheader|Wren-math}}
{{libheader|Wren-fmt}}
<lang ecmascript>import "/math" for Int
import "/fmt" for Fmt
var isSquare = Fn.new { |n|
var s = n.sqrt.floor
return s*s == n
}
var primes = Int.primeSieve(15999)
System.print("Quadrat special primes under 16,000:")
System.print(" Prime1 Prime2 Gap Sqrt")
var lastQuadSpecial = 3
var gap = 1
var count = 1
Fmt.print("$,7d $,7d $,6d $4d", 2, 3, 1, 1)
for (p in primes.skip(2)) {
gap = p - lastQuadSpecial
if (isSquare.call(gap)) {
Fmt.print("$,7d $,7d $,6d $4d", lastQuadSpecial, p, gap, gap.sqrt)
lastQuadSpecial = p
count = count + 1
}
}
System.print("\n%(count+1) such primes found.")</lang>
{{out}}
<pre>
Quadrat special primes under 16,000:
Prime1 Prime2 Gap Sqrt
2 3 1 1
3 7 4 2
7 11 4 2
11 47 36 6
47 83 36 6
83 227 144 12
227 263 36 6
263 587 324 18
587 911 324 18
911 947 36 6
947 983 36 6
983 1,019 36 6
1,019 1,163 144 12
1,163 1,307 144 12
1,307 1,451 144 12
1,451 1,487 36 6
1,487 1,523 36 6
1,523 1,559 36 6
1,559 2,459 900 30
2,459 3,359 900 30
3,359 4,259 900 30
4,259 4,583 324 18
4,583 5,483 900 30
5,483 5,519 36 6
5,519 5,843 324 18
5,843 5,879 36 6
5,879 6,203 324 18
6,203 6,779 576 24
6,779 7,103 324 18
7,103 7,247 144 12
7,247 7,283 36 6
7,283 7,607 324 18
7,607 7,643 36 6
7,643 8,219 576 24
8,219 8,363 144 12
8,363 10,667 2,304 48
10,667 11,243 576 24
11,243 11,279 36 6
11,279 11,423 144 12
11,423 12,323 900 30
12,323 12,647 324 18
12,647 12,791 144 12
12,791 13,367 576 24
13,367 13,691 324 18
13,691 14,591 900 30
14,591 14,627 36 6
14,627 14,771 144 12
14,771 15,671 900 30
49 such primes found.
</pre>
| |||