Iterated Linear Function

链接：http://acm.hust.edu.cn/vjudge/problem/401745/origin

题目：Consider a linear function f(x) = Ax + B. Let's define g⁽⁰⁾(x) = x and g⁽ⁿ⁾(x) = f(g^(n - 1)(x)) for n > 0. For the given integer values A, B, n and x find the value of g⁽ⁿ⁾(x) modulo 10⁹ + 7.

题意：求把abnx带入上述公式的结果对1e9+7取模的结果。

分析：非常多的乘方运算，可以推出通项公式用快速幂求，也可以直接用矩阵描述这个公式，在用矩阵快速幂直接得出结果。（这是我第一次手写矩阵快速幂，非常的丑陋。。）

题解：

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <queue>
#include <stack>
#include <vector>
#include <map>
#include <cstring>
#include <functional>
#include <cmath>
using namespace std;

__int64 mn;
__int64 a,b,n,x;
__int64 m[2][2],mt[2][2],mans[2][2],manst[2][2],ans[2];

void s()
{
	m[0][0]=((mt[0][0]*mt[0][0])%mn+(mt[0][1]*mt[1][0])%mn)%mn;
	m[0][1]=((mt[0][0]*mt[0][1])%mn+(mt[0][1]*mt[1][1])%mn)%mn;
	m[1][0]=((mt[1][0]*mt[0][0])%mn+(mt[1][1]*mt[1][0])%mn)%mn;
	m[1][1]=((mt[1][0]*mt[0][1])%mn+(mt[1][1]*mt[1][1])%mn)%mn;
	mt[0][0]=m[0][0];
	mt[0][1]=m[0][1];
	mt[1][0]=m[1][0];
	mt[1][1]=m[1][1];
}

void t()
{
	ans[0]=((mans[0][0]*ans[0])%mn+mans[0][1])%mn;
}

int main()
{
	//freopen("in.txt","r",stdin);
	mn=1000000007;
	while(~scanf("%I64d %I64d %I64d %I64d",&a,&b,&n,&x))
	{
		mt[0][0]=m[0][0]=a;
		mt[0][1]=m[0][1]=b;
		mt[1][0]=m[1][0]=0;
		mt[1][1]=m[1][1]=1;
		ans[0]=x;
		ans[1]=1;
		manst[0][0]=mans[0][0]=1;
		manst[0][1]=mans[0][1]=0;
		manst[1][0]=mans[1][0]=0;
		manst[1][1]=mans[1][1]=1;
		while(n>0)
		{
			if(n&1) 
			{
				mans[0][0]=((manst[0][0]*mt[0][0])%mn+(manst[0][1]*mt[1][0])%mn)%mn;
				mans[0][1]=((manst[0][0]*mt[0][1])%mn+(manst[0][1]*mt[1][1])%mn)%mn;
				mans[1][0]=((manst[1][0]*mt[0][0])%mn+(manst[1][1]*mt[1][0])%mn)%mn;
				mans[1][1]=((manst[1][0]*mt[0][1])%mn+(manst[1][1]*mt[1][1])%mn)%mn;
				manst[0][0]=mans[0][0];
				manst[0][1]=mans[0][1];
				manst[1][0]=mans[1][0];
				manst[1][1]=mans[1][1];
			}
			s();
			n/=2;
		}
		t();
		printf("%I64d\n",ans[0]);
	}
	return 0;
}