Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A369952
a(n) is the number of distinct values of the permanent of an n X n Hankel matrix using the first 2*n - 1 prime numbers.
3
1, 1, 2, 59, 2493, 180932, 19939272
OFFSET
0,3
MATHEMATICA
a[n_] := CountDistinct[Table[Permanent[HankelMatrix[Join[Drop[per = Part[Permutations[Prime[Range[2 n - 1]]], i], n], {Part[per, n]}], Join[{Part[per, n]}, Drop[per, - n]]]], {i, (2 n - 1) !}]]; Join[{1}, Array[a, 5]]
PROG
(PARI) a(n) = my(v=[1..2*n-1], list=List()); forperm(v, p, listput(list, matpermanent(matrix(n, n, i, j, prime(p[i+j-1])))); ); #Set(list); \\ Michel Marcus, Feb 08 2024
(Python)
from itertools import permutations
from sympy import primerange, prime, Matrix
def A369952(n): return len({Matrix([p[i:i+n] for i in range(n)]).per() for p in permutations(primerange(prime((n<<1)-1)+1))}) if n else 1 # Chai Wah Wu, Feb 12 2024
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Stefano Spezia, Feb 06 2024
EXTENSIONS
a(6) from Michel Marcus, Feb 08 2024
STATUS
approved