|
%I #19 Oct 11 2019 04:05:35
%S 4,16,27,308,440,528,594,5445,5808,6534,8085,9702,11550,13860,14784,
%T 16500,16632,18711,19800,21120,22275,23760,25344,26730,28512,32076,
%U 37268,53240,63888,94864,135520,152460,162624,181500,182952,193600
%N Composite numbers divisible by the palindromic sum of their palindromic prime factors (counted with multiplicity).
%C The subsequence of numbers k of A046355 such that A262049(k) divides k. - _R. J. Mathar_, Sep 09 2015
%H Giovanni Resta, <a href="/A046366/b046366.txt">Table of n, a(n) for n = 1..1387</a> (terms < 10^10)
%e 1041714 = 2 * 3^3 * 101 * 191 -> Sum of factors is 303 -> 1041714 / 303 = 3438 exactly.
%p isA046366 := proc(n)
%p local sofpp ;
%p if isA046355(n) then
%p sofpp := A262049(n) ;
%p if modp(n,sofpp) = 0 then
%p true;
%p else
%p false;
%p end if;
%p else
%p false;
%p end if;
%p end proc:
%p for n from 2 to 1000 do
%p if isA046366(n) then
%p printf("%d,",n) ;
%p end if;
%p end do: # _R. J. Mathar_, Sep 09 2015
%t palQ[n_]:=Reverse[x=IntegerDigits[n]]==x; Select[Range[4,10^5],!PrimeQ[#]&&And@@palQ/@Join[{y=Total[Times@@@(x=FactorInteger[#])]},First/@x]&&IntegerQ[#/y]&](* _Jayanta Basu_, Jun 05 2013 *)
%Y Cf. A046367.
%K nonn,base
%O 1,1
%A _Patrick De Geest_, Jun 15 1998
|
