login
A372693
Numbers k such that A372692(k) = A372692(k+1) > 1.
2
7380, 18755, 24804, 25631, 26299, 27467, 32799, 44891, 49196, 49725, 50940, 53603, 59652, 64386, 71027, 79739, 85788, 89300, 94275, 103212, 105056, 105875, 124992, 129348, 132011, 138060, 141899, 147100, 149435, 155484, 158147, 164196, 170324, 175571, 181620, 184283
OFFSET
1,1
COMMENTS
The numbers k such that A372692(k) = A372692(k+1) = 1 are in A372690.
LINKS
MATHEMATICA
f[p_, e_] := p^(2^(-1 + Position[Reverse@ IntegerDigits[e, 2], _?(# == 0 &)])); s[1] = 1; s[n_] := s[n] = Times @@ (Flatten@ (f @@@ FactorInteger[n]) + 1);
Select[Range[10^5], (s1 = s[#]) > 1 && s1 == s[# + 1] &]
PROG
(PARI) s(n) = {my(f = factor(n), k); prod(i = 1, #f~, k = apply(x -> 1 - x, binary(f[i, 2])); prod(j = 1, #k, if(k[j], f[i, 1]^(2^(#k-j)) + 1, 1))); }
lista(kmax) = {my(s1 = s(1), s2); for(k = 2, kmax, s2 = s(k); if(s1 > 1 && s1 == s2, print1(k - 1, ", ")); s1 = s2); }
CROSSREFS
Subsequence of A068781.
A372694 is a subsequence.
Similar sequences: Cf. A002961, A064125, A293183, A306985, A343819, A348346.
Sequence in context: A284925 A128478 A257498 * A043589 A043806 A043824
KEYWORD
nonn
AUTHOR
Amiram Eldar, May 10 2024
STATUS
approved