To execute what we have so far, run

$ ./generate-pddl.py N

where N is the number you want to factor. This will produce two files: domain.pddl and problem.pddl.

This code was only tested with python3.10 on Ubuntu 22.04. Moreover ,the PDDL files produced were tested with Fast Downward and Powerlifted -- also note that Fast Downward is very slow with N > 10.