cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A380233 Odd abundant numbers not divisible by 5 that are also doublets (cf. A020338).

This page as a plain text file.
%I A380233 #40 Jan 18 2025 09:35:38
%S A380233 153153,171171,189189,207207,243243,261261,279279,297297,351351,
%T A380233 459459,513513,567567,621621,729729,783783,837837,891891,999999,
%U A380233 1392313923,1556115561,1719917199,1883718837,2034920349,2211322113,2375123751,2538925389,2702727027,3194131941,4176941769,4668346683
%N A380233 Odd abundant numbers not divisible by 5 that are also doublets (cf. A020338).
%C A380233 There are 26 odd abundant numbers not divisible by 5 less than 10^6. The surprising fact is that 18 of them are doublets.
%C A380233 Another interesting fact is that here there are no terms with 8 digits.
%p A380233 doublet:= n -> n * (10^(1+ilog10(n))+1):
%p A380233 select(t -> numtheory:-sigma(t) > 2*t, [seq(seq(doublet(10*x+i),i=[1,3,7,9]),x=1..10000); # _Robert Israel_, Jan 17 2025
%t A380233 Select[Table[FromDigits[Join[#, #] &@ IntegerDigits[n]], {n, Select[Range[50000], CoprimeQ[#, 10] &]}], DivisorSigma[-1, #] > 2 &] (* _Amiram Eldar_, Jan 17 2025 *)
%o A380233 (PARI) select(x->((x%5) && (sigma(x)>2*x)), vector(50000, n, eval(Str(2*n-1, 2*n-1)))) \\ _Michel Marcus_, Jan 17 2025
%Y A380233 Intersection of A064001 and A020338.
%Y A380233 Cf. A005101, A005231, A047201.
%K A380233 nonn,base
%O A380233 1,1
%A A380233 _Omar E. Pol_, Jan 17 2025
%E A380233 More terms from _Michel Marcus_, Jan 17 2025