A333432 A(n,k) is the n-th number m that divides k^m - 1 (or 0 if m does not exist); square array A(n,k), n>=1, k>=1, read by antidiagonals.

1, 1, 2, 1, 0, 3, 1, 2, 0, 4, 1, 3, 4, 0, 5, 1, 2, 9, 8, 0, 6, 1, 5, 4, 21, 16, 0, 7, 1, 2, 25, 6, 27, 20, 0, 8, 1, 7, 3, 125, 8, 63, 32, 0, 9, 1, 2, 49, 4, 625, 12, 81, 40, 0, 10, 1, 3, 4, 343, 6, 1555, 16, 147, 64, 0, 11, 1, 2, 9, 8, 889, 8, 3125, 18, 171, 80, 0, 12

Offset

1,3

Links

Seiichi Manyama, Antidiagonals n = 1..25, flattened

OEIS Wiki, 2^n mod n

Example

Square array A(n,k) begins:
  1, 1,  1,   1,  1,     1,  1,     1,  1, ...
  2, 0,  2,   3,  2,     5,  2,     7,  2, ...
  3, 0,  4,   9,  4,    25,  3,    49,  4, ...
  4, 0,  8,  21,  6,   125,  4,   343,  8, ...
  5, 0, 16,  27,  8,   625,  6,   889, 10, ...
  6, 0, 20,  63, 12,  1555,  8,  2359, 16, ...
  7, 0, 32,  81, 16,  3125,  9,  2401, 20, ...
  8, 0, 40, 147, 18,  7775, 12,  6223, 32, ...
  9, 0, 64, 171, 24, 15625, 16, 16513, 40, ...

Maple

A:= proc() local h, p; p:= proc() [1] end;
      proc(n, k) if k=2 then `if`(n=1, 1, 0) else
        while nops(p(k))<n do for h from 1+p(k)[-1]
          while k&^h mod h <> 1 do od;
          p(k):= [p(k)[], h]
        od; p(k)[n] fi
      end
    end():
seq(seq(A(n, 1+d-n), n=1..d), d=1..12);  # Alois P. Heinz, Mar 24 2020

Mathematica

A[n_, k_] := Module[{h, p}, p[_] = {1}; If[k == 2, If[n == 1, 1, 0], While[ Length[p[k]] < n, For[h = 1 + p[k][[-1]], Mod[k^h, h] != 1, h++]; p[k] = Append[p[k], h]]; p[k][[n]]]];
Table[A[n, 1+d-n], {d, 1, 12}, {n, 1, d}] // Flatten (* Jean-François Alcover, Nov 01 2020, after Alois P. Heinz *)

Crossrefs

Columns k=1-20 give: A000027, A063524, A067945, A014945, A067946, A014946, A067947, A014949, A068382, A014950, A068383, A014951, A116621, A177805, A014957, A177807, A128358, A333506, A128360.

Main diagonal gives A333433.

Cf. A333429, A333500.

Sequence in context: A130106 A127093 A141543 * A333500 A182720 A347316

Adjacent sequences: A333429 A333430 A333431 * A333433 A333434 A333435

Keyword

nonn,tabl

Author

Seiichi Manyama, Mar 21 2020

Status

approved