HDU 1121 Complete the Sequence 数学题

题目大意：

就是给出一个长度为n的数列， 是一个多项式通项的数列，要求填补后面的m项

大致思路：

布吉岛为什么这题被HDU分在动态规划里面=_=

我用数学方法解的

首先如果给出的数列有n个数，那么这个数列至多只能是一个 (n - 1) 次的多项式 (拉格朗日插值公式可以看出,不过不用在意这个细节,,知道一定是 (n - 1)次以下即可)

那么对于给出的 n 个数， 每次都相邻两个数作差，后面的一个数减去前面的一个，这样当得到下一个数列时都会比原来少一项

如果得到的新数列不是每项都相等就继续作差直到相等为止， 这样就是利用了两个齐次多项式 f ( n ) - f ( n - 1 )一定会降次， 同时也就是n阶等差数列的性质

这样由于输入的数最多100个，我需要最差的次数在10^4以内得到最后相同的数列之后，在后面加上C个( 需要填的个数 ) 得到的相同数，然后往回加过去就得到了原来需要补充的C个后续的数

代码如下：

Result  :  Accepted     Memory  :  312 KB     Time  :  109 ms

/*
 * Author: Gatevin
 * Created Time:  2014/8/15 19:58:24
 * File Name: haha.cpp
 */
#include<iostream>
#include<sstream>
#include<fstream>
#include<vector>
#include<list>
#include<deque>
#include<queue>
#include<stack>
#include<map>
#include<set>
#include<bitset>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cctype>
#include<cmath>
#include<ctime>
#include<iomanip>
using namespace std;
const double eps(1e-8);
typedef long long lint;

int S,C;
int a[110];
int now;
vector <int> v[110];
int last;
bool check(int x)
{
    for(unsigned int i = 0; i <= v[x].size() - 1; i++)
    {
        if(v[x][i] != v[x][0]) return true;
    }
    return false;
}

int main()
{
    int t;
    scanf("%d", &t);
    while(t--)
    {
        scanf("%d %d", &S, &C);
        for(int i = 1; i <= S; i++)
        {
            scanf("%d", a + i);
            v[0].push_back(a[i]);
        }
        now = 0;
        while(check(now))//切克闹~
        {
            for(unsigned int i = 1; i <= v[now].size() - 1; i++)
            {
                v[now + 1].push_back(v[now][i] - v[now][i - 1]);
            }
            now++;
        }
        last = v[now].size() - 1;
        for(int i = 1; i <= C; i++)
        {
            v[now].push_back(v[now][0]);
        }
        int tmp = 0;
        for(int i = now - 1; i >= 0; i--)
        {
            tmp = v[i].size() - 1;
            for(int j = 1; j <= C; j++)
            {
                v[i].push_back(v[i][tmp + j - 1] + v[i + 1][last + j]);
            }
            last = tmp;
        }
        for(int i = 0; i < C; i++)
        {
            printf("%d%c",v[0][S + i], i == C - 1? '\n' : ' ');
        }
        for(int i = 0; i <= now; i++)
        {
            if(!v[i].empty()) v[i].clear();
        }
    }
    return 0;
}