快速幂（别人那看来的）

链接：https://ac.nowcoder.com/acm/contest/322/A
来源：牛客网
 

时间限制：C/C++ 1秒，其他语言2秒
空间限制：C/C++ 32768K，其他语言65536K
64bit IO Format: %lld

题目描述

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

聪明的你可以帮助他们吗？

输入描述:

第一行有一个整数n，n <= 1e5

接下来一行有n个数，每个数的大小不超过1e16

输出描述:

输出取模之后的和

示例1

输入

复制

4
1 6 9 12

输出

复制

21502

 

#include<iostream>
#include<algorithm>
#include<vector>
#include<cstring>
#include<queue>
#include<stack>
using namespace std;
typedef long long ll;
const ll maxn = 10000019;
ll pow1(ll a,ll b)//快速幂 
{
    ll r=1;
    if(b==0)
        return 1;
    if(b==1)
        return a;
    while(b>1)
    {
        if(b%2!=0)
        {
        	r*=a;
        	r%=maxn; 
		}
        a*=a;
        a%=maxn; 
        b/=2;
    }
    return a*r%maxn;
}
int main(){
    int n;
    ll a,sum=0,b,j;
	cin>>n;
	for(ll i=1;i<=n;i++){
		cin>>a;
		sum+=pow1(a,i);
		sum%=maxn;
	} 
	cout<<sum;
    return 0;
}

 

a^b

求a的b次方对p取模的值(1<=a,b,p<=1e9)

int power(int a,int b,int p){
    int ans = 1 % p;
    for(; b; b >>= 1){
        if(b & 1) ans = (long long)ans * a % p;
        a = (long long)a * a %p;
    }
    return ans;
}

a*b(64位整数乘法)

求a乘b对p取模的值，其中1<=a,b,p<=1e18

long long mul(long long a,long long b,long long p){
    long long ans = 0;
    for(; b; b >>= 1){
        if(b & 1) ans = (ans + a) % p;
        a = a * 2 %p;
    }
    return ans;
}