|
|
A072394
|
|
Numbers n such that sigma(n)=reversal(n)-n.
|
|
2
|
|
|
1563, 1633, 18673, 32207, 1405313, 1567563, 1656833, 193613415, 325933027, 376491249, 2287850446, 2432416646, 13823276223, 14055445313, 19087920283, 23804849568, 36303512827, 148868530953
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
COMMENTS
|
If (58*1000^n-169)/111 is prime then (58*1000^n-169)/37 is in the sequence (the proof is easy). Next term is greater than 12*10^8. - Farideh Firoozbakht, Jan 29 2006
If p = 156/101*(10^(4n)-1)-1 is prime then 91*p is in the sequence (the proof is easy).
A178321 gives numbers n such that (58*1000^n-169)/111 = 58/111*(10^(3n)-1)-1 is prime and A178322 gives numbers n such that 156/101*(10^(4n)-1)-1 is prime. (End)
|
|
LINKS
|
|
|
EXAMPLE
|
reverse(1563) - 1563 = 3651 - 1563 = 2088 = sigma(1563), so 1563 is a term of the sequence.
376491249 is in the sequence because sigma(376491249)=565703424 =942194673-376491249=reversal(376491249)-376491249.
|
|
MATHEMATICA
|
Select[Range[10^6], FromDigits[Reverse[IntegerDigits[n]]] - # == DivisorSigma[1, # ] &]
Do[If[DivisorSigma[1, n]==FromDigits[Reverse[IntegerDigits[n]]]- n, Print[n]], {n, 1200000000}] (* Farideh Firoozbakht *)
|
|
CROSSREFS
|
|
|
KEYWORD
|
base,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|