【nowcoder 214267】广义肥波

最新推荐文章于 2023-04-12 07:00:00 发布

SSL_TJH

最新推荐文章于 2023-04-12 07:00:00 发布

阅读量148

点赞数

分类专栏： # 数学或数论牛客题目

要转发的话记得附上链接哦~（不过也不会有人转发的啦）

本文链接：https://blog.csdn.net/weixin_43346722/article/details/112726718

版权

广义斐波那契递推取模优化模运算枚举算法

关键词由CSDN通过智能技术生成

数学或数论同时被 2 个专栏收录

285 篇文章 3 订阅

订阅专栏

牛客题目

89 篇文章 0 订阅

订阅专栏

广义肥波

题目链接：nowcoder 214267

到牛客看：

——>点我跳转<——

题目大意

定义广义斐波那契数列：
$\left\{\begin{matrix} f_1=1\\ f_2=1\\ f_i=a\times f_{i-1}+b\times f_{i-2}(n\geq3,n\in \mathbb{Z}) \end{matrix}\right.$
然后让你求 $m^{f_n} \%(1e9+7)$

思路

看到题目的范围，我们可以看出应该是枚举一个 $n$ ，然后其他数的 log，大概这样的复杂度。

然后很明显枚举 $n$ 可以得出 $f_n$ ，但是你要求的是 $m^{f_n}$ ， $f_n$ 不能取模，你总不能搞高精吧。

那我们考虑不求 $f_n$ ，枚举 $n$ 求出 $m^{f_i}$ 的值 $g_i$ 。
那我们可以看出转移就是： $g_i=g_{i-1}^a\times g_{i-2}^b$ 。（就是全部运算往上诺一级）

这样就可以在取模的意义下进行而不出问题了。

代码

#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define ll long long
#define mo 1000000007

using namespace std;

ll a, b, m, re, tmp, l, r;
int n;

ll ksm(ll x, ll y) {
	re = 1ll;
	while (y) {
		if (y & 1) re = (re * x) % mo;
		x = (x * x) % mo;
		y >>= 1;
	}
	return re;
}

int main() {
	scanf("%lld %lld %lld %d", &a, &b, &m, &n);
	
	l = r = m;
	for (int i = 3; i <= n; i++) {
		tmp = (ksm(l, b) % mo * ksm(r, a)) % mo;
		l = r;
		r = tmp;
	}
	
	printf("%lld", r);
	
	return 0;
}