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

 

Logo

Many excellent designs for a new banner were submitted. We will use the best of them in rotation.

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 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 [From Stefan Steinerberger, Aug 29 2008]

PROG

(PARI code from M. F. Hasler, Aug 25 2008) 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", "))

CROSSREFS

Cf. A063655, A142591.

Sequence in context: A117287 A121615 A097605 * A187731 A182140 A190223

Adjacent sequences:  A143575 A143576 A143577 * A143579 A143580 A143581

KEYWORD

nonn,changed

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 | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement .

Last modified April 19 05:22 EDT 2014. Contains 240738 sequences.