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.

A335939 Bi-unitary weird numbers (A292986) that are not exponentially odd numbers (A268335).

This page as a plain text file.
%I A335939 #12 Mar 25 2023 05:29:29
%S A335939 5390,7400,11830,20230,24010,25270,37030,58870,67270,95830,117670,
%T A335939 129430,154630,196630,243670,260470,314230,352870,373030,436870,
%U A335939 482230,491744,507232,554470,658630,714070,742630,801430,831670,893830,1129030,1201270,1313830,1352470
%N A335939 Bi-unitary weird numbers (A292986) that are not exponentially odd numbers (A268335).
%C A335939 Weird numbers (A006037) that are exponentially odd (A268335) are also bi-unitary weird numbers (A292986), since all of their divisors are bi-unitary.
%H A335939 Amiram Eldar, <a href="/A335939/b335939.txt">Table of n, a(n) for n = 1..10000</a>
%t A335939 fun[p_, e_] := If[OddQ[e], (p^(e+1)-1)/(p-1), (p^(e+1)-1)/(p-1)-p^(e/2)]; bsigma[1] = 1; bsigma[n_] := Times @@ (fun @@@ FactorInteger[n]); biuabQ[n_] := bsigma[n] > 2*n; f[n_] := Select[Divisors[n], Function[d, CoprimeQ[d, n/d]]]; bdiv[m_] := Select[Divisors[m], Last@Intersection[f@#, f[m/#]] == 1 &]; bweirdQ[n_] := biuabQ[n] && Module[{d = Most @ bdiv[n], x}, SeriesCoefficient[Series[Product[1 + x^d[[i]], {i, Length[d]}], {x, 0, n}], n] == 0]; expOddQ[n_] := AllTrue[Last /@ FactorInteger[n], OddQ]; Select[Range[1000], !expOddQ[#] && bweirdQ[#] &]
%Y A335939 Intersection of A072587 and A292986.
%Y A335939 Cf. A006037, A222266, A268335, A335936, A335938.
%K A335939 nonn
%O A335939 1,1
%A A335939 _Amiram Eldar_, Jun 30 2020
%E A335939 More terms from _Amiram Eldar_, Mar 25 2023