[seqfan] Re: request for help with A036952

Richard Mathar mathar at strw.leidenuniv.nl
Fri Mar 12 19:50:43 CET 2010


One could define A036952 as "Numbers n such that A007088(n) is prime."
or "..such that their binary representation interpreted as a base-10 number is a prime" or...
It seems that this is "Union of A156059 and A065720"

A007088 := proc(n)
	dgs := convert(n,base,2) ;
	add(op(i,dgs)*10^(i-1),i=1..nops(dgs)) ;
end proc:
isA036952 := proc(n)
	isprime( A007088(n)) :
end proc:
A036952 := proc(n)
	if n =1 then
		3;
	else
		for a from procname(n-1)+1 do
			if isA036952(a) then
				return a ;
			end if;
		end do:
	end if;
end proc:
seq(A036952(n),n=1..80) ;





More information about the SeqFan mailing list