#include <stdio.h>
#include <stdbool.h>

#define MAXSIZE 9

typedef struct {
    int r[MAXSIZE+1]; // first index used as tmp, not real data
    int len;
}SqList;

void swap(SqList *L, int i, int j) {
    int tmp = L->r[i];
    L->r[i] = L->r[j];
    L->r[j] = tmp;
}

void merge(int sr[], int tr[], int s, int m, int t) {
    // 本函數(shù)的任務(wù)就是比對sr中兩個分組(s..m, m+1..t)的元素大小并歸并到tr
    int j,k,l;

    j = m + 1; // 第2分組的起始位置
    k = s; // k用于tr數(shù)組中的游標與sr中的起始位置對應(yīng)起來
    while (s<=m && j<=t) {
        if (sr[s] < sr[j]) {
            tr[k++] = sr[s++];
        }
        else {
            tr[k++] = sr[j++];
        }
    }

    // 只要是合并，就肯定至少是2個序列合并，肯定會在比對后剩下1個未消耗完元素的序列分組
    while (s<=m) {
        tr[k++] = sr[s++];
    }

    while (j<=t) {
        tr[k++] = sr[j++];
    }
}

void msort(int sr[], int tr[], int s, int t) {
    /*
     * 把sr進行歸并排序并有序保存到(歸并到)tr中
     */
    int m;
    int tmpr[MAXSIZE+1]; // 每層遞歸的臨時數(shù)組，存放本次被調(diào)用時s到t歸并后的下標值(位置與首次傳入的L->r相同)

    if (s == t) {
        tr[s] = sr[s]; // 歸并的思想，1個元素的分組為有序
    }
    else { // 不是1個元素的分組，繼續(xù)分組
        m = (s+t)/2;
        msort(sr, tmpr, s, m);
        msort(sr, tmpr, m+1, t);
        // 合并tmpr到tr完成本層的排序任務(wù)
        merge(tmpr, tr, s, m, t);
    }
}

void merge_sort(SqList *L) {
    msort(L->r, L->r, 1, L->len); // 因為在msort中第1個參數(shù)sr數(shù)組只是讀取，所以這里這樣傳遞沒有問題
}


int main(void) {
    SqList list = {
        {999,50,10,90,30,70,40,80,60,20},
        MAXSIZE
    };

    merge_sort(&list);
    printf("after merge_sort:\n");
    for (int i=0; i<=MAXSIZE; i++) {
        printf("index: %d, value: %d\n",i,list.r[i]);
    }
    return 0;
}

#include <stdio.h>
#include <stdbool.h>
#include <stdlib.h>

#define MAXSIZE 9

typedef struct {
    int r[MAXSIZE+1]; // first index used as tmp, not real data
    int len;
}SqList;

void merge(int sr[], int tr[], int s, int m, int t) {
    // 本函數(shù)的任務(wù)就是比對sr中兩個分組(s..m, m+1..t)的元素大小并歸并到tr
    int j,k,l;

    j = m + 1; // 第2分組的起始位置
    k = s; // k用于tr數(shù)組中的游標與sr中的起始位置對應(yīng)起來
    while (s<=m && j<=t) {
        if (sr[s] < sr[j]) {
            tr[k++] = sr[s++];
        }
        else {
            tr[k++] = sr[j++];
        }
    }

    // 只要是合并，就肯定至少是2個序列合并，肯定會在比對后剩下1個未消耗完元素的序列分組
    while (s<=m) {
        tr[k++] = sr[s++];
    }

    while (j<=t) {
        tr[k++] = sr[j++];
    }
}


void merge_pass(int sr[], int tr[], int k, int len) {
    int i=1; // 合并時的游標
    while (i < len-2*k+1) { // 也就是每次循環(huán)后，當前所剩余的是否還夠2個完整子序列
        merge(sr, tr, i, i+k-1, i+2*k-1); //合并本輪掃描到的2個子序列
        i+=2*k; // 賦值后的i為下一輪2個子序列的起始位置
    }

    // 下面是掃尾工作，**可能會**出現(xiàn)2種情況，a. 剩余1～2個子序列之間的情況, b. 剩余<=1個子序列的情況
    if (i < len-k+1) {
        merge(sr, tr, i, i+k-1, len);
    }
    else { // 這里加else也可以 如果上面i正好把序列消耗完，則循環(huán)不會執(zhí)行
        while (i<len) {
            tr[i] = sr[i];
            i++;
        }
    }
}

void merge_sort(SqList *L) {
    int *tr = (int *)malloc(L->len*sizeof(int));
    int k=1;
    /*
     * 循環(huán)k為序列長度，與遞歸的方式相比，正好反過來，非遞歸方式直接從序列為1開始合并，直到序列不小于待排序的數(shù)組長度為止
     * 每次循環(huán)都是子序列*4變長的過程
     */ 
    while (k<L->len) {
        merge_pass(L->r, tr, k, L->len); // 序列1變2
        k++;
        merge_pass(tr, L->r, k, L->len); // 序列2變4
        k++;
    }
}


int main(void) {
    SqList list = {
        {999,50,10,90,30,70,40,80,60,20},
        MAXSIZE
    };

    merge_sort(&list);
    printf("after merge_sort2:\n");
    for (int i=0; i<=MAXSIZE; i++) {
        printf("index: %d, value: %d\n",i,list.r[i]);
    }
    return 0;
}