《算法竞赛入门经典》第十章 数学概念与方法 UVa11582(幂取模，循环规律，模算术）

#include<iostream>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<algorithm>
#define N 1005
using namespace std;
typedef unsigned long long ull;
ull a,b;
int n;
vector<int> table[N];
int rem[N];//用于记录每个n所求得的循环数 
ull pow_mod(ull a,ull b,int n)//快速幂取模，分治二分思想 
{
	if(b==0) return 1;//结束条件，开始"归” 
	ull x=pow_mod(a,b/2,n);//分治开始 ,开始"递” 
	x=(x*x)%n;
	if(b%2) x=(x*a)%n;//如果是奇数的话，就额外乘一个a 
	return x;
} 
void make_table()//制表，以不同的n作为行数，每行记录一个循环内的全部数，并用rem数组来记录循环的长度 
{
	for(int i=2;i<N; i++ )
	{
	    table[i].push_back(0);
		table[i].push_back(1);//将每一行的前两个初始化为1和0 
		for(int j=2;;j++)
		{
			table[i].push_back((table[i][j-1]+table[i][j-2])%i);
			if(table[i][j]==1&&table[i][j-1]==0)//检验是否出现了循环 
	        {
  				rem[i]=j-1;//记录循环的长度 
  				break;
  			}
		}	
	}
}

int main()
{
   //freopen("E:\\in.txt","r",stdin); 
   //freopen("E:\\out.txt","w",stdout);
   make_table();//事先制表 
   int t;
   cin>>t;
   while(t--)
   {
   	  cin>>a>>b>>n;
   	  if(a==0) cout<<0<<endl;
   	  else if(rem[n]==0) cout<<0<<endl;//对两种特殊情况进行处理 
   	  else{
		 int ans=pow_mod(a%rem[n],b,rem[n]);//得到a^b，在一个循环中所处的位置 
   	     cout<<table[n][ans]<<endl;
   	  }
   }
   return 0;
}