# Code for Pari/GP # \\K(i,j): Kimberling Expulsion Recurrence # K(i,j) = # { # my(i1,j1); # i1=i; # j1=j; # while(j1<(2*i1-3), # if(j1%2, # j1=i1+((j1-1)/2), # j1=i1-((j1+2)/2) # ); # i1--; # ); # return(i1+j1-1); # } # addhelp(K,"K(i,j): Kimberling Expulsion Array"); # n=0; # A038807 = 2; # max_iter = 75; # while(n<max_iter, print(n" "A038807); write("b175785.txt",n" "A038807); n++; A038807=K(A038807, A038807)) 0 2 1 3 2 5 3 10 4 9 5 20 6 46 7 83 8 12 9 24 10 23 11 36 12 79 13 124 14 172 15 56 16 119 17 61 18 169 19 17 20 42 21 84 22 232 23 285 24 596 25 1186 26 3190 27 6857 28 14225 29 12495 30 30482 31 45827 32 79090 33 144112 34 423486 35 1087497 36 2443796 37 628733 38 871389 39 1199242 40 2787410 41 7975876 42 14345987 43 576445 44 1425599 45 1529574 46 2834470 47 4361479 48 2859775 49 5101756 50 1444243 51 3091522 52 9240735 53 2837036 54 6761707 55 10649036 56 17944639 57 34251976 58 101333993 59 125551569 60 216748684 61 528889352 62 64112048 63 19597024 64 47920943 65 106574679 66 232784460 67 238739142 68 510423893 69 252813731 70 380987480 71 274339926 72 257662655 73 104541876 74 253842952