login
Write 2 - Zeta(s) in the form 1/Product_{n > 1}(1 + a(n)/n^s).
2

%I #5 Dec 19 2017 02:38:00

%S 1,1,2,1,2,1,2,2,2,1,4,1,2,2,6,1,4,1,4,2,2,1,8,2,2,2,4,1,6,1,6,2,2,2,

%T 12,1,2,2,8,1,6,1,4,4,2,1,16,2,4,2,4,1,8,2,8,2,2,1,16,1,2,4,10,2,6,1,

%U 4,2,6,1,24,1,2,4,4,2,6,1,16,6,2,1,16,2,2

%N Write 2 - Zeta(s) in the form 1/Product_{n > 1}(1 + a(n)/n^s).

%t nn=100;

%t facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];

%t -Solve[Table[-1==Sum[Times@@a/@f,{f,facs[n]}],{n,2,nn}],Table[a[n],{n,2,nn}]][[1,All,2]]

%Y Cf. A001055, A045778, A050376, A220418, A220420, A273866, A273873, A289501, A290261, A290262, A290971, A290973, A295279, A295632, A295636.

%K nonn

%O 2,3

%A _Gus Wiseman_, Nov 24 2017