本题的算法中涉及的三个函数：

double bbp(int n,int k,int l) 其中n为十六进制位第n位，k取值范围为0到n+7，用来计算16nS1，16nS2，16nS3，16nS4小数部分的每一项。返回每一项的小数部分。

void pi(int m,int n,int p[]) 计算从n位开始的连续m位的十六进制数字。其中p为存储十六进制数字的数组。

void div(int p[])

void add(int a[],int b[]) 这两个函数都是为最后把十六进制数字转换为十进制数字服务的。

最后把1000个数字分别存储在整型数组r[]中，输出就是按顺序输出该数组。