Evaluate binomial coefficients: Difference between revisions

Content added Content deleted
Line 900: Line 900:
write(*,fmt='(i0,a)',advance='no') n_C_r(j,i),' '
write(*,fmt='(i0,a)',advance='no') n_C_r(j,i),' '
end do
end do
write(*,*)
write(*,'(a,i0)') ' 60C30 = ',n_C_r(60,30)
end do
end do
stop
stop
Line 907: Line 907:


pure function n_C_r(n, r) result(bin)
pure function n_C_r(n, r) result(bin)
use, intrinsic :: iso_fortran_env, only : int64
integer(16) :: bin
integer(int64) :: bin
integer, intent(in) :: n
integer, intent(in) :: n
integer, intent(in) :: r
integer, intent(in) :: r
integer(int64) :: num
integer(int64) :: den
integer(16) :: num
integer :: i, k
integer(16) :: den
integer, parameter :: primes(*) = [2,3,5,7,11]
integer :: i
integer :: k
integer, parameter :: primes(*) = [2,3,5,7,11,13,17,19]
num = 1
num = 1
den = 1
den = 1
Line 962: Line 963:
24Cr = 1 24 276 2024 10626 42504 134596 346104 735471 1307504 1961256 2496144 2704156 2496144 1961256 1307504 735471 346104 134596 42504 10626 2024 276 24 1
24Cr = 1 24 276 2024 10626 42504 134596 346104 735471 1307504 1961256 2496144 2704156 2496144 1961256 1307504 735471 346104 134596 42504 10626 2024 276 24 1
25Cr = 1 25 300 2300 12650 53130 177100 480700 1081575 2042975 3268760 4457400 5200300 5200300 4457400 3268760 2042975 1081575 480700 177100 53130 12650 2300 300 25 1
25Cr = 1 25 300 2300 12650 53130 177100 480700 1081575 2042975 3268760 4457400 5200300 5200300 4457400 3268760 2042975 1081575 480700 177100 53130 12650 2300 300 25 1
60C30 = 118264581564861424


=={{header|FreeBASIC}}==
=={{header|FreeBASIC}}==