這個程式使用分治法計算最大子段，結果為最大子段之和，用遞迴程式實現。

原始資料使用隨機函式生成。

採用結構化程式設計，可以很容易改為從標準輸入或檔案讀入資料，只需要修改函式getData即可。

/*
 * 最大子段演算法程式
 */

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

#define N 13

void getData(int [], int);

int maxsumfun(int a[], int left, int right);

int main(void)
{
    int a[N], i;

    getData(a, N); /* 獲得資料放入陣列a中 */

    printf("datas: ");
    for (i = 0; i < N; i++)
        printf("%d ", a[i]);
    printf("\n");

    int ms;
    ms = maxsumfun(a, 0, N - 1);

    printf("maxsum=%d\n", ms);

    return 0;
}

int maxsumfun(int a[], int left, int right)
{
    int maxsum = 0;

    if(left == right)
    {
        if(a[left] > 0)
            maxsum = a[left];
        else
            maxsum = 0;
    }
    else
    {
        int mid, leftsum, rightsum, midsum;

        mid = (left + right) / 2;
        leftsum = maxsumfun(a, left, mid);
        rightsum = maxsumfun(a, mid+1, right);

        int leftsum1 = 0;
        int lefts = 0;
        int i;
        for(i=mid; i>=left; i--)
        {
            lefts += a[i];
            if(lefts > leftsum1)
                leftsum1 = lefts;
        }

        int rightsum1 = 0;
        int rights = 0;
        for(i=mid+1; i<=right; i++)
        {
            rights += a[i];
            if(rights > rightsum1)
                rightsum1 = rights;
        }

        midsum = leftsum1 + rightsum1;
        if(midsum < leftsum)
            maxsum = leftsum;
        else
            maxsum = midsum;
        if(maxsum < rightsum)
            maxsum = rightsum;
    }
    return maxsum;
}

void getData(int d[], int n)
{
    time_t t;
    srand((unsigned) time(&t));  /* 設定隨機數起始值 */

    int i;
    for(i=0; i < n; i++)
        d[i] = rand() % 30 - 10; /* 獲得-10--20之間的整數值 */
}
 

// 遞迴計算最大子段
int maxsumfun(int a[], int left, int right)
{
    int maxsum = 0;

    if(left == right)
    {
        if(a[left] > 0)
            maxsum = a[left];
        else
            maxsum = 0;
    }
    else
    {
        int mid, leftsum, rightsum, midsum;

        mid = (left + right) / 2;
        leftsum = maxsumfun(a, left, mid);
        rightsum = maxsumfun(a, mid+1, right);

        int leftsum1 = 0;
        int lefts = 0;
        int i;
        for(i=mid; i>=left; i--)
        {
            lefts += a[i];
            if(lefts > leftsum1)
                leftsum1 = lefts;
        }

        int rightsum1 = 0;
        int rights = 0;
        for(i=mid+1; i<=right; i++)
        {
            rights += a[i];
            if(rights > rightsum1)
                rightsum1 = rights;
        }

        midsum = leftsum1 + rightsum1;
        if(midsum < leftsum)
            maxsum = leftsum;
        else
            maxsum = midsum;
        if(maxsum < rightsum)
            maxsum = rightsum;
    }
    return maxsum;
}