login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A216894 n - (sum of prime factors of n) is a positive square. 1
1, 6, 22, 54, 68, 164, 166, 336, 388, 454, 588, 854, 886, 1086, 1122, 1124, 1636, 1710, 1828, 2182, 2356, 2468, 2702, 2960, 3046, 3048, 3708, 3748, 3770, 4036, 4054, 4655, 5106, 5394, 5636, 6502, 7108, 7368, 7956, 8324, 9170, 9188, 9412, 9438, 9471, 9726 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
Contains 4p for odd prime p if 3p-2 is a square, in particular if p is in A122430. - Robert Israel, Apr 13 2014
LINKS
EXAMPLE
54 = 2*3^3 and 54 -(2+3) = 49 is a square, hence 54 is in the sequence.
MAPLE
filter:= proc(x) local y;
y:= x - add(i, i=numtheory:-factorset(x));
y > 0 and issqr(y)
end proc;
N:= 10000; # to get all entries <= N
A216894:= select(filter, [$1..N]); # Robert Israel, Apr 13 2014
MATHEMATICA
nspfQ[n_]:=Module[{c=n-Total[Transpose[FactorInteger[n]][[1]]]}, c>0 && IntegerQ[ Sqrt[c]]]; Join[{1}, Select[Range[10000], nspfQ]] (* Harvey P. Dale, Feb 12 2015 *)
CROSSREFS
Cf. A008472.
Sequence in context: A202803 A366917 A092185 * A212692 A034134 A081441
KEYWORD
nonn
AUTHOR
Michel Lagneau, Sep 19 2012
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 15 11:27 EDT 2024. Contains 374332 sequences. (Running on oeis4.)