小米OJ-dreamstart 的催促-01

题意：有一天集训队的学弟们正在计算一堆数，但是 dreamstart 感觉他们算的太慢了，就让他们坐在一起想出一个快速计算的方法，但是由于他们一时想不出来，想让你帮助他们。他们说现在有一个数列，要算出第 i 个数的 i^i次幂并且把每个数计算出来的值加到一起，最后答案模 10000019。聪明的你可以帮助他们吗？

输入要求：第 1 个数字为整数n，1<=n<=10^5；余下 n 个数，每个数的大小不超过10^15。

举例：有1,6,9,12一堆数，现要计算(1的次幂+6的次幂+9的次幂+12的次幂)%10000019 。

思路：快速幂取模运算、直接用for循环会产生超时；另外涉及到大数操作，若求幂后直接相加运算，即使用long long也会出现溢出。需计算完一个数的幂次后，然后取余，结果作为下一次的被加数。循环操作求得结果。

计算幂次：

long pow(long a, long b, long m) {  //计算幂次
 	    long r = 1;
        a %= m;
        while (b > 0) {
            if ((b & 1) == 1) {  //判断奇偶性：相与结果为1则为奇数。
                r = (r * a) % m;
            }
            a = (a * a) % m;
            b >>= 1;  //右移一位(b/2)
        }
        return r;
}

代码：

/**
  MIOJ dreamstart的催促
*/
#include <bits/stdc++.h>
#define M 10000019L
using namespace std;
long pow(long a, long b, long m) {  //计算幂次
 	  long r = 1;
      a %= m;
      while (b > 0) {
         if ((b & 1) == 1) {  //判断奇偶性
            r = (r * a) % m;
         }
         a = (a * a) % m;
         b >>= 1;  //右移一位(b/2)
      }
      return r;
}
int main()
{
    int n;
    long res = 0,z;
    scanf("%d",&n);
    for(int i = 1; i<=n; ++i) {
        scanf("%ld",&z);
        long x = pow(i,i,M-1);        //​​​​​​​i的i次幂
        long y = pow(z,x,M);         //第i个数(z)的i的i次幂
        res = (res + y)%M;
    }
    cout<<res;
    return 0;
}