List Commands Library for 50g

[John Keith]

There are number theory uses as well. I\->BL works like the Mathematica command IntegerDigits, which is very useful. I just used I\->BL in a "slow and dirty" program for A000966 based on the first Mathematica program in the linked article.

Edit: shorter and faster thanks to the 2nd formula in A027868

Code:


\<< LSEQ 1.
   \<<  DUP 5 I\->BL LSUM - 4 / 
   \>> DOLIST LGRP DUP 0 SWAP Last LSEQR SWAP Diff
\>>

Given an integer N as input, this program returns a list of all the missing numbers up to N!. For an input of 1000, it returns a list of the first 49 terms in 53s (approx. mode) or 72s (exact mode). Requires both ListExt and GoferLists.