Add XS routine for ([p1,p2,p3],[e1,e2,e3]) return from factor
Lucas sequence should return a bigint if arg is a bigint