Experimental trainning #0173 (Jan 4th 2013)

/*
 *

Instructions:

  1) Download the b-files for the following sequences: A173952, A003952, and A048441.

  2) Delete the offsets column in each b-file: In this step are generated the files

      b173952woOffstes.txt,
      b003952woOffstes.txt, and
      b048441woOffsets.txt

      Suggested method (Unix/Linux):

    Use the "cut" program (GNU file utils / core utils).

    Examples:

      a) In a command line shell try
      ------------------------------------------------------------------------
        seqfan@arithmoz$ cut -d' ' -f2 b173952.txt > b173952woOffsets.txt
      ------------------------------------------------------------------------

      b) Or also
      ------------------------------------------------------------------------
        seqfan@arithmoz$ cut -f2 b173952.txt > b173952woOffsets.txt
      ------------------------------------------------------------------------

    The proper invokation depends on which kind of separators are used inside the b-file
    (either tabs or whitespaces).

  3) Place the files generated in the step-2 at the same directory, and run
  this script there.

*
*/

luckyUpperBound=10; /* Just it, a bound. */

E=readvec("b048441woOffsets.txt");  /* All available terms for A048441 in OEIS */
conservativeUpperBound=min(20,#E);
A=readvec("b173952woOffsets.txt");  /* All the known terms for A173952 */
B0=readvec("b003952woOffsets.txt");
B=vector(#B0-1,k,B0[k+1]);          /* Terms in A003952 excepting the first. */

C=vector(#A-1,k,A[k+1]-sum(j=1,k-1,A[j]));  /* A particular kind of "First differences" for A173952... */
C/=2^5;                                     /* ... divided by 2^5. */

D=vector(min(min(#B,#C),#E),k,B[k]-C[k]);   /* More differences: This time between A003952 and the previous definition for A173952 */

/*==================
      Sumary
==================*/

for(i=1,luckyUpperBound,print(D[i]));
print("\nWith the exception of the two zeroes at the beginning, these \nnumbers listed previously are the first ",luckyUpperBound-2," terms of A048441.");
printf("\n\nNow these other numbers liste below are the differences against the next terms of A048441.\n");
for(j=luckyUpperBound+1,conservativeUpperBound,print(E[j-2]-D[j]));
quit();