Hdu 1576

給出n=A%9973,B,求出(A/B)%9973;

設k=(A/B)%9973;
(A/B)=k+9973y(y未知)
  A=kB+9973yB
在等式兩邊%9973
n=A%9973
n=kB%9973
kB-9973y=n
gcd(B,9973)=1;
所以套擴展歐幾里得再乘上n即可;

#include <cstdio>
#include <iostream>
using namespace std;

typedef long long LL;
LL A,B,k,n,t;
LL x,y;
const int m=9973;

void Exgcd(LL a,LL b,LL &x,LL &y)
{
    if(!b)  {x=1,y=0;return ;}
    else Exgcd(b,a%b,y,x),y-=(a/b)*x;
}

int main()
{
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld%lld",&n,&B);
        Exgcd(B,m,x,y);
        x=((x%m)+m)%m;
        printf("%lld\n",(((x*n)%m)+m)%m);
    }
    return 0;
}
發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章