vijos1781 同余方程

描述

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

格式

输入格式

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

输出格式

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

样例1

样例输入1

 

样例输出1

 

限制

每个测试点1s

提示

对于40%的数据，2 ≤b≤ 1,000； 
对于60%的数据，2 ≤b≤ 50,000,000； 
对于100%的数据，2 ≤a, b≤ 2,000,000,000。

来源

Noip2012提高组复赛Day2T1

 

又是复习题……练练exgcd

感觉这么水的题发上去显得我很没有水平

算了……反正我也确实没水平

#include<cstdio>

#include<iostream>

#include<cstring>

#include<cstdlib>

#include<algorithm>

#include<cmath>

#include<queue>

#include<deque>

#include<set>

#include<map>

#include<ctime>

#define LL long long

#define inf 0x7ffffff

#define pa pair<int,int>

#define pi 3.1415926535897932384626433832795028841971

using namespace std;

inline LL read()

{

    LL x=0,f=1;char ch=getchar();

    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}

    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}

    return x*f;

}

int a,b,x,y;

inline int gcd(int a,int b)

{

	if (!b)return a;

	return gcd(b,a%b);

}

inline void exgcd(int a,int b,int &x,int &y)

{

	if (!b){x=1;y=0;return;}

	exgcd(b,a%b,x,y);

	int t=x;x=y;y=t-a/b*y;

}

int main()

{

	a=read();b=read();

	exgcd(a,b,x,y);

	x=(x%b+b)%b;

	printf("%d\n",x);

}