弱校联盟-ZZOIOJ-2180(乘法逆元详解)

题目

GJJ每天都要膜拜一发数学大佬,因为GJJ的数学太差了。这不,GJJ又遇到难题了,他想求助WJJ,但是WJJ这几天忙于追妹子,哪有时间给他讲题, 于是GJJ求助于热爱ACM的你,Acmer们能帮帮他吗?问题是求: k^0 + k^1 +…+ k^(n) mod p (0 < k < 100, 0 <= n <= 10^9, p = 1000000007)

例如:6^0 + 6^1 +…+ 6^(10) mod 1000000007 (其中k = 6, n = 10, p = 1000000007)

Input

输入测试数据有多组,每组输入两个整数k, n

Output

每组测试数据输出:Case #: 计算结果

Sample Input

2 1
6 10

Sample Output

Case 1: 3
Case 2: 72559411

HINT

Source

这道题当初没写出来,也没想到是等比数列,还是思维发散不行,还有一点就是在取模的时候有个除数为了精确就可以求乘法逆元转化为乘法来做

若对于数字A,C 存在X,使A * X = 1 (mod C) ,那么称X为 A 对C的乘法逆元。

乘法逆元怎么说呢,例如 :

12 / 4 mod 7 = ? 很显然结果是3
我们现在对于数对 (4,7), 可以知道 X = 2是 4 对7的乘法逆元即2*4=1(mod 7)
那么我们有(12 / 4) * (4 * 2 ) = (?) * (1) (mod 7)
除法被完美地转化为了乘法从而避免精度丢失

这里只说当 a p互质时(就是 a p最大公约数为1)
求法 一 费马小定理

如果a p互质 a^(p-1) mod p = 1

a的乘法逆元就是 a^(p-2)

算法时间复杂度O(log2 n) 就是个快速幂

求法二 扩展欧几里得

ax= 1mod p;
可以写成
ax - py =1
这就看出来了吧 时间复杂度 O(ln n)

那这道题就好写了

一 费马小定理

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<set>
#include<map>
#include<stack>
#include<queue>
#include<string>
#include<cstring>
#include<cmath>
#include<sstream>
#include<cstdlib>
using namespace std;
#define inf 0x3f3f3f3f
long long mod = 1000000007; 
long long poww(long long a,long long b)
{
    long long ans=1;
    while(b)
    {
        if(b&1) ans=ans*a%mod;
        b>>=1;
        a=a*a%mod;

    }
    return ans;
}

int main()
{
    int k,n;
    int tt=1;
    while(cin>>k>>n)
    {   cout<<"Case "<<tt++<<": ";
        if(k==1) 
        {
            cout<<(n+1)%mod<<endl;continue;
        }
        long long kk=poww(k,n+1)-1;
        long long nn=k-1;
        int ans=poww(nn,mod-2);  费马小定理 
        cout<<ans*kk%mod<<endl;


    }



    return 0;
}

二 扩展欧几里得

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<set>
#include<map>
#include<stack>
#include<queue>
#include<string>
#include<cstring>
#include<cmath>
#include<sstream>
#include<cstdlib>
using namespace std;
#define inf 0x3f3f3f3f
long long mod = 1000000007; 

void exgcd(int a,int b,long long &gg,long long &x,long long &y)
{
    if(!b) gg=a,x=1,y=0;
    else
    {
        exgcd(b,a%b,gg,y,x);
        y-=x*(a/b);
    }

}
int main()
{
    int k,n;
    int tt=1;
    while(cin>>k>>n)
    {   cout<<"Case "<<tt++<<": ";
        if(k==1) 
        {
            cout<<(n+1)%mod<<endl;continue;
        }
        long long kk=poww(k,n+1)-1;
        long long nn=k-1;
        long long gg,x,y;  //扩展欧几里得 gg  公约数 ax+by=z; 
        exgcd(nn,mod,gg,x,y);

         long long  ans=(x%mod+mod)%mod;
        cout<<ans*kk%mod<<endl;


    }



    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
系统根据B/S,即所谓的电脑浏览器/网络服务器方式,运用Java技术性,挑选MySQL作为后台系统。系统主要包含对客服聊天管理、字典表管理、公告信息管理、金融工具管理、金融工具收藏管理、金融工具银行卡管理、借款管理、理财产品管理、理财产品收藏管理、理财产品银行卡管理、理财银行卡信息管理、银行卡管理、存款管理、银行卡记录管理、取款管理、转账管理、用户管理、员工管理等功能模块。 文中重点介绍了银行管理的专业技术发展背景和发展状况,随后遵照软件传统式研发流程,最先挑选适用思维和语言软件开发平台,依据需求分析报告模块和设计数据库结构,再根据系统功能模块的设计制作系统功能模块图、流程表和E-R图。随后设计架构以及编写代码,并实现系统能模块。最终基本完成系统检测和功能测试。结果显示,该系统能够实现所需要的作用,工作状态没有明显缺陷。 系统登录功能是程序必不可少的功能,在登录页面必填的数据有两项,一项就是账号,另一项数据就是密码,当管理员正确填写并提交这二者数据之后,管理员就可以进入系统后台功能操作区。进入银行卡列表,管理员可以进行查看列表、模糊搜索以及相关维护等操作。用户进入系统可以查看公告和模糊搜索公告信息、也可以进行公告维护操作。理财产品管理页面,管理员可以进行查看列表、模糊搜索以及相关维护等操作。产品类型管理页面,此页面提供给管理员的功能有:新增产品类型,修改产品类型,删除产品类型。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值