I wrote a program to produce the exact cover for the puzzle and I did some experiments with my exact cover solver to find solutions and generate an estimate of the number of solutions. I am not very sure on the estimate. One run returned a number of 10^50 and another one (with a slightly different, but more successful search strategy) returned a number of 10^148.
The generated exact cover is about 5 megabytes using a notation the positions for the ones in the vector are listed.
The second run has now found about 200 solutions in about six hours, which is far less than the ten thousands that they have found [2].
I do not know of any other programs that can estimate the number of solutions of an exact cover.
[1] https://www.iwriteiam.nl/D2510.html#16
[2] https://news.ycombinator.com/item?id=45559402