Re: FindInstance for sum of primes

*To*: mathgroup at smc.vnet.net*Subject*: [mg115586] Re: FindInstance for sum of primes*From*: leigh pascoe <leigh at evry.inserm.fr>*Date*: Sat, 15 Jan 2011 04:46:59 -0500 (EST)

Thank you to all who replied to my question. It seems I should have been more explicit with my question. I was originally sent the observation that 2011 = 157 + 163 + 167 + 173 + 179 + 181 + 191 + 193 + 197 + 199 + 211 =Prime[37]+.... Prime[47] which I took to be of interest because the year (ignoring the 1000's) is the sum of that same number of consecutive primes. I wondered whether this was a common phenomenon and tried to answer the question with Mathematica i.e. if year (mod 1000)=N where N can be expressed as the sum of N consecutive primes. So the search was not restricted to 11 consecutive primes but to cases that followed the same rule. The fastest way to find examples seems to be by inspection with a little help from Ma. For the first 1000 primes, generate the lists list = Total /@ Partition[Prime[Range[1000]], #, 1] & /@ Range[100]; which gives lists of the sums of 1,2,.. 100 primes then inspect the lists for cases where the sum equals a year (mod 1000) In[8]:= MemberQ[Mod[list[[#]], 1000], #] & /@ Range[100] Out[8]= {True, True, True, True, False, True, True, True, True, True, \ True, True, True, True, True, True, True, True, True, True, True, \ True, True, True, True, True, True, True, True, True, True, True, \ True, True, True, True, True, True, False, True, True, False, True, \ True, True, True, True, True, True, True, True, True, True, False, \ False, True, True, True, True, True, True, True, True, True, True, \ True, True, True, True, False, True, True, True, True, True, True, \ True, False, True, True, True, True, False, False, True, True, True, \ False, True, True, True, True, True, True, True, True, True, True, \ True, True} and we see that such years are quite common. We can use the function f[n_, m_] := Sum[Prime[i], {i, n, n + m}]; to find particular examples by trial and error Some examples would be {3001, 4001, 7001...} are the sums of 1 (sic) consecutive primes {1002,4002,10002} are the sums of 2 consecutive primes etc. also near in time 5014 is the sum of 14 consecutive primes 6023 is the sum of 23 consecutive primes No solution was observed for 5 consecutive primes in the range scanned. However a solution for 5 consecutive primes exists when the number of primes examined is raised to 10000. I guess that responds to my curiosity even if it is not a very elegant way of tackling the problem.. In passing I noticed that the year 1212 was the sum of 12 consecutive primes. Are there any other years of the form NN that are the sum of N consecutive primes? Leigh Le 14/01/2011 12:16, DrMajorBob a =C3=A9crit : > Here are the first 100 sums of 11 consecutive primes: > > Clear[f]; > f[1] = Total@Prime@Range@11 > f[n_] := f[n] = f[n - 1] + Prime[n + 10] - Prime[n - 1] > Array[f, {100}] // Timing > > 160 > > {0.0009, {160, 195, 233, 271, 311, 353, 399, 443, 491, 539, 583, 631, > 677, 725, 779, 833, 883, 931, 979, 1025, 1081, 1139, 1197, 1253, > 1313, 1367, 1423, 1483, 1543, 1607, 1673, 1727, 1787, 1843, 1901, > 1951, 2011, 2077, 2141, 2203, 2263, 2323, 2383, 2443, 2507, 2573, > 2643, 2703, 2757, 2811, 2865, 2925, 2993, 3063, 3125, 3185, 3253, > 3321, 3397, 3469, 3541, 3617, 3691, 3757, 3825, 3895, 3967, 4033, > 4097, 4159, 4229, 4297, 4369, 4435, 4501, 4565, 4631, 4699, 4763, > 4825, 4883, 4943, 5009, 5069, 5135, 5199, 5265, 5337, 5403, 5483, > 5567, 5657, 5741, 5823, 5903, 5981, 6065, 6149, 6227, 6305}} > > Bobby > > On Thu, 13 Jan 2011 02:27:33 -0600, Peter Pein<petsie at dordos.net> wrote: > >> On 12.01.2011 01:25, leigh pascoe wrote: >>> Dear Mathgroup, >>> >>> 2011 is the sum of 11 consecutive primes. I want to check if any other >>> years have a similar property. >>> >>> Now define >>> >>> f[n_, m_] := Sum[Prime[i], {i, n, n + m}]; >>> eq = Mod[f[n, m] - year, 1000] == 0 >>> >>> and we see that >>> >>> In[64]:= Mod[f[37, 10] - 2011, 1000] >>> >>> Out[64]= 0 >>> >>> but >>> >>> In[65]:= FindInstance[eq, {n, m, year}, Integers] >>> >>> During evaluation of In[65]:= FindInstance::exvar: The system contains a >>> nonconstant expression i independent of variables {n,m,year}.>> >>> >>> Out[65]= FindInstance[Mod[-year + \!\(\*UnderoverscriptBox[\(\[Sum]\), >>> \(i = n\), \(m + n\)]\(Prime[i]\)\), 1000] == 0, {n, m, year}, Integers] >>> >>> Apparently FindInstance doesn't like the dummy variable "i". How can we >>> perform this search in Mathematica?? >>> >>> LP >>> >> >> Hello Leigh, >> >> I can't imagine how to solve your problem using Reduce/FindInstance, but >> you can construct a table containing sums of k conecutive primes: >> >> In[1]:= pSumTable = >> Mod[ >> With[{pr = Prime[Range[1000]]}, >> Table[Total /@ Partition[pr, k, 1], {k, 100}] >> ], >> 10000]; >> >> say, you want to find 2011. Then use >> >> In[2]:= LengthAndStart = Position[pSumTable, 2011, 2] >> >> Out[2]= {{1, 305}, {3, 121}, {3, 551}, {9, 494}, {11, 37}, {11, >> 737}, {29, 170}, {35, 675}, {49, 342}, {49, 395}, {59, 16}, {61, >> 614}, {83, 241}, {91, 6}, {95, 606}} >> >> {11, 37} is the case you mentioned in your post. >> >> to verify the results (without the help of pSumTable): >> >> In[3]:= Total[Prime[Range[#2, #2 + #1 - 1]]]& @@@ LengthAndStart >> >> Out[3]= {2011, 2011, 12011, 32011, 2011, 62011, 32011, 182011, \ >> 122011, 142011, 12011, 292011, 152011, 22011, 462011} >> >> >> Cheers, >> Peter >> >

**Re: Caution on BibTex references with Mathematica**

**Re: NSolve unable to find all possible roots**

**Re: FindInstance for sum of primes**

**Cubic equations again...**