hdu 2793 矩阵快速幂
2014-08-24 14:53
183 查看
构造矩阵
AC代码如下:
a[1][1] = 0;a[1][2] = 1;a[1][3] = 0;a[1][4] = 0; a[2][1] = 0;a[2][2] = 0;a[2][3] = 1;a[2][4] = 0; a[3][1] = 1;a[3][2] = 1;a[3][3] = 1;a[3][4] = 0; a[4][1] = 0;a[4][2] = 0;a[4][3] = 1;a[4][4] = 1;
AC代码如下:
#include <iostream> #include <cstring> #include <cstdio> #include <algorithm> using namespace std; const __int64 MAX_N = 5; const __int64 MOD = 1000000007; __int64 N; void multipy( __int64 a[MAX_N][MAX_N], __int64 b[MAX_N][MAX_N], __int64 c[MAX_N][MAX_N] ){ for( __int64 i = 1; i <= 4; i++ ){ for( __int64 j = 1; j <= 4; j++ ){ c[i][j] = 0; for( __int64 k = 1; k <= 4; k++ ){ c[i][j] = ( c[i][j] + a[i][k] * b[k][j] % MOD ) % MOD; } } } } void get_matrix_pow( __int64 a[MAX_N][MAX_N], __int64 n ){ __int64 ans[MAX_N][MAX_N] = {0}; __int64 temp[MAX_N][MAX_N]; for( __int64 i = 1; i <= 4; i++ ) ans[i][i] = 1; while( n ){ if( n % 2 == 1 ){ multipy( ans, a, temp ); memcpy( ans, temp, sizeof( __int64 ) * MAX_N * MAX_N ); } multipy( a, a, temp ); memcpy( a, temp, sizeof( __int64 ) * MAX_N * MAX_N ); n /= 2; } memcpy( a, ans, sizeof( __int64 ) * MAX_N * MAX_N ); } __int64 solve( __int64 n ){ __int64 a[MAX_N][MAX_N] = {0}; if( n < 0 ){ return 0; } if( n == 0 ){ return 1; }else if( n == 1 ){ return 2; } a[1][1] = 0;a[1][2] = 1;a[1][3] = 0;a[1][4] = 0; a[2][1] = 0;a[2][2] = 0;a[2][3] = 1;a[2][4] = 0; a[3][1] = 1;a[3][2] = 1;a[3][3] = 1;a[3][4] = 0; a[4][1] = 0;a[4][2] = 0;a[4][3] = 1;a[4][4] = 1; get_matrix_pow( a, n - 1 ); __int64 ans = 0; __int64 b[MAX_N]; b[1] = b[2] = b[3] = 1; b[4] = 2; for( __int64 i = 1; i <= 4; i++ ){ ans = ( ans + a[4][i] * b[i] % MOD ) % MOD; } return ans; } int main(){ __int64 A, B; while( scanf( "%I64d%I64d", &A, &B ) != EOF ){ printf( "%I64d\n", ( solve( B ) - solve( A - 1 ) + MOD ) % MOD ); } return 0; }
相关文章推荐
- 【矩阵快速幂】hdu 1757
- HDU 2807 The Shortest Path【矩阵的快速比较】
- hdu 2157 (K步路+矩阵快速幂)
- HDU 4565——So Easy!(矩阵快速幂乘法)
- HDU4549M-斐波那契数列(矩阵快速幂,二分幂)
- 【矩阵快速幂】hdu 1575
- hdu-4471-Homework-矩阵快速幂+优化加速
- HDU 1005 矩阵的快速二分幂
- hdu 3411 推公式+矩阵快速幂
- hdu 1575 (矩阵快速幂)
- hdu 2232 机器人的舞蹈 矩阵快速幂
- 数和矩阵快速幂hdu 3493
- HDU 4291:A Short problem_成都赛区网络赛—矩阵快速幂
- hdu 4291 矩阵快速幂+循环节
- HDU 2807 The Shortest Path(最短路+矩阵快速比较)
- hdu 1575 矩阵快速幂
- hdu 2254 奥运(矩阵快速幂)
- hdu 1575 Tr A(矩阵乘法快速求幂)
- HDU 4549 M斐波那契数列(矩阵快速幂+欧拉定理)
- hdu-4471-Homework-矩阵快速幂+优化加速