A061712 extension requested

Max relf at unn.ac.ru
Wed Aug 3 11:09:50 CEST 2005


with(combstruct);
a:=proc(n)
  local m,is,s,t,r;
  if n=1 then return 2 fi;
  r:=+infinity;
  for m from 0 to 100 do
    is := iterstructs(Combination(n-2+m),size=n-2);
    while not finished(is) do
      s := nextstruct(is);
      t := 2^(n-1+m)+1+add(2^i,i=s);
      # print(s,t);
      if isprime(t) then r:=min(t,r) fi;
    od;
    if r<+infinity then return r fi;
   od;
   return 0;
end;
seq(a(n),n=1..60);

2, 3, 7, 23, 31, 311, 127, 383, 991, 2039, 3583, 6143, 8191, 73727, 63487, 129023, 131071, 522239, 524287, 1966079, 4128767, 16250879, 14680063, 33546239, 67108351, 201064447, 260046847, 536739839, 1073479679, 5335154687, 2147483647, 8581545983, 16911433727, 32212254719, 51539607551, 206141652991, 266287972351, 1030791102463, 824633720831, 2196875771903, 4260607557631, 15384572854271, 17591112302591, 26388279066623, 70368744046591, 140735340871679, 246290604621823, 560750930165759, 1125762467889151, 2251799813685119, 4362862139015167, 9005000231485439, 17996806323437567, 63050394783186911, 72053195991416831, 108086391056891903, 504403149675560959, 576460752169205759, 1152917106560335871, 2305834213120671743

Max

David Wilson wrote:
> Could some kind soul please compute A061712 for me out to a(60)?
> 
> --------------------------------
> - David Wilson
> 
> 





More information about the SeqFan mailing list