关于快速幂取余运算问题

p次方求和

时间限制： 1000 ms  |  内存限制： 65535 KB

难度：3

描述

一个很简单的问题，求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

# include <stdio.h>

int pow(int x, int y, int z);

int main()
{
    int t;
    scanf("%d", &t);
    while (t--)
    {
        int n, p, i;
        int sum = 0;
        scanf("%d %d", &n, &p);

        for (i=1; i<=n; i++)
          sum += pow(i, p, 10003);

        printf("%d\n", sum%10003);
    }

    return 0;
}

int pow(int x, int y, int z)  //x底数  y指数  z取余的数
{
	int ans = 1;
	x = x % z;        //之前一定要对x取余，因为这样不会影响结果，而且减小运算量
	while (y > 0)
	{
	    if (y % 2 == 1)             //如果要是奇数次方的话，先让ans乘上一个x
                  ans = ( ans * x ) % z;    

           x = (x * x) % z;           //让x自乘，y减半  然后迭代
	    y /= 2;

	}

	return ans;              //可能你会有疑问，为什么返回的是ans  这是因为无论次幂是奇数还是偶数
                                 //对它一直除以2的话，最后一定是一个1，然后就会乘到ans上边
}