【PAT B1049】 数列的片段和

给定一个正数数列，我们可以从中截取任意的连续的几个数，称为片段。例如，给定数列 { 0.1, 0.2, 0.3, 0.4 }，我们有 (0.1) (0.1, 0.2) (0.1, 0.2, 0.3) (0.1, 0.2, 0.3, 0.4) (0.2) (0.2, 0.3) (0.2, 0.3, 0.4) (0.3) (0.3, 0.4) (0.4) 这 10 个片段。

给定正整数数列，求出全部片段包含的所有的数之和。如本例中 10 个片段总和是 0.1 + 0.3 + 0.6 + 1.0 + 0.2 + 0.5 + 0.9 + 0.3 + 0.7 + 0.4 = 5.0。

输入格式：
输入第一行给出一个不超过 10^5 的正整数 N，表示数列中数的个数，第二行给出 N 个不超过 1.0 的正数，是数列中的数，其间以空格分隔。

输出格式：
在一行中输出该序列所有片段包含的数之和，精确到小数点后 2 位。

输入样例：

4
0.1 0.2 0.3 0.4

输出样例：

5.00

思路一：一开始直接还是想的暴力。虽然知道数据这么大肯定过不了。（但是能骗分）。

#include <cstdio>

const int maxn = 100010;
double a[maxn];

int main() {
    int N;
    scanf("%d", &N);
    for (int i = 0; i < N; i++) {
        scanf("%lf", &a[i]);
    }
    double sum = 0;
    for (int i = 0; i < N; i++) {
        for (int j = i; j < N; j++) {
            for (int k = i; k <= j; k++) {
                sum += a[k];
            }
        }
    }
    printf("%.2f", sum);
}

思路二：后面就在IPAD上面比划，发现了数学规律。如下图，对于数列中第i个数，想办法找到其所在的片段规律，发现片段首肯定在a1到ai之间，一共有i种选择方法，片段尾同理，有n-i+1种选择方法。由组合可知，ai总共会出现在i * (n-i+1)个片段中。由此可得解：

#include <cstdio>

int main() {
    int N;
    scanf("%d", &N);
    double sum = 0, temp = 0;
    for (int i = 1; i <= N; i++) {
        scanf("%lf", &temp);
        sum += temp * i * (N - i + 1);
    }
    printf("%.2f", sum);
}

然而很不幸，测试点3仍无法通过。

思路三：去柳神博客查找解决方法，发现原来是double类型在计算机中表示精确度的问题。也找到了另外一位大神的解答。
暂时的解决方法是：尝试把输入的double类型的值扩大1000倍后转为long long整型累加，同时使用long long类型保存sum的值，输出时除以1000.0转为浮点型再输出（相当于把小数点向后移动3位后再计算，避免double类型的小数部分存储不精确，多次累加后对结果产生影响）

#include <cstdio>

int main() {
    int N;
    scanf("%d", &N);
    long long sum = 0;
    double temp = 0;
    for (int i = 1; i <= N; i++) {
        scanf("%lf", &temp);
        sum += (long long) (1000 * temp) * i * (N - i + 1);
    }
    printf("%.2f", sum / 1000.0);
}

参考链接：
1.柳神博客 link
2.Ruihan Zheng 大神博客 link

在大量或高精度计算中要避免使用 double 类型。