GNU/R: How to compute blocks of binary relations?

Question

Let S be a binary relation on a set K (say). A \'block\' B of S is a maximal subset (of K) satisfying the property B^2 \\subset S. Problem: Find all blocks of a given relation.<