cp's OEIS Frontend

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.

Showing 1-4 of 4 results.

A039623 a(n) = n^2*(n^2+3)/4.

Original entry on oeis.org

1, 7, 27, 76, 175, 351, 637, 1072, 1701, 2575, 3751, 5292, 7267, 9751, 12825, 16576, 21097, 26487, 32851, 40300, 48951, 58927, 70357, 83376, 98125, 114751, 133407, 154252, 177451, 203175, 231601, 262912, 297297, 334951, 376075, 420876, 469567
Offset: 1

Views

Author

Christian Meland (christian.meland(AT)pfi.no)

Keywords

Comments

Previous definition was: Consider a figure like this <> (a squashed square, symmetric about both axes); each side is given 1 of n colors; a(n) = number of possibilities, allowing turning over.
Also number of 2 X 2 matrices with entries mod n, up to row and column permutation. Number of k X l matrices with entries mod n, up to row and column permutation is Z(S_k X S_l; n,n,...) where Z(S_k X S_l; x_1,x_2,...) is cycle index of Cartesian product of symmetric groups S_k and S_l of degree k and l, respectively. - Vladeta Jovovic, Nov 04 2000
Also, if a 2-set Y and a 3-set Z are disjoint subsets of an n-set X then a(n-5) is the number of 6-subsets of X intersecting both Y and Z. - Milan Janjic, Sep 08 2007

Crossrefs

Programs

  • Magma
    [n^2*(n^2+3)/4 : n in [1..50]]; // Wesley Ivan Hurt, Dec 26 2016
  • Maple
    A039623:=n->n^2*(n^2+3)/4: seq(A039623(n), n=1..50); # Wesley Ivan Hurt, Dec 26 2016
  • Mathematica
    Table[(n^2 (n^2+3))/4,{n,40}] (* or *) LinearRecurrence[{5,-10,10,-5,1},{1,7,27,76,175},40] (* Harvey P. Dale, Oct 01 2011 *)
  • PARI
    Vec((-1-2*x-2*x^2-x^3)/(x-1)^5 + O(x^50)) \\ Michel Marcus, Aug 23 2015
    
  • PARI
    a(n) = (1/4)*n^2*(n^2+3); \\ Altug Alkan, Apr 16 2016
    

Formula

From Harvey P. Dale, Oct 01 2011: (Start)
G.f.: (1 + 2*x + 2*x^2 + x^3)/(1 - x)^5.
a(1)=1, a(2)=7, a(3)=27, a(4)=76, a(5)=175; for n>5, a(n) = 5*a(n-1) - 10*a(n-2) + 10*a(n-3) - 5*a(n-4) + a(n-5). (End)
E.g.f.: x*(4 + 10*x + 6*x^2 + x^3)*exp(x)/4. - Ilya Gutkovskiy, Apr 16 2016
a(n) = t(n-1)*t(n) + t(n-1) + t(n) where t=A000217. - J. M. Bergot, Apr 16 2016
a(n) = A000217(n)^2 - n*A000217(n-1). - Bruno Berselli, Feb 14 2017
a(n) = T(T(n-1)) + T(T(n)) where T(n) = A000217(n). - Charlie Marion, Feb 09 2023
Sum_{n>=1} 1/a(n) = 2*(1 + Pi^2 - sqrt(3)*Pi*coth(sqrt(3)*Pi))/9. - Amiram Eldar, Feb 13 2023
a(n) = binomial(n,2)*binomial(n+1,2) + n^2 = A006011(n) + A000290(n). - Detlef Meya, Nov 23 2023

Extensions

More terms from Sam Alexander
Simplified the definition. - N. J. A. Sloane, Apr 20 2016

A020478 Number of singular 2 X 2 matrices over Z(n) (i.e., with determinant = 0).

Original entry on oeis.org

1, 10, 33, 88, 145, 330, 385, 736, 945, 1450, 1441, 2904, 2353, 3850, 4785, 6016, 5185, 9450, 7201, 12760, 12705, 14410, 12673, 24288, 18625, 23530, 26001, 33880, 25201, 47850, 30721, 48640, 47553, 51850, 55825, 83160, 51985, 72010, 77649, 106720
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := p^(2*e - 1)*(p^(e + 1) + p^e - 1); a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Oct 22 2020 *)
  • PARI
    a(n)=if(n<1, 0, direuler(p=2, n, (1-p*X)/((1-p^2*X)*(1-p^3*X)))[n])
    
  • PARI
    a(n)=local(c=0); forvec(x=vector(4,k,[1,n]),c+=((x[1]*x[2]-x[3]*x[4])%n==0)); c

Formula

From Vladeta Jovovic, Apr 22 2002: (Start)
a(n) = n^4 - A005353(n).
Multiplicative with a(p^e) = p^(2*e - 1)*(p^(e+1) + p^e - 1). (End)
Dirichlet g.f.: zeta(s-2)*zeta(s-3)/zeta(s-1).
A102631(n) | a(n). - R. J. Mathar, Mar 30 2011
Sum_{k=1..n} a(k) ~ Pi^2 * n^4 / (24*Zeta(3)). - Vaclav Kotesovec, Jan 31 2019
From Piotr Rysinski, Sep 11 2020: (Start)
a(n) = n * A069097(n).
Proof: a(n) is multiplicative with a(p^e) = p^(2*e - 1)*(p^(e+1) + p^e - 1), A069097(n) is multiplicative with A069097(p^e) = p^(e-1)*(p^e*(p+1)-1), so a(p^e) = p^e*A069097(p^e). (End)

A115077 Number of 2 X 2 symmetric matrices over Z(n) having nonzero determinant.

Original entry on oeis.org

0, 4, 18, 44, 100, 180, 294, 432, 630, 900, 1210, 1548, 2028, 2548, 3150, 3744, 4624, 5436, 6498, 7500, 8820, 10164, 11638, 13104, 14900, 16900, 18792, 20972, 23548, 26100, 28830, 31360, 34848, 38148, 41650, 44676, 49284, 53428, 57798, 62000
Offset: 1

Views

Author

T. D. Noe, Jan 12 2006

Keywords

Crossrefs

Cf. A005353 (number of 2 X 2 matrices over Z(n) having nonzero determinant), A115075.

Programs

  • Mathematica
    Table[cnt=0; Do[m={{a, b}, {b, c}}; If[Det[m, Modulus->n]>0, cnt++ ], {a, 0, n-1}, {b, 0, n-1}, {c, 0, n-1}]; cnt, {n, 50}]
    f[p_, e_] := p^e*(p^e + p^(e-1) - p^(Ceiling[e/2] - 1)); a[1] = 0; a[n_] := n^3 - Times @@ f @@@ FactorInteger[n]; Array[a, 100]  (* Amiram Eldar, Oct 31 2023 *)
  • PARI
    a(n) = {my(f = factor(n), p, e); n^3 - prod(i = 1, #f~, p = f[i, 1]; e = f[i, 2]; p^e*(p^e + p^(e-1) - p^((e+1)\2 - 1)));} \\ Amiram Eldar, Oct 31 2023

Formula

a(n) = n^3 - A115075(n).
For squarefree n, a(n) = (n-1)*n^2.

A181107 Triangle read by rows: T(n,k) is the number of 2 X 2 matrices over Z(n) having determinant congruent to k mod n, 1 <= n, 0 <= k <= n-1.

Original entry on oeis.org

1, 10, 6, 33, 24, 24, 88, 48, 72, 48, 145, 120, 120, 120, 120, 330, 144, 240, 198, 240, 144, 385, 336, 336, 336, 336, 336, 336, 736, 384, 576, 384, 672, 384, 576, 384, 945, 648, 648, 864, 648, 648, 864, 648, 648, 1450, 720, 1200, 720, 1200, 870, 1200, 720, 1200, 720
Offset: 1

Views

Author

Erdos Pal, Oct 03 2010

Keywords

Comments

The n-th row is {T(n,0),T(n,1),...,T(n,n-1)}.
Let m denote the prime power p^e.
T(m,0) = A020478(m) = (p^(e+1) + p^e-1)*p^(2*e-1).
T(m,1) = A000056(m) = (p^2-1)*p^(3*e-2).
T(prime(n),1) = A127917(n).
Sum_{k=1..n-1} T(n,k) = A005353(n).
T(n,1) = n*A007434(n) for n>=1 because A000056(n) = n*Jordan_Function_J_2(n).
T(2^n,1) = A083233(n) = A164640(2n) for n>=1. Proof: a(n):=T(2^n,1); a(1)=6, a(n)=8*a(n-1); A083233(1)=6 and A083233(n) is a geometric series with ratio 8 (because of its g.f.), too; A164640 = {b(1)=1, b(2)=6, b(n)=8*b(n-2)}.
T(2^n,0) = A165148(n) for n>=0, because 2*T(2^n,0) = (3*2^n-1)*4^n.
T(2^e,2) = A003951(e) for 2 <= e. Proof: T(2^e,2) = 9*8^(e-1) is a series with ratio 8 and initial term 72, as A003951(2...inf) is.
Working with consecutive powers of a prime p, we need a definition (0 <= i < e):
N(p^e,i):=#{k: 0 < k < p^e, gcd(k,p^e) = p^i} = (p-1)*p^(e-1-i). We say that these k's belong to i (respect to p^e). Note that N(p^e,0) = EulerPhi(p^e), and if 0 < k < p^e then gcd(k,p^e) = gcd(k,p^(e+1)). Let T(p^e,[i]) denote the common value of T(p^e,k)'s, where k's belong to i (q.v.PROGRAM); for example, T(p^e,[0]) = T(p^e,1). The number of the 2 X 2 matrices over Z(p^e), T(p^e,0) + Sum_{i=0..e-1} T(p^e,[i])*N(p^e,i) = p^(4e) will be useful.
On the hexagon property: Let prime p be given and let T(p^e,[0]), T(p^e,[1]), T(p^e,[2]), ..., T(p^e,[e-2]), T(p^e,[e-1]) form the e-th row of a Pascal-like triangle, e>=1. Let denote X(r,s) an element of the triangle and its value T(p^r,[s]). Let positive integers a and b given, so that the entries A(m-a,n-b), B(m-a,n), C(m,n+a), D(m+b,n+a), E(m+b,n), F(m,n-b) of the triangle form a hexagon spaced around T(p^m,[n]); if a=b=1 then they surround it. If A*C*E = B*D*F, then we say that the triangle T(.,.) has the "hexagon property". (In the case of binomial coefficients X(r,s) = COMB(r,s), the "hexagon property" holds (see [Gupta]) and moreover gcd(A,C,E) = gcd(B,D,F) (see [Hitotumatu & Sato]).)
Corollary 2.2 in [Brent & McKay] says that, for the d X d matrices over Z(p^e), (mutatis mutandis) T_d(p^e,0) = K*(1-P(d+e-1)/P(e-1)) and T_d(p^e,[i]) = K*(q^e)*((1-q^d)/(1-q))*P(d+i-1)/P(i), where q=1/p, K=(p^e)^(d^2), P(t) = Product_{j=1..t} (1-q^j), P(0):=1. (For the case d=2, we have T(p^e,[i]) = (p+1)*(p^(i+1)-1)*p^(3*e-i-2).) Due to [Brent & McKay], it can be simply proved that for d X d matrices the "hexagon property" is true. The formulation implies an obvious generalization: For the entries A(r,u), B(r,v), C(s,w), D(t,w), E(t,v), F(s,u) of the T_d(.,.)-triangle, a hexagon-like property A*C*E = B*D*F holds. This is false in general for the COMB(.,.)-triangle.
Another (rotated-hexagon-like) property: for the entries A(m-b1,n), B(m-a1,n+c2), C(m+a2,n+c2), D(m+b2,n), E(m+a2,n-c1), F(m-a1,n-c1) of the T_d(.,.)-triangle, the property A*C*E = B*D*F holds, if and only if 2*(a1 + a2) = b1 + b2. This is also in general false for COMB(.,.)-triangle.

Examples

			From _Andrew Howroyd_, Jul 16 2018: (Start)
Triangle begins:
    1;
   10,   6;
   33,  24,  24;
   88,  48,  72,  48;
  145, 120, 120, 120, 120;
  330, 144, 240, 198, 240, 144;
  385, 336, 336, 336, 336, 336, 336;
  736, 384, 576, 384, 672, 384, 576, 384;
  945, 648, 648, 864, 648, 648, 864, 648, 648;
  ... (End)
		

Crossrefs

Column k=0 is A020478.
Column k=1 is A000056.
Row sums are A005353.

Programs

  • Other
      (* computing T(p^e,k) ; p=prime, 1<=e, 0<=k
    				
  • PARI
    S(p,e)={my(u=vector(p^e)); my(t=(p-1)*p^(e-1)); u[1] = p^e + e*t; for(j=1, p^e-1, u[j+1] = t*(1+valuation(j, p))); vector(#u, k, sum(j=0, #u-1, u[j + 1]*u[(j+k-1) % #u + 1]))}
    T(n)={my(f=factor(n), v=vector(n,i,1)); for(i=1, #f~, my(r=S(f[i,1], f[i,2])); for(j=0, #v-1, v[j + 1] *= r[j % #r + 1])); v}
    for(n=1, 10, print(T(n))); \\ Andrew Howroyd, Jul 16 2018

Formula

T(a*b,k) = T(a,(k mod a))*T(b,(k mod b)) if gcd(a,b) = 1.
Sum_{k=1..n-1, gcd(k,n)=1} T(n,k) = A000252(n). - Andrew Howroyd, Jul 16 2018

Extensions

Terms a(24)-a(55) from b-file by Andrew Howroyd, Jul 16 2018
Showing 1-4 of 4 results.