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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A247824 Least positive integer m such that m + n divides prime(m) + prime(n). 23
1, 5, 5, 5, 2, 2, 38, 16, 40, 12, 13, 1, 11, 1, 11, 4, 35, 38, 35, 35, 38, 35, 35, 36, 31, 31, 33, 33, 36, 36, 25, 25, 2, 25, 4, 3, 4, 6, 6, 8, 222, 8, 95, 223, 99, 98, 95, 88, 222, 94, 93, 94, 95, 92, 226, 88, 83, 92, 225, 92 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,2

COMMENTS

Conjecture: a(n) exists for any n > 0. Moreover, a(n) < n*(n-1) for all n > 2.  - Zhi-Wei Sun, Sep 25 2014

A247869(n) = (prime(a(n)) + prime(n)) / (a(n) + n). - Reinhard Zumkeller, Sep 27 2014

I have verified the conjecture for n up to 10^5, and noted that max{a(n): n=1..10^5} = a(79276) = 3141281384 > 3*10^9. - Zhi-Wei Sun, Oct 08 2014

I would like to offer 500 US dollars as the prize for the first proof of the above conjecture. - Zhi-Wei Sun, Feb 24 2018

LINKS

Zhi-Wei Sun, Table of n, a(n) for n = 1..100000

Zhi-Wei Sun, A new theorem on the prime-counting function, arXiv:1409.5685, 2014.

Zhi-Wei Sun, m+n divides prime(m)+prime(n) for some n>0, a message to Number Theory List, Sept. 27, 2014.

Zhi-Wei Sun, A new theorem on the prime-counting function, Ramanujan J. 42(2017), no.1, 59-67.

EXAMPLE

a(2) = 5 since 5 + 2 = 7 divides prime(5) + prime(2) = 11 + 3 = 14.

a(10409) = 69804276 since 69804276 + 10409 = 69814685 divides prime(10409) + prime(69804276) = 109481 + 1396184219 = 1396293700 = 20*69814685.

a(35980) = 180302246 since 35980 + 180302246 = 180338226 divides prime(35980) + prime(180302246) = 427727 + 3786675019 = 3787102746 = 21*180338226.

a(79276) = 3141281384 since 79276 + 3141281384 = 3141360660 divides prime(79276) + prime(3141281384) = 1010431 + 75391645409 = 75392655840 = 24*3141360660.

MATHEMATICA

Do[m=1; Label[aa]; If[Mod[Prime[m]+Prime[n], m+n]==0, Print[n, " ", m]; Goto[bb]]; m=m+1; Goto[aa]; Label[bb]; Continue, {n, 1, 60}]

lpi[n_]:=Module[{k=1, p=Prime[n]}, While[!Divisible[p+Prime[k], k+n], k++]; k]; Array[lpi, 60] (* Harvey P. Dale, Apr 23 2015 *)

PROG

(PARI) a(n) = {m = 1; while ((prime(m) + prime(n)) % (m + n), m++); m; } \\ Michel Marcus, Sep 25 2014

(PARI) a(n)=my(p=prime(n), m); forprime(q=2, , if((p+q)%(n+m++)==0, return(m))) \\ Charles R Greathouse IV, Sep 25 2014

(Haskell)

import Data.List (genericIndex)

a247824 n = genericIndex a247824_list (n - 1)

a247824_list = f ips where

   f ((x, p) : xps) = head

     [y | (y, q) <- ips, (p + q) `mod` (x + y) == 0] : f xps

   ips = zip [1..] a000040_list

-- Reinhard Zumkeller, Sep 27 2014

CROSSREFS

Cf. A000040, A247600, A247793.

Sequence in context: A299784 A021647 A181668 * A094851 A173602 A054245

Adjacent sequences:  A247821 A247822 A247823 * A247825 A247826 A247827

KEYWORD

nonn

AUTHOR

Zhi-Wei Sun, Sep 24 2014

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 November 20 06:06 EST 2018. Contains 317385 sequences. (Running on oeis4.)