HDU——T 1576 A/B
阿新 • • 發佈:2017-08-09
auth namespace ace key using u+ href chm esp
每組數據有兩個數n(0 <= n < 9973)和B(1 <= B <= 10^9)。
Output
對應每組數據輸出(A/B)%9973。
http://acm.hdu.edu.cn/showproblem.php?pid=1576
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 6418 Accepted Submission(s): 5075
Input 數據的第一行是一個T,表示有T組數據。
每組數據有兩個數n(0 <= n < 9973)和B(1 <= B <= 10^9)。
Sample Input 2 1000 53 87 123456789
Sample Output 7922 6060
Author xhd
Source HDU 2007-1 Programming Contest 膜的性質+逆元
(a/b)mod p=?
定義 c為b在mod p意義下的逆元
=a*c mod p = (a mod p * c mod p)mod p
① exgcd求逆元
1 #include <algorithm> 2#include <cstdio> 3 4 using namespace std; 5 6 const int mod(9973); 7 int T,n,b,x,y; 8 9 int exgcd(int a,int b,int &x,int &y) 10 { 11 if(!b) 12 { 13 x=1; y=0; 14 return a; 15 } 16 int ret=exgcd(b,a%b,x,y); 17 int tmp=x; x=y;18 y=tmp-a/b*y; 19 return ret; 20 } 21 22 int main() 23 { 24 scanf("%d",&T); 25 for(;T--;) 26 { 27 scanf("%d%d",&n,&b); 28 int gcd=exgcd(b,mod,x,y); 29 x*=gcd; 30 x=(x%mod+mod)%mod; 31 printf("%d\n",(x*n)%mod); 32 } 33 return 0; 34 }
② 歐拉定理求逆元
1 #include <algorithm> 2 #include <cstdio> 3 4 using namespace std; 5 6 #define LL long long 7 const int mod(9973); 8 int T,n,b; 9 10 LL phi(LL x) 11 { 12 LL ret=1; 13 for(LL i=2;i*i<=x;i++) 14 if(x%i==0) 15 { 16 x/=i; 17 ret*=i-1; 18 for(;x%i==0;) 19 ret*=i,x/=i; 20 } 21 if(x>1) ret*=x-1; 22 return ret; 23 } 24 LL Q_pow(LL a,LL b) 25 { 26 LL ret=1,base=a; 27 for(;b;b>>=1) 28 { 29 if(1&b) ret=(ret*base)%mod; 30 base=(base*base)%mod; 31 } 32 return ret; 33 } 34 35 int main() 36 { 37 scanf("%d",&T); 38 for(;T--;) 39 { 40 scanf("%d%d",&n,&b); 41 LL x=Q_pow((LL)b,(LL)phi((LL)mod)-1); 42 x=(x%mod+mod)%mod; 43 printf("%d\n",(x*n)%mod); 44 } 45 return 0; 46 }
HDU——T 1576 A/B