PAT B1049. 数列的片段和 (20)

1049 数列的片段和 (20 分)
给定一个正数数列，我们可以从中截取任意的连续的几个数，称为片段。例如，给定数列 { 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>
int main(){
	int n;
	scanf("%d",&n);
	double a[100010];
	for(int i=0;i<n;i++){
		scanf("%lf",&a[i]);
	}
	double sum=0;
	for(int i=0;i<n;i++){//计算第i个数字打头的和sum； 
		double sum1=0;
		for(int j=i;j<n;j++){
			sum1+=a[j];
			sum+=sum1; 
		}
	}
	printf("%.2f",sum);
	return 0;
}

思路2：

找规律：
第i个数出现的总次数等于i*（n+1-i），于是总和就等于Σ（i*次数）。

参考代码2：

#include <cstdio>
int main(){
	int N;
	scanf("%d",&N);
	double a;
    long long ans=0;
	for(int i=1;i<=N;i++){
		scanf("%lf",&a);
		ans+=(long long) (a*1000)*i*(N+1-i);//第i位的总出现次数为a*i*(N+1-i); 
	}
	printf("%.2f\n",ans/1000.0);
	return 0;
}

点评：

①此题属于数学问题型，关键是要透过表象看见本质。
②思路1暴力法会超时，思路2原代码会有一个样例点错误。现找到其他同学的评论供参考。

以下是柳神博客原文 “N比较大时，double类型的值多次累加导致的精度误差，因为输入为十进制小数，存储到double中时，计算机内部使用二进制表示，且计算机的字长有限，有的十进制浮点数使用二进制无法精确表示只能无限接近，在字长的限制下不可避免会产生舍入误差，这些细微的误差在N较大时多次累加会产生较大误差，所以建议不要使用double类型进行多次累加的精确计算，而是转为能够精确存储的整型。尝试把输入的double类型的值扩大1000倍后转为long long整型累加，同时使用long long类型保存sum的值，输出时除以1000.0转为浮点型再输出（相当于把小数点向后移动3位后再计算，避免double类型的小数部分存储不精确，多次累加后对结果产生影响） 但我觉得乘以1000也未必严谨，可能测试样例最小只有小数点后三位，如果测试样例变成小数点后四位、五位、六位，乘以1000相当于直接在小数点后三位处截断，而原本第四五六位经过多次累加进位后依然可能会引起精度问题，但如果乘以10000就会超出long long的值，我认为最精确的应该是截取到所有小数中最大的位数的那一位。。可能我的想法有疏漏，经过测试，测试样例确实没有超过小数点后三位，虽然修改为乘以1000后代码已经AC，但如果对测试样例稍加修改，可能又会导致不AC了…所以这道题先打个问号吧，我猜可能将来题目样例还会被修改…”