cp's OEIS Frontend

A034094 (-1)sigma perfect numbers: (-1)sigma(a) = m*a for some integer m, where if a = Product p(i)^r(i) then (-1)sigma(a) = Product_{i} (-1 + Sum_{s=1..r(i)} p(i)^s).

%I A034094 #22 Jul 07 2022 08:11:46
%S A034094 1,20,312,9744,29280,53352,1666224,5006880,106798080,133301760,
%T A034094 980733600,9099742080,18262471680,22794600960,1556055895680,
%U A034094 3577201689600,4464942451200,380428773854896765462278360268800000
%N A034094 (-1)sigma perfect numbers: (-1)sigma(a) = m*a for some integer m, where if a = Product p(i)^r(i) then (-1)sigma(a) = Product_{i} (-1 + Sum_{s=1..r(i)} p(i)^s).
%C A034094 The indices of some terms are 1, so these numbers are fixed points of (-1)sigma where (-1)sigma is A049060.
%e A034094 Factorizations 2^2*5, 2^3*3*13, 2^4*3*7*29, 2^5*3*5*61, 2^3*3^3*13*19, 2^4*3^3*7*19*29, 2^5*3^3*5*19*61, 2^10*3*5*17*409, 2^9*3*5*17*1021, 2^5*3^2*5^2*7*11*29*61, 2^7*3*5*11^2*13*23*131, 2^10*3^3*5*17*19*409, 2^9*3^3*5*17*19*1021, 2^7*3^3*5*11^2*13*19*23*131, 2^10*3^2*5^2*7*11*17*29*409, 2^9*3^2*5^2*7*11*17*29*1021, 2^24*3^3*5^5*7^2*11*17*19*29*61*233*239*467*479*70051.
%t A034094 f[p_, e_] := (p^(e+1)-2*p+1)/(p-1); s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[10^5], Divisible[s[#], #] &] (* _Amiram Eldar_, Jul 07 2022 *)
%o A034094 (PARI) msig(n) = {f = factor(n); for (i=1, #f~, f[i, 1] = (f[i,1]^(f[i,2]+1)-2*f[i,1]+1)/(f[i,1]-1); f[i, 2] = 1;); factorback(f);}
%o A034094 isok(n) = denominator(msig(n)/n) == 1; \\ _Michel Marcus_, Jun 02 2016
%Y A034094 Cf. A034095, A049060.
%K A034094 nonn,more
%O A034094 1,2
%A A034094 _Yasutoshi Kohmoto_
%E A034094 a(1)=1 prepended by _Michel Marcus_, Jun 02 2016
%E A034094 a(10) and a(11) switched and missing term a(13) inserted by _Amiram Eldar_, Jul 07 2022