费马小定理: 假如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)相比效率有了极大的提高。
原理
编辑
把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);
}