UVA11582 Colossal Fibonacci Numbers!

题目地址

题目链接

题解

get一个新操作，序列的循环节

这据说这一种模数很小的序列的话都会有循环节？对这玩意有点萌币。

如果\(Fib\)序列中在\(mod\space n\)的意义下，出现了连续的0和1，就说明出现了循环节，对这个循环节进行取膜就行了

本来还打了矩阵快速幂的结果没用上

#include <bits/stdc++.h>
using namespace std;
#define ll unsigned long long

ll a, b, mod;
ll p[1010];
vector<ll>f[1010];

ll power(ll a, ll b, ll mod) {
    ll ans = 1;
    while(b) {
        if(b & 1) ans = ans * a % mod;
        a = a * a % mod;
        b >>= 1;
    }
    return ans % mod;
}

struct mat {
    ll m[2][2];
    mat() { memset(m, 0, sizeof(m)); }
    mat operator * (mat &x) const {
        mat c;
        for(int k = 0; k < 2; ++k) {
            for(int i = 0; i < 2; ++i) {
                for(int j = 0; j < 2; ++j) {
                    c.m[i][j] += m[i][k] * x.m[k][j];
                    c.m[i][j] %= mod;
                }
            }
        }
        return c;
    }
};

mat mat_pow(ll b) {
    mat ans, a;
    ans.m[0][0] = ans.m[1][1] = 1;
    a.m[0][0] = a.m[1][0] = a.m[0][1] = 1;
    a.m[1][1] = 0;
    while(b) {
        if(b & 1) ans = ans * a;
        a = a * a;
        b >>= 1;
    }
    return ans;
}

int solve() {
    if(!a || mod == 1) return 0;
    ll t = power(a % p[mod], b, p[mod]);
    return f[mod][t];
}

int main() {
    for(int n = 2; n <= 1000; ++n) {
        f[n].push_back(0); f[n].push_back(1);
        for(int i = 2; ; ++i) {
            f[n].push_back((f[n][i - 1] + f[n][i - 2]) % n);
            if(f[n][i] == 1 && f[n][i - 1] == 0) {
                p[n] = i - 1;
                break;
            }
        }
    }
    int T;
    scanf("%d", &T);
    while(T--) {
        scanf("%llu%llu%llu", &a, &b, &mod);
        printf("%llu\n", solve());
    }
    return 0;
}

转载于:https://www.cnblogs.com/henry-1202/p/10286321.html