OFFSET
1,2
COMMENTS
All divisors of numbers from this sequence are in this sequence. Likewise, the product of any terms in this sequence is a number that is also in this sequence.
Subsequence of A017173 (numbers of form 9n+1). - Krizek
For prime numbers, it is enough to verify that the number itself is congruent to 1 mod 9. The first composite term is 361, which is the square of the first prime in this sequence. - Alonso del Arte, May 02 2012
FORMULA
a(n) = 9*k(n) + 1 for k(n) = A211823(n).
EXAMPLE
Number 703 with divisors 1, 19, 37, 703 is in sequence because all divisors have additive digital root of 1.
MATHEMATICA
digitalRoot[n_, b_:10] := FixedPoint[Plus@@IntegerDigits[#, b] &, n]; A211821 = Select[Range[1, 1999, 9], Union[digitalRoot[Divisors[#]]] == {1} &] (* Alonso del Arte, May 02 2012 *)
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Jaroslav Krizek, Apr 26 2012
STATUS
approved