NOIp2012TG/Luogu P1082 同余方程 解题报告

原创 2016年08月28日 17:54:28

这是一道数论题,是扩展欧几里得算法的裸题
博主为了让自己记住这个算法,特地来写一篇博文

下面来看看一下题:

题目描述

求关于 x 的同余方程 ax1 (mod b)的最小正整数解。

输入输出格式

输入格式:
输入只有一行,包含两个正整数 a, b,用一个空格隔开。

输出格式:
输出只有一行,包含一个正整数 x0,即最小正整数解。输入数据保证一定有解。

输入输出样例

输入样例#1:
3 10
输出样例#1:
7

说明

【数据范围】

对于 40%的数据,2 ≤b≤ 1,000;

对于 60%的数据,2 ≤b≤ 50,000,000;

对于 100%的数据,2 ≤a, b≤ 2,000,000,000。

NOIP 2012 提高组 第二天 第一题


解题思路

扩展欧几里得算法

分析:
ax1(modb)可以通过移项变为:ax10(modb) .
然后再次变为ax1=bt 。即:axbt1
所以用扩展欧几里德方法解就可以了。

下面是C++代码:

#include <iostream>
#include <cstdio>
using namespace std;
long long x,y,a,b,ans;
long long gcd(long long a,long long b)//扩展欧几里得算法
{
    long long t,ret;
    if(b==0)//欧几里得算法
    {
        x=1;
        y=0;
        return a;
    }
    ret=gcd(b,a%b);
    t=y;//扩展欧几里得
    y=x-(a/b)*y;
    x=t;
    return ret;
}
int main()
{
    cin>>a>>b;
    ans=gcd(a,b);
    while(x>b)x-=b; //输出最小正整数解
    while(x<0)x+=b;
    cout<<x<<endl;
    return 0;
}

写给自己:如果无法理解,建议背过代码

感谢惠读

版权声明:本文为博主原创文章,未经博主允许不得转载。 举报

相关文章推荐

NOIp1998PJ/Luogu P1010 幂次方 解题报告

关于此题,我要感谢CoolTeam的题解,帮助我完成了此题 题目描述任何一个正整数都可以用2的幂次方表示。例如137=2^7+2^3+2^0 同时约定方次用括号来表示,即a^b 可表示为...

Luogu P1569 KC与龙珠 解题报告

不费话了,直接切入正题。下面看一下题目: 题目背景KC与龙珠题目描述ZKC大神和LSH大神省队集训时住在同一个房间,他们有时候会一起去楼下的85℃喝奶茶,同时在电脑上玩《龙珠》。《龙珠》中有许多不同的...

精选:深入理解 Docker 内部原理及网络配置

网络绝对是任何系统的核心,对于容器而言也是如此。Docker 作为目前最火的轻量级容器技术,有很多令人称道的功能,如 Docker 的镜像管理。然而,Docker的网络一直以来都比较薄弱,所以我们有必要深入了解Docker的网络知识,以满足更高的网络需求。

POJ 1061 青蛙(蛤蟆)的约会

青蛙的约会 Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 113194 Accepted: 23080...

The 14th Zhejiang University Programming Contest(未完工)

The 14th Zhejiang University Programming Contest
  • RaynX
  • RaynX
  • 2014-05-10 23:02
  • 528

HDU3605 匈牙利匹配

这题的时间 卡得很彻底....各种循环优化,外加一个前驱记录,匈牙利匹配的运用。 #include #include #include using namespace std; int n, ...

hdoj-1008-Elevator(解题报告)

题目来源:http://acm.hdu.edu.cn/showproblem.php?pid=1008 问题描述: Elevator Time Limit: 2000/1000 MS (Java/...

Matrix Multiplication解题报告

题目摘要:Johnny and John are good friends.Johnny is going to take the entrance exams for postgraduate sc...

Elevator Trouble解题报告

题目摘要:You are on your way to your firstjob interview as a program tester, and you are already late. T...

1165 纪念品分组(要认真看题目!!!!!!!!!!)

纪念品分组 时间限制(普通/Java):1000MS/3000MS          运行内存限制:65536KByte 总提交:35            测试通过:19 描述 A...

hdoj-1283-最简单的计算机(解题报告)

题目来源:http://acm.hdu.edu.cn/showproblem.php?pid=1283 问题描述: 最简单的计算机 Time Limit: 2000/1000 MS (Java/O...
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

(最多只允许输入30个字)