psh-prng/ran

77 lines
2.8 KiB
Plaintext
Raw Normal View History

2021-10-29 17:59:05 +00:00
#!/bin/sh
[ -c /dev/urandom ] && {
read -r seed < /dev/urandom # read a single line of urandom
seed=${#seed}
} || { # otherwise generate a seed by counting stuff
n=0;
[ -x /bin/ ] && { # if we can access /bin count the files in it
for i in /bin/*; do
[ -f "$i" ] && : $((n+=1))
done
}
[ -x /tmp/ ] && {
for i in /tmp/*; do
[ -e "$i" ] && : $((n+=1))
done
}
[ -x /var/ ] && {
for i in /var/*; do
[ -e "$i" ] && : $((n+=1))
done
}
: $((n+=${#PWD})); : $((n+=${#0}))
seed="${n:-100}"
}
[ "$1" ] && seed="$1" # allow a seed to be set via $1
2021-10-29 17:59:05 +00:00
# $seed should now contain a number
# the idea here is the $seed should be /slightly/ random or at least dicated by the contents of the current systen
# now we must remove all possible leading 0's
until [ ${seed#0} = $seed ]; do
seed=${seed#0}
2021-10-29 17:59:05 +00:00
done
# echo "being LCG"
2021-10-29 17:59:05 +00:00
# LCG imp follows
# generate a modulus (m) from the length of $seed and then detemine if its a prime
m="$(( ${#seed} * 10))"; [ "$((m%2))" -eq 0 -o "$((m%3))" -eq 0 -a "$m" -gt 2 -a "$m" -gt 3 ] || : $((m+=1)) # if prime
# ^ if not a multiple of 2 or 3 expect the number to be prime and increase it by 1
2021-10-29 17:59:05 +00:00
p=0 # $p is the next prime and can be determined using a while loop that adds 1 to it until it becomes a prime
#echo "determine modulus and its factors"
if [ $((m%2)) -eq 0 ]; then # now determine the prime factors of $m
2021-10-29 17:59:05 +00:00
f=$((m/2)) # next prime is 3; 5
p=2
until [ $((f % 2)) -gt 0 -a $(( f % 3)) -gt 0 ]; do
[ $((f%p)) -eq 0 ] && {
f=$((f/p))
} || {
: $((p+=1))
until [ $((p%2)) -gt 0 -a $((p%3)) -gt 0 ] || [ "$p" -eq 2 -o "$p" -eq 3 ]; do
: $((p+=1))
done
}
2021-10-29 17:59:05 +00:00
done
elif [ $((m%3)) -eq 0 ]; then
2021-10-29 17:59:05 +00:00
f=$((m/3)) # next prime is 5; 7
p=3
until [ $((f % 2)) -gt 0 -a $((f % 3)) -gt 0 ]; do
[ $((f%p)) -eq 0 ] && {
f=$((f/p))
} || {
: $((p+=1)); until [ $((p%2)) -gt 0 -a $((p%3)) -gt 0 ] || [ "$p" -eq 2 -o "$p" -eq 3 ]; do
: $((p+=1))
done
}
2021-10-29 17:59:05 +00:00
done # until $f is prime do the above
fi # $f is the prime factor of $m; $f*2 = $m
if [ "$seed" -gt 75 ]; then
a="$((seed-75))" # here the multiplier (a) is seed (z) - 75 # the 75 here is taken from the ZX81's multiplier (a) from its own LCG
elif [ "$seed" -gt 6 ]; then # the 6 here for fallback is literally random
a="$((seed-6))"
fi
until [ "$(( (a-1)%f ))" -eq 0 ]; do # the multiplier (a) -1 MUST be a multiple of $f
: $((a+=1)) # add to a until a-1 is a multiple of $f
done
# increment (c) and the modulus (m) must be co-prime; and two random primes should be co-prime therefore set $c to $p; see above
# the final digit used by LCG will be the seed (z) # z & c will not be defined to save on memory
seed=$(( ((a*seed)+p)%m )) # finally calculate a need seed via the formula: (a*z + c)%m and return it
echo "$seed"