OFFSET
1,3
COMMENTS
Row n begins with 1 and has prime(n)-1 terms. The first differences of each row are symmetric. For k > p^2, the solutions are just shifted by m*p^2 for m > 0. An open question is whether every integer appears in this sequence. For instance, 2 does not appear until the prime 1093 and 5 does not appear until the prime 20771.
For row n > 1, the sum of the terms in row n is (p-1)*p^2*(p+1)/2, which is A138416. - T. D. Noe, Aug 24 2008, corrected by Robert Israel, Sep 27 2016
Max Alekseyev points out that there is a much faster method of computing these numbers. Let p=prime(n) and let r be a primitive root of p (see A001918 and A060749). Then the terms in row n are r^(k*p) (mod p^2) for k=0..p-2. - T. D. Noe, Aug 26 2008
The numbers in this sequence are the bases to Euler pseudoprimes q, which are squares of prime numbers, such that n^((q-1)/2) == +-1 mod q. An exception is the first number 9 = 3*3, which is, following the strict definition in Crandall and Pomerance, no Fermat pseudoprime and hence no Euler pseudoprime. - Karsten Meyer, Jan 08 2011
For row n > 1, the sum is zero modulo p^2 (rows are antisymmetric due to Binomial Theorem). - Peter A. Lawrence, Sep 11 2016
REFERENCES
R. Crandall and C. Pomerance, Prime Numbers: A Computational Perspective, Springer, NY, 2005
LINKS
T. D. Noe, Rows n=1..50 of triangle, flattened
EXAMPLE
(2) 1,
(3) 1, 8,
(5) 1, 7, 18, 24,
(7) 1, 18, 19, 30, 31, 48,
(11) 1, 3, 9, 27, 40, 81, 94, 112, 118, 120,
(13) 1, 19, 22, 23, 70, 80, 89, 99, 146, 147, 150, 168,
(17) 1, 38, 40, 65, 75, 110, 131, 134, 155, 158, 179, 214, 224, 249, 251, 288,
MAPLE
f:= proc(n) local p, j, x;
p:= ithprime(n);
x:= numtheory:-primroot(p);
op(sort([seq(x^(i*p) mod p^2, i=0..p-2)]))
end proc:
map(f, [$1..20]); # Robert Israel, Sep 27 2016
MATHEMATICA
CROSSREFS
KEYWORD
nonn,tabf
AUTHOR
T. D. Noe, Aug 24 2008
STATUS
approved