A049763 Triangular array T, read by rows: T(n,k) = n^4 mod k, for k = 1..n and n >= 1.

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

Author

Clark Kimberling

Keywords

• nonn
• tabl

Examples

			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;
  ...
		

Links

• G. C. Greubel, Rows n = 1..100 of triangle, flattened

Crossrefs

Row sums are in A049764.

Cf. A049759, A049760, A049761, A049762.

Programs

GAP

Flat(List([1..15], n-> List([1..n], k-> PowerMod(n,4,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

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 *)

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

Sage

[[power_mod(n,4,k) for k in (1..n)] for n in (1..15)] # G. C. Greubel, Dec 13 2019