login
A049763
Triangular array T, read by rows: T(n,k) = n^4 mod k, for k = 1..n and n >= 1.
4
0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 1, 4, 1, 0, 0, 1, 0, 1, 1, 3, 2, 1, 0, 0, 0, 1, 0, 0, 4, 4, 0, 1, 0, 0, 1, 1, 1, 1, 1, 4, 1, 7, 1, 0, 0, 0, 0, 0, 1, 0, 2, 0, 0, 6, 1, 0, 0, 1, 1, 1, 1, 1
OFFSET
1,34
EXAMPLE
Triangle T(n,k) (with rows n >= 1 and columns k >= 1) begins as follows:
0;
0, 0;
0, 1, 0;
0, 0, 1, 0;
0, 1, 1, 1, 0;
0, 0, 0, 0, 1, 0;
0, 1, 1, 1, 1, 1, 0;
0, 0, 1, 0, 1, 4, 1, 0;
0, 1, 0, 1, 1, 3, 2, 1, 0;
0, 0, 1, 0, 0, 4, 4, 0, 1, 0;
...
MAPLE
seq(seq( `mod`(n^4, k), k = 1..n), n = 1..20); # G. C. Greubel, Dec 13 2019
MATHEMATICA
Flatten[Table[PowerMod[n, 4, k], {n, 20}, {k, n}]] (* Harvey P. Dale, Jan 19 2015 *)
PROG
(PARI) T(n, k) = lift(Mod(n, k)^4);
for(n=1, 15, for(k=1, n, print1(T(n, k), ", "))) \\ G. C. Greubel, Dec 13 2019
(Magma) [[Modexp(n, 4, k): k in [1..n]]: n in [1..15]]; // G. C. Greubel, Dec 13 2019
(Sage) [[power_mod(n, 4, k) for k in (1..n)] for n in (1..15)] # G. C. Greubel, Dec 13 2019
(GAP) Flat(List([1..15], n-> List([1..n], k-> PowerMod(n, 4, k) ))); # G. C. Greubel, Dec 13 2019
CROSSREFS
Row sums are in A049764.
Sequence in context: A228623 A036875 A036877 * A328290 A182878 A221971
KEYWORD
nonn,tabl
STATUS
approved