Quadrat special primes: Difference between revisions

Content added Content deleted

Inline

@@ Line 108: / Line 108: @@
 11243  11279  11423  12323  12647  12791  13367  13691  14591  14627  14771
+</pre>
+=={{header|AWK}}==
+<lang AWK>
+# syntax: GAWK -f QUADRAT_SPECIAL_PRIMES.AWK
+# converted from FreeBASIC
+BEGIN {
+    stop = 15999
+    p = 2
+    j = 1
+    printf("%5d ",p)
+    count++
+    while (1) {
+      while (1) {
+        if (is_prime(p+j*j)) { break }
+        j++
+      }
+      p += j*j
+      if (p > stop) { break }
+      printf("%5d%1s",p,++count%10?"":"\n")
+      j = 1
+    }
+    printf("\nQuadrat special primes 1-%d: %d\n",stop,count)
+    exit(0)
+}
+function is_prime(x,  i) {
+    if (x <= 1) {
+      return(0)
+    }
+    for (i=2; i<=int(sqrt(x)); i++) {
+      if (x % i == 0) {
+        return(0)
+      }
+    }
+    return(1)
+}
+</lang>
+{{out}}
+<pre>
+3     7    11    47    83   227   263   587   911
+983  1019  1163  1307  1451  1487  1523  1559  2459
+4259  4583  5483  5519  5843  5879  6203  6779  7103
+7283  7607  7643  8219  8363 10667 11243 11279 11423
+12647 12791 13367 13691 14591 14627 14771 15671
+Quadrat special primes 1-15999: 49
 </pre>