青藤训练营 **Pell数列**求解

原题：

	**题目描述**
	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;	
	}
}