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!)
A278847 a(n) = permanent M_n where M_n is the n X n matrix m(i,j) = i^2 + j^2. 10
1, 2, 41, 3176, 620964, 246796680, 174252885732, 199381727959680, 345875291854507584, 864860593764292790400, 2996169331694350840741440, 13929521390709644084719495680, 84659009841182126038701730464000, 658043094413184868424932006273344000 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,2

COMMENTS

From Zhi-Wei Sun, Aug 19 2021: (Start)

I have proved that a(n) == (-1)^(n-1)*2*n! (mod 2n+1) whenever 2n+1 is prime.

Conjecture 1: If 2n+1 is composite, then a(n) == 0 (mod 2n+1).

Conjecture 2: If p = 4n+1 is prime, then the sum of those Product_{j=1..2n}(j^2-f(j)^2)^{-1} with f over all the derangements of {1,...,2n} is congruent to 1/(n!)^2 modulo p. (End)

LINKS

Vaclav Kotesovec, Table of n, a(n) for n = 0..36

Zhi-Wei Sun, Arithmetic properties of some permanents, arXiv:2108.07723 [math.GM], 2021.

FORMULA

a(n) ~ c * d^n * (n!)^3 / n, where d = 3.809076776112918119... and c = 1.07739642254738...

MAPLE

with(LinearAlgebra):

a:= n-> `if`(n=0, 1, Permanent(Matrix(n, (i, j)-> i^2+j^2))):

seq(a(n), n=0..16);  # after Alois P. Heinz

MATHEMATICA

Flatten[{1, Table[Permanent[Table[i^2+j^2, {i, 1, n}, {j, 1, n}]], {n, 1, 15}]}]

PROG

(PARI) a(n)={matpermanent(matrix(n, n, i, j, i^2 + j^2))} \\ Andrew Howroyd, Aug 21 2018

CROSSREFS

Cf. A005249, A085750, A085807, A204249, A278845, A278925, A278926, A346934, A346949.

Sequence in context: A054742 A113634 A098634 * A297007 A193270 A350859

Adjacent sequences:  A278844 A278845 A278846 * A278848 A278849 A278850

KEYWORD

nonn

AUTHOR

Vaclav Kotesovec, Nov 29 2016

STATUS

approved

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 July 7 11:47 EDT 2022. Contains 355148 sequences. (Running on oeis4.)