|
%I #19 Jan 23 2024 20:34:44
%S 49,98,103,107,161,187,196,197,206,214,237,239,249,253,322,347,374,
%T 389,392,394,412,417,425,428,443,474,478,479,491,498,499,501,503,506,
%U 509,561,569,644,685,691,694,725,729,735,737,748,753,765,778,779,784,788,789,797,809,817,823,824,829,833,834,837,841,849,850
%N Representable positive integers n that are not the inverse of their inverse in binary64 (double precision) IEEE 754 floating-point arithmetic (Version where 1 and n*(1/n) are unequal).
%H Andrey Zabolotskiy, <a href="/A309541/b309541.txt">Table of n, a(n) for n = 1..10000</a>
%H Efstratios Gallopoulos, <a href="http://scgroup.hpclab.ceid.upatras.gr/class/SC/Notes/book08.pdf">Scientific Computation I</a>, five terms provided (in Greek) [broken link]
%H Mark, comment on blog post <a href="http://www.wolter.biz/2008/08/rundungsfehler-im-flash-player/">"Rundungsfehler im Flash Player?"</a>, eight terms (in German)
%H koDoz, comment in thread <a href="https://scratch.mit.edu/discuss/topic/43137/">"Welche Grenzwerte gelten für Zahlen, Strings usw?"</a>, nine terms (in German)
%o (Python 2 or Python 3)
%o for i in range(1,1000):
%o if i*(1./i) != 1: print(i)
%o (Fortran)
%o doubleprecision one, r
%o integer i
%o parameter (one=1.0D0)
%o do 10 i = 1, 500
%o R = one / dble(i)
%o if ( R * dble(i) .ne. one) write (*,1000) i
%o 1000 format (i0)
%o 10 continue
%o end
%o C _Hugo Pfoertner_, Jan 18 2024
%Y See A275419 for the n != 1/(1/n) version.
%K nonn,fini
%O 1,1
%A _Adam M. Scherlis_, Aug 06 2019
|
