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”).
%I #36 Nov 21 2024 21:25:46
%S 1,2,3,4,5,6,7,9,10,11,12,13,15,17,19,23,25,29,30,31,35,36,37,40,41,
%T 43,47,48,49,50,53,59,60,61,65,67,71,73,76,79,83,89,97,101,103,107,
%U 109,113,120,121,127,130,131,132,136,137,139,140,149,150,151,157,163,167,169,173,175,179,180
%N Integers m such that m divides sigma_2(m) - k where k is some divisor of m.
%e 2 is in this sequence because 2 divides A001157(2) - 1 = 5 - 1 = 4.
%t Select[Range@ 180, Function[n, Total@ Boole@ Map[Divisible[ DivisorSigma[2, n] - #, n] &, Divisors@ n] > 0]] (* _Michael De Vlieger_, Mar 19 2017 *)
%o (Magma) [[n: k in [1..n] | Denominator(n/k) eq 1 and
%o Denominator(((DivisorSigma(2, n))-k)/n) eq 1]: n in [1..100]];
%o (PARI) isok(n) = fordiv(n, d, if (!((sigma(n, 2) - d) % n), return (1))); \\ _Michel Marcus_, Mar 18 2017
%Y Supersequence of A166684.
%Y Cf. A001157 (sigma_2(n): sum of squares of divisors of n), A205523 (integers n such that n divides sigma_1(n) - i where i is some divisor of n), A284082.
%K nonn
%O 1,2
%A _Juri-Stepan Gerasimov_, Mar 18 2017