|
|
A265578
|
|
LCM-transform of number of divisors function (A000005).
|
|
2
|
|
|
1, 2, 1, 3, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
Terms larger than one occur at n = 2, 4, 6, 16, 24, 36, 64, 120, 840, 900, 1296, 7560, 44100, 46656, 83160, ... - Antti Karttunen, Nov 06 2018
|
|
LINKS
|
|
|
MAPLE
|
LCMXfm:=proc(a) local L, i, n, g, b;
L:=nops(a);
g:=Array(1..L, 0); b:=Array(1..L, 0);
b[1]:=a[1]; g[1]:=a[1];
for n from 2 to L do g[n]:=ilcm(g[n-1], a[n]); b[n]:=g[n]/g[n-1]; od;
lprint([seq(b[i], i=1..L)]);
end;
with(numtheory);
t1:=[seq(tau(n), n=1..100)];
LCMXfm(t1);
|
|
MATHEMATICA
|
LCMXfm[a_List] := Module[{L = Length[a], b, g}, b[1] = g[1] = a[[1]]; b[_] = 0; g[_] = 0; Do[g[n] = LCM[g[n-1], a[[n]]]; b[n] = g[n]/g[n-1], {n, 2, L}]; Array[b, L]];
|
|
PROG
|
(PARI)
up_to = 16384;
LCMtransform(v) = { my(len = length(v), b = vector(len), g = vector(len)); b[1] = g[1] = 1; for(n=2, len, g[n] = lcm(g[n-1], v[n]); b[n] = g[n]/g[n-1]); (b); };
v265578 = LCMtransform(vector(up_to, i, numdiv(i)));
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|