http://acm.hdu.edu.cn/showproblem.php?pid=5451
通过一道名为So easy!的共轭构造矩阵幂取模的问题来联想到这个题,不难发现应该用同样的套路来进行共轭构造。
但是问题来了,指数过分强大,取模数超小,于是一开始考虑欧拉降幂再矩阵快速幂。数据一大,发现答案全错!后请教得知无理数无法欧拉降幂。
那么就只有一种可能了——有规律可循!
找循环节,其实欧拉降幂不就是在找循环节嘛(只是这个题用不了),瞬间变简单。
#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<iomanip>
using namespace std;
const int maxn = 1e5 + 12;
#define ll long long int
#define Build build
#define Node node
#define charmax(a,b) a=max(a,b)
#define charmin(a,b) a=min(a,b)
#define clr(a,b) memset(a,b,sizeof a)
bool vis[maxn];
int prim[maxn];
int a[maxn];
//I'm very sorry for my friends and our team "goodbye".That's not my time to say goodbye.I'm back now! ——Irish_Moonshine
int qkm(int x, int mod) {
ll ans = 1, res = 2;
while (x) {
if (x % 2) {
ans = ans * res%mod;
}
res = res * res%mod;
x >>= 1;
}
return ans;
}
int work(int m) {
a[0] = 10 % m;
a[1] = 98 % m;
for (int j = 2;; j++) {
a[j] = (a[j - 1] * 10 - a[j - 2] + m) % m;
if (a[j - 1] == a[0] && a[j] == a[1]) return j - 1;
}
}
int main()
{
int T;
scanf("%d", &T);
int za = 0;
while (T--) {
za++;
int m; ll x;
scanf("%lld %d", &x, &m);
int cs = work(m);
int zy = qkm(x, cs);
printf("Case #%d: ", za);
printf("%d\n", (a[zy] - 1 + m) % m);
}
return 0;
}