[問題] 這支快速排序法的副程式怎麼使用
找到副程式如下
(來源：https://gist.github.com/t-nissie/479f0f16966925fa29ea)
---------------------------------------------------
recursive subroutine quicksort(a, first, last)
implicit none
real*8 a(*), x, t
integer first, last
integer i, j
x = a( (first+last) / 2 )
i = first
j = last
do
do while (a(i) < x)
i=i+1
end do
do while (x < a(j))
j=j-1
end do
if (i >= j) exit
t = a(i); a(i) = a(j); a(j) = t
i=i+1
j=j-1
end do
if (first < i-1) call quicksort(a, first, i-1)
if (j+1 < last) call quicksort(a, j+1, last)
end subroutine quicksort
---------------------------------------------------
自己試寫了一下主程式如下
---------------------------------------------------
program quick_sort_test
real*8 a(*)
integer first, last
data a/12.0,54.0,65.0,2.0,3.0,40.0,91.0,7.0,321.0,50.0/
first=0
last=22
call quicksort(a, first, last)
end program
---------------------------------------------------
因為不太懂
1. real*8 a(*) 的星號
2. first, last 的意義
結果一直顯示array boundary有問題
調了 first 和 last 也還是一樣
主程式應該要怎麼修改呢?
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 42.76.232.247 (臺灣)
※ 文章網址: https://www.ptt.cc/bbs/Fortran/M.1590375785.A.145.html
※ 編輯: BanPeeBan (42.76.232.247 臺灣), 05/25/2020 11:04:07
推
05/25 11:52,
3年前
, 1F
05/25 11:52, 1F
→
05/25 11:52,
3年前
, 2F
05/25 11:52, 2F
→
05/25 11:53,
3年前
, 3F
05/25 11:53, 3F
推
05/25 11:59,
3年前
, 4F
05/25 11:59, 4F
→
05/25 11:59,
3年前
, 5F
05/25 11:59, 5F
推
05/26 20:46,
3年前
, 6F
05/26 20:46, 6F
→
05/26 20:48,
3年前
, 7F
05/26 20:48, 7F
→
05/26 20:51,
3年前
, 8F
05/26 20:51, 8F
→
05/26 20:53,
3年前
, 9F
05/26 20:53, 9F
→
05/26 20:55,
3年前
, 10F
05/26 20:55, 10F
推
05/26 21:03,
3年前
, 11F
05/26 21:03, 11F
→
05/26 21:04,
3年前
, 12F
05/26 21:04, 12F
Fortran 近期熱門文章
PTT數位生活區 即時熱門文章
