hdu 1576 A/B 拓展欧几里得算法

A/B

Time Limit: 1000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others)
Total Submission(s): 2390    Accepted Submission(s): 1731

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)。

Output

对应每组数据输出(A/B)%9973。

Sample Input

2
1000 53
87 123456789

Sample Output

7922
6060

Author

xhd

  对于拓欧我用的一点也不熟练,特别是限制解必须为正数时,而本题规定了b,9973互质,直接取模至正数,还变简单了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long qword;
qword ext_gcd(qword a,qword b,qword &x,qword &y)
{
        if (a%b==0)
        {
                //a*x+b*y==b
                x=0;y=1;
                return y;
        }
        qword ret=ext_gcd(b,a%b,x,y);
        qword tx=x,ty=y;
        x=ty;
        y=tx-a/b*ty;
        return ret;
}
int main()
{
        //freopen("input.txt","r",stdin);
        //A=9973*x+n
        //(9973*x+n)=y*B
        //9973*x-B*y==-n
        qword n,a,b,x,y,yy,xx;
        int nn;
        scanf("%d",&nn);
        qword g;
        while (nn--)
        {
                scanf("%I64d%I64d",&n,&b);
                g=ext_gcd(9973,b,x,y);
                x*=-n;y*=n;
            //    cout<<9973*x-b*y<<endl;
                yy=(y%9973+9973)%9973;
                xx=x-(y-yy)/9973*b;
            //    cout<<9973*xx-b*yy<<endl;
                cout<<yy<<endl;
        }
}
时间: 2024-10-14 10:28:59

hdu 1576 A/B 拓展欧几里得算法的相关文章

HDU - 1576 A/B(扩展欧几里得算法)

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=1576 题意:要求(A/B)%9973,但由于A很大,我们只给出n(n=A%9973)(我们给定的A必能被B整除,且gcd(B,9973) = 1). 普通版欧几里得算法(辗转相除): 1 typedef long long LL; 2 LL gcd(LL a,LL b){ 3 return (b==0) ? a : gcd(b,a%b); 4 } 扩展欧几里得算法(理论):对于不完全为0的非负整数,

HDU - 1356 The Balance(拓展欧几里得算法的解空间结构)

题目: Description Ms. Iyo Kiffa-Australis has a balance and only two kinds of weights to measure a dose of medicine. For example, to measure 200mg of aspirin using 300mg weights and 700mg weights, she can put one 700mg weight on the side of the medicin

HDU 1576 A/B(拓展欧几里得)

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). Output 对应每组数据输出(A/B)%9973. Sample Input 2 1000 53 87 123456789 Sample Output 7922 60

[数论]拓展欧几里得算法

欧几里得算法(辗转相除法) 用来求解最大公约数 1 int gcd(int a,int b){ 2 return b ? gcd(b,a%b) : a; 3 } 在 #include<algorithm> 中也可以直接调用 __gcd(a,b) 拓展欧几里得算法 求解不定方程: 引理:存在 x , y 使得 ax+by=gcd(a,b) 设a,b,c为任意整数,若方程ax+by=c的一组解是(x0,y0),则它的任意整数解都可以写成(x0+k*b/gcd(a,b),y0-k*a/gcd(a,b

拓展欧几里得算法及代码实现

扩展欧几里得算法就是求: ax + by = gcd(a, b) 的一组整数解(x, y) 一.非递归的实现: 首先看a = 60, b = 22的情况: 表格左边是欧几里得算法,右边等式计算ax + by = gcd(a, b)的解 a = 2 × b + 16 16 = a - 2b b = 1 × 16 + 6 6 = b - 1 × 16 = b - 1 × (a - 2b) = -a + 3b 16 = 2 × 6 + 4 4 = 16 - 2 × 6 = (a - 2b) - 2 ×

数学问题——拓展欧几里得算法

一.拓展欧几里得算法 该算法用来解决这样一个问题:给定两个非零整数 a 和 b,求一组整数解 (x,y) ,使得 ax + by = gcd(a,b) 成立,其中 gcd(a,b) 表示 a 和 b 的最大公约数. 递归边界:当 b 为 0 时,此时的 a 就等于 gcd,显然有 a*1+b*0=gcd 成立,此时 x=1,y=0: 递推公式:设当计算 gcd(a,b) 时,有 ax1 + by1 = gcd 成立:而在下一步计算gcd(b,a%b) 时,又有 bx2 + (a%b)y2 = g

如何使用拓展欧几里得算法求解模线性方程组(详解)

式子a≡b(mod n)称为a和b关于模n同余,它的充要条件是a-b是n的整数倍,即a-b=zn(其中z取整数). 而模线性方程组ax≡b(mod n)可以写成ax-b=zn(其中z取整数),移项可得 ax-zn=b,也即二元一次方程ax+by=c的形式,利用拓展欧几里得算法(extgcd)可以求解该方程是否有解及其一组解,并可根据该组解写出解系,进而求出一个特解,比如最小正整数解. 下面给出拓展欧几里得算法的程序. 1 typedef long long LL; 2 void extgcd(L

Sumdiv|同余|约数|拓展欧几里得算法

目录 Sumdiv|同余|约数|拓展欧几里得算法 Problem 分析 约数个数定理部分 约数和定理部分 等比数列部分 题目分析 扩展欧几里得算法部分 Code 呕,我吐了. Sumdiv|同余|约数|拓展欧几里得算法 Problem \[ 求A^{B}的所有约数之和 \ mod \ 9901\left(1\leqslant A,B \leqslant 5*10^{7}\right) \] 分析 约数个数定理部分 定理内容: 对于一个大于1的正整数n可以分解质因数: 则n的正约数个数为: 定理证

对拓展欧几里得算法的一点理解

首先需要明确的一点是:这是一种算法,而非一个证明题. 算法的需求与数学证明题是不一样的,数学证明题要求严谨完整,而算法只需要证明我用到的某个的性质成立即可,相当于是“恰好发现了这一点”. 于是对于拓展欧几里得,我们是从欧几里得算法中发现了一个递推的性质,从而受到启发,产生猜想:可不可以利用递推求出二元一次方程的解? 我们把猜想建立在欧几里得算法之上,利用该算法的递推过程,贯穿该过程来得到想要解决问题的答案. 也就是说,我们仅仅需要证明我们得到的答案是正确的,而并非深刻挖掘欧几里得算法的过程. 拓