2015 ACM/ICPC Asia Regional Shenyang Online & HDU 5451 Best Solver【循环节+数论||共轭矩阵幂+广义斐波那契降幂】

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;
}