STATUS
reviewed
approved
reviewed
approved
proposed
reviewed
editing
proposed
This is different from A001422 is a finite subsequence.
proposed
editing
editing
proposed
2229 is in the sequence because A033461(2229) = 51267 = 23 * 2229.
proposed
editing
editing
proposed
Numbers k such that k divides A033461(k) is divisible by k.
max = 100; p = ConstantArray[0, max^2 + 1]; p[[1]] = 1; p[[2]] = 1; Do[Do[p[[j + 1]] += p[[j - k^2 + 1]], {j, max^2, k^2, -1}]; , {k, 2, max}]; Select[Range[1, max^2], ModDivisible[p[[# + 1]], #] == 0 &]