P = Cases[Import["https://oeis.org/
A001035/b001035.txt", "Table"], {_, _}][[All, 2]];
a[n_] := Sum[P[[k+1]] Sum[Binomial[n, s] StirlingS2[n-s, k-s], {s, 0, k}], {k, 0, n}];
transitive[r_]:=Catch[Do[If[a[[2]]==b[[1]]&&FreeQ[r, {a[[1]], b[[2]]}], Throw[False]], {a, r}, {b, r}]; True];