login
A195375
Numbers whose cube expanded in base 14 uses only prime digits.
0
13, 21, 27, 693, 3587, 75285, 435365, 437155, 929509, 1731765, 1902599, 1926903, 104916853, 605497891, 1914409189, 4634485093, 12882557131, 112438392323
OFFSET
1,1
COMMENTS
14 is believed to be the first base in which there are infinitely many such numbers.
EXAMPLE
27 is in the sequence because 27^3 = 13 + 5*14 + 2*14^2 + 7*14^3 uses only prime coefficients, namely 13, 5, 2 and 7.
MAPLE
isA195375 := proc(n)
dset := convert(convert(n^3, base, 14), set) ;
for d in dset do
if not isprime(d) then
return false;
end if;
end do;
return true;
end proc:
for n from 1 do
if isA195375(n) then
print(n);
end if;
end do: # R. J. Mathar, Sep 20 2011
PROG
(PARI) isA195375(n, b=14)={n=divrem(n^3, b); until(!n=divrem(n[1], b), isprime(n[2])|return); 1}
/* The following version allows us to vary the set of allowed digits: */
is_A195375(n, a=10412/*=2^2+2^3+2^5+2^7+2^11+2^13*/, b=14)={
n=divrem(n^3, b); until(!n=divrem(n[1], b), bittest(a, n[2])|return); 1
} \\ M. F. Hasler, Sep 21 2011
CROSSREFS
Sequence in context: A304006 A304948 A087683 * A128819 A165955 A331116
KEYWORD
base,nonn
AUTHOR
Kausthub Gudipati, Sep 17 2011
EXTENSIONS
a(13)-a(18) from D. S. McNeil, Sep 20 2011
STATUS
approved