HDU-1163 Eddy's digital Roots

前置芝士: 同余定理，快速幂

这个题目是杭电上的一道数论相关的题目
如果暴力的话是一定会超时的，这个时候我们就要用到同余定理来优化。

同余定理是数论中的重要概念，下面贴上百度百科的定义：同余定理

给定一个正整数m，如果两个整数a和b满足a-b能够被m整除，即(a-b)/m得到一个整数，那么就称整数a与b对模m同余，记作a≡b(mod m)。对模m同余是整数的一个等价关系。

在这里我们简单地证明一下：
对于给定的十进制数$$x=1000\ast a+100\ast b+10\ast c+d$$我们将这个式子进行转换
$$x=9(99\ast a+9\ast b+c)+(a+b+c+d)$$因为9一定能整除 3 和 9 ，所以只要 $(a+b+c+d)$ 能除尽 9 。那么 $x$ 也一定能除尽3和9。
所以我们可以认为$x\equiv (a+b+c+d)(mod9)$
那么根据这个定理我们就可以很容易得出题目的解法

假设 $N^N=1234567$ , 那么 $$1234567mod9=1$$ $1+2+3+4+5+6+7=28$$$28mod9=1$$$2+8=10$$$10mod9=1$$$1+0=1$$$1mod9=1$$
所以 $N^N$ 的数字根 = $N^N\%9$

注意 当 $N^N\%9=9$ 时,对应的数字根为 $9$!

本题因为$N$的范围过大,所以要配合快速幂来使用

代码如下:

#include<bits/stdc++.h> 
using namespace std;
int qPow(long long A , long long x)  //A^x
{
    if (x == 0) return 1;
    long long ans = 1;

    while (x)
    {
        if ( x & 1) 	//如果n的尾数为1; 
        {
            ans = ans * A % 9 ;  //  ans = ans * A % 10
        }
        A = A * A % 9; 	// A = A * A % 10
        x >>= 1;  
    }
    return ans;
}
int main()
{
	int n,ans;
	
	while(~scanf("%d",&n) && n)
	{
		ans=qPow(n,n);
		printf("%d\n",ans==0? 9:ans);
	}
	return 0;
}

参考博客：
https://www.cnblogs.com/zzqc/p/6684794.html
https://blog.csdn.net/Techmonster/article/details/50113789