n! = a^b

Wouter Meeussen w.meeussen.vdmcc at vandemoortele.be
Thu Sep 3 16:23:09 CEST 1998


hi,

would it make sense to define 2 sequences on the basis of:

a[n] is the best approximation to (n-1)! of the form a[n]^b with a[n]<=n 
b[n] is the best approximation to (n-1)! of the form a^b[n] with a   <=n

with a[n]^b[n] implied, since 
Log[ (n-1)!]/Log[a[n]]=~     b[n]   and
Log[ (n-1)!]/    b[n] =~ Log[a[n]]

of course, the relative error is soon les than a percent, but the absolute error
gets as humongous as n! does.

*****
I really don't know if 'best approximations' have a place in EIS.
do you think they belong there?
*****

a variation is the analogous sequence, where we allow the base b to be
*twice* n:
c[n] is the best approximation to (n-1)! of the form c[n]^d with c[n]<=2n
d[n] is the best approximation to (n-1)! of the form c^d[n] with c   <=2n
(it gives a better fit, at the detriment of 'pertinence').

a[n]  =
{2, 2, 4, 5, 5, 3, 8, 6, 5, 4, 7, 12, 5, 7, 16, 9, 11, 16, 6, 12, 17, 4, 3, 
21, 25, 25, 12, 4, 14, 9, 18, 23, 30, 5, 24, 8,  
  12, 25, 28, 11, 26, 19, 18, 27, 40, 18, 5, 9, 34, 45, 32, 26, 50, 17, 55, 
55, 12, 43, 6, 47, 38, 16, 32, 16, 56, 16, 30,  
  61, 66, 66, 37, 29, 28, 26, 15, 6, 24, 67, 40, 26, 37, 60, 68, 5, 49, 52, 
43, 19, 53, 83, 66, 16, 94, 94, 79, 10, 51, 89, 
  24, 25, 69, 100, 58, 77, 54, 33, 57, 86, 33, 39, 106, 15, 96, 101, 45, 77, 
112, 112, 69, 93, 107, 10, 98, 94, 118, 54, 11, 
  113, 124, 124, 48, 62, 74, 130, 30, 136, 35, 43, 106, 74, 38, 111, 107, 50, 
63, 16, 126, 17, 117, 105, 60, 35, 132, 55, 82, 
  123, 88, 128, 6, 108, 59, 161, 161, 161, 139, 37, 130, 44, 78, 58, 115, 119,
 54, 156, 103, 128, 110, 162, 114, 142, 168,  
  138, 134, 69, 71, 139, 120, 35, 37, 121, 103, 109, 193, 193, 15, 130, 100, 
160, 151}

b[n]=
{0, 1, 1, 2, 3, 6, 4, 6, 8, 11, 9, 8, 14, 13, 10, 14, 14, 13, 22, 17, 16, 35, 
47, 18, 18, 19, 26, 49, 27, 34, 27, 26, 25, 55,  
  29, 46, 40, 32, 32, 46, 35, 40, 42, 38, 35, 46, 85, 64, 41, 39, 44, 48, 41, 
58, 42, 43, 71, 48, 103, 49, 53, 71, 58, 74,  
  52, 77, 64, 54, 54, 55, 65, 71, 73, 76, 93, 143, 82, 63, 73, 84, 77, 69, 68,
 181, 76, 76, 81, 105, 79, 72, 77, 118, 73, 74,  
  78, 150, 89, 79, 113, 113, 87, 81, 93, 88, 97, 112, 98, 90, 116, 112, 89, 
155, 93, 93, 114, 101, 94, 95, 107, 101, 99, 203,  
  103, 105, 101, 122, 205, 105, 104, 105, 132, 125, 121, 108, 156, 109, 152, 
145, 118, 129, 154, 120, 122, 147, 140, 211,  
  122, 210, 126, 130, 149, 173, 127, 156, 143, 132, 143, 133, 363, 140, 162, 
131, 132, 133, 138, 190, 142, 184, 161, 174,  
  150, 150, 181, 144, 158, 152, 158, 147, 159, 153, 149, 156, 158, 184, 184, 
160, 166, 225, 223, 169, 176, 175, 157, 158,  
  309, 173, 184, 168, 171}

and 

c[n]=
 {2, 2, 6, 5, 11, 9, 16, 14, 13, 21, 7, 12, 25, 23, 22, 30, 21, 38, 36, 34, 
44, 32, 40, 50, 48, 46, 36, 54, 14, 35, 61, 59,  
   30, 56, 55, 46, 75, 73, 28, 11, 26, 19, 66, 57, 74, 18, 5, 71, 34, 90, 69, 
77, 61, 96, 55, 83, 92, 91, 81, 47, 66, 16, 32,  
   16, 86, 94, 85, 61, 123, 66, 37, 119, 77, 98, 127, 6, 150, 67, 40, 112, 
103, 142, 95, 5, 49, 149, 126, 19, 53, 83, 115,  
   166, 107, 114, 79, 100, 51, 113, 24, 106, 69, 127, 58, 188, 112, 33, 197, 
160, 194, 149, 132, 15, 96, 101, 45, 225, 124,  
   163, 194, 93, 107, 137, 190, 94, 130, 151, 11, 208, 207, 124, 205, 62, 157,
 253, 30, 164, 35, 261, 106, 163, 38, 229, 107,  
   252, 195, 16, 194, 247, 117, 105, 232, 255, 132, 219, 82, 263, 238, 161, 
216, 108, 59, 245, 214, 174, 139, 319, 231, 44,  
   78, 220, 297, 238, 54, 156, 187, 256, 110, 162, 114, 225, 325, 138, 134, 
69, 271, 139, 200, 249, 328, 248, 267, 109, 193,  
   193, 229, 237, 100, 273, 151},

d[n]=
 {0, 1, 1, 2, 2, 3, 3, 4, 5, 5, 9, 8, 7, 8, 
9, 9, 11, 10, 11, 12, 12, 14, 14, 14, 15, 16,  
   18, 17, 27, 21, 19, 20, 25, 22, 23, 25, 23, 24, 32, 46, 35, 40, 29, 31, 30,
 46, 85, 33, 41, 33, 36, 36, 39, 36, 42, 39,  
   39, 40, 42, 49, 46, 71, 58, 74, 47, 47, 49, 54, 47, 55, 65, 50, 56, 54, 52,
 143, 52, 63, 73, 58, 60, 57, 63, 181, 76, 60,  
   63, 105, 79, 72, 68, 64, 71, 71, 78, 75, 89, 75, 113, 78, 87, 77, 93, 73, 
82, 112, 75, 79, 77, 82, 85, 155, 93, 93, 114,  
   81, 92, 88, 86, 101, 99, 95, 90, 105, 99, 97, 205, 93, 94, 105, 96, 125, 
103, 95, 156, 105, 152, 98, 118, 109, 154, 104,  
   122, 104, 110, 211, 112, 108, 126, 130, 112, 111, 127, 116, 143, 114, 117, 
127, 121, 140, 162, 121, 125, 131, 138, 119,  
   127, 184, 161, 131, 125, 131, 181, 144, 140, 133, 158, 147, 159, 140, 132, 
156, 158, 184, 140, 160, 150, 145, 139, 147,  
   146, 175, 157, 158, 154, 154, 184, 152, 171}

wouter.

NV Vandemoortele Coordination Center
Oils & Fats Applied Research
Prins Albertlaan 79
Postbus 40
B-8870 Izegem (Belgium)
Tel: +/32/51/33 21 11
Fax: +/32/51/33 21 75
vdmcc at vandemoortele.be






More information about the SeqFan mailing list