HDU2086简单找通项的题。。。

A1 = ?

Time Limit: 5000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 4389    Accepted Submission(s): 2816

Problem Description

有如下方程：A _i = (A _i-1 + A _i+1)/2 - C _i (i = 1, 2, 3, .... n).
若给出A ₀, A _n+1, 和 C ₁, C ₂, .....C _n.
请编程计算A ₁ = ?

 

Input

输入包括多个测试实例。
对于每个实例，首先是一个正整数n,(n <= 3000); 然后是2个数a ₀, a _n+1.接下来的n行每行有一个数c _i(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

题意很清楚了，不过做这道题还是花了蛮久时间的- -（智商低没有找到正确的方法），开始想通过写N=2 3 4时的A1 表达式来找规律 发现不靠谱，后面换成了写相邻两项的
表达式，然后累加发现也错了，最后才走上正轨发现要充分利用好A0和AN+1这两个数，然后只要求A1那么递推的式子里面尽可能多的出现A1就好。然后试着去写了一下大概是
这样的：
A1+A2=A0+A2-2*(C1+C2);
A1+A3=A0+A4-2*(C1+C2+C3);// 实际上是：2*(A1+A3)=A0+A2+A2+A4-2*(C1+C3) ,又2A2=A1+A3-2*C2代换可以得到左边的式子后面以此类推。。。
A1+A4=A0+A5-2*(C1+C2+C3+C4);
...
A1+AN=A0+AN+1-2*(C1+C2+C3+...CN)；
发现累加还是有点不太好表达，然后去给它凑了一项
A1+A1=A0+A2-2*C1；
再次累加可以得到：
(N+1)A1+A2+A3+A4+...AN=N*A0+A2+A3+A4+...AN+AN+1-2*(N*C1+(N-1)*C2+(N-2)*C3+....CN)；
消元然后可以得到：A1=（N*A0+AN+1-2*(N*C1+(N-1)*C2+...CN)）/(N+1)；

代码很简单。。不过我竟然用了两张草稿纸- -。。。 这智商低的- -。。。 不能忍啊TVT。。。



   
   
#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
using namespace std;

int main()
{
    int n;
    double  c[3500];
    while(scanf("%d",&n)!=EOF)
    {
        double  a,b;
        scanf("%lf %lf ",&a,&b);
        int i;
        for(i=1;i<=n;i++)
        {
            scanf("%lf",&c[i]);
        }
        double  sum=0;
        for(i=1;i<=n;i++)
        {
            sum+=(n-i+1)*c[i];
        }
        sum*=2;
        double ans;
        ans=(n*a+b-sum)/(n+1);
        printf("%.2lf\n",ans);



    }



    return  0;
}