HDU4291A Short problem
来源:互联网 发布:cupid美瞳淘宝 编辑:程序博客网 时间:2024/05/20 23:35
随意选了一题,就作为博客生涯的开端吧^^
#include<stdio.h>typedef __int64 LL;/*g(n) % MOD1 = g(n % MOD2) % MOD1;//MOD2为循环节,以下如上g(n) % MOD2 = g(n % MOD3) % MOD2; g(g(g(n))) % MOD1= g(g(g(n)) % MOD2) % MOD1= g(g(g(n) % MOD3) % MOD2) % MOD1;n = g(n) % MOD3;n = g(n) % MOD2;n = g(n) % MOD1;*/const LL MOD1 = 1000000007;const LL MOD2 = 222222224;const LL MOD3 = 183120;const LL base[2][2] = {3, 1, 1, 0};LL ans[2][2], tmp[2][2];/*int main()//求循环节{ LL n0 = 0, n1 = 1, tt; LL mod = MOD2;//改变该值 for(int i = 1; ;i++) { tt = (3 * n1 + n0) % mod; n0 = n1; n1 = tt; if(n0 == 0 && n1 == 1) { printf("%d\n", i); break; } }}*/inline void Mult(LL pro[][2], LL mul1[][2], LL mul2[][2], LL MOD)//ans = mul1 * mul2;{ LL m1[2][2], m2[2][2]; for(int i = 0; i < 2; i++) for(int j = 0; j < 2; j++) m1[i][j] = mul1[i][j], m2[i][j] = mul2[i][j]; for(int i = 0; i < 2; i++) for(int j = 0; j < 2; j++) { pro[i][j] = 0; for(int k = 0; k < 2; k++) pro[i][j] += m1[i][k] * m2[k][j]; pro[i][j] %= MOD; }}inline void FastPow(LL n, LL MOD)//矩阵的base^n{ ans[0][0] = 1, ans[0][1] = 0; ans[1][0] = 0, ans[1][1] = 1;//单位矩阵 for(int i = 0; i < 2; i++) for(int j = 0; j < 2; j++) tmp[i][j] = base[i][j]; while(n) { if(n & 1) { Mult(ans, ans, tmp, MOD); } Mult(tmp, tmp, tmp, MOD); n >>= 1; }}int main(){ LL m; while(scanf("%I64d", &m) != EOF) { FastPow(m, MOD3); m = ans[1][0]; FastPow(m, MOD2); m = ans[1][0]; FastPow(m, MOD1); m = ans[1][0]; printf("%I64d\n", m); } return 0;}
- HDU4291A Short problem
- hdu4291A Short problem 矩阵快速幂
- hdu4291 A Short problem
- HDU4291-- A Short problem
- A Short problem HDU4291
- short-path problem (Dijkstra)
- short-path problem (Floyd)
- short-path problem (Spfa)
- hdu4291 A Short problem
- hdu 4291A Short problem
- hdu 4291 A Short problem
- hdu 4291 A Short problem
- HDU 4291 A Short problem
- hdoj-4921-A Short problem
- FOJ 2013 A short problem 解题报告
- FZU 2013 A short problem(DP)
- hdu 4291 A Short problem 矩阵+循环
- hdu 5107 K-short Problem(线段树)
- 程序员成长之路(一)_附教程下载
- 设置 Eclipse/MyEclipse Alt+/ 快速提示快捷键
- 关于页缓存的几点总结
- C# StringHelper
- JAVA---JDK环境变量的配置
- HDU4291A Short problem
- Mysql索引的使用
- 程序员成长之路(二)_附教程下载
- 10047 - The Monocycle
- C# 文件操作类1
- 文件操作类2
- Eclipse常用快捷键
- 弹出消息对话框类
- java.lang.OutOfMemoryError: PermGen space