1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77
|
#include <string> #include <cstring> #include <cstdio> using namespace std;
typedef long long ll;
const int N=2; ll n,p,m; struct { ll g[3][3]; }f,res;
void build(node &x) { for(int i=1;i<=N;i++) for(int j=1;j<=N;j++) if(i==j) x.g[i][j]=1LL; else x.g[i][j]=0LL; }
void cheng(node &x,node &y,node &z) { memset(z.g,0,sizeof(z.g)); for(int i=1;i<=N;i++) for(int j=1;j<=N;j++) if(x.g[i][j]) { for(int k=1;k<=N;k++) { z.g[i][k]+=x.g[i][j]*y.g[j][k]; if(z.g[i][k]>=m) z.g[i][k]%=m; } } }
void kasumi(int k) { build(res); node tmp=f,t; while(k) { if(k&1) { cheng(res,tmp,t); res=t; } cheng(tmp,tmp,t); tmp=t; k>>=1; } }
ll solve() { if(n<=2) return 1ll; kasumi(n-2); ll ret=res.g[1][1]+res.g[2][1]; if(ret>=m) ret-=m; return ret; } int main() { cin>>n>>p; m=p; f.g[1][1]=f.g[1][2]=f.g[2][1]=1; f.g[2][2]=0; int res=(int)solve(); cout<<res<<endl; return 0; }
|
近期评论