library(sets)

a=set(“x”,“y”,“z”)
b=set(“y”,“q”)

axa=set_cartesian(a,a)
p1=set(“x”,“y”)

p_complement=axa-p1