1,2

These numbers are (apart from 1) the numbers of elements in finite fields. - Franz Vrabec (franz.vrabec(AT)planetuniqa.at), Aug 11 2004

Numbers whose divisors form a geometrical progresion. The divisors of p^k are 1, p, p^2, p^3, ...p^k. - Amarnath Murthy (amarnath_murthy(AT)yahoo.com), Jan 09 2002

a(n) = A025473(n)^A025474(n). - David Wasserman (wasserma(AT)spawar.navy.mil), Feb 16 2006

a(n) = A117331(A117333(n)). - Reinhard Zumkeller (reinhard.zumkeller(AT)lhsystems.com), Mar 08 2006

These are also precisely the orders of those finite affine planes that are known to exist as of today. (The order of a finite affine plane is the number of points in an arbitrarily chosen line of that plane. This number is unique for all lines comprise the same number of points.) - Peter C. Heinig (algorithms(AT)gmx.de), Aug 09 2006

Except for first term, the index of the second number divisible by n in A002378, if the index equals n. - Mats Granvik (mgranvik(AT)abo.fi), Nov 18 2007

These are precisely the numbers such that lcm(1,...,m-1)<lcm(1,...,m) (=A003418(m) for m>0; here for m=1, the l.h.s. is taken to be 0). We have a(n+1)=a(n)+1 iff n is a Mersenne prime or n+1 is a Fermat prime. - M. F. Hasler (Maximilian.Hasler(AT)gmail.com), Jan 18 2007

The sequence is A000015 without repetitions, or more formally, A000961=Union[A000015]. - Zak Seidov (zakseidov(AT)yahoo.com), Feb 06 2008

Except for a(1)=1, indices for which the cyclotomic polynomial Phi[k] yields a prime at x=1, cf. A020500. - M. F. Hasler (www.univ-ag.fr/~mhasler), Apr 04 2008

Also, {A138929(k) ; k>1} = {2*A000961(k) ; k>1} = {4,6,8,10,14,16,18,22,26,32,34,38,46,50,54,58,62,64,74,82,86,94,98,...} are exactly the indices for which Phi[k](-1) is prime. - M. F. Hasler (www.univ-ag.fr/~mhasler), Apr 04 2008

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 870.

R. Lidl and H. Niederreiter, Introduction to Finite Fields and Their Applications, Cambridge 1986, Theorem 2.5, p. 45.

M. Koecher and A. Krieg, Ebene Geometrie, Springer, 1993

T. D. Noe, Table of n, a(n) for n = 1..10000

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, December 1972 [alternative scanned copy].

Eric Weisstein's World of Mathematics, Link to a section of The World of Mathematics.

Eric Weisstein's World of Mathematics, Link to a section of The World of Mathematics.

Index entries for "core" sequences

m=a(n) for some n <=> lcm(1,...,m-1)<lcm(1,...,m), where lcm(1...0):=0 as to include a(1)=1. a(n+1)=a(n)+1 <=> a(n+1)=A019434(k) or a(n)=A000668(k) for some k (by Catalan's conjecture). - M. F. Hasler (Maximilian.Hasler(AT)gmail.com), Jan 18 2007

readlib(ifactors): for n from 1 to 250 do if nops(ifactors(n)[2])=1 then printf(`%d, `, n) fi: od:

Select[ Range[ 2, 250 ], Mod[ #, # - EulerPhi[ # ] ] == 0 & ]

Select[ Range[ 2, 250 ], Length[FactorInteger[ # ] ] == 1 & ]

max = 0; a = {}; Do[m = FactorInteger[n]; w = Sum[m[[k]][[1]]^m[[k]][[2]], {k, 1, Length[m]}]; If[w > max, AppendTo[a, n]; max = w], {n, 1, 1000}]; a(*Artur Jasinski*)

(MAGMA) [ n : n in [1..1000] | IsPrimePower(n) ];

(PARI) A000961(n, l=-1, k=0)=until(n--<1, until(l<lcm(l, k++), ); l=lcm(l, k)); k print_A000961(lim=999, l=-1)=for(k=1, lim, l==lcm(l, k)&next; l=lcm(l, k); print1(k, ", ")) - M. F. Hasler (Maximilian.Hasler(AT)gmail.com), Jan 18 2007

Cf. A010055, A065515, A095874, A025473.

Cf. indices of record values of A003418; A000668 and A019434 give a member of twin pairs a(n+1)=a(n)+1.

A138929(n) = 2*a(n).

Sequence in context: A087441 A059046 A036116 this_sequence A128603 A096165 A115919

Adjacent sequences: A000958 A000959 A000960 this_sequence A000962 A000963 A000964

nonn,easy,core,nice

njas