A320632 Gus Wiseman OEIS

nn = 390;
TableForm[
  M = Table[
    Table[If[Mod[n, k] == 0, If[n == k, 0, 1], 0], {k, 1, nn}], {n, 1,
      nn}]];
Table[MatrixPower[M, p][[All, 1]], {p, 1, 16}];
Transpose[
  Table[Table[(+1)^(k - 1)*Binomial[PrimeOmega[n] - 1, k]*
     Binomial[PrimeNu[n] - 1 + k, k], {k, 0, 15}], {n, 1, nn}]];
Sign[Total[%% - %]]
Flatten[Position[%, 1]]