|
|
A309520
|
|
Primes p for which h1(p)/G(p) has a record value.
|
|
0
|
|
|
3, 5, 7, 11, 23, 73, 89, 179, 233, 761, 1451, 2741, 4391, 5231, 42611, 198221, 305741, 6766811, 1326662801, 1979990861, 4735703723, 9697282541, 35285447111, 45169368641, 169684421321, 187946428721
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
h1 is given by A000927, and G(p) = 2*p*(p/(4*Pi^2))^((p-1)/4).
a(1)-a(12) from Fung et al., a(13)-a(14) from Shokrollahi, a(15)-a(17) from Broadhurst, a(18) from Languasco et al. and Broadhurst, a(19)-a(26) from Broadhurst.
|
|
LINKS
|
|
|
PROG
|
(PARI) h1(p) = if (p<5, 1, abs( matdet(matrix((p-1)/2-2, (p-1)/2-2, i, j, ((i+2)*(j+2))\p - ((i+1)*(j+2))\p)) )); \\ A000927
G(p) = 2*p*(p/(4*Pi^2))^((p-1)/4);
lista(nn) = {my(m = 0, nm); for (n=2, nn, p = prime(n); if ((nm = h1(p)/G(p)) > m, print1(p, ", "); m = nm); ); }
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,more
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|