login
Indices n, not squarefree, where A055231(n) = A055231(n-A055231(n)).
0

%I #16 Jan 09 2013 15:18:13

%S 9,45,63,99,117,153,171,207,261,279,289,315,333,369,387,423,477,495,

%T 531,549,585,603,639,657,676,693,711,747,765,801,819,855,873,909,927,

%U 963,981,1017,1035,1071,1143,1179,1197,1233,1251,1287,1305,1341,1359,1395

%N Indices n, not squarefree, where A055231(n) = A055231(n-A055231(n)).

%C A055231(n) is the powerfree part of n. This sequence is infinite because all numbers of the form n = 9p, where p is a prime > 3, are in the sequence : A055231(9p) = p and A055231(9p - p) = A055231(8p) = p. The positive numbers of A055792 are also in the sequence because A055792(n) are squares and A055792(n)-1 are also squares.

%e 63 is in the sequence because A055231(63) = A055231(7*3^2) = 7, A055231(63 - 7) = A055231(56) = A055231(7*2^3) = 7.

%p isA013929 := proc(n)

%p n>3 and not numtheory[issqrfree](n) ;

%p end proc:

%p isA207359 := proc(n)

%p isA013929(n) and (A055231(n)- A055231(n- A055231(n))=0);

%p end proc:

%p for n from 1 to 5000 do

%p if isA207359(n) then

%p printf(`%d, `,n);

%p end if;

%p end do: # adapted from A140394.

%Y Cf. A055231, A140394.

%K nonn

%O 1,1

%A _Michel Lagneau_, Feb 17 2012