解析:
Fn =
[
f
n
f
n
+
1
]
\begin{bmatrix}fn & fn+1 \end{bmatrix}
[fnfn+1]
Fn+1 =
[
f
n
+
1
f
n
+
2
]
\begin{bmatrix}fn+1 & fn+2\end{bmatrix}
[fn+1fn+2]
[ f n f n + 1 ] \begin{bmatrix}fn & fn+1\end{bmatrix} [fnfn+1] * [ 0 c 1 b ] \begin{bmatrix}0 & c\\ 1 & b\\\end{bmatrix} [01cb] = [ f n + 1 f n + 2 ] \begin{bmatrix}fn+1 & fn+2 \end{bmatrix} [fn+1fn+2]
[ 0 1 ] \begin{bmatrix}0 & 1\end{bmatrix} [01] * [ 0 c 1 b ] \begin{bmatrix}0 & c\\ 1 & b\\\end{bmatrix} [01cb] n = [ f n + 1 f n + 2 ] \begin{bmatrix}fn+1 & fn+2 \end{bmatrix} [fn+1fn+2]
class Solution {
public:
long long MOD=1000000007;
struct lxw
{
long long res[2][2];
}node,base;
lxw multi(lxw a,lxw b)
{
lxw tmp;
memset(tmp.res,0,sizeof tmp.res);
for(int i=0;i<2;i++)
for(int j=0;j<2;j++)
for(int k=0;k<2;k++)
tmp.res[i][j]=(tmp.res[i][j]+a.res[i][k]*b.res[k][j])%MOD;
return tmp;
}
long long nthElement(long long n, long long b, long long c) {
memset(node.res,0,sizeof node.res);
memset(base.res,0,sizeof base.res);
node.res[0][0]=0;node.res[0][1]=1;
base.res[0][1]=c;base.res[1][0]=1;base.res[1][1]=b;
while(n){
if(n&1) node=multi(node,base);
base=multi(base,base);
n>>=1;
}
return node.res[0][0]%MOD;
}
};