A sequence unique including differences

franktaw at netscape.net franktaw at netscape.net
Fri May 30 02:12:48 CEST 2008


I'm surprised that the following was not already in the OEIS:

%I A140778
%S A140778 
1,3,7,12,18,8,17,28,13,27,43,19,39,60,22,45,70,26,55,85,31,63,96,34,69,10
5,37,77,
118,42,88,135,48,97,147,52,103,156,56,113,171,59,120,184,65,131,198,71,14
3,216,74,
149,227,79,159,240,82,165,249,86,175,265,91,183,276,94,192,291,101,203,30
7,106,
213,321,109,219,330,115,229,345,117,238,360,123,247,372,126,254,383,130,2
62,395,
134,270,407,138,277,417,144,285
%N A140778 a(n) is the smallest positive integer such that no number 
occurs
twice in the sequence and its absolute first differences.
%C A140778 This sequence and its first differences include every 
positive
integer (exactly once).
%e A140778 For a(5), the sequence to that point is [1,3,7,12], with 
absolute
differences [2,4,5].  The next number cannot be 6, because then 6 would 
be in
both the sequence and the first differences.  Since all values smaller 
than 6
are taken, the difference must be positive and at least 6.  A 
difference of 6
works, a(5) = 18.
%o A140778 (PARI) IsInList(v, k) = 
for(i=1,#v,if(v[i]==k,return(1)));return(0)
IsInDiff(v, k) = for(i=2,#v,if(abs(v[i]-v[i-1])==k,return(1)));return(0)
NextA140778(v)={
 local(i,d);
 if(#v==0,return(1));
 i=2;
 while(1,
  d=abs(i-v[#v]);
  if(!(i==d || IsInList(v,i) || IsInDiff(v,i) || IsInList(v,d) ||
IsInDiff(v,d)), return(i));
  i++)
}
v=[];for(i=1,100,v=concat(v,NextA140778(v)));v
%Y A140778 Cf. A140779, A081145.
%O A140778 1
%K A140778 ,easy,nonn,
%A A140778 Franklin T. Adams-Watters (FrankTAW at Netscape.net), May 29 
2008

Question: can an upper bound for a(n)/n be established?  If so, what is 
the least upper bound?  Up to n=100, the maximum is n=74, a(n) = 321, 
a(n)/n ~ 4.338.

The same question can be asked for A081145; the maximum for the values 
in the database is n=51, a(n) = 115, a(n)/n ~ 2.255.

Franklin T. Adams-Watters





More information about the SeqFan mailing list