• 乘法逆元整理


    乘法逆元

    1.费马小定理

    如果 $ a,p $ 互质,那么 $ a^{p-1} ≡ 1(mod p) $
    又由逆元方程知 $ a*x ≡ 1 (mod p) $ ,得到 $ ax ≡ a^{p-1}(mod p) $
    所以逆元 $ x $ 为 $ a^{p - 2} mod p $,快速幂求解即可。

        #include<iostream>
        #include<cstdio>
        #include<cstring>
        #include<algorithm>
        #include<cmath>
        
        using namespace std;
        
        #define LL long long
        
        LL n,p;
        
        inline LL fast_pow(LL a,LL b,LL p) {
        	LL ans = 1;
        	a %= p;
        	while(b) {
        		if(b&1) ans = ans * a % p;
        		a = a * a % p;
        		b >>= 1;
        	}
        	ans %= p;
        	return ans;
        }
        
        int main() {
        	scanf("%lld%lld",&n,&p);
        	for(int i = 1 ; i <= n ; i++) 
        		printf("%lld 
    ",fast_pow(i,p-2,p));
        	return 0;
        }
    

    2.线性求逆元

    这个不会证。
    不过不会证问题也不大,可以直接当一个结论用。
    $ inv[i] = p - frac{p}{i} * inv[p % i] % p$
    不过要注意 $ inv[1] = 1,inv[0] = tan90^o = 0 ( ) for $ 循环从2开始,时间复杂度 $ O(n) $ 。

    
      
      using namespace std;
      
      #define LL long long
      const int N = 3e6 + 10;
      
      LL n,p,inv[N];
      
      int main() {
      	scanf("%lld%lld",&n,&p);
      	inv[1] = 1;inv[0] = 0;
      	for(int i = 2; i <= n ; i++) 
      		inv[i] = p - (p / i) * inv[p % i] % p;
      	for(int i = 1 ; i <= n ; i++)
      		printf("%lld 
    ",inv[i]);
      	return 0;
      }
    
    
  • 相关阅读:
    SQL SERVER 2005 行转列
    为什么jQuery.get、jQuery.getJSON、jQuery.post无法返回JSON
    Silverlight 布局控件
    python socket connection
    linux下缓存的查看/修改
    openstackflat 网络问题
    Fail to start neutronserver
    copy module
    python ConfigParser
    itertools 介绍
  • 原文地址:https://www.cnblogs.com/Repulser/p/11207140.html
Copyright © 2020-2023  润新知