This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A320161 #27 Oct 20 2018 10:20:47 %S A320161 0,1,3,0,1,8,0,1,24,0,1,10,39,48,0,1,27,36,37,84,85,94,120,0,1,6,29, %T A320161 34,61,108,135,140,163,168,0,1,25,45,56,75,82,132,157,207,214,233,244, %U A320161 264,288,0,1,42,43,73,88,106,120,161,200,241,255,273,288,318,319,360 %N A320161 Irregular triangle read by rows: row n lists 0 <= k < p^2 such that p^2 divides A316269(k, p-Kronecker(k^2-4, p)), p = prime(n). %C A320161 p always divides A316269(k, p-Kronecker(k^2-4, p)), so it's interesting to see when p^2 also divides A316269(k, p-Kronecker(k^2-4, p)). %C A320161 In the following comments, let p = prime(n). Note that A316269(0, m) and A316269(1, m) is not defined, so here k must be understood as a remainder modulo p^2. because A316269(k+s*p^2, m) == A316269(k, m) (mod p^2). %C A320161 Let p = prime(n). Every row contains 0, 1 and p^2 - 1. For n >= 3, the n-th row contains p - 2 numbers, whose remainders modulo p form a permutation of {0, 1, 3, 4, ..., p - 3, p - 1}. %C A320161 Every row is antisymmetric, that is, k is a member iff p^2 - k is, k > 0. As a result, the sum of the n-th row is prime(n)^2*(prime(n) - 3)/2. %C A320161 Equivalently, for n >= 2, row n lists 0 <= k < p^2 such that p^2 divides A316269(k, (p-Kronecker(k^2-4, p))/2), p = prime(n). %H A320161 Jianing Song, <a href="/A320161/b320161.txt">Table of n, a(n) for n = 1..29083</a> (primes below 600) %H A320161 Jianing Song, <a href="/A320161/a320161.txt">Table of n, a(n) for n = 1..75796</a> (primes below 1000) %H A320161 Wikipedia, <a href="https://en.wikipedia.org/wiki/Wall-Sun-Sun_prime">Wall-Sun-Sun prime</a> %e A320161 Table starts %e A320161 p = 2: 0, 1, 3, %e A320161 p = 3: 0, 1, 8, %e A320161 p = 5: 0, 1, 24, %e A320161 p = 7: 0, 1, 10, 39, 48, %e A320161 p = 11: 0, 1, 27, 36, 37, 84, 85, 94, 120, %e A320161 p = 13: 0, 1, 6, 29, 34, 61, 108, 135, 140, 163, 168, %e A320161 p = 17: 0, 1, 25, 45, 56, 75, 82, 132, 157, 207, 214, 233, 244, 264, 288, %e A320161 p = 19: 0, 1, 42, 43, 73, 88, 106, 120, 161, 200, 241, 255, 273, 288, 318, 319, 360, %e A320161 p = 23: 0, 1, 12, 15, 60, 86, 105, 141, 142, 156, 223, 306, 373, 387, 388, 424, 443, 469, 514, 517, 528, %e A320161 p = 29: 0, 1, 42, 46, 80, 101, 107, 120, 226, 227, 327, 330, 358, 409, 432, 483, 511, 514, 614, 615, 721, 734, 740, 761, 795, 799, 840, %e A320161 ... %o A320161 (PARI) %o A320161 B(k, p) = (([k, -1; 1, 0]^(p-kronecker(k^2-4,p)))[1,2])%(p^2) %o A320161 forprime(p=2, 50, for(k=0, p^2-1, if(!B(k, p), print1(k, ", ")));print) %Y A320161 Cf. A143548, A316269, A320162 (discriminant k^2+4, a more studied case). %Y A320161 Cf. A238490 (primes p such that 4 occurs in the corresponding row), A238736 (primes p such that 6 occurs in the corresponding row). %K A320161 nonn,tabf %O A320161 1,3 %A A320161 _Jianing Song_, Oct 06 2018