A/B
Time Limit: 1000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 5121 Accepted Submission(s): 3972
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
Author
xhd
Source
Recommend
linle
#include <iostream>
using namespace std;
typedef long long int ll;
/* (A/B)%9973=k 等效于 A/B=9973*x+k */
/* 由于A%9973=n 可得 k*B%9973=n */
/* 即(k/n)*B+(-y/n)*9973=gcd(B,9973)=1 */
/* 即 k*B+(-y)*9973=n */
ll exgcd(ll a,ll b,ll &x,ll &y)
{
ll temp,ans;
if(b==0)
{
x=1;
y=0;
return a;
}
ans=exgcd(b,a%b,x,y);
temp=x;
x=y;
y=temp-a/b*y;
return ans;
}
ll cal(ll a,ll b,ll c)
{
ll x,y;
ll gcd=exgcd(a,b,x,y);
if(c%gcd!=0) return -1;
x*=c/gcd;
b/=gcd;
if(b<0) b=-b;
ll ans=x%b;
if(ans<=0) ans+=b;
return ans;
}
int main()
{
ll T,n,B;
while(cin>>T)
{
while(T--)
{
cin>>n>>B;
ll ans=cal(B,9973,n);
cout<<ans<<endl;
}
}
return 0;
}