login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A276517
Indices k such that A276516(k) = 0.
10
2, 3, 6, 7, 8, 11, 12, 15, 18, 19, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 41, 43, 44, 45, 46, 47, 48, 53, 54, 60, 61, 67, 70, 72, 74, 76, 79, 82, 84, 87, 90, 92, 93, 96, 105, 106, 107, 108, 111, 112, 114, 117, 122, 128, 133, 135, 139, 141, 148, 159
OFFSET
1,1
COMMENTS
This is different from A001422, first difference: a(14) = 25, A001422(14) = 27.
Conjecture: for k > 7169 there are no more terms in this sequence (tested for k < 10000000).
LINKS
EXAMPLE
3 is in the sequence because A276516(3) = 0
4 is not in the sequence because A276516(4) = -1
4222 is in the sequence because A276516(4222) = 0
7169 is in the sequence because A276516(7169) = 0
MATHEMATICA
nn = 100; A276516 = Rest[CoefficientList[Series[Product[(1-x^(k^2)), {k, nn}], {x, 0, nn^2}], x]]; Select[Range[nn^2], A276516[[#]]==0&]
nmax = 10000; nn = Floor[Sqrt[nmax]]+1; poly = ConstantArray[0, nn^2 + 1]; poly[[1]] = 1; poly[[2]] = -1; poly[[3]] = 0; Do[Do[poly[[j + 1]] -= poly[[j - k^2 + 1]], {j, nn^2, k^2, -1}]; , {k, 2, nn}]; A276516 = Take[poly, {2, nmax+1}]; Select[Range[nmax], A276516[[#]]==0&]
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Dec 12 2016
STATUS
approved