This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A134018 #15 Dec 01 2017 09:04:51 %S A134018 0,1,3,10,45,226,1113,5230,23565,102826,438273,1836550,7601685, %T A134018 31183426,127084233,515429470,2083077405,8396552026,33779262993, %U A134018 135696871990,544528258725,2183337968626,8749031918553,35043178292110,140313885993645,561679104393226 %N A134018 Let P(A) be the power set of an n-element set A. Then a(n) = the number of pairs of elements {x,y} of P(A) for which either 0) x and y are disjoint and for which either x is a subset of y or y is a subset of x, or 1) x and y are intersecting but for which x is not a subset of y and y is not a subset of x. %H A134018 Ross La Haye, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL12/LaHaye/lahaye5.html">Binary Relations on the Power Set of an n-Element Set</a>, Journal of Integer Sequences, Vol. 12 (2009), Article 09.2.6. [_Ross La Haye_, Feb 22 2009] %H A134018 <a href="/index/Rec#order_04">Index entries for linear recurrences with constant coefficients</a>, signature (10,-35,50,-24). %F A134018 a(n) = (1/2)(4^n - 3^(n+1) + 5*2^n - 3) = 3*StirlingS2(n+1,4) + StirlingS2(n+1,2). %F A134018 G.f.: x*(1-7*x+15*x^2)/((1-x)*(1-2*x)*(1-3*x)*(1-4*x)). [_Colin Barker_, Jul 29 2012] %e A134018 a(3) = 10 because for P(A) = {{},{1},{2},{3},{1,2},{1,3},{2,3},{1,2,3}} we have for case 0 {{},{1}}, {{},{2}}, {{},{3}}, {{},{1,2}}, {{},{1,3}}, {{},{2,3}}, {{},{1,2,3}} and we have for case 1 {{1,2},{1,3}}, {{1,2},{2,3}}, {{1,3},{2,3}}. %t A134018 LinearRecurrence[{10,-35,50,-24},{0,1,3,10},30] (* _Harvey P. Dale_, Dec 01 2017 *) %Y A134018 Cf. A000225, A032263. %K A134018 nonn,easy %O A134018 0,3 %A A134018 _Ross La Haye_, Jan 10 2008 %E A134018 More terms from _Harvey P. Dale_, Dec 01 2017