|
%I #14 Sep 03 2018 22:58:50
%S 1,1,3,1,1,3,1,1,9,1,11,3,1,1,3,1,1,9,1,1,3,11,1,3,1,1,9,1,1,3,1,1,33,
%T 1,1,9,1,1,3,1,1,3,1,11,9,1,1,3,1,1,3,1,1,9,11,1,3,1,1,3,1,1,9,1,1,33,
%U 1,1,3,1,1,9,1,1,3,1,11,3,1,1,9,1,1,3,1,1,3,11,1,9,1,1,3,1,1,3,1,1,99,1,1,3
%N Highest common factor of n and its 9's complement.
%H Harvey P. Dale, <a href="/A109848/b109848.txt">Table of n, a(n) for n = 1..1000</a>
%e a(55)=11 because the 9's complement of 55 is 99 - 55 = 44 and the greatest common divisor of 44 and 55 is 11.
%p digs := proc(inp::integer) local resul,shiftinp : resul := 1 : shiftinp := iquo(inp,10) : while shiftinp > 0 do resul := resul+1 : shiftinp := iquo(shiftinp,10) : od : RETURN(resul) : end: nin := proc(inp::integer) RETURN(10^digs(inp)-1-inp) : end : for n from 1 to 80 do comp := nin(n) ; #print(n,comp,gcd(n,comp)) ; printf("%d,",gcd(n,comp)) ; od : # _R. J. Mathar_, Mar 27 2006
%t hcf[n_]:=Module[{idn=IntegerDigits[n],c9},c9=FromDigits[PadRight[{}, Length[idn],9]-idn];GCD[n,c9]]; Array[hcf,110] (* _Harvey P. Dale_, Dec 18 2012 *)
%Y Cf. A109847, A084021.
%K base,easy,nonn
%O 1,3
%A _Amarnath Murthy_, Jul 06 2005
%E Corrected and extended by _R. J. Mathar_, Mar 27 2006
%E More terms from _Joshua Zucker_, May 03 2006
|
