最大子序列和/积 寻找主要元素

《数据结构与算法分析——C语言描述》第二章看书笔记

第二章

练习题主要就是各种最大子序列和/积和寻找主要元素的问题。

正文：

其实在之前在cf上遇到过一个类似于最大子序列和的问题，当时心里就想着肯定是划定的那个区域中间慢慢移动，当时的题没有那么复杂，所以我当时的算法并没有那么精妙也通过了。书里让自己思考为什么这种线性算法是正确可行的：

int MaxSubsequenceSum(const int A[], int N) {
    int ThisSum, MaxSum, j;

    ThisSum = MaxSum = 0;
    for (j = 0; j < N; j++) {
        ThisSum += A[j];

        if (ThisSum > MaxSum)
            MaxSum = ThisSum;
        else if (ThisSum < 0)
            ThisSum = 0;
    }
    return MaxSum;
}

我是觉得它保证了每一次对ThisSum的累加都是值得的，ThisSum为负表示现有的这样的子序列不是最好的（而存在MaxSum里的也不一定是最好的 # 不一定而已 也有可能就是最好的了）。那么既然前面的子序列不够好，就继续向后寻找。MaxSum已经存好了现在已经遍历过的序列的最大子序列和了，所以不用担心会因为直接放弃ThisSum前面的一大堆而影响结果（这也是我第一次看这个算法时最大的质疑）。这样，在线性时间内，就完成了求最大子序列和。（对了，它默认了如果整个序列全是正数那么最大子序列和为0）

另外还讲了折半查找法、欧几里得算法（辗转相除法）、高效幂运算（分奇偶进行递归处理）。这些都比较基础，而且好理解。

练习

有些练习很基础（甚至我觉得第一章的练习其实挺难的，就算做出来我也做的并不完善），这里只分享我觉得有一些障碍的（作为菜鸡如我）。

2.6 算法运算时间分析
以前从来没有做过时间分析…不过粗略的估计也不是很难，练习过后应该会比较熟练了。
比较有趣的是最后两个例子。

Sum = 0;
for (i = 0; i < N; i++)
    for (j = 0; j < i * i; j++)
        for (k = 0; k < j; k++)
            Sum++;

Sum = 0;
for (i = 0; i < N; i++)
    for (j = 0; j < i * i; j++)
        if (j % i == 0)
            for (k = 0; k < j; k++)
                Sum++;

唯一的区别就是那一条if语句。通过简单的证明，可以发现，在最坏的情况下(i ==N)，j的取值也仅仅是：$N, 2N, 3N, ... , N^2$ 所以最内层的语句其实也只执行了$N^4$的时间，也就是$O(N^4)$.(不知道这么说准不准确，反正大概意思表达出来了…)

2.10 Horner法则…其实就是秦九韶算法嘛…高中数学讲这些的时候我记得我听的还很起劲…

2.12
(a) 和最大子序列和完全类似，反几个符号就可以了，自己写了一个：
（同样我也默认了如果整个序列每个数都为正，那么最小子序列和为0）

#include <stdio.h>

#define MAXN 10000

int MinSubsequenceSum(const int A[], int N);
int main(void) {
    int N;
    printf("Please enter N: (N <= 10000)");
    scanf("%d", &N);

    int A[MAXN];
    printf("Please enter N numbers:");
    int i;
    for (i = 0; i < N; i++)
        scanf("%d", &A[i]);

    int min_sum = MinSubsequenceSum(A, N);
    printf("min_part_sum is : %d\n", min_sum);

    return 0;
}

int MinSubsequenceSum(const int A[], int N) {
    int ThisSum, MinSum, j;
    ThisSum = MinSum = 0;

    for (j = 0; j < N; j++) {
        ThisSum += A[j];
        if (ThisSum < MinSum)
            MinSum = ThisSum;
        else if (ThisSum > 0)
            ThisSum = 0;
    }
    return MinSum;
}

(b) 求最小正子序列和
这道题有点想不明白…正子序列和最小的话，岂不就是整个序列中最小的那个正数？还是我断句有问题？难道题的意思是说： 正(子序列和) 吗？（应该是这样了…有点木）
尝试着写了一个，要是有问题望不吝赐教…学编程没多久…

#include <stdio.h>

#define MAXN 10000

int MinSubsequenceSum(const int A[], int N);
int main(void) {
    int N;
    printf("Please enter N: (N <= 10000)");
    scanf("%d", &N);

    int A[MAXN];
    printf("Please enter N numbers:");
    int i;
    for (i = 0; i < N; i++)
        scanf("%d", &A[i]);

    int min_sum = MinSubsequenceSum(A, N);
    printf("min_part_sum is : %d\n", min_sum);

    return 0;
}

int MinSubsequenceSum(const int A[], int N) {
    int ThisSum, MinSum, j;
    ThisSum = MinSum = 0;

    for (j = 0; j < N; j++) {
        ThisSum += A[j];
        if (ThisSum <= 0 && A[j+1] + ThisSum < 0)
            ThisSum = 0;
        else if ((ThisSum < MinSum && ThisSum > 0) || MinSum == 0)
            MinSum = ThisSum;
    }
    return MinSum;
}

(c) 求最大子序列乘积
似乎我只能想到最倒霉的时候$O(N^2)$的辣鸡算法…
我假设在只输入一个负数的情况下，最大子序列乘积为0（改掉这个默认也非常简单 在这里我懒得改了）。代码经过我修改之后在LintCode上跑过，所以除了上面说的那么默认以外应该还是没太大问题（要改的话把FinalMax 的默认值改成INT_MIN后再稍微加点东西就行了） 代码如下：

#include <stdio.h>

#define MAXN 1000

int MaxSubsequencePro(const int A[], int N);
int main(void) {
    int N;
    printf("Please enter N: (N <= 1000)");
    scanf("%d", &N);

    int A[MAXN] = { 0 };
    printf("Please enter N numbers:");
    int i;
    for (i = 0; i < N; i++)
        scanf("%d", &A[i]);

    int MaxPro = MaxSubsequencePro(A, N);
    printf("MaxPro is : %d\n", MaxPro);

    return 0;
}

int MaxSubsequencePro(const int A[], int N) {
    int PartMax, FinalMax, j, k, f1, f2, part1, part2;
    FinalMax = 0;
    PartMax = 1;
    part1 = part2 = 1;
    f1 = f2 = 0; // f标志每一个区间的第一个负号已经读完 f2标志这个数组里给的数字不全为0
    for (j = 0; j <= N; j++) {

        if (!A[j]) {
            if (PartMax < 0) {
                for (k = j - 1; A[k] != 0, k >= 0; k--) {
                    if (A[k] < 0) {
                        part2 *= A[k];
                        break;
                    }
                    else 
                        part2 *= A[k];
                }
                PartMax /= (part1 > part2) ? part1 : part2;
            }

            if (PartMax > FinalMax)
                FinalMax = PartMax;

            PartMax = part1 = part2 = 1;
            f1 = 0;
            continue;
        }

        f2 = 1;
        if (!f1) {
            part1 *= A[j];
            if (A[j] < 0)
                f1 = 1;
        }
        PartMax *= A[j];

    }

    if (!f2)
        FinalMax = 0;

    return FinalMax;
}

2.14 Eratosthenes筛…这个用C和py都看过也写过很多次了…虽然有一个*号但也不写代码了（python里面结合着做筛选的filter来做 非常的舒服）
2.15 看了答案 $X^2, X^4, X^8, X^{10}, X^{20}, X^{40}, X^{60}, X^{62}$
2.16 又看了答案，就是把N转换成二进制之后，从$X$一直乘到$X^{2^{[log^{N}]}}$ 就那样吧
2.19 寻找主要元素
这里其实我有一个问题，当N是奇数的时候，答案里说应该看前面有没有出现那个候选元素，如果有就直接丢掉最后那个数，如果没有就把它加进去…我不是很懂怎么看有没有出现候选元素…所以我直接把它加进去了，按道理来说不会出问题，只是效率低一丢丢而已。时间复杂度是$O(N)$

#include <stdio.h>

#define MAXN 1000

int * find_candidates(int[], int N);

int main(void) {
    int N, i;
    int A[MAXN];
    printf("Enter N:");
    scanf("%d", &N);
    printf("Enter N numbers: ");
    for (i = 0; i < N; i++)
        scanf("%d", &A[i]);

    int main_element = find_candidates(A, N);
    printf("The main element is : %d\n", main_element);

    return 0;
}

int * find_candidates(int A[], int N) {
    int j, k, end;
    end = (N & 1) ? N - 1 : N;
    for (j = 0, k = 0; j < end - 1; j += 2) {
        if (A[j] == A[j + 1])
            A[k++] = A[j];
    }
    // 看不懂答案里说的意思 奇数的时候我偷懒直接把最后一个数给加进去了......不知道会不会有问题orz
    if (N & 1)
        A[k] = A[N - 1];
    if (k + 1 <= 2)
        return A[0];
    else
        return find_candidates(A, k + 1);
}

2.25 立方最大子序列和
这本书真是跟玄学一样，终于快肝完第二章了…怀着激动的心情思考最后一道题…
诶这几个公式不就是高中数学求前n项和的公式么…（貌似是）
不过…本辣鸡还是不懂这道题的意思…

唉 好好修炼……不然心里的想法都没法实现了啊…