login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A268983 Numbers n such that Sum_{d|n} sigma(n/d)^d/d is an integer, where sigma is A000203. 2
1, 39, 793, 25930965, 76147347 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
Conjecture: This sequence contains only the three numbers listed.
This conjecture is disproven. The number 2408321608150261253061174553 belongs to this sequence, but it may not be a(4). - Gevorg Hmayakyan, May 11 2017
a(6) > 10^8. - Hiroaki Yamanouchi, Aug 10 2017
Numbers n such that n divides Sum_{d|n} d*sigma(d)^(n/d). - Chai Wah Wu, Oct 02 2017
LINKS
MAPLE
a := proc (n) options operator, arrow; add(numtheory:-sigma(n/d)^d/d, d in numtheory:-divisors(n)) end proc:
seq(`if`(type(a(n), 'integer'), n, NULL), n = 1 .. 1000)
MATHEMATICA
Select[Range[10^5], IntegerQ@ Sum[DivisorSigma[1, #/d]^d/d, {d, Divisors@ #}] &] (* Michael De Vlieger, Mar 30 2016 *)
PROG
(PARI) isok(n) = denominator(sumdiv(n, d, sigma(n/d)^d/d)) == 1; \\ Michel Marcus, Feb 17 2016
(Python)
from __future__ import division
from sympy import divisors, divisor_sigma
A268983_list = [n for n in range(1, 10000) if not sum(d*pow(int(divisor_sigma(d)), n//d, n) % n for d in divisors(n, generator=True)) % n] # Chai Wah Wu, Oct 02 2017
CROSSREFS
Sequence in context: A232354 A232067 A247222 * A293918 A264458 A004330
KEYWORD
nonn,more
AUTHOR
Gevorg Hmayakyan, Feb 16 2016
EXTENSIONS
a(4)-a(5) from Hiroaki Yamanouchi, Aug 10 2017
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 23 16:40 EDT 2024. Contains 371916 sequences. (Running on oeis4.)