// C++ source code that prints sequence A083116 and A083117 // Bo Gyu Jeong, 12~13 Jun 2012 #include int f[30010][2]; const int mm=1000; int main(void){ FILE *out1=fopen("b083116.txt","w"), *out2=fopen("b083117.txt","w"); int t,i,j,k; for(i=1;i<=mm;i++) f[i][0]=f[i][1]=999999; // marks the 0 terms for(i=1;i*10<=mm;i++) f[i*10][0]=f[i*10][1]=0; for(i=1;i*16<=mm;i++) f[i*16][0]=f[i*16][1]=0; for(i=1;i*25<=mm;i++) f[i*25][0]=f[i*25][1]=0; for(t=1;t<=1100;t++){ // 1100 is somehow arbitrary for(i=1;i<=mm;i+=2){ // checks if i is a factor of R_t for(k=0,j=1;j<=t;j++) k=(k*10+1)%i; if(k) continue; for(j=1;j<10;j++){ // R_t ~ 9*R_t if(j*i>mm) break; if(f[j*i][0]=i) ch=1; if(ch){ fprintf(out2,"%d",p/i); p%=i; } } fprintf(out2,"\n"); } return 0; }