关于逆元的概念、用途和可行性的思考(附51nod 1013 和 51nod 1256)
【逆元的概念】
逆元和单位元这个概念在群中的解释是: 逆元是指数学领域群G中任意一个元素a,都在G中有唯一的逆元a',具有性质a×a'=a'×a=e,其中e为该群的单位元。
群的概念是: 如果独异点(幺半群)中每一个元素都有逆元,那么这个独异点(幺半群)叫做群。
独异点(幺半群): 有单位元的半群。
半群: 可结合的代数系统。即
,有
。
代数系统:我的理解是代数系统包含一个数的集合A和至少一个运算规则,所有的运算都是封闭的,不会产生不在A集合中的数。
我们知道的实数集合R和加减乘除等一系列运算规则就组成了一个代数系统。根据上面的概念我们当然知道这是一个群。
简单来说:对于任意群中元素a,b,如果a*b = 1 ,那么a就是b的左逆元,b就是a的右逆元。(如果这个群满足交换律,这个群就是交换群,那么a和b互为逆元。)
这里有一个例题,就是求逆元的:

当然这是一道单纯求逆元的题。
(K*M)% N = 1
看到这个我们想把%消掉,看起来就会简单了。
==> (K*M-1)%N = 0
==> K*M-1 = S*N (S为未知数)
现在我们成功的消掉了%,这个等式只有K和S两个未知数。
如果还没看出来的话,我们把K换成x,S换成y,再移项看看:
==> M*x - N*y = 1
是不是很熟悉,对,就是拓展欧几里得。
ll gcd(ll a,ll b,ll &x,ll &y){
if (b==){
x=,y=;
return a;
}
ll q=gcd(b,a%b,y,x);
y-=a/b*x;
return q;
}
这样能够解出x,y的一对解,再把它移到适合的范围内就得到了我们的结果。
这题的代码如下:
#include <iostream>
#include <algorithm>
#include <string>
#include <stdio.h>
#include <string.h>
#include <math.h>
#define rep(u,i,n) for(int u = i;u <= n; u++)
typedef long long ll;
using namespace std;
ll gcd(ll a,ll b,ll &x,ll &y){
if (b==){
x=,y=;
return a;
}
ll q=gcd(b,a%b,y,x);
y-=a/b*x;
return q;
}
int main() {
ll n,m,x = ,y = ;
while(cin >> m >> n){
gcd(n,m,x,y);
if(y > ) cout << y << endl;
else cout << n+y << endl;
}
return ;
}
【逆元的用途】 除法取模
我们知道 (a*b)%n = c --> ((a%n)*(b%n))%n = c;
但是(a/b)%n 该怎么求呢?
如果n = 11, a = 3, b = 10 的话,直接算会导致结果错误(3/10)%11 = 0。
我们知道3/10是有值的,但是直接算结果会变成0,肯定出了某种错误。
这个错误我们暂时不做讨论,着重解决问题。
这时乘法逆元就派上用场了我们知道(3/10)%11 ==> (3*(1/10))%11
而1/10在乘法上是10的逆元(mod n = 11),意思就是我们用10的逆元取代1/10的位置就可以解决了。
(3/(10的逆元))%11就是我们要的结果。
于是我们成功的解决了除法时取模的问题。
这里有一个例题:

这个逆元是手动求的,懒得写求逆元代码。
代码如下:
#include <bits\stdc++.h>
typedef long long ll;
using namespace std;
const int mod = ;
ll mod_pow(ll x,ll n)
{
ll ans = ;
while(n > ) {
if(n % == ){
ans = ans * x % mod;
}
n /= ;
x = x*x % mod;
}
return ans;
}
int main()
{
ll n,ans;
cin >> n;
n++;
ans = (mod_pow(,n)-)*%mod; // 500000004是2对mod的逆元 ,逆元在除以后取模时使用
cout << ans << endl;
return ;
}
关于逆元的概念、用途和可行性的思考(附51nod 1013 和 51nod 1256)的更多相关文章
- 51Nod 1013 3的幂的和 快速幂 | 乘法逆元 | 递归求和公式
1.乘法逆元 直接使用等比数列求和公式,注意使用乘法逆元 ---严谨,失细节毁所有 #include "bits/stdc++.h" using namespace std; #d ...
- 51nod 1013【快速幂+逆元】
等比式子: Sn=(a1-an*q)/(1-q) n很大,搞一发快速幂,除法不适用于取膜,逆元一下(利用费马小定理) 假如p是质数,且gcd(a,p)=1,那么 a^(p-1)≡1(mod p).刚好 ...
- 51Nod 1013 3的幂的和(快速幂+逆元)
#include <iostream> #include <algorithm> #include <string> #define MOD 1000000007 ...
- POJ1845:Sumdiv(求因子和+逆元+质因子分解)好题
题目链接:http://poj.org/problem?id=1845 定义: 满足a*k≡1 (mod p)的k值就是a关于p的乘法逆元. 为什么要有乘法逆元呢? 当我们要求(a/b) mod p的 ...
- luogu3811 【模板】乘法逆元
题目大意:给出n,求1~n所有数的乘法逆元. 乘法逆元的概念是:如果b*rev(b)≡1 (mod p),p与b互质,则rev(b)就是b的模p乘法逆元.乘法逆元往往用于除法取模. 具体操作详见htt ...
- poj 1845 【数论:逆元,二分(乘法),拓展欧几里得,费马小定理】
POJ 1845 题意不说了,网上一大堆.此题做了一天,必须要整理一下了. 刚开始用费马小定理做,WA.(poj敢说我代码WA???)(以下代码其实都不严谨,按照数据要求A是可以等于0的,那么结果自然 ...
- 了解dto概念,什么是DTO
了解dto概念 此博文收集整理了一些主流的文章对于DTO模式的解读,他们大体相似而又各有所不同.对于设计模式的解读也是一个仁者见仁智者见智的事情,不过设计模式往往都是前辈们在遇到一类特定的问题下而 ...
- 001_C/C++笔试题_考察C/C++语言基础概念
(一)文章来自:C/C++笔试题-主要考察C/C++语言基础概念.算法及编程,附参考答案 (二)基础概念 2. 头文件中的ifndef/define/endif的作用? 答:防止该头文件被重复引用. ...
- PMP用语集
AC actual cost 实际成本 ACWP actual cost of work performed 已完工作实际成本 BAC budget at completion 完工预算 BCWP b ...
随机推荐
- Java基础---其他对象
第一讲 System类 一.概述 1.System是描述系统一些信息的类,类中的属性和方法都是静态的.不能被实例化,没有提供构造函数. 2.字段摘要 out:标准输出流.默认是控制台. in: ...
- javascript实现jsp页面的打印预览
1.加载WebBrowser打印预览控件 <OBJECT classid="CLSID:8856F961-340A-11D0-A96B-00C04FD705A2" heigh ...
- MySQL(二)之服务管理与初始化文件修改和连接MySQL
上一篇给大家介绍了怎么在linux和windows中安装mysql,本来是可以放在首页的,但是博客园说“安装配置类文件”不让放在首页.接下来给大家介绍一下在linux和windows下MySQL的一下 ...
- C++类静态成员与类静态成员函数
当将类的某个数据成员声明为static时,该静态数据成员只能被定义一次,而且要被同类的所有对象共享.各个对象都拥有类中每一个普通数据成员的副本,但静态数据成员只有一个实例存在,与定义了多少类对象 ...
- 配置文件http://www.cnblogs.com/Jevon-ran/p/7112007.html
author:JevonWei 版权声明:原创作品 /etc/centos-release 记录系统版本信息 /etc/issue 系统登录前提示信息 /etc/motd 系统登陆后提示信息 /etc ...
- 关联本地文件夹到github项目
git init git remote add origin https://自己的仓库url地址 git status git add . git commit -m '[提交内容的描述]' 先 p ...
- SourceTree使用方法介绍
SourceTree比命令行更容易操作,能更直观看到发生了什么.但是没有哪一家git图形化软件能完成git的所有操作,封装后的使用也隐藏了git的一些细节,在图形化工具出现一些非常罕见的情况时,还是需 ...
- JS查错小工具-三生有幸【推荐】
H5和CSS语言在开发者官网上都有在线查错工具,同样的,更加复杂的JavaScript也需要一个查错工具,(别指望DreamWeaver了,debug功能做的太垃圾,还不如Firefox自带的强..) ...
- (4)ES6解构赋值-字符串篇
字符串的解构赋值 let [a,b,c,d,e] = 'Apple'; console.log(a); //A console.log(b); //p console.log(c); //p cons ...
- 团队作业8——第二次项目冲刺(Beta阶段)--第六天
一.Daily Scrum Meeting照片 二.燃尽图 三.项目进展 学号 成员 贡献比 201421123001 廖婷婷 16% 201421123002 翁珊 16% 201421123004 ...