login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A126161
Number of augmented amicable pairs (m,n) with m<n and for which m<=10^k, k=1,2,3,...
0
0, 0, 0, 1, 4, 9, 36, 84, 188, 420, 930, 1931
OFFSET
1,5
FORMULA
An augmented amicable pair (m,n) is a pair of integers m, n with m<n and for which sigma(m)=sigma(n)=m+n-1
EXAMPLE
a(6)=9 because there are 9 augmented amicable pairs with m<=10^6
MATHEMATICA
s[n_]:=DivisorSigma[1, n]-n; AugmentedAmicableNumberQ[n_]:=If[s[s[n]+1]==n-1 && !DivisorSigma[1, n]==2n-1, True, False]; AugmentedAmicablePairList[ k_]:=(bnlist=Select[Range[k], AugmentedAmicableNumberQ[ # ]&]; newprlist= Table[Sort[{bnlist[[n]], s[bnlist[[n]]]+1}], {n, 1, Length[bnlist]}]; augamprlist=Union[newprlist, newprlist]); data=AugmentedAmicablePairList[10^7]; Table[Length[Select[data, First[ # ]<10^k &]], {k, 1, 7}]
CROSSREFS
Sequence in context: A115700 A029806 A133125 * A179934 A239213 A346537
KEYWORD
hard,nonn
AUTHOR
Ant King, Dec 20 2006
STATUS
approved