A372693 Numbers k such that A372692(k) = A372692(k+1) > 1.
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
Keywords
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
Crossrefs
Programs
-
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] &]
-
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);}
Comments