A070019 Duplicate of A056075.
1, 4, 56, 7192, 7232, 7912, 10792, 17272, 30592, 114256, 2154584, 3428368, 44375136
Offset: 1
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.
Fixed points: 56, 7192, 7232, 7912, 10792, ... sigma(66) = 144, tau(66) = 8 and 144 - 8 - 66 = 70. sigma(70) = 144, tau(70) = 8 and 144 - 8 - 70 = 66.
with(numtheory); P:=proc(q)local a,n; for n from 1 to q do a:=sigma(n)-tau(n)-n; if sigma(a)-tau(a)-a=n then print(n); fi; od; end: P(10^6);
with(numtheory): A055681:=n->`if`(sigma(n)-phi(n) mod n=0,n,NULL): seq(A055681(n), n=1..10^5); # Wesley Ivan Hurt, Sep 13 2014
Do[If[Mod[DivisorSigma[1, n]-EulerPhi[n], n]==0, Print[n]], {n, 1, 10^9}]
for(n=1,10^8,if((sigma(n)-eulerphi(n))%n==0,print1(n,", "))) \\ Derek Orr, Sep 13 2014
11 is a term as its abundance is -10 and its number of divisors is 2, the former number being divisible by the latter.
Select[Range[180], Divisible[DivisorSigma[1,#]-2#, DivisorSigma[0,#]]&]
for(n=1, 180, ((sigma(n)-2*n)%numdiv(n)==0) && print1(n ", "))
isok(n) = !((sigma(n)-2*n)%numdiv(n)); \\ Michel Marcus, Apr 09 2018
4 is a term since 2^4 - (2*4+1) = 16 - 9 = 7 is a prime. 7 is a term since 2^7 - (2*7+1) = 128 - 15 = 113 is a prime.
Select[Range[2400], PrimeQ[2^# - 2*# - 1] &]
Comments