A007669 Duplicate of A034343.
1, 2, 4, 8, 16, 36, 80, 194, 506, 1449, 4631, 17106, 74820, 404283, 2815595
Offset: 1
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.
Triangle T(n,k) (with rows n >= 1 and columns k >= 1) begins as follows: 1; 1, 2; 1, 3, 4; 1, 4, 7, 8; 1, 5, 11, 15, 16; 1, 7, 19, 30, 35, 36; 1, 8, 29, 56, 73, 79, 80; 1, 10, 44, 107, 161, 186, 193, 194; ...
# We illustrate how to get a g.f. for column k in Maple when k is small. with(GroupTheory); G := ProjectiveGeneralLinearGroup(4, 2); GroupOrder(G); # We get that the order is 20160. f:=CycleIndexPolynomial(G, [x||(1..20160)]); # We get # 1/20160*x1^15 + 1/192*x1^7*x2^4 + 1/96*x1^3*x2^6 + 1/16*x1^3*x2^2*x4^2 + # 1/18*x1^3*x3^4 + 1/6*x1*x2*x3^2*x6 + 1/8*x1*x2*x4^3 + 1/180*x3^5 + 2/7*x1*x7^2 + # 1/12*x3*x6^2 + 1/15*x5^3 + 2/15*x15 # The only dummy variables that appear are x1, x2, x3, x4, x5, x6, x7, and x15. g:=subs(x1 = 1/(1 - y), subs(x2 = 1/(-y^2 + 1), subs(x3 = 1/(-y^3 + 1), subs(x4 = 1/(-y^4 + 1), subs(x5 = 1/(-y^5 + 1), subs(x6 = 1/(-y^6 + 1), subs(x7 = 1/(-y^7 + 1), subs(x15 = 1/(-y^15 + 1), f)))))))); # Then we take a Taylor expansion of the above g.f. taylor(g,y=0,50); # We get a Taylor expansion for column k = 4 (i.e., A034338). # Petros Hadjicostas, Sep 30 2019
# Fripertinger's method to find the g.f. of column k for small k: def A076832col(k, length): G = PSL(k, GF(2)) D = G.cycle_index() f = sum(i[1]*prod(1/(1-x^j) for j in i[0]) for i in D) return f.taylor(x, 0, length).list() # For instance the Taylor expansion for column k = 4 gives A034338: print(A076832col(4, 30)) # Petros Hadjicostas, Sep 30 2019
Comments