51 Nod 1013 3的幂的和 矩阵链乘法||逆元+快速幂

这道题我写了两种写法

一种利用逆元 a/b%mod=a*c%mod; (c是b的逆元）易得2的逆元就是5~~~04；

一种是矩阵快速幂 利用递推式得出结论

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int mod=;
int read(){
    int ans=,f=,c=getchar();
    while(c<''||c>''){if(c=='-') f=-; c=getchar();}
    while(c>=''&&c<=''){ans=ans*+(c-''); c=getchar();}
    return ans*f;
}
int n;
int qmod(long long a,int b){
    long long ans=; a%=mod;
    while(b){
        if(b&) ans=ans*a%mod;
        b>>=;
        a=a*a%mod;
    }
    return ans;
}
int main()
{
    n=read()+;
    printf("%lld\n",(qmod(,n)-)*(long long)%mod);
    return ;
}

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
using namespace std;
const int mod=;
int read(){
    int ans=,f=,c=getchar();
    while(c<''||c>''){if(c=='-') f=-; c=getchar();}
    while(c>=''&&c<=''){ans=ans*+(c-''); c=getchar();}
    return ans*f;
}
typedef LL mat[][];
int n;
void quickmod(mat a,mat b){
    mat c={};
    for(int i=;i<;i++)
     for(int k=;k<;k++)
      for(int j=;j<;j++)
       c[i][j]=(c[i][j]+a[i][k]*b[k][j])%mod;
    for(int i=;i<;i++)
     for(int j=;j<;j++)
      a[i][j]=c[i][j];
}
void fastpow(int n){
    mat a={,,,},b={,,,};
    while(n){
        if(n&) quickmod(a,b);
        n>>=;
        quickmod(b,b);
    }
    printf("%d\n",(a[][]+a[][])%mod);
}
int main()
{
    n=read(); fastpow(n);
    return ;
}