模取余

很久没写了。有次看到csdn里模取余这个题，看了好多次，查阅了一些资料，今天就发一下对此题的见解。

题目是这样的：

大家对指数运算都非常熟悉，定义f(n,k) = n ^ k （n的k次方），
例如f(2, 3) = 8, f(3,4) = 81。 我们定义f(0,0) = 1。 我们的目标是给定整数n1,k1,n2,k2,n，求f(f(n1,k1),f(n2,k2)) % n。 （％是取余数运算）。
 输入格式 多组数据，每组数据就一行包含5个整数n1,k1,n2,k2,n。 (0<=n1,k1,n2,k2<=1000000000, 1 <= n <=10000000)
输出格式 每组数据一行，表示最终结果。

 

其实这里就是求a^b mod c   觉得这公式够明确的了，就直接写代码了。 

#include "math.h"

void main(){
	int n1,k1,n2,k2,n;
	n1=1;
	k1=2;
	n2=3;
	k2=4;
	n=5;
	int p=pow(n1,k1);
	int q=pow(n2,k2);
	int result=(int)pow(p,q)%n;
	printf("%d",result);
	} 

很明显 只有很小的数据才能算出结果。这个时候，我就觉得a^b mod c  是可以继续化简的。

首先是化简指数b，a^b mod c=a*a^(b-1) mod c （其中a^(b-1)可以继续化简，这样可将指数化简为乘积）    这样可以看做A*B mod c   继续化简A*B    有个公式  A*B mod c = ((A mod c )*B) mod c
最终改进代码如下：

#include "math.h"
void modExp(int a,int b,int c){
 int d,i;
 
 d=1,i=1;
 for(;i<=b;i++) {
 d=d*a%c;
 } 
 printf("%d",d);
 }
	
void main(){

	modExp(45,67,89);
} 

这样a^b mod c 就搞定了，关于该题的代码稍加修改就ok啦。这里就省略了。