HDOJ 2086 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位小数).

-----------------------------------------------------------------------------------------------------

纯粹数学题，找规律:

An = (1/2)An-1 + (1/2)An+1 - Cn

An-1 = (2/3)An-2 + (1/3)An+1 - (2/3)Cn - (4/3)Cn-1

An-2 = (3/4)An-3 + (1/4)An+1 - (1/2)Cn - Cn-1 - (3/2)Cn-2

An-3 = (4/5)An-4 + (1/5)An+1 - (2/5)Cn - (4/5)Cn-1 - (6/5)Cn-2 - (8/5)Cn-3

......

（是不是有点感觉了呢）

接着:

A1 = (n/(n+1))A0 + (1/(n+1))An+1 - (2/(n+1))Cn - (4/(n+1))Cn-1 - ... -(2n/(n+1))C1

     = [ nA0 + An+1 - 2(Cn + 2Cn-1 + 3Cn-2 + ... + nC1) ]/(n+1)

那么，下面就可以直接套公式了,我的AC代码如下：

C++语言 : Codee#7546

//HDOJ 2086
#include <iostream>
using namespace std ;
int main ()
{
    const int MaxSize = 3001 ;
    int n;
    double a1 , a0 , end ;
    double c [ MaxSize ];
    while ( cin >> n)
    {
        cin >> a0 >> end ;
        for ( int i = 0 ; i < n; ++ i )
            cin >> c [ i ];
        a1 = n * a0 + end ;
        for ( int i = 0 ; i < n; ++ i )
            a1 -= 2 * ( i + 1 ) * c [ n - 1 - i ];
        a1 = a1 / ( 1 + n);
        printf ( "%.2f /n " , a1 );
    }
    return 0 ;
}