hdu1576-A/B-(同余定理+乘法逆元+费马小定理+快速幂)
Posted shoulinniao
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了hdu1576-A/B-(同余定理+乘法逆元+费马小定理+快速幂)相关的知识,希望对你有一定的参考价值。
A/B
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 10383 Accepted Submission(s):
8302
Problem Description
要求(A/B)%9973,但由于A很大,我们只给出n(n=A%9973)(我们给定的A必能被B整除,且gcd(B,9973)
= 1)。
Input
数据的第一行是一个T,表示有T组数据。
每组数据有两个数n(0 <= n < 9973)和B(1 <= B <= 10^9)。
每组数据有两个数n(0 <= n < 9973)和B(1 <= B <= 10^9)。
Output
对应每组数据输出(A/B)%9973。
Sample Input
2
1000 53
87 123456789
Sample Output
7922
6060
逆元概念:p是素数,在%p的情况下,如果 ( a*inv(a) ) ≡ 1 (%p),则inv(a)是a的乘法逆元
解题过程:
(A/B)%P
= (A%P) * (1/B)%P 同余定理,并且(A%P)=n,已知
= n * (1/B)%P
= (A%P) * (1/B)%P 同余定理,并且(A%P)=n,已知
= n * (1/B)%P
先不管n,
只看(1/B)%P
找B的逆元,按概念就是 ( B*inv(B) ) % P = 1%P = 1;
代入( 1/B )%P
=( B*inv(B)/B )%P
=inv(B) %P
费马小定理概念:如果p是素数,a和p互质,即gcd(a,p)=1,则 a^(p-1)≡1( %p )
逆元和费马小定理联立可得 a * a^(p-2) ≡ 1( %p ), inv(a) = a^(p-2)
代入inv(B) %P
=B^(P-2) %P
则原式 = n * B^(P-2) % P
#include<stdio.h> #include<iostream> #include<algorithm> #include<cstring> #include<math.h> #include<string> #define ll long long #define inf 0x3f3f3f3f using namespace std; const ll p=9973; ll power(ll a,ll b)//快速幂 { ll res=1; while(b) { if(b%2==1) res=res*a%p; b=b/2; a=a*a%p; } return res%p; } int main() { ll t,n,b; scanf("%lld",&t); while(t--) { scanf("%lld%lld",&n,&b); ll ans; ans=n*power(b,p-2)%p; printf("%lld ",ans); } return 0; }
以上是关于hdu1576-A/B-(同余定理+乘法逆元+费马小定理+快速幂)的主要内容,如果未能解决你的问题,请参考以下文章