|
1, 1, 1, 3, 1, 3, 1, 3, 1, 30, 1, 3, 1, 3, 10, 21, 1, 3, 1, 30, 1, 126, 1, 21, 10, 3, 1, 315, 1, 30, 1, 21, 42, 66, 10, 3, 1, 3, 1, 11550, 1, 315, 1, 126, 10, 990, 1, 21, 1, 30, 22, 693, 1, 3, 420, 2205, 1, 2310, 1, 1650, 1, 3, 1, 273, 10, 126, 1, 66, 330, 245700, 1, 21, 1, 3, 10, 585, 42, 693, 1, 11550, 1, 546, 1, 315, 220, 3, 770
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,4
|
|
LINKS
|
|
|
FORMULA
|
a(n) = Product_{d|n, d<n} A019565(d)^[2 == d mod 3].
For all n >= 1:
|
|
PROG
|
(PARI)
A019565(n) = {my(j, v); factorback(Mat(vector(if(n, #n=vecextract(binary(n), "-1..1")), j, [prime(j), n[j]])~))}; \\ From A019565
A319992(n) = { my(m=1); fordiv(n, d, if((d<n)&&(2==(d%3)), m *= A019565(d))); m; };
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|