Forbidden numbers: Difference between revisions

Line 521:

83331 = forbidden number count up to five hundred thousand

833329 = forbidden number count up to five million

</pre>

=={{header|Phix}}==

function forbidden(integer n)

while n>7 and remainder(n,4)=0 do

n /= 4

end while

return remainder(n,8) = 7

end function

function forbidden_le_count(atom lmt)

atom result = 0, Power4 = 1, n = 7

while n<=lmt do

atom delta := Power4*8

result += floor((lmt-n)/delta)+1;

Power4 *= 4;

n = Power4*7

end while

return result

end function

sequence f50 = {}

integer i = 0

while length(f50)<50 do

if forbidden(i) then f50 &= i end if

i += 1

end while

printf(1,"The first 50 forbidden numbers are:\n%s\n",

{join_by(f50,1,10," ",fmt:="%3d")})

for t=2 to iff(machine_bits()=32?16:17) do

atom lmt = 5*power(10,t), count = forbidden_le_count(lmt)

printf(1,"Forbidden numbers up to %,d: %,d\n",{lmt,count})

end for

<pre>

The first 50 forbidden numbers are:

7 15 23 28 31 39 47 55 60 63

71 79 87 92 95 103 111 112 119 124

127 135 143 151 156 159 167 175 183 188

191 199 207 215 220 223 231 239 240 247

252 255 263 271 279 284 287 295 303 311

Forbidden numbers up to 500: 82

Forbidden numbers up to 5,000: 831

Forbidden numbers up to 50,000: 8,330

Forbidden numbers up to 500,000: 83,331

Forbidden numbers up to 5,000,000: 833,329

Forbidden numbers up to 50,000,000: 8,333,330

Forbidden numbers up to 500,000,000: 83,333,328

Forbidden numbers up to 5,000,000,000: 833,333,330

Forbidden numbers up to 50,000,000,000: 8,333,333,327

Forbidden numbers up to 500,000,000,000: 83,333,333,328

Forbidden numbers up to 5,000,000,000,000: 833,333,333,327

Forbidden numbers up to 50,000,000,000,000: 8,333,333,333,327

Forbidden numbers up to 500,000,000,000,000: 83,333,333,333,326

Forbidden numbers up to 5,000,000,000,000,000: 833,333,333,333,327

Forbidden numbers up to 50,000,000,000,000,000: 8,333,333,333,333,325

Forbidden numbers up to 500,000,000,000,000,000: 83,333,333,333,333,323

</pre>