UVa11582 - Colossal Fibonacci Numbersl (打表找规律)

UVa11582 - Colossal Fibonacci Numbers

解题思路

因为所需的结果是取模的，所以可以在计算序列的时候直接计算取模后的序列，即让$f(i) = f(i)\ mod \ n$，因为n的取值范围是[1, 1000]，所以整个序列元素的取值最多也只有1000种情况。
根据斐波那契数列的性质，可以知道，当$(\ f_i\ $ ,$\ f_{i+1}\ )$开始循环的时候，整个序列就开始循环，因此最多$n^2$时序列就开始循环，（就此题来说，在打表时发现，最长的序列也只在3000左右）。
计算时，枚举n列出所有可能的序列，然后算出$f(a^b)$等于循环中的那一项就行。
注意事项

• 当n=1时，整个序列为（0，0，$\dots$，0），当a=0时，结果都为0；
• 计算$f(a^b)$处于序列那一项时，不能直接计算$a^b\ \%\ mod$，可能会溢出。可以计算$(a\%mod)^b\ \%\ mod$。

代码

#include <iostream>
#include <cstdio>
#include <vector>

using namespace std;
typedef unsigned long long ull;
typedef long long ll;
const int maxn = 1000 + 5;
vector<int> ans[maxn];

void init()
{
    for(int i=2; i<=1000; i++)
    {
        ans[i].push_back(0);
        ans[i].push_back(1);
        for(int j=2; ; j++)
        {
            int temp = (ans[i][j-1]+ans[i][j-2]) % i;
            if(temp==1 && ans[i][j-1]==0) break;
            ans[i].push_back(temp);
        }
    }
}

int pow_mod(ull a, ull b, int n)
{
    if(b==0) return 1;

    ll re = pow_mod(a, b>>1, n);
    re = re * re % n;
    if(b&1) re = re * a % n;

    return re;
}

int main()
{
    init();

    int n, t;
    ull a, b;

    cin>>t;
    while(t--)
    {
        cin>>a>>b>>n;
        if(a==0||n==1)
        {
            cout<<0<<"\n";
        }
        else
        {
            int i = pow_mod(a%(ans[n].size()-1), b, (ans[n].size()-1));
            cout<<ans[n][i]<<"\n";
        }
    }

    return 0;
}