[seqfan] Re: A067266

[Robert G. Wilson v]

Et al,

	Without dissecting Harvey's Mmca, I wrote mine from scratch. It is:
m = 0; k = 1; lst = {}; While[k < 10001, m = m + MoebiusMu[k]; If[m == PrimeOmega[k], AppendTo[lst, k]; Print[k]]; k++]; lst

	His output and mine are identical.

Bob.

-----Original Message-----
From: SeqFan [mailto:seqfan-bounces at list.seqfan.eu] On Behalf Of Harvey P. Dale
Sent: Monday, July 14, 2014 9:47 AM
To: Sequence Fanatics Discussion list
Subject: [seqfan] A067266

	My program (below) generates terms that are significantly different (albeit with many terms the same) from those in the above sequence.  Is my program wrong or should I correct the terms of the sequence?

	The program:

With[{nn=1000},Flatten[Position[Thread[{Accumulate[Array[MoebiusMu,nn]],PrimeOmega[Range[nn]]}],_?(First[#]==Last[#]&),{1},Heads->False]]]

	Best,

	Harvey
 

_______________________________________________

Seqfan Mailing list - http://list.seqfan.eu/