[openjudge]035:Pell数列[题解及讲解]

最新推荐文章于 2025-03-22 13:45:54 发布

xu_haiao

最新推荐文章于 2025-03-22 13:45:54 发布

阅读量8.5k

点赞数 14

分类专栏： Openjudge题解文章标签： c++

本文链接：https://blog.csdn.net/xu_haiao/article/details/119565275

版权

Openjudge题解专栏收录该内容

1 篇文章

订阅专栏

博客探讨了Pell数列的定义及快速计算方法。通过优化算法，避免了大数值计算导致的时间消耗，利用模运算保持计算结果在合理范围内。代码示例展示了如何在C++中实现这一优化，从而在限制时间内求得Pell数列的指定项并输出。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

035: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行，每行输出对应一个输入。输出应是一个非负整数。

样例输入

2
1
8

样例输出

1
408

代码如下：（思路借鉴于035:Pell数列.md · 徒步天下/程序设计与算法一OpenJudge题解 - Gitee.com）

#include <iostream>
#include <cstdio>
using namespace std;

int pell(int N){
	if (N == 1){
		return 1;
	}
	else if (N == 2){
		return 2;
	}
	else{
		int f1 = 1,f2 = 2,f3;
		for (int i = 0;i < (N-2);i++){
			f3 = (2*f2 + f1) % 32767;
			f1 = f2;
			f2 = f3;
		}
		return f2;
	}
}

int main(){
	int n;
	int data;
	cin >> n;
	for (int i = 0;i < n;i++){
		cin >> data;
		cout << pell(data) << endl;
	}
	return 0;
}

当初自己做题的时候用的递归，思路很简单但是一旦数字超过30速度会明显变慢，因为pell数列的值上升非常快：

1 ----- 1
2 ----- 2
3 ----- 5
4 ----- 12
5 ----- 29
6 ----- 70
7 ----- 169
8 ----- 408
9 ----- 985
10 ----- 2378
11 ----- 5741
12 ----- 13860
13 ----- 33461
14 ----- 80782
15 ----- 195025
16 ----- 470832
17 ----- 1136689
18 ----- 2744210
19 ----- 6625109
20 ----- 15994428
21 ----- 38613965
22 ----- 93222358
23 ----- 225058681
24 ----- 543339720