D - 3的幂的和 （扩展欧几里得求解逆元、费马小定理求解逆元）

D - 3的幂的和

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

Input

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

Output

输出：计算结果

Sample Input

3

Sample Output

40

解析：

扩展欧几里得求解逆元模板：

int x,y;
int gcd(int a,int b)//模版题 
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    int r=gcd(b,a%b);
    int t=y;
    y=x-(a/b)*y;
    x=t;
    return r;
}

费马小定理求解逆元模板：

long long quickpow ( long long a, long long b) 
{
    if (b<0) return 0;
    long long ret=1;
    a%=p;
    while(b) 
    {
      if(b&1) ret=(ret*a)%p ;
      b>>=1;
      a=(a*a)%p;
    }
    return ret ;
}
/*long long inv(long long a) 
{
   return quickpow (a, p - 2);//没看懂？？？
}*/

 

 

代码：

#include<iostream>
#include<math.h>
#include<stdio.h>
using namespace std;
int p=1000000007;
int x,y;
int gcd(int a,int b)//模版题 
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    int r=gcd(b,a%b);
    int t=y;
    y=x-(a/b)*y;
    x=t;
    return r;
}
long long quickpow ( long long a, long long b) 
{
    if (b<0) return 0;
    long long ret=1;
    a%=p;
    while(b) 
    {
      if(b&1) ret=(ret*a)%p ;
      b>>=1;
      a=(a*a)%p;
    }
    return ret ;
}
/*long long inv(long long a) 
{
   return quickpow (a, p - 2);
}*/
int main()
{
    int n;
    cin>>n;
    gcd(2,p);
    if(x<0) x+=p;
    long long sum=(quickpow(3,n+1)-1)*x%p; 
    cout<<sum<<endl;
return 0;
}