原题:
**题目描述**
Pell数列a1,a2,a3...的定义是这样的,a1=1,a2=2,...an=2an−1+an−2(n>2)。
给出一个正整数 k,要求Pell数列的第 k 项模上 32767 是多少。
**输入格式**
第1行是测试数据的组数 ,后面跟着 行输入。每组测试数据占 1 行,包括一个正整数(1≤k<1000000)。
**输出格式**
n 行,每行输出对应一个输入。输出应是一个非负整数。
**样例**
样例输入
2
1
8
样例输出
1
408
解题思路:
最简单的思路,按照递归算法算出pell数列的第n项,然后输出第n项模32767的值。但是我们很快会发现,这个方法会导致运行超时,第n项的值太大,导致数据溢出,那么我们可以稍加改进,在算出一项后,只保留其模32767的值,代码如下
#include <bits.stdc++.h>
using namespace std;
int k=32767;
int f(int x)
{
if(x==1) return 1;
if(x==2) return 2;
int a=(2*(f(x-1))%k+(f(x-2))%k)%k;
return a;
}
int main()
{
int n,a;
cin>>n;
for(int i=1;i<=n;i++)
{
cin>>a;
cout<<f(a)<<endl;
}
return 0;
}
但是提交后,其运行状态是这样的:
没有想到具体的原因,但是发现网上在解这道题目时,大部分使用了数组,数组房间数是超过了所要求的的个数的,代码如下
#include <iostream>
#include <cstdio>
#define k 32767
#define Max 1000000
using namespace std;
int a[Max+5];
int f(int x)
{
if(a[x]!=0) return a[x];
if(x==1) return 1;
if(x==2) return 2;
return a[x]=(2*(f(x-1))%k+(f(x-2))%k)%k;
}
int main()
{
int n,a;
cin>>n;
for(int i=1;i<=n;i++)
{
cin>>a;
cout<<f(a)<<endl;
}
return 0;
————————————————
版权声明:本文为CSDN博主「走马观花_MO」的原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/qq_32447387/article/details/80106633
这篇代码一开始我并没有完全理解,毕竟还是刚刚学习递归的小白,询问了一圈后,才认识到这种新的方法——打表法。
根据递归的原理,假如我要求第12项的值。即a12,那么此时我将要求a11以及a10的值,而a11的值又需要去求a10的值与a9的值,以此类推。我原本的方法,在求a12时,求取了a11以及a10的值,而求取a11的值又需要去求a10的值与a9的值,那么此时求取a[10]的值就进行了两次,之后的a9,a8,a7每次都需要求两次,而这,才仅仅只是求取一个数a12,那如果说我要求很多项,此时计算机所要做的工作重复量是非常大的,所以,作者使用了打表法,每求出一个数,就将其存储在数组中,当需要用到该项时,直接调用即可。
以下附上非递归解法:
#include<bits/stdc++.h>
using namespace std;
#define m 1000000
#define mod 32767
int a[m+1];
int main()
{
int n;
cin>>n;
a[1]=1;
a[2]=2;
for(int i=3;i<=1000010;i++)
{
a[i]=(a[i-1]*2+a[i-2])%mod;
}
int k;
for(int i=1;i<=n;i++)
{
cin>>k;
cout<<a[k]<<endl;
}
}