permcount[v_] := Module[{m = 1, s = 0, k = 0, t}, For[i = 1, i <= Length[v], i++, t = v[[i]]; k = If[i > 1 && t == v[[i - 1]], k + 1, 1]; m *= t*k; s += t]; s!/m];
edges[v_, k_] := Product[SeriesCoefficient[Product[g = GCD[v[[i]], v[[j]]]; (1 + x^(v[[j]]/g) + O[x]^(k + 1))^g, {j, 1, Length[v]}], {x, 0, k}], {i, 1, Length[v]}];
a[n_] := Module[{s = 0}, Do[s += permcount[p]*edges[p, 3], {p, IntegerPartitions[n]}]; s/n!];
