cp's OEIS Frontend

A353933 a(n) is the permanent of the n X n symmetric matrix M(n) whose generic element M[i,j] is equal to the digital root of i*j.

%I A353933 #14 Sep 28 2022 08:13:10
%S A353933 1,1,8,216,7344,168183,7226091,295259094,11801772252,1673511251940,
%T A353933 65568867621336,2710049208604776,202103867012027328,
%U A353933 12881755844526953376,736186737257150962752,70484099228399057425344,5507570249593121504026368,434305172863416192470350848,122043063804581668929348667392
%N A353933 a(n) is the permanent of the n X n symmetric matrix M(n) whose generic element M[i,j] is equal to the digital root of i*j.
%C A353933 The matrix M(n) is nonsingular only for n = 1, 5 and 6 with determinant equal respectively to 1, 6561 and 59049.
%C A353933 The rank of M(n) is 1 for 1 <= n <= 3, 3 for n = 4, 5 for n = 5, 6 for 6 <= n <= 8, and 7 for n >= 9. - _Jianing Song_, Sep 28 2022
%F A353933 Sum_{i=1..n} M[n-i+1,i] = A353128(n+1).
%e A353933 a(7) = 7226091:
%e A353933      1, 2, 3, 4, 5, 6, 7
%e A353933      2, 4, 6, 8, 1, 3, 5
%e A353933      3, 6, 9, 3, 6, 9, 3
%e A353933      4, 8, 3, 7, 2, 6, 1
%e A353933      5, 1, 6, 2, 7, 3, 8
%e A353933      6, 3, 9, 6, 3, 9, 6
%e A353933      7, 5, 3, 1, 8, 6, 4
%t A353933 M[i_, j_]:=If[i*j==0, 0, 1+Mod[i*j-1, 9]]; Join[{1},Table[Permanent[Table[M[i, j], {i,  n}, {j, n}]],{n,18}]]
%o A353933 (PARI) a(n) = matpermanent(matrix(n, n, i, j, (i*j-1)%9+1)); \\ _Michel Marcus_, May 12 2022
%Y A353933 Cf. A003991, A010888, A353109, A353128, A353974 (trace of the matrix M(n)).
%K A353933 nonn,base
%O A353933 0,3
%A A353933 _Stefano Spezia_, May 11 2022