HDU2086 A1 = ?（组合数学 + 思维）

A1 = ?

Problem Description
有如下方程：Ai = (Ai-1 + Ai+1)/2 - Ci (i = 1, 2, 3, … n).
若给出A0, An+1, 和 C1, C2, …Cn.
请编程计算A1 = ?

Input
输入包括多个测试实例。
对于每个实例，首先是一个正整数n,(n <= 3000); 然后是2个数a0, an+1.接下来的n行每行有一个数ci(i = 1, …n);输入以文件结束符结束。

Output
对于每个测试实例，用一行输出所求得的a1(保留2位小数).

Sample Input
1
50.00
25.00
10.00
2
50.00
25.00
10.00
20.00

Sample Output
27.50
15.00

题目分析

题目让计算某个序列中的一个数字
那么我想应该分为两种情况思考

1. 需要求出完整的序列，然后再求出答案
2. 直接根据已知条件，构造一个表达式，求出答案

现在题目让我们求a₁，就第一种情况而言，不能直接求出a₁，那应该能够求出a_n倒推出a₁，然而通过阅读题目发现a_n和a₁的求法一样，是一个对称的序列，所以先算出a_n再算a₁不可行，然后我想能不能通过求a_n/2对分递归求解，因为我想到 a_n = (a_n-1 + a_n+1) / 2 中 a_n/2 = (a₀ + a_n) / 2，希望直接求出a_n/2，然而在这道题似乎不大行或许是我当时没找到

直接算的话，先写几项
a₁ = a₀/2 + a₂/2 - c₁
a₂ = a₁/2 + a₃/2 - c₂
a₃ = a₂/2 + a₄/2 - c₃
a₄ = a₃/2 + a₅/2 - c₄

然后把上面几项都加起来 消去公共项得，
a₁ + a₄ = a₀ + a₅ - 2 * (c₁ + c₂ + c₃ + c₄)
可以从中抽象出a₁ + a_n = a₀ + a_n+1 - 2 * (c₁ + … + c_n)
在n = 2 的时候 有
a₂ = a₀ + a₃ - a₁ - 2 * (c₁ + c₂)
同时，
a₁ = a₀/2 + a₂/2 - c₁
所以，
2a₁ = a₀ + a₀ + a₃ - a₁ - 2 * (c₁ + c₂) - 2 * c₁

3a₁ = 2a₀ + a₃ - 4 * c₁ - 2 * c₂

a₁ = (2a₀ + a₃ - 4 * c₁ - 2 * c₂) / 3

到这里发现a₁已经只和已知条件有关
同理，我们尝试写更多项，就能够推导到到他的公式
在n = n的时候，

a₁ = n * a₀ + a_n - 2 * (n * c₁ + (n - 1) * c₂ + … + c_n) / （n + 1）

推导到这个之后就算OK了

完整代码如下，

#include <bits/stdc++.h>
using namespace std;

int main() {
	int n;
	while (~scanf("%d",&n)) {
		double a0, an1;
		scanf("%lf%lf", &a0, &an1);
		double sum = 0;
		for (int i = 0; i < n; ++i) {
			double t;
			scanf("%lf", &t);
			sum += (n - i) * t;
		}
		double a1 = (n * a0 + an1 - 2 * sum) / (n + 1);
		printf("%.2lf\n", a1);
	}
	return 0;
}