Number of distinct exponentials

Nick Hobson nickh at qbyte.org
Fri Dec 22 10:42:58 CET 2006


Hi Seqfans,

What do you think of this sequence: 1, 3, 7, 11, 19, 28, 40, 50, 60, 76,  
96, 115, 139, 163, 189, 207, 239, 270, 306, 340, 378, 417, 461, 503, 539,  
585, 621, 670, 726, 779, 839, 881, 941, ... ?

It is the number of distinct terms i^j for 1 <= i,j <= n.  (Cf. A027424).   
I don't have a formula for the sequence, or any nice bounds, such as those  
for A027424.

Similarly, the corresponding sequence for 2 <= i,j <= n may be of  
interest?  (1, 4, 8, 15, 23, 34, 44, 54, 69, 88, 106, 129, 152, 177, 195,  
226, 256, 291, 324, 361, 399, 442, 483, 519, 564, 600, 648, 703, 755, 814,  
856, 915, 976, ... .)  A complementary sequence is n^2 - a(n): the number  
of duplicates (0, 0, 1, 1, 2, 2, 5, 10, 12, 12, ...), but is it worth  
adding both?

PARI script for 1 <= i,j <= n:
g(lim)={local(z,x,c,f); z=listcreate(lim^2); vector(lim, m, for(i=1, m,  
for(j=1, m, x=factor(i); x[,2]*=j; c=Str(x); f=setsearch(z, c, 1); if(f,  
listinsert(z, c, f)))); #z)};
g(100)

Nick






More information about the SeqFan mailing list