%I #24 Jan 26 2022 17:57:40
%S 1,2,3,4,5,6,7,8,9,1,1,3,1,1,3,1,1,9,1,2,3,2,1,6,1,2,9,2,1,3,1,1,3,1,
%T 1,9,1,1,3,4,1,6,1,4,9,2,1,12,1,5,3,1,1,9,5,1,3,1,1,6,1,2,9,2,1,6,1,2,
%U 3,7,1,9,1,1,3,1,7,3,1,8,9,2,1,12,1,2,3,8,1,9
%N Greatest common divisor of n and its digit sum.
%C a(n) = 1 iff n is in A339076. - _Bernard Schott_, Jan 17 2022
%H Harry J. Smith, <a href="/A066750/b066750.txt">Table of n, a(n) for n = 1..1000</a>
%t Table[GCD[n,Total[IntegerDigits[n]]],{n,100}] (* _Harvey P. Dale_, Aug 13 2015 *)
%o (PARI) a(n) = gcd(n, sumdigits(n)); \\ _Michel Marcus_, Jan 17 2022
%o (Python)
%o from math import gcd
%o def A066750(n): return gcd(n,sum(int(d) for d in str(n))) # _Chai Wah Wu_, Jan 26 2022
%Y Cf. A007953, A339076.
%K nonn,base
%O 1,2
%A _Reinhard Zumkeller_, Jan 16 2002