[seqfan] Re: A119911

[Harvey P. Dale]

Peter:
	Many thanks.  That did it, and I've now added the corrected Mma program plus some additional terms.
	Thanks again.
	Best,
	Harvey
 

-----Original Message-----
From: SeqFan [mailto:seqfan-bounces at list.seqfan.eu] On Behalf Of Moses, Peter J. C.
Sent: Sunday, May 11, 2014 9:52 AM
To: Sequence Fanatics Discussion list
Subject: [seqfan] Re: A119911

Hi Harvey,

I find that this:

LucasQ[n_]:=IntegerQ[Sqrt[5 n^2+20]]||IntegerQ[Sqrt[5 n^2-20]] Select[Range[102],!LucasQ[DivisorSigma[0,#]]&]

gives the same numbers as A119911.

To make yours give the same numbers, make lnos=LucasL[Range[0,20]].  This will then include 2 as a Lucas number.

Best regards,
Peter.


Peter. J. C. Moses.
Moparmatic Co.
1154 Evesham Road.
Astwood Bank.
Redditch.
Worcs.
B96 6DT.
England.
Tel & Fax: +44 1527 892413.
52°15.713' N. 001°56.256' W

-----Original Message-----
From: Harvey P. Dale
Sent: Sunday, May 11, 2014 2:23 PM
To: Sequence Fanatics Discussion list
Subject: [seqfan] A119911

Are the terms of this sequence correct?  The sequence lists only 38 terms less than 103, whereas I get 64 such terms.
Here is the Mma program I used: 
With[{lnos=LucasL[Range[20]]},Select[Range[300],!MemberQ[lnos,DivisorSigma[0,#]]&]]
Best,
Harvey


_______________________________________________

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


---
This email is free from viruses and malware because avast! Antivirus protection is active.
http://www.avast.com


_______________________________________________

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