次方求模

描述

一个很简单的问题，求1^p+2^p+3^p+……+n^p的和。

输入

第一行单独一个数字t表示测试数据组数。接下来会有t行数字，每行包括两个数字n,p，
输入保证0<n<=1000,0<=p<=1000。

输出

输出1^p+2^p+3^p+……+n^p对10003取余的结果，每个结果单独占一行。

样例输入

2

10 1

10 2

样例输出

55

385

</pre></dd><dd style="font-family:Arial"><pre id="sample_output" style="white-space:pre-wrap; word-wrap:break-word"><pre code_snippet_id="599904" snippet_file_name="blog_20150207_2_290467" name="code" class="cpp">//快速幂模板。/* 快速幂模板， 用到位运算
思路： n为偶数时： x^n = x^(n / 2 + n / 2) = x ^ (n / 2) *ｘ＾(n / 2) 
       n为奇数时： x^n = x^(n / 2 + n / 2 + 1) = x ^ (n / 2) *ｘ＾(n / 2) * x 

<pre name="code" class="cpp">#include <iostream>
#include <stdio.h>
using namespace std;
int f(int base, int pow)
{
    if(pow == 0)
        return 1%10003;
    long long t = f(base, pow >> 1);
    long long ans = t * t % 10003;
    if(pow & 1)
        ans = ans * base % 10003;
    return ans;

}
int main()
{
    int t, n, p, sum, i;
    scanf("%d", &t);
    while(t--)
    {
        sum = 0;
        scanf("%d %d", &n, &p);
        /*if(n == 1) //这个要不得。。否则会wa。。至今不解。。。
            sum = 1;*/
         if(p == 0)
            sum = n % 10003;
        else
        {
            for(i = 1 ; i <= n ; i++)
                sum = (sum+ f(i, p))%10003;
        }
        printf("%d\n", sum);
    }
    return 0;
}