UVa 11582 (快速幂取模) Colossal Fibonacci Numbers!

题意：

斐波那契数列f(0) = 0, f(1) = 1, f(n+2) = f(n+1) + f(n) (n ≥ 0)

输入a、b、n，求f(a^b)%n

分析：

构造一个新数列F(i) = f(i) % n，则所求为F(a^b)

如果新数列中相邻两项重复出现的话，则根据递推关系这个数列是循环的。

相邻两项所有可能组合最多就n²中，所以根据抽屉原理得到这个数列一定是循环的。

求出数列的周期，然后快速幂取模即可。

 1 #include <cstdio>
 2 #include <iostream>
 3 using namespace std;
 4
 5 const int maxn = 1000 + 10;
 6 typedef unsigned long long ULL;
 7 int f[maxn][maxn * 6], period[maxn];
 8
 9 int pow_mod(ULL a, ULL b, int n)
10 {
11     if(b == 0) return 1;
12     int k = pow_mod(a, b/2, n);
13     k = k * k % n;
14     if(b % 2 == 1) k = k * a % n;
15     return k;
16 }
17
18 int solve(ULL a, ULL b, int n)
19 {
20     if(a == 0 || n == 1) return 0;
21     int p = pow_mod(a % period[n], b, period[n]);
22     return f[n][p];
23 }
24
25 int main(void)
26 {
27     freopen("11582in.txt", "r", stdin);
28     for(int n = 2; n < maxn; ++n)
29     {
30         f[n][0] = 0, f[n][1] = 1;
31         for(int i = 2; ; ++i)
32         {
33             f[n][i] = (f[n][i-2] + f[n][i-1]) % n;
34             if(f[n][i-1] == 0 && f[n][i] == 1)
35             {
36                 period[n] = i - 1;
37                 break;
38             }
39         }
40     }
41
42     ULL a, b;
43     int n, T;
44     scanf("%d", &T);
45     while(T--)
46     {
47         cin >> a >> b >> n;
48         printf("%d\n", solve(a, b, n));
49     }
50
51     return 0;
52 }

代码君