费马小定理

定理解释

费马小定理：如果质数p，和满足gcd(a,p) == 1的任意数字a，那么就一定满足下面的式子：

证明等这里就不叙述了，想要了解的话可以查看：百度百科：费马小定理

应用

我们这里就引用一个百度百科中的例子：计算2^100（mod 13）的数值

在这个例子中我们知道p = 13，并且p是质数，那么我们就可以去找最大的m值满足：（2^(p-1)）^m

我们求解到最大的m值为： 100/（p-1）向下取整为8

那么我们的式子就可以转化为：

那么根据费马小定理得前面余数为1，那么最终答案为：

结果为：3

这样我们就可以不必去计算2^100这样一个大数了

延伸

还有与费马小定理的推广：欧拉定理，对于任何两个互质的正整数a,g (gcd(a,g) = 1)，m >= 2 有 a^(φ(m)) = 1 (mod m)  

这里的φ(m)表示的是小于或等于m的正整数中与m互质的数的数目

具体内容可以参见：维基百科：欧拉函数

栗子

ZOJ-3785

这个题中我们的p也就是一周的天数7，最开始的形式：

我们对他们进行求余操作：

然后利用费马小定理：

然后我们发现这个是42一循环的，那么我们就可以打表然后进行求解了。

#include <iostream>
#include <algorithm>
#include <cstring>
#include <cstdio>
#include <cmath>
#define INF 0x3f3f3f3f
#define ll long long

using namespace std;
char s[10][10] = {"Sunday","Monday","Tuesday","Wednesday","Thursday","Friday","Saturday"};
int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}
ll quick_pow(ll a,ll b)
{
    ll ans = 1;
    while(b)
    {
        if(b & 1) ans *= a;
        a *= a;
        b >>= 1;
    }
    return ans;
}
int main()
{
    for(int i = 1;i < ;)
    ll ans[50] = {0,1,5,4,1,4,5,5,6,0, 4, 6 ,0 ,6, 6, 0, 2, 0, 1, 6, 0, 0, 1, 5, 6, 3, 0, 6, 6, 0, 1, 4, 6, 5, 6, 6, 0, 2, 4, 5, 0, 6, 6};
//    int cnt = 1;           打表操作
//    for(int i = 1;i <= 6;i ++)
//    {
//        int num = 0;
//        for(int j = i;;j ++)
//        {
//            ans[cnt+7*num] = quick_pow(i,j>=7?j-6:j);
//            ans[cnt+7*num] %= 7;
//            num ++;
//            if(num == 6) break;
//        }
//        cnt++;
//    }
//    for(int i = 1;i <= 42;i ++)
//        ans[i] += ans[i-1],ans[i] %= 7;
//   这里我们得到的表就是上面ans数组
    ll T,n;
    scanf("%I64d",&T);
    while(T--)
    {
        scanf("%I64d",&n);
        int num = (((n/42%7) * ans[42]) % 7 + ans[n%42] % 7 )% 7;
        printf("%s\n",s[(num+6) % 7]);
    }
    return 0;
}

参考博客

https://blog.csdn.net/zcy_2016/article/details/55054146