博客
关于我
hdu 1757 A Simple Math Problem(矩阵快速幂)
阅读量:136 次
发布时间:2019-02-27

本文共 1427 字,大约阅读时间需要 4 分钟。

分析

f的递推关系

在这里插入图片描述

代码

#include
using namespace std;#define mod(x) ((x)%m)int k, m;struct mat{ int d[10][10]; mat operator*(const mat x){ mat ret; int tmp; for(int i = 0; i < 10; i++){ for(int j = 0; j < 10; j++){ tmp = 0; for(int k = 0; k < 10; k++){ tmp = mod(tmp + d[i][k]* x.d[k][j]); } ret.d[i][j] = tmp; } } return ret; } void init_unit(){ for(int i = 0; i < 10; i++) for(int j = 0; j < 10; j++) d[i][j] = i == j ? 1 : 0; } void init(){ for(int i = 0; i < 10; i++) for(int j = 0; j < 10; j++) d[i][j] = i == j+1 ? 1 : 0; } void set(int i, int v){ d[9-i][9] = v; }}a;mat fastPow(mat base, int pow){ mat res; res.init_unit(); while(pow){ if(pow & 1) res = res * base; base = base * base; pow >>= 1; } return res;}int main(){ int tmp; while(scanf("%d%d", &k, &m) == 2){ if(k < 10){ // 猜测:本题无此测试数据 printf("%d\n", k%m); continue; } a.init(); for(int i = 0; i < 10; i++) scanf("%d", &tmp), a.set(i, tmp); a = fastPow(a, k-9); int ans = 0; for(int i = 0; i < 10; i++) ans = mod(ans+a.d[i][9]*i); printf("%d\n", ans); } return 0;}

转载地址:http://dwib.baihongyu.com/

你可能感兴趣的文章
mysql generic安装_MySQL 5.6 Generic Binary安装与配置_MySQL
查看>>
Mysql group by
查看>>
MySQL I 有福啦,窗口函数大大提高了取数的效率!
查看>>
mysql id自动增长 初始值 Mysql重置auto_increment初始值
查看>>
MySQL in 太多过慢的 3 种解决方案
查看>>
Mysql Innodb 锁机制
查看>>
MySQL InnoDB中意向锁的作用及原理探
查看>>
MySQL InnoDB事务隔离级别与锁机制深入解析
查看>>
Mysql InnoDB存储引擎 —— 数据页
查看>>
Mysql InnoDB存储引擎中的checkpoint技术
查看>>
Mysql InnoDB存储引擎中缓冲池Buffer Pool、Redo Log、Bin Log、Undo Log、Channge Buffer
查看>>
MySQL InnoDB引擎的锁机制详解
查看>>
Mysql INNODB引擎行锁的3种算法 Record Lock Next-Key Lock Grap Lock
查看>>
mysql InnoDB数据存储引擎 的B+树索引原理
查看>>
mysql interval显示条件值_MySQL INTERVAL关键字可以使用哪些不同的单位值?
查看>>
mysql problems
查看>>
MySQL replace函数替换字符串语句的用法(mysql字符串替换)
查看>>
mysql workbench6.3.5_MySQL Workbench
查看>>
MySQL Workbench安装教程以及菜单汉化
查看>>
MySQL Xtrabackup 安装、备份、恢复
查看>>