扩展欧几里得,逆元初识(poj 1061+codeforce 7C line+hdu 1576 A/B)

2024-03-27 23:58

本文主要是介绍扩展欧几里得,逆元初识(poj 1061+codeforce 7C line+hdu 1576 A/B),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!




poj 1061 青蛙的约会:

#include <iostream>
#include<cstdio>
#define LL long long
using namespace std;
LL gcd(LL a, LL b){return b?gcd(b,a%b):a;
}
void extend_Euclid(LL a,LL b,LL &x,LL &y){
if(b == 0) {x = 1;y = 0;return ;
}
extend_Euclid(b,a%b,x,y);
LL tmp = x;  
x = y;
y = tmp - (a / b) * y;
}
LL res(int a,int b){while(a<0){if(b<0)a-=b;else a+=b;}return a;
}
int main()
{LL x,y,m,n,L;LL a,b,c,xi,yi;while(~scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&L)){a=m-n; b=L; c=y-x;LL d=gcd(a,b);if(c%d!=0){cout<<"Impossible\n";continue;}a/=d; b/=d; c/=d;extend_Euclid(a,b,xi,yi);xi=(xi*c)%b; //通解x不会比L长,所以要提前mod(b).xi=res(xi,b);cout<<xi<<endl;}return 0;
}

codeforce 7C line:

题意:判断一个线性方程组Ax+By+C=0是否有整数解,如果有,则输出整数解,否则输出-1.

并非是需要x一定要>0,利用扩展欧几里得求得x0,y0后再乘以相关倍数(-C/d)即可  [ d=gcd(A,B) ].

#include <iostream>
#include<cstdio>
#include<cmath>
#define LL long long
#include<cmath>
using namespace std;
void extend_Euclid(LL a,LL b,LL &f,LL &x,LL &y){
if(b == 0) {x = 1;y = 0;f=a;return ;
}
extend_Euclid(b,a%b,f,x,y);
LL tmp = x;
x = y;
y = tmp - (a / b) * y;
}
int main()
{LL A,B,C,f,x,y;while(cin>>A>>B>>C){extend_Euclid(A,B,f,x,y);if(C%f!=0){cout<<"-1\n";continue;}cout<<x*(-C/f)<<" "<<y*(-C/f)<<endl;}return 0;
}

hdu 1576 A/B

题意:要求(A/B)%9973,但由于A很大,我们只给出n(n=A%9973)(我们给定的A必能被B整除,且gcd(B,9973) = 1)。数据的第一行是一个T,表示有T组数据。
每组数据有两个数n(0 <= n < 9973)和B(1 <= B <= 10^9)。对应每组数据输出(A/B)%9973。

已知n=A%9973,设x=A/B,那么n=A-A/9973*9973=Bx-A/9973*9973,求出x大问题就解决了,最后为了防止负数产生,x=(x%9973+9973)%9973.

#include <iostream>
#include<cstdio>
#include<cmath>
#define LL long long
#include<cmath>
using namespace std;
void extend_Euclid(LL a,LL b,LL &f,LL &x,LL &y){
if(b == 0) {x = 1;y = 0;f=a;return ;
}
extend_Euclid(b,a%b,f,x,y);
LL tmp = x;
x = y;
y = tmp - (a / b) * y;
}
int main()
{//freopen("cin.txt","r",stdin);LL n,B,N;while(cin>>N){while(N--){LL a,b,c,f,x,y;scanf("%lld%lld",&n,&B);c=n;  b=9973;  a=B;extend_Euclid(a,b,f,x,y);x=x*(c/f);cout<<(x%9973+9973)%9973<<endl;}}return 0;
}


这篇关于扩展欧几里得,逆元初识(poj 1061+codeforce 7C line+hdu 1576 A/B)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/853837

相关文章

hdu4828(卡特兰数+逆元)

这题的前几个数据分别为1,2,5,14,32......................然后确定这是个卡特兰数列 下面来介绍下卡特兰数,它的递推式为f[i+1] = f[i]*(4*n - 6)/n,其中f[2] = f[3] =1;f[4] = 2;f[5] = 14;f[6] = 32.................................. 但是这题的n太大了,所以要用到逆元,

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

hdu4869(逆元+求组合数)

//输入n,m,n表示翻牌的次数,m表示牌的数目,求经过n次操作后共有几种状态#include<iostream>#include<algorithm>#include<cstring>#include<stack>#include<queue>#include<set>#include<map>#include<stdio.h>#include<stdlib.h>#includ

usaco 1.3 Mixing Milk (结构体排序 qsort) and hdu 2020(sort)

到了这题学会了结构体排序 于是回去修改了 1.2 milking cows 的算法~ 结构体排序核心: 1.结构体定义 struct Milk{int price;int milks;}milk[5000]; 2.自定义的比较函数,若返回值为正,qsort 函数判定a>b ;为负,a<b;为0,a==b; int milkcmp(const void *va,c

poj 3974 and hdu 3068 最长回文串的O(n)解法(Manacher算法)

求一段字符串中的最长回文串。 因为数据量比较大,用原来的O(n^2)会爆。 小白上的O(n^2)解法代码:TLE啦~ #include<stdio.h>#include<string.h>const int Maxn = 1000000;char s[Maxn];int main(){char e[] = {"END"};while(scanf("%s", s) != EO

hdu 2093 考试排名(sscanf)

模拟题。 直接从教程里拉解析。 因为表格里的数据格式不统一。有时候有"()",有时候又没有。而它也不会给我们提示。 这种情况下,就只能它它们统一看作字符串来处理了。现在就请出我们的主角sscanf()! sscanf 语法: #include int sscanf( const char *buffer, const char *format, ... ); 函数sscanf()和

hdu 2602 and poj 3624(01背包)

01背包的模板题。 hdu2602代码: #include<stdio.h>#include<string.h>const int MaxN = 1001;int max(int a, int b){return a > b ? a : b;}int w[MaxN];int v[MaxN];int dp[MaxN];int main(){int T;int N, V;s

poj 1511 Invitation Cards(spfa最短路)

题意是给你点与点之间的距离,求来回到点1的最短路中的边权和。 因为边很大,不能用原来的dijkstra什么的,所以用spfa来做。并且注意要用long long int 来存储。 稍微改了一下学长的模板。 stack stl 实现代码: #include<stdio.h>#include<stack>using namespace std;const int M

poj 3259 uva 558 Wormholes(bellman最短路负权回路判断)

poj 3259: 题意:John的农场里n块地,m条路连接两块地,w个虫洞,虫洞是一条单向路,不但会把你传送到目的地,而且时间会倒退Ts。 任务是求你会不会在从某块地出发后又回来,看到了离开之前的自己。 判断树中是否存在负权回路就ok了。 bellman代码: #include<stdio.h>const int MaxN = 501;//农场数const int

poj 1258 Agri-Net(最小生成树模板代码)

感觉用这题来当模板更适合。 题意就是给你邻接矩阵求最小生成树啦。~ prim代码:效率很高。172k...0ms。 #include<stdio.h>#include<algorithm>using namespace std;const int MaxN = 101;const int INF = 0x3f3f3f3f;int g[MaxN][MaxN];int n