experimental-projects/oeis-research/Paolo P. Lava/Squarefree difference of semiprimes/filter_primes.sf

#!/usr/bin/ruby

var n = 7136306
var P = primes(20000)

say (n, " -> ", P.grep {|p|
    var t = (p * (p + n))
    t.is_semiprime && !t.is_square
}.grep {|p|

    var w = (p * (n + p))

    P.any {|q|
        var t = (q * (q + w))
        t.is_semiprime && !t.is_square
    }
})