OFFSET
1,4
LINKS
Pedersen J. M., Known amicable pairs.
FORMULA
The number of integer pairs (m,n) with m<n for which usigma(m)=usigma(n)=m+n.
EXAMPLE
a(7)=78 because there are 78 unitary amicable pairs with m<=10^7
MATHEMATICA
UnitaryDivisors[n_Integer?Positive] := Select[Divisors[n], GCD[ #, n/# ] == 1 \ &]; sstar[n_] := Plus @@ UnitaryDivisors[ n] - n; UnitaryAmicableNumberQ[k_] := If[Nest[sstar, k, 2] == k && ! sstar[k] == k, True, False]; UnitaryAmicablePairList[ k_] := (anlist = Select[Range[k], UnitaryAmicableNumberQ[ # ] &]; prlist = \ Table[Sort[{anlist[[n]], sstar[anlist[[ n]]]}], {n, 1, Length[anlist]}]; data = Union[prlist, prlist]); UnitaryAmicablePairList[10^7]; Table[Length[Select[data, First[ # ] < \ 10^k &]], {k, 1, 7}]
CROSSREFS
KEYWORD
hard,nonn
AUTHOR
Ant King, Dec 20 2006
STATUS
approved