A conjecture by the late Romanian mathematician Alexandru Lupas.
Posted in [sci.math][1] in 2005, but no proof was found.
Physicist [Alan Sokal][2] just reminded me of it, saying it was related to something he is working on.

Let $P_n(z)$ be the Legendre polynomials, defined by the generating function
$$
\big(1-2tz+t^2\big)^{-1/2} = \sum_{k=0}^\infty t^k P_k(z) .
$$
Let $g(\alpha,\beta) = 4\cos(2\alpha)+8\sin(\beta)\cos(\alpha)+5$ be defined for 
$(\alpha,\beta) \in (-\pi,\pi)\times (-\pi,\pi)$ .  Let $A_n$ be these [Apéry numbers][3]
$$
A_n = \sum_{k=0}^n \binom{n}{k}\binom{n+k}{k}\sum_{j=0}^k \binom{k}{j}^3 .
$$
  
**examples**  
$$
P_0(z)=1;\qquad
P_1(z)=z;\qquad
P_2(z)=\frac{3}{2}\;z^2-\frac{1}{2};\qquad
P_3(z)=\frac{5}{2}\;z^3-\frac{3}{2}\;z;
\\
A_0 = 1;\qquad A_1=5;\qquad A_2=73;\qquad A_3=1445;\qquad A_4=33001 .
$$

**Prove or disprove:**  
$A_n$ is the average of $P_n\circ g$.  More explicitly:  for all natural numbers $n$,
$$
A_n = \frac{1}{4\pi^2}\int_{-\pi}^\pi \int_{\pi}^\pi 
P_n\big(g(\alpha,\beta)\big)\;d\beta\;d\alpha
$$
  
This is surely true (Sokal says he has checked it through $n=123$).  But can **you** prove it?

**additional notes**  

Also
$$
A_n = \sum_{k=0}^n \binom{k}{n}^2\binom{n+k}{k}^2
$$
The conjecture should be equivalent to
$$
\frac{1}{4\pi^2}\int_{-\pi}^\pi\int_{-\pi}^\pi\frac{d\alpha\;d\beta}{\sqrt{
t^2-2t(4\cos(2\alpha)+8\sin\beta\cos\alpha+5)+1}\;}
=\sum_{k=0}^\infty A_k t^k
$$
In the integral we can change variables to get
$$
\frac{1}{\pi^2}\int_{-1}^1\int_{-1}^1\frac{dp\;dq}{\sqrt{1-p^2}
\sqrt{1-q^2}\sqrt{1-2t(8pq+8q^2+1)+t^2}\;}
$$

  [1]: http://mathforum.org/kb/message.jspa?messageID=3704516
  [2]: http://www.physics.nyu.edu/faculty/sokal/
  [3]: http://oeis.org/A005259