51Nod_1126 求递推序列的第N项

最新推荐文章于 2019-09-14 10:31:24 发布

原创最新推荐文章于 2019-09-14 10:31:24 发布 · 304 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

算法专栏收录该内容

356 篇文章

订阅专栏

博客围绕51Nod_1126题目展开，该题要求根据给定的递推式f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7，以及A、B和N的值，求出f(n)。输入为A、B、N，输出为f(n)的值。分析指出可通过找循环节求解，还提及会给出C++程序。

51Nod_1126 求递推序列的第N项

http://www.51nod.com/Challenge/Problem.html#!#problemId=1126

题目

有一个序列是这样定义的：f(1) = 1, f(2) = 1, f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7。给出A，B和N，求f(n)的值。

输入

输入3个数：A,B,N。数字之间用空格分割。(-10000 <= A, B <= 10000, 1 <= N <= 10^9)

输出

输出f(n)的值。

样例输入

3 -1 5

样例输出

分析

由递推式 f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7可知，f(n)由后两项决定，又因为f(1)=f(2)=1，因此只要找到下一个相邻两项都为1的位置就可以找到循环节

C++程序

#include<iostream>

using namespace std;

typedef long long ll;

const int MOD=7;
const int N=10000000;

int f[N]={1,1,1};

ll fun(ll a,ll b,ll n)
{
	if(n==1||n==2)
	  return 1;
	int len=0;
	for(ll i=3;i<=n;i++){
		f[i]=(a*f[i-1]+b*f[i-2])%MOD;
		if(f[i]==1&&f[i-1]==1){
			len=i-2;
			f[0]=f[i-2];
			return f[n%len];
		} 
	}
	return f[n];
}


int main()
{
	ll a,b,n;
	cin>>a>>b>>n;
	cout<<fun(a,b,n)<<endl;
	return 0;
}