login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A053191 a(n) = n^2 * phi(n). 27

%I #51 Sep 08 2022 08:45:00

%S 1,4,18,32,100,72,294,256,486,400,1210,576,2028,1176,1800,2048,4624,

%T 1944,6498,3200,5292,4840,11638,4608,12500,8112,13122,9408,23548,7200,

%U 28830,16384,21780,18496,29400,15552,49284,25992,36504,25600,67240

%N a(n) = n^2 * phi(n).

%C Number of invertible 2 X 2 symmetric matrices over Z(n). - _T. D. Noe_, Jan 13 2006

%C Note that A115077 gives the number of 2 X 2 symmetric matrices having nonzero determinant. However, for composite n, a nonzero determinant is not sufficient for the matrix to be invertible; the determinant must also be relatively prime to n. - _T. D. Noe_, Jan 13 2006

%C Also Euler phi function of n^3.

%C For n^k, EulerPhi(n^k) = n^(k-1)*EulerPhi(n). The same holds if Phi is replaced by the cototient function.

%C Also, the sum of the degrees of the irreducible representations of the group GL(2,Z_n) (sequence A000252). - Sharon Sela (sharonsela(AT)hotmail.com), Feb 06 2002

%H Seiichi Manyama, <a href="/A053191/b053191.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = n^2 * phi(n) = A000010(n^3).

%F Dirichlet g.f.: zeta(s-3)/zeta(s-2). - _R. J. Mathar_, Feb 09 2011

%F The n-th term of the Dirichlet inverse is n^2 * A023900(n) = (-1)^omega(n) * a(n) / A003557(n), where omega = A001221. - _Álvar Ibeas_, Nov 24 2017

%F Sum_{n>=1} 1/a(n) = Product_{p prime} (1 + p/(p^4 - p^3 - p + 1)) = 1.38097852211302096879... - _Amiram Eldar_, Dec 06 2020

%e n=5: n^3 = 125, EulerPhi(125) = 125 - 25 = 100.

%p with(numtheory):a:=n->phi(n^3): seq(a(n), n=1..41); # _Zerinvary Lajos_, Oct 07 2007

%t Table[cnt=0; Do[m={{a, b}, {b, c}}; If[Det[m, Modulus->n]>0 && MatrixQ[Inverse[m, Modulus->n]], cnt++ ], {a, 0, n-1}, {b, 0, n-1}, {c, 0, n-1}]; cnt, {n, 2, 50}] (* _T. D. Noe_, Jan 13 2006 *)

%t Table[n^2*EulerPhi[n],{n,1,40}] (* _Vladimir Joseph Stephan Orlovsky_, Nov 10 2009 *)

%o (Sage) [n^2*euler_phi(n) for n in range(1, 42)] # _Zerinvary Lajos_, Jun 06 2009

%o (Magma) [ n^2*EulerPhi(n) : n in [1..100] ]; // _Vincenzo Librandi_, Apr 21 2011

%o (PARI) a(n) = n^2*eulerphi(n); \\ _Michel Marcus_, Oct 31 2017

%Y Cf. A000252 (number of invertible 2 X 2 matrices over Z(n)), A115075, A115076, A115077.

%Y Cf. A000010, A051953, A002618, A053650, A053192, A001248, A319087.

%K nonn,mult

%O 1,2

%A _Labos Elemer_, Mar 02 2000

%E Edited by _N. J. A. Sloane_ at the suggestion of _Andrew S. Plewe_, Jun 05 2007

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 25 01:35 EDT 2024. Contains 371964 sequences. (Running on oeis4.)