|
%I M0017 #32 Jul 04 2024 19:59:59
%S 0,0,1,2,0,1,2,2,1,2,4,3,4,4,5,5,4,5,8,6,8,7,8,8,9,9,10,10,15,11,12,
%T 12,11,12,16,13,16,14,15,15,17,16,17,17,19,18,19,19,20,20,21,21,23,22,
%U 23,23,25,24,25,25,27,26,27,27,26,27,31,28,32,29,30,30,35,31,32,32,34,33,34,34,36,35,36,36,38,37,38,38
%N 2-part of number of graphs on n nodes.
%D N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
%H Chai Wah Wu, <a href="/A007149/b007149.txt">Table of n, a(n) for n = 0..97</a>
%H Steven C. Cater and Robert W. Robinson, <a href="http://cobweb.cs.uga.edu/~rwr/publications/techtournament.pdf">Exponents of 2 in the numbers of unlabeled graphs and tournaments</a>, Congressus Numerantium, 82 (1991), pp. 139-155.
%H Steven C. Cater and Robert W. Robinson, <a href="/A007149/a007149.pdf">Exponents of 2 in the numbers of unlabeled graphs and tournaments</a>, Preprint. (Annotated scanned copy)
%F a(n) = A007814(A000088(n)). - _Michel Marcus_, Jan 06 2020
%t A000088 = Cases[Import["https://oeis.org/A000088/b000088.txt", "Table"], {_, _}][[All, 2]];
%t IntegerExponent[#, 2]& /@ A000088 (* _Jean-François Alcover_, Jan 06 2020 *)
%o (Python)
%o from itertools import combinations
%o from math import prod, factorial, gcd
%o from fractions import Fraction
%o from sympy.utilities.iterables import partitions
%o def A007149(n): return (~(m:=int(sum(Fraction(1<<sum(p[r]*p[s]*gcd(r,s) for r,s in combinations(p.keys(),2))+sum((q>>1)*r+(q*r*(r-1)>>1) for q, r in p.items()),prod(q**r*factorial(r) for q, r in p.items())) for p in partitions(n)))) & m-1).bit_length() # _Chai Wah Wu_, Jul 02 2024
%Y Power of 2 dividing A000088. Cf. A007814.
%K nonn
%O 0,4
%A _N. J. A. Sloane_
%E More terms from _Alois P. Heinz_, Aug 15 2019
|
