51nod-1013 3的幂的和(乘法逆元)

最新推荐文章于 2019-10-23 18:24:35 发布

z岁月无声

最新推荐文章于 2019-10-23 18:24:35 发布

阅读量163

点赞数

分类专栏： 51Nod 数论文章标签：算法

本文链接：https://blog.csdn.net/C_13579/article/details/79763898

版权

51Nod 同时被 2 个专栏收录

90 篇文章 1 订阅

订阅专栏

数论

44 篇文章 0 订阅

订阅专栏

1013 3的幂的和

基准时间限制：1 秒空间限制：131072 KB 分值: 20 难度：3级算法题

关注

求：3^0 + 3^1 +...+ 3^(N) mod 1000000007

Input

输入一个数N(0 <= N <= 10^9)

Output

输出：计算结果

Input示例

Output示例

思路：列举前几个可以看出 ans= (f(n)+(f(n)-1)/2 )%MOD --f(n)=3^n

但是除法对于取余并没有封闭性，因此要把除法转换掉，就要用到乘法逆元.

乘法逆元：若数 x 是 a*x%p=1 则x为 a对于p的乘法逆元，记为 a^(-1)

即 a*a^(-1) = 1 (mod p)

乘法逆元的一大应用是模意义下的除法，除法在模意义下并不是封闭的，但我们可以根据上述公式，将其转化为乘法

x/y=x/y*(y*y^(-1))=x*y^(-1) (mod p)

利用费马小定理 a^(p-1)=1 (mod p) $p$ 为素数。

则 a^(p-1)=1 (mod p) --> a*a^(p-2)=1 (mod p)

由乘法逆元的定义 $a^(p-2$ ) 即为 $a$ 的乘法逆元。

使用快速幂计算 $a^(p-2),$ 时间复杂度为 $O (lo g a)$ 。

Code:

#include<iostream>
using namespace std;
typedef long long LL; 

const LL MOD=1000000007;
int n;

int pow(int n,LL p){
    LL ans=1,num=n;
    while(p){
    	if(p&1)	ans=ans*num%MOD;
    	num=num*num%MOD;
		p>>=1;
	}
    return ans;
}
int main()
{
	ios::sync_with_stdio(false);
	while(cin>>n){
		LL a=3,ans=1;
		while(n>0){
			if(n&1)	ans=ans*a%MOD;
			a=a*a%MOD;
			n>>=1;
		}
		LL bi=pow(2,MOD-2);
		ans=(3*ans+2*MOD-1)*bi%MOD;
		cout<<ans<<endl;
	}
	return 0;
}