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

 

Logo

Annual Appeal: Today, Nov 11 2014, is the 4th anniversary of the launch of the new OEIS web site. 70,000 sequences have been added in these four years, all edited by volunteers. Please make a donation (tax deductible in the US) to help keep the OEIS running.

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 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 | 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 December 22 14:10 EST 2014. Contains 252364 sequences.