标签:Operations Binary const Matrix GukiZ int ++ ull define
显然我们要拆位, 因为每位都独立, 然后问题就变成能用dp求的东西,然后用矩阵快速幂优化一下。
注意mod为1的情况。
#include<bits/stdc++.h> #define LL long long #define fi first #define se second #define mk make_pair #define PLL pair<LL, LL> #define PLI pair<LL, int> #define PII pair<int, int> #define SZ(x) ((int)x.size()) #define ull unsigned long long using namespace std; const int N = 1e6 + 7; const int inf = 0x3f3f3f3f; const LL INF = 0x3f3f3f3f3f3f3f3f; const int mod = 1e9 + 7; const double eps = 1e-8; ull n, k, l, m; ull op[] = {1, 0, 1, 0}; struct Matrix { ull a[4][4]; Matrix() { memset(a, 0, sizeof(a)); } void init() { for(int i = 0; i < 4; i++) a[i][i] = 1; } Matrix operator * (const Matrix &B) const { Matrix C; for(int i = 0; i < 4; i++) for(int j = 0; j < 4; j++) for(int k = 0; k < 4; k++) C.a[i][j] = (C.a[i][j] + a[i][k] * B.a[k][j]) % m; return C; } Matrix operator ^ (ull b) { Matrix C; C.init(); Matrix A = (*this); while(b) { if(b & 1) C = C * A; A = A * A; b >>= 1; } return C; } } M; int main() { cin >> n >> k >> l >> m; ull ans = 1, ret0 = 0, ret1 = 0; M.a[0][0] = 1, M.a[0][1] = 0, M.a[0][2] = 1, M.a[0][3] = 0; M.a[1][0] = 0, M.a[1][1] = 1, M.a[1][2] = 0, M.a[1][3] = 1; M.a[2][0] = 1, M.a[2][1] = 0, M.a[2][2] = 0, M.a[2][3] = 0; M.a[3][0] = 0, M.a[3][1] = 1, M.a[3][2] = 1, M.a[3][3] = 1; Matrix tmp = M ^ (n - 1); for(int i = 0; i < 4; i += 2) for(int j = 0; j < 4; j++) ret0 = (ret0 + tmp.a[i][j] * op[j]) % m; for(int i = 1; i < 4; i += 2) for(int j = 0; j < 4; j++) ret1 = (ret1 + tmp.a[i][j] * op[j]) % m; for(int i = 0; i < l; i++) { if(k >> i & 1) ans = ans * ret1 % m, k ^= (1ll << i); else ans = ans * ret0 % m; } if(k) ans = 0; cout << ans % m << "\n"; return 0; } /* */
标签:Operations,Binary,const,Matrix,GukiZ,int,++,ull,define 来源: https://www.cnblogs.com/CJLHY/p/10356348.html
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。