login
This site is supported by donations to The OEIS Foundation.

 

Logo

Annual Appeal: Please make a donation (tax deductible in USA) to keep the OEIS running. Over 5000 articles have referenced us, often saying "we discovered this result with the help of the OEIS".

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A143578 A positive integer n is included if j+n/j divides k+n/k for every divisor k of n, where j is the largest divisor of n that is <= sqrt(n). 2
1, 2, 3, 5, 7, 11, 13, 15, 17, 19, 23, 29, 31, 35, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 95, 97, 101, 103, 107, 109, 113, 119, 127, 131, 137, 139, 143, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 199, 209, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263, 269, 271, 277, 281, 283, 287, 293 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

This sequence trivially contains all the primes.

There is no term <= 5*10^7 with bigomega(n)>2, i.e., with more than 2 prime factors. - M. F. Hasler, Aug 25 2008. Compare A142591.

If it is always true that the terms have <= 2 prime divisors, then this sequence is equal to {1} U primes U {pq: p, q prime, p+q | p^2-1}. - David W. Wilson, Aug 25 2008

LINKS

Table of n, a(n) for n=1..70.

EXAMPLE

The divisors of 35 are 1,5,7,35. The sum of the two middle divisors is 5+7 = 12. 12 divides 7 + 35/7 = 5+35/5 = 12, of course. And 12 divides 1 + 35/1 = 35 +35/35 = 36. So 35 is in the sequence.

MATHEMATICA

a = {}; For[n = 1, n < 200, n++, b = Max[Select[Divisors[n], # <= Sqrt[n] &]]; If[ Length[Union[Mod[Divisors[n] + n/Divisors[n], b + n/b]]] == 1, AppendTo[a, n]]]; a (* Stefan Steinerberger, Aug 29 2008 *)

PROG

(PARI) isA143578(n)={ local( d=divisors(n), j=(1+#d)\2, r=d[ j ]+d[ 1+#d-j ]); for( k=1, j, ( d[k]+d[ #d+1-k] ) % r & return ); 1 }

for(n=1, 300, isA143578(n) && print1(n", ")) \\ M. F. Hasler, Aug 25 2008

CROSSREFS

Cf. A063655, A142591.

Sequence in context: A117287 A121615 A097605 * A187731 A182140 A240960

Adjacent sequences:  A143575 A143576 A143577 * A143579 A143580 A143581

KEYWORD

nonn

AUTHOR

Leroy Quet, Aug 24 2008

EXTENSIONS

More terms from M. F. Hasler, Aug 25 2008 and Stefan Steinerberger, Aug 29 2008

STATUS

approved

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

License Agreements, Terms of Use, Privacy Policy .

Last modified December 3 21:14 EST 2016. Contains 278745 sequences.