let strongset graph q =
  if G.mem_vertex graph q then
    G.fold_succ (fun p acc -> S.add p acc) graph q S.empty
  else S.empty