一.题目链接:
generator 1
二.题目大意:
给出
有递推式:
求第 n 项. ()
三.分析:
这题如果直接用二进制倍增的话会 TLE. (要将 n 转换为 2 进制)
但如果十进制倍增的话,总共需要操作 次. (貌似还不错的样子)
又学会了一个新知识点(从来没听说过。。。)
出题人:这题我小时候做过(orzzzzzzzz)
四.代码实现:
#include <set>
#include <map>
#include <ctime>
#include <queue>
#include <cmath>
#include <stack>
#include <bitset>
#include <vector>
#include <cstdio>
#include <sstream>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#define eps 1e-8
#define lc k * 2
#define rc k * 2 + 1
#define pi acos(-1.0)
#define ll long long int
using namespace std;
const int M = (int)1e6;
//const ll mod = (ll)1e9 + 7;
const ll inf = 0x3f3f3f3f3f;
struct node
{
ll D[2][2];
}a, b1, b2;;
ll mod;
char s[M + 5];
node mul(struct node a, struct node b)
{
struct node c;
memset(c.D, 0, sizeof(c.D));
for(int i = 0; i < 2; ++i)
{
for(int j = 0; j < 2; ++j)
{
for(int k = 0; k < 2; ++k)
{
c.D[i][j] = (c.D[i][j] + a.D[i][k] * b.D[k][j] % mod) % mod;
}
}
}
return c;
}
int main()
{
scanf("%lld %lld %lld %lld", &a.D[0][0], &a.D[0][1], &b1.D[1][1], &b1.D[0][1]);
b1.D[1][0] = 1;
scanf("%s", s + 1);
scanf("%lld", &mod);
int len = strlen(s + 1);
for(int i = len; i; --i)
{
b2 = b1;
for(int j = 1; j <= 9; ++j)
{
if(s[i] - '0' == j)
a = mul(a, b2);
b2 = mul(b2, b1);
}
b1 = b2;
}
printf("%lld\n", a.D[0][0]);
return 0;
}