login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A093321 a(1) = 1; for m >= 2, a(m) =sum_{primes p<m which do not divide m} a(pi(p)), where pi(p) is the order of p among the primes = the number of primes <= p. 2

%I

%S 1,0,1,0,1,1,2,1,2,0,2,2,3,3,3,3,4,5,6,5,7,5,7,8,8,7,9,8,9,7,9,10,10,

%T 8,10,10,11,11,12,11,13,15,16,17,18,16,19,21,22,20,20,20,22,24,23,24,

%U 24,24,25,27,29,31,34,33,32,32,34,37,38,38,40,44,45,49,51,50,51,50,52,55

%N a(1) = 1; for m >= 2, a(m) =sum_{primes p<m which do not divide m} a(pi(p)), where pi(p) is the order of p among the primes = the number of primes <= p.

%t PrimeFactors[n_Integer] := Flatten[ Table[ # [[1]], {1}] & /@ FactorInteger[n]]; a[1] = 1; a[n_] := a[n] = (Plus @@ (a[ # ] & /@ Complement[ Range[ PrimePi[n]], PrimePi[ PrimeFactors[n]]])); Table[ a[n], {n, 80}] (* _Robert G. Wilson v_, May 04 2004 *)

%Y Cf. A093320, A094163.

%K nonn,easy

%O 1,7

%A _Leroy Quet_, Apr 26 2004

%E More terms from _Robert G. Wilson v_, May 04 2004

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 18 07:05 EDT 2019. Contains 325134 sequences. (Running on oeis4.)