hdu5451 Best Solver （广义斐波那契数列求通项+矩阵快速幂）

写在最前面，方便复习：
已知数列 { a n } \{a_n\} {an​} 其满足递推式 $a_{n+2}=c_1{a}_{n+1}+c_2{a}_n$ 其中 $c_1,c_2$ 为常数且 $c_2\ne 0$
记 $x^2=c_{1}x+c_2$ 的两根为 $x_1,x_2$ ，此时有：

1. $x1\ne x_2$ 时， $a_n=b_1{x}_1^n+b_2{x}_2^n$ ，$b_1,b_2$ 满足$\{\begin{array}{c}{a}_1=b_1{x}_1+b_2{x}_2\\ a_2=b_1{x}_1^2+b_2{x}_2^2\end{array}$
2. $x1=x_2=x$ 时， $a_n=b_1{x}_1^n+b_2{x}_2^n$ ，$b_1,b_2$ 满足$\{\begin{array}{c}{a}_1=(b_1+b_2)x\\ a_2=(b_1+b_2)x^2\end{array}$

题目链接：点击这里

题目大意：
$t$ 组输入，每组输入给定正整数 $x(1\le x<2^{32})$ ，和质数 $p(p\le 46337)$ 求：
$$\lfloor y=(5+2\sqrt{6}{)}^{1+2x}\rfloor \mathrm{m}\mathrm{o}\mathrm{d}p$$

题目分析：
设 $x_1=5+2\sqrt{6},x_2=5-2\sqrt{6}$ ，则这两根显然是方程 $x^2-10x+1=0$ 的两根
故其递推方程为 $a_{n+2}=10a_{n+1}-a_n$ ，其通项为 $a_n=(5+2\sqrt{6}{)}^n+(5-2\sqrt{6}{)}^n$
因为 $5-2\sqrt{6}<1$ ，所以 $(5-2\sqrt{6}{)}^n<1$ ，所以题目所求就是 $a_n-1$
因为此题模数不到 $1e5$ 所以循环节不会很长可以直接暴力算出，但是如果模数很大则需要构造矩阵来求循环节，具体方法点击这里
之所以要找循环节是因为 $2^x$ 很大 对于每一个模数都存在循环节，这样就可以避免直接进行幂次运算

具体细节见代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<set>
#include<map>
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
ll read()
{
	ll res = 0,flag = 1;
	char ch = getchar();
	while(ch<'0' || ch>'9')
	{
		if(ch == '-') flag = -1;
		ch = getchar();
	}
	while(ch>='0' && ch<='9')
	{
		res = (res<<3)+(res<<1)+(ch^48);//res*10+ch-'0';
		ch = getchar();
	}
	return res*flag;
}
const int maxn = 1e5+5;
ll mod;
const double pi = acos(-1);
const double eps = 1e-8;
ll n,cnt,a[maxn],loop[maxn];
ll get_loop()
{
	a[0] = 2%mod,a[1] = 10%mod;
	for(int i = 2;i < maxn;i++)
	{
		a[i] = (10*a[i-1]%mod-a[i-2]+mod)%mod;
		if(a[i] == a[1] && a[i-1] == a[0])
			return i-1;
	}
}
ll qpow(ll a,ll b,ll mod)
{	
	ll res = 1;
	while(b)
	{
		if(b&1) res = res*a%mod;
		a = a*a%mod;
		b >>= 1;
	}
	return res;
}
int main()
{
	ll t = read();
	while(t--)
	{
		n = read(),mod = read();
		if(!loop[mod]) loop[mod] = get_loop(); //这里要记忆化一下不然就t掉了
		ll pos = (1+qpow(2,n,loop[mod]))%loop[mod]; 
		for(int i = 2;i <= pos;i++) //因为记忆化过所以上次get_loop 不一定枚举到了pos，所以需要重新扫一遍
			a[i] = (10*a[i-1]%mod-a[i-2]+mod)%mod;
		printf("Case #%lld: %lld\n",++cnt,(a[pos]-1+mod)%mod);
	}
	return 0;
}