OFFSET
1,2
COMMENTS
The integer sums d(k)/d(1) + d(k-1)/d(2) + ... + d(k)/d(1) are given by A229999. - Clark Kimberling, Jun 16 2018
Also numbers m such that the sum of the squares of the unitary divisors of m is divisible by m (the unitary version of A046762) - Amiram Eldar, Jun 16 2018
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..333 (terms below 10^10)
EXAMPLE
The first 10 sums: 1, 5/2, 10/3, 17/4, 26/5, 25/3, 50/7, 65/8, 82/9, 13, so that a(1) = 1 and a(10) = 13.
MATHEMATICA
z = 1000; r[n_] := Select[Divisors[n], GCD[#, n/#] == 1 &]; k[n_] := Length[r[n]];
t[n_] := Table[r[n][[k[n] + 1 - i]]/r[n][[k[1] + i - 1]], {i, 1, k[n]}];
s = Table[Plus @@ t[n], {n, 1, z}]; a[n_] := If[IntegerQ[s[[n]]], 1, 0]; u = Table[a[n], {n, 1, z}]; Flatten[Position[u, 1]] (* A229996 *)
usigma2[n_] := If[n == 1, 1, Times @@ (1 + Power @@@ FactorInteger[n]^2)]; seqQ[n_] := Divisible[usigma2[n], n]; Select[Range[10^6], seqQ] (* Amiram Eldar, Jun 16 2018 *
PROG
(PARI) is(n) = {my(f = factor(n)); !(prod(i = 1, #f~, f[i, 1]^(2*f[i, 2]) + 1) % n); } \\ Amiram Eldar, Jun 16 2024
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Oct 31 2013
EXTENSIONS
Definition corrected by Clark Kimberling, Jun 16 2018
STATUS
approved