OpenJudge Pell数列

来源:互联网 发布:软件维护工程师 编辑:程序博客网 时间:2024/06/05 01:19

Pell数列





总时间限制: 
3000ms 
内存限制: 
65536kB
描述
Pell数列a1, a2, a3, ...的定义是这样的,a1 = 1, a2 = 2, ... , an = 2 * an − 1 + an - 2 (n > 2)。
给出一个正整数k,要求Pell数列的第k项模上32767是多少。
输入
第1行是测试数据的组数n,后面跟着n行输入。每组测试数据占1行,包括一个正整数k (1 ≤ k < 1000000)。
输出
n行,每行输出对应一个输入。输出应是一个非负整数。
样例输入
218
样例输出
1408
源代码:
#include<iostream>#include<cstdio>#include<cstdlib>using namespace std;long pell[1000000];long Pell(long n){    if (n == 1)    {        return 1;    }    if (n == 2)    {        return 2;    }    return 2 * pell[n - 1] + pell[n - 2];}int main(){    int n, m, i, j;    long s[10000];    cin >> n;    for (i = 0; i < n; i++)    {        cin >> s[i];    }    for (i = 1; i < 1000000; i++)    {        pell[i] = Pell(i);        pell[i] %= 32767;    }    for (i = 0; i < n; i++)    {        cout << pell[s[i]] << endl;    }    return 0;}


原创粉丝点击