let add c p1 p2 =
    PTbl.set c p1 (PSet.add p2 (PTbl.get c p1));
    PTbl.set c p2 (PSet.add p1 (PTbl.get c p2))