OFFSET
1,1
COMMENTS
Is the sequence infinite?
Odd terms in the sequence are: a(18) = 1575, a(25) = 2457, a(26) = 2925, a(34) = 3915, a(35) = 4125, a(47) = 6345, a(50) = 6669, ...
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..1000
EXAMPLE
a(1) = 40 = 2^3*5^1, with a = 11 and b = 13.
a(2) = 48 = 2^4*3^1, with a = 11 and b = 19.
Notice that a and b are both prime and not equal.
MAPLE
a:= proc(n) option remember; local an, bn, k, l;
for k from 1 +`if`(n=1, 0, a(n-1)) do
l:= ifactors(k)[2];
an:= add( i[1] * i[2], i=l);
bn:= add( i[1] ^ i[2], i=l);
if isprime(an) and isprime(bn) and an<>bn then break fi
od; k
end:
seq(a(n), n=1..50); # Alois P. Heinz, Jan 20 2011
MATHEMATICA
a[n_] := a[n] = Module[{an, bn, k, p, e}, For[k = 1 + If[n==1, 0, a[n-1]], True, k++, {p, e} = Transpose[FactorInteger[k]]; an = p.e; bn = Total[p^e]; If[PrimeQ[an] && PrimeQ[bn] && an != bn, Break[]]]; k];
Array[a, 50] (* Jean-François Alcover, Nov 20 2020 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Bobby Browning and Rohan Hemasinha (rhemasin(AT)uwf.edu), Jan 07 2011
EXTENSIONS
More terms from Alois P. Heinz, Jan 20 2011
STATUS
approved