多次乘方运算 超级数

 

一、超级数

【题目描述】：

众所周知，多次加法运算是乘法，多次乘法运算是乘方，那么多次乘方运算呢？

 

 

但是，这个数光是想想就远远超过了整个宇宙的大小，所以你只要求出它的值MOD P就可以了

 

【输入】：

第一行两个数N，P，如题目描述。

第二行N个数，表示a_1~n。

 

【输出】：

仅一个数，表示它的值mod p的结果。

 

【样例输入】：

5 13

2 2 2 2 2

 

【样例输出】：

3

 

【说明】：

20%的数据N≤3，ai<1000；

40%的数据ai，p为质数或1；

100%的数据N≤20，ai，p≤Maxlongint。

 

算法简介：

如果X < PHI(C) 则 A^X MOD C = A^(X MOD PHI(C)) MOD C 

 否则 A^X MOD C = A^(X MOD PHI(C)+ PHI(C)) MOD C 

 

 

PHI（X）函数表示的是X前有多少个数与他互质......

 

 

利用欧拉函数来做，不懂的可以去找去找我初中老师，证明可见：http://hi.baidu.com/aekdycoin/blog/item/b6f1762565bb403fc8955908.html

 

 

 

关于欧拉函数：http://wenku.baidu.com/view/2cdefb62caaedd3383c4d35c.html

 

program super; var i,j,n,p:longint; ans:int64; a:array[1..21]of int64; function phi(x:longint):longint; var i:longint; begin phi:=x; for i:=2 to x do if x mod i=0 then begin while x mod i=0 do x:=x div i; phi:=round(phi*(1-1/i)); end; if x>1 then phi:=1; end; function mi(a,b,p:int64):int64;inline; var sum:int64; begin mi:=1; while b>0 do begin if b and 1=1 then mi:=mi*a; b:=b >> 1; a:=a*a; if a>=p then a:=a mod p; if mi>=p then mi:=mi mod p; end; end; function cimi(x,p:int64):int64;inline; var i,j,q:int64; begin q:=phi(p); if x=n-1 then begin if a[x+1]>=q then cimi:=mi(a[x],a[x+1] mod q+q,p)mod p else cimi:=mi(a[x],a[x+1],p) mod p; end else begin i:=cimi(x+1,phi(p)); if i=0 then cimi:=mi(a[x],i+q,p)mod p else cimi:=mi(a[x],i,p) mod p; end; end; begin assign(input,'super.in'); reset(input); assign(output,'super.out'); rewrite(output); read(n,p); for i:=1 to n do read(a[i]); if n=0 then begin close(input);close(output);exit;end; if p=1 then begin writeln(0);close(input);close(output);exit;end; if n=1 then begin writeln(a[1] mod p);close(input);close(output);exit;end; ans:=cimi(1,p)mod p; writeln(ans); close(input); close(output); end.

 

这个是新的，话说大家可以把我的程序翻出来，把PHI中的TRUNC改成ROUND，就对了