File: primes.sh
   1 #!/bin/sh
   2 
   3 # The MIT License (MIT)
   4 #
   5 # Copyright © 2020-2025 pacman64
   6 #
   7 # Permission is hereby granted, free of charge, to any person obtaining a copy
   8 # of this software and associated documentation files (the “Software”), to deal
   9 # in the Software without restriction, including without limitation the rights
  10 # to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
  11 # copies of the Software, and to permit persons to whom the Software is
  12 # furnished to do so, subject to the following conditions:
  13 #
  14 # The above copyright notice and this permission notice shall be included in
  15 # all copies or substantial portions of the Software.
  16 #
  17 # THE SOFTWARE IS PROVIDED “AS IS”, WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
  18 # IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  19 # FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
  20 # AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
  21 # LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
  22 # OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
  23 # SOFTWARE.
  24 
  25 
  26 # primes [options...] [count...]
  27 #
  28 # Find the first n prime numbers, emitting each one on its own line, with a
  29 # default of 100 thousand primes when not given a count.
  30 #
  31 # The help option is `-h`, `--h`, `-help`, or `--help`.
  32 
  33 
  34 # handle leading options
  35 case "$1" in
  36     -h|--h|-help|--help)
  37         awk '/^# +primes /, /^$/ { gsub(/^# ?/, ""); print }' "$0"
  38         exit 0
  39     ;;
  40 esac
  41 
  42 awk '
  43     BEGIN {
  44         left = 100 * 1000
  45         if (ARGV[1] ~ /^0+/ || ARGV[1] + 0 != 0) {
  46             left = ARGV[1] + 0
  47             delete ARGV[1]
  48         }
  49 
  50         if (left < 0) left = 0
  51 
  52         printf "showing first %d prime numbers\n", left > "/dev/stderr"
  53 
  54         if (left > 0) {
  55             print 2
  56             left--
  57         }
  58 
  59         for (n = 3; left > 0; n += 2) {
  60             if (positive_odd_is_prime(n)) {
  61                 print n
  62                 left--
  63             }
  64         }
  65 
  66         exit
  67     }
  68 
  69     function positive_odd_is_prime(n, max, div) {
  70         max = sqrt(n)
  71         for (div = 3; div <= max; div += 2) {
  72             if (n % div == 0) {
  73                 return 0
  74             }
  75         }
  76         return 1
  77     }
  78 ' "$@"