【思维+快速幂】2018年东北农业大学春季校赛-K wyh的数列

wyh的数列

时间限制：C/C++ 1秒，其他语言2秒
空间限制：C/C++ 262144K，其他语言524288K
64bit IO Format: %lld

题目描述

wyh学长特别喜欢斐波那契数列，F（0）=0，F（1）=1，F（n)=F(n-1)+F(n-2)(n>=2)

一天他突发奇想，想求F（a^b）%c

输入描述:

输入第一行一个整数T(1<=T<=100)，代表测试组数
接下来T行，每行三个数 a,b,c (a,b<=2^64) (1

输出描述:

输出第a^b项斐波那契数对c取余的结果

示例1

输入

3
1 1 2
2 3 1000
32122142412412142 124124124412124 123

输出

1
21
3

思路

快速幂这一点大家应该都可以想到，但是对于这极大的数据可能大家会无从下手。
我们不妨再来找找规律
假设c为3

我们不难发现，在 f(x) mod c当中
当x = 8时，其的值与x= 0的相等
当x = 9时，其的值与x= 1的相等
那么按照规律，下面应该就是前8个mod值反复循环的过程。
那么我们无疑就可以找到规律：
只要我们知道它是从第几项开始循环，那么我们只要找到开始循环之前的所有f(x) mod c值第一遍出现的值，看看我们所要算的数据是对应其循环的第几个就好。
而且我们不难发现最多其是从n²开始循环或者更小，那么我们只要找到f[i] = 1&&f[i-1]==0,那么其不难就是每i-1次循环一次。
而阶乘我们快速幂就可以实现。

坑点

1、想不到
2、没有紫书
3、ull和llu写法和用法忘记了。

AC代码

#include<bits/stdc++.h>
using namespace std;
int mod;
typedef unsigned long long ull;

int f[1000006];

int mod_pow(ull x,ull n)
{
   ull res = 1;
   while(n>0){
    if(n&1) res = res * x %mod;
    x = x * x % mod;
    n>>=1;
   }
   return res;
}

void solve(void)
{
    int t;
    cin>>t;
    while(t--){
        ull a,b,c;
        scanf("%llu%llu%llu",&a,&b,&c);
        f[0] = 0;
        f[1] = 1;
        for(int i = 2 ; i < 1000006 ; i++){
            f[i] = (f[i-1] + f[i-2]) % c;
            if(f[i]== 1 &&f[i-1]==0){
                mod = i - 1;
                break;
            }
        }
        printf("%d\n",f[mod_pow(a%mod,b)]);
    }
}

int main(void)
{
    solve();
    return 0;
}