Polyanalytic boundary value problems for planar domains with harmonic Green function