扩展欧几里得,逆元初识(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

相关文章

css中的 vertical-align与line-height作用详解

《css中的vertical-align与line-height作用详解》:本文主要介绍了CSS中的`vertical-align`和`line-height`属性,包括它们的作用、适用元素、属性值、常见使用场景、常见问题及解决方案,详细内容请阅读本文,希望能对你有所帮助... 目录vertical-ali

Java常用注解扩展对比举例详解

《Java常用注解扩展对比举例详解》:本文主要介绍Java常用注解扩展对比的相关资料,提供了丰富的代码示例,并总结了最佳实践建议,帮助开发者更好地理解和应用这些注解,需要的朋友可以参考下... 目录一、@Controller 与 @RestController 对比二、使用 @Data 与 不使用 @Dat

Spring组件初始化扩展点BeanPostProcessor的作用详解

《Spring组件初始化扩展点BeanPostProcessor的作用详解》本文通过实战案例和常见应用场景详细介绍了BeanPostProcessor的使用,并强调了其在Spring扩展中的重要性,感... 目录一、概述二、BeanPostProcessor的作用三、核心方法解析1、postProcessB

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