HDU 2086 A=？ 数学题

题目描述：有一个公式，A_i = (A_i-1 + A_i+1)/2 - C_i (i = 1, 2, 3, .... n).，如果给出A₀, A_n+1, 和 C₁, C₂, .....Cn要你计算出A1是多少。

解题报告：一个简单的数学题，我的做法先将公式变形得到：A[n] = 2*A[n-1]+2*C[n-1]-A[n-2],是分别定义两个数组，设A1 = x，然后两个数组里面一个存放第n项的系数，另一个存放第n项的常数，这样一直递推到第n+1项，然后直接解一元一次方程就可以了。代码附上：

#include<cstdio>
#include<cmath>
const int MAX = 3000+5;
int n;
double xishu[MAX],changshu[MAX],C[MAX],a0,an1;
void dabiao() {
    xishu[1] = 1;
    changshu[1] = 0;
    xishu[2] = 2;
    changshu[2] = 2.0*C[1]-a0;
    for(int i = 3;i<=n+1;++i) {
        xishu[i] = 2.0*xishu[i-1]-xishu[i-2];
        changshu[i] = 2.0*changshu[i-1]+2*C[i-1]-changshu[i-2];
    }
}
int main() {
     while(scanf("%d",&n)!=EOF) {
         scanf("%lf%lf",&a0,&an1);
         for(int i = 1;i<=n;++i)
         scanf("%lf",&C[i]);
         dabiao();
         double ans = (an1-changshu[n+1])/xishu[n+1];
         printf("%.2lf
",ans);
     }
     return 0;
}