作者BanPeeBan (踢屁屁)
看板Fortran
标题[问题] 这支快速排序法的副程式怎麽使用
时间Mon May 25 11:03:02 2020
找到副程式如下
(来源：
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://webptt.com/cn.aspx?n=bbs/Fortran/M.1590375785.A.145.html
※ 编辑: BanPeeBan (42.76.232.247 台湾), 05/25/2020 11:04:07
1F：推 fragmentwing: 如果我没弄错 你这样参考点会直接选到a(11)吧 可是 05/25 11:52
2F：→ fragmentwing: 你的数列是不是只有10个值 05/25 11:52
3F：→ fragmentwing: 先试着其他别动，把last改成20看看? 05/25 11:53
4F：推 fragmentwing: 痾不 改成10好了 他应该是想抓无序下排在中间的值当 05/25 11:59
5F：→ fragmentwing: 参考值吧 05/25 11:59
6F：推 espresso1: a(*)是在副程式里假设阵列大小的用法 05/26 20:46
7F：→ espresso1: 不能用在主程式，主程式的阵列要有明确的上下限 05/26 20:48
8F：→ espresso1: data有10个a就宣告10或以上 05/26 20:51
9F：→ espresso1: first, last就a是要传递给副程式的上下界 05/26 20:53
10F：→ espresso1: 这里就是1跟10，表示传递a(1)到a(10) 05/26 20:55
11F：推 espresso1: 因此a(*)改为a(10) 05/26 21:03
12F：→ espresso1: first=1, last=10 05/26 21:04