题意:有一天集训队的学弟们正在计算一堆数,但是 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;
}