费马小定理 + 费马大定理 + 勾股数的求解 + 快速幂 + 矩阵快速幂 【模板】

 

 

 

费马小定理： 假如p是质数，且gcd(a,p)=1，那么 a^(p-1)≡1（mod p） 两边都mod p；
即：假如a是整数，p是质数，且a,p互质(即两者只有一个公约数1)，那么a的(p-1)次方除以p的余数恒等于1。

延伸：1. n*a^(p-1) ≡ n (mod p)
　　    2. a^(p-2) ≡ a^(-1)（mod p）

 正序逆序要灵活运用！

 

例题：

C~K的难题：费马小定理+快速幂

Problem Description

众所周知 C~K 喜欢数学，但是他最近被一个题给难住了，题目是这样的。
要求 (A/B)%10007，但由于 A 很大，我们只给出 n (n = A%10007)（我们给定的A必能被B整除，且 gcd(B,10007) = 1）。
你能帮助他解答吗？他会很感谢你的。

Input

数据的第一行是一个 T，表示有 T 组数据。
每组数据有两个数 n (0 <= n < 10007) 和 B (1 <= B <= 10^9)。

Output

对应每组数据输出 (A/B)%10007。

Sample Input

2
1000 53
87 123456789

Sample Output

8893
7424

 

/*
 
(A/B)%mod => A%mod/B => n/B => n*B^(-1) => 延伸2 => n*(B^(mod-2))%mod
 
*/
#include<iostream>
#include<algorithm>
using namespace std;
#define ll long long
 
int qmi(ll a, ll b, ll mod) {   //快速幂
    ll flag = 1;
    while (b) {
        if (b & 1)flag = (flag*a) % mod;
        a = (a*a) % mod;
        b = b >> 1;
    }
    return flag%mod;
}
 
int main() {
    int t, n, b,mod=10007;
    cin >> t;
    while (t--) {
        cin >> n >> b;
 
        cout << n*qmi(b, mod-2, mod) % mod << "\n";
    }
    return 0;
}

 

Find Integer
Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/65536 K (Java/Others)
Total Submission(s): 6597    Accepted Submission(s): 1999
Special Judge
Problem Description

people in USSS love math very much, and there is a famous math problem .

give you two integers n,a,you are required to find 2 integers b,c such that a^n+b^n=c^n.

Input

one line contains one integer T;(1≤T≤1000000)

next T lines contains two integers n,a;(0≤n≤1000,000,000,3≤a≤40000)

 Output

print two integers b,c if b,c exits;(1≤b,c≤1000,000,000);

else print two integers -1 -1 instead.

 Sample Input

1 2 3

Sample Output

4 5

费马大定理（基础知识）：
 费马大定理，又被称为“费马最后的定理”，由17世纪法国数学家皮耶·德·费玛提出。

它断言当整数n >2时，关于x, y, z的方程 x^n + y^n = z^n 没有正整数解。

历经三百多年的历史，最终该定理在1995年被英国数学家安德鲁·怀尔斯彻底证明。

另附费马大定理的简单证明过程（来自百度）有兴趣的可以看一下哦

勾股数的两种计算方法：
1）当a为大于1的奇数2n+1时，b=2n^2+2n, c=2n^2+2n+1；当a为大于4的偶数2n时，b=n^2-1, c=n^2+1。

2）

题目解释：
给出a和n，请找到符合a^n+b^n=c^n的b和c，如果存在请输出一组，否则输出-1 -1

解题思路：

对n分n>2, 2，1，0，四种情况讨论

n=2时应用勾股数的计算方法

ac代码：
 

#include <iostream>
#define ll long long int
using namespace std;
int main()
{
    ll t;
    scanf("%lld",&t);
    while(t--)
    {
        ll n,a;
        scanf("%lld%lld",&n,&a);
        if(n==2)
        {
            ll m=a/2;
            if(a%2==1 && a>1)
            {
 
                 printf("%lld %lld\n",2*m*m+2*m,2*m*m+2*m+1);
            }
            else if(a%2==0 && a>4)
            {
                printf("%lld %lld\n",m*m-1,m*m+1);
            }
            if(a<3)
                printf("-1 -1\n");
            continue;
 
        }
        if(n>2||n==0)
        {
            printf("-1 -1\n");
            continue;
        }
        if(n==1)
        {
            printf("%lld %lld\n",a+1,a+a+1);
            continue;
        }
 
    }
    return 0;
}
//当a为大于1的奇数2n+1时，b=2n^2+2n, c=2n^2+2n+1。
//当a为大于4的偶数2n时，b=n^2-1, c=n^2+1

 

 

 

快速幂

顾名思义，快速幂就是快速算底数的n次幂。其时间复杂度为 O(log₂N)， 与朴素的O(N)相比效率有了极大的提高。

 

原理

编辑

以下以求a的b次方来介绍 [1] 

把b转换成二进制数。

该二进制数第i位的权为

 

例如

11的二进制是1011

11 = 2³×1 + 2²×0 + 2¹×1 + 2º×1

因此，我们将a¹¹转化为算

 

 

快速幂可以用位运算来实现

b and 1{也就是取b的二进制最低位(即第0位)　判断b是否为奇数，是则为1} 

b shr 1{就是去掉b的二进制最低位(即第0位)} 

C++实现为

b & 1//取b二进制的最低位，判断和1是否相同，相同返回1，否则返回0，可用于判断奇偶 

b>>1//把b的二进制右移一位，即去掉其二进制位的最低位

 

O(logn):

#include<bits/stdc++.h>
using namespace std;
int pow_mod(int a, int n, int m)
{
    if(n == 0) return 1;
    int x = pow_mod(a, n/2, m);
    long long ans = (long long)x * x % m;
    if(n % 2 == 1) ans = ans *a % m;
    return (int)ans;
}
int main()
{
    int a, n, m;
    cin >> a >> n >> m;
    cout << pow_mod(a, n, m);
}

#include<bits/stdc++.h>
using namespace std;
int pow_mod(int a, int n, int m)
{
    long long ans = 1;
    while(n){
        if(n&1){
            ans = (ans * a) % m;
        }
        a = (a * a) % m;
        n >>= 1;
    }
    return ans;
}
int main()
{
    int a, n, m;
    cin >> a >> n >> m;
    cout << pow_mod(a, n, m);
}

矩阵快速幂（原理+模板）

https://blog.csdn.net/wust_zzwh/article/details/52058209