login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A337518
Number of non-isomorphic graphs on n unlabeled nodes modulo 3.
1
1, 1, 2, 1, 2, 1, 0, 0, 1, 0, 2, 2, 0, 0, 1, 2, 0, 1, 0, 2, 1, 2, 0, 0, 2, 0, 0, 1, 1, 1, 0, 0, 2, 0, 1, 0, 2, 1, 0, 0, 2, 2, 0, 2, 0, 0, 2, 1, 2, 1, 0, 1, 0, 1, 1, 1, 0, 1, 2, 2, 1, 0, 1, 2, 0, 1, 2, 1, 1, 1, 1, 0, 0, 2, 1, 1, 2, 0, 2, 0, 0, 0, 2, 1, 2, 2, 1
OFFSET
0,3
COMMENTS
For the mod-2 case, the sequence is eventually constant, because there are an even number of graphs on n vertices for n>4. (In fact, the number of factors of 2 in A000088(n) is asymptotically n/2; see Cater and Robinson in the Links section.)
LINKS
Steven C. Cater and Robert W. Robinson, Exponents of 2 in the numbers of unlabeled graphs and tournaments, Congressus Numerantium, 82 (1991), pp. 139-155.
FORMULA
a(n) = A000088(n) mod 3.
EXAMPLE
For n = 4, there are 11 graphs on 4 nodes up to isomorphism, so a(4) = 2 = 11 mod 3.
PROG
(Python)
from itertools import combinations
from math import prod, factorial, gcd
from fractions import Fraction
from sympy.utilities.iterables import partitions
def A337518(n): return 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()))%3 for p in partitions(n))) % 3 # Chai Wah Wu, Jul 02 2024
CROSSREFS
Sequence in context: A024356 A143947 A226518 * A073781 A048622 A105661
KEYWORD
nonn
AUTHOR
Drake Thomas, Nov 21 2020
STATUS
approved