[seqfan] Re: A154676
Peter Pein
petsie at dordos.net
Fri Jan 16 02:24:53 CET 2009
zak seidov schrieb:
> Terms in A154676 are of the form 3708 n^2,
> and 1262 values of n<10^5 are in
> http://zak08.livejournal.com/4070.html
>
> In general A154676 is rather raw and needs some editing -
> sorry Vlad!
>
>
>
>
>
>
>
> _______________________________________________
>
> Seqfan Mailing list - http://list.seqfan.eu/
>
Dear group,
could someone please check, if the following Mathematica-code gives correct
results. It seems slightly more efficient than that published in OEIS.
sorry for the unformatted output, but it is nearly 2:30am and I want to go to bed.
Have a nice day/night/time
Peter
Timing[lst = {};
Do[If[PrimeQ[n - 1] &&
PrimeQ[n + 1],
If[IntegerQ[Sqrt[n*103]],
AppendTo[lst, n]]],
{n, 9!, 2*11!, 6}]; lst]
Out[1]=
{139.704731*Second, {2317500, 12047292, 26163648, 43250112, 47347452, 61704828}}
In[2]:=
Timing[lst = {}; Do[n = k^2/103;
If[PrimeQ[n - 1] && PrimeQ[
n + 1], AppendTo[lst, n]],
{k, 103, 10300000, 103}]; lst]
Out[5]=
{1.096068*Second, {2317500, 12047292,
26163648, 43250112, 47347452,
61704828, 168228252, 333720000,
351755712, 426127068, 513127872,
840143808, 979638768, 998790588,
1089276912, 1330434108, 1357220700,
1388809152, 1694467008, 1927570428,
1986835392, 2035992348, 2136108348,
2858437872, 3070594800, 3241626300,
3903322608, 5136799932, 5250898800,
5303982528, 5812304832, 5830889328,
6598593648, 7298860572, 7371874800,
8232130800, 8566963200, 9779390208,
10640610288, 11097673200,
11136195612, 11174784768,
11316849372, 11708859132,
13428184572, 14027456700,
15265246428, 15719995548,
15919126272, 16445855088,
18012040128, 18422979228,
18821841372, 18955744668,
19496104092, 19547152128,
20095817472, 20407778928,
20774770812, 21198194748,
21429332928, 23603381532,
24186972528, 24452865792,
24720212592, 24950523888,
25085365308, 25667372988,
26572747932, 27211829148,
31356568512, 32028669468,
33039114300, 34879261212,
37780307712, 38877059952,
38949132348, 42663002112,
43904503548, 45034494300,
45215589312, 46152675072,
46336002300, 46362221568,
50872336128, 51423196608,
52227213372, 54193639932,
55764181872, 57707696700,
58647689532, 59061324348,
61936833852, 62361889308,
64539905472, 64570848732,
67702625532, 70293727728,
72670603932, 74620518768,
74886901488, 76326840252,
80486889948, 81665915292,
81700722288, 84226315248,
84650985072, 88629806268,
96293852928, 96596425728,
97316938332, 100573062912,
101075782428, 101812780800,
108726810300, 110379506688,
110500925148, 110703437568,
111068427132, 112943310588,
114256832508, 114503963292,
118201087728, 118368630000,
119839686300, 120980708352,
124866191772, 125167618800,
125771563008, 126723554928,
126853650108, 127114040700,
129732630300, 134334762588,
135453002688, 139608039168,
144936464832, 146377338012,
149468742108, 150506388828,
152640027648, 160055449200,
160201651932, 160689476412,
163927717308, 164865926592,
165311275932, 169246231488,
172720957500, 175006105200,
176282120700, 189085310748,
189244217088, 197599320000,
199118057472, 202392073152,
203379483888, 204314244732,
206854343388, 209800434672,
212655019932, 214681634748,
215472358332, 219676043772,
222079984668, 230129173872,
247686720732, 251458854300,
255628811772, 259398687168,
260392075200, 278726281200,
281885723388, 285974026992,
286886536128, 289829445948,
292655590848, 295826553792,
296356708800, 297750623868,
305643917628, 305711251200,
309426018300, 310103840700,
314186348700, 323400343068,
331694556912, 343797012828,
345584546928, 347233171932,
360417158748, 363788264700,
364817368188, 364890931200,
376980824412, 395370231408,
396596422512, 400671740700,
401520171888, 405465532092,
407795342652, 410756017500,
413727401052, 417574729152,
419307892848, 420649614108,
421914368412, 423339486300,
426515017500, 430261428672,
439981820892, 443543429052,
453493891548, 458840438208,
464052377052, 477673830000,
484516057788, 487402602300,
489189750768, 489360125952,
498260512512, 504815366592,
515602086912, 516214414908,
518228908272, 523590694812,
525089880000, 529688927232,
530398252800, 532795730652,
548556213852, 562078455552,
563997512412, 564454912752,
564820966512, 571154071068,
575118646128, 581509061532,
588962438172, 606854428092,
608373343548, 610845630300,
636621362172, 649316367612,
654922162800, 655316456688,
663624037872, 685226667888,
692605784412, 703697476608,
719621797788, 721586070000,
723242270448, 729990107388,
743579779068, 755070793200,
761965170300, 769959154800,
786828289212, 803552674428,
805081794300, 806502992832,
809239589532, 820563254208,
823544608572, 827860553712,
835301460348, 841099437468,
845237131632, 849946276800,
876878533872, 887860858608,
890272352700, 900297668592,
907476401088, 913984920000,
922621133808, 938716987068,
944862804252, 949010702832,
954833022972, 959957074800,
966171920112, 975411503388,
981676095228}}
More information about the Seqfan
mailing list