|
|
A372902
|
|
Numbers k that divide the k-th large Schröder number.
|
|
1
|
|
|
1, 2, 6, 33, 42, 154, 198, 258, 270, 342, 850, 1170, 1666, 1806, 2295, 2574, 3262, 3366, 3834, 4070, 4654, 4970, 5439, 6006, 6118, 6162, 6699, 7095, 7254, 7497, 7595, 10241, 11475, 12642, 14014, 15345, 17470, 17670, 18018, 19845, 22446, 23994, 24570, 24651, 25245, 25974, 26334
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Numbers k such that k | A006318(k).
|
|
LINKS
|
|
|
EXAMPLE
|
2 is a term since A001850(2) = 6 = 2 * 3 is divisible by 2.
6 is a term since A001850(6) = 1806 = 6 * 301 is divisible by 6.
|
|
MATHEMATICA
|
seq[kmax_] := Module[{sc0 = 1, sc1 = 2, sc2, s = {1}}, Do[sc2 = ((6*k-3)*sc1 - (k-2)*sc0)/(k+1); If[Divisible[sc2, k], AppendTo[s, k]]; sc0 = sc1; sc1 = sc2, {k, 2, kmax}]; s]; seq[27000]
|
|
PROG
|
(PARI) lista(kmax) = {my(sc0 = 1, sc1 = 2, sc2); print1(1, ", "); for(k = 2, kmax, sc2 = ((6*k-3)*sc1 - (k-2)*sc0)/(k+1); if(!(sc2 % k), print1(k, ", ")); sc0 = sc1; sc1 = sc2); }
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|