博客
关于我
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执行计划【explain】,看这一篇就够啦!
查看>>
Mysql执行计划字段解释
查看>>
mysql执行计划怎么看
查看>>
MySQL执行计划解读
查看>>
mysql执行顺序与索引算法
查看>>
mysql批量update优化_Mysql中,21个写SQL的好习惯,你值得拥有呀
查看>>
mysql批量update操作时出现锁表
查看>>
MYSQL批量UPDATE的两种方式
查看>>
mysql批量修改字段名(列名)
查看>>
MySQL批量插入数据遇到错误1213的解决方法
查看>>
mysql技能梳理
查看>>
MySQL报Got an error reading communication packets错
查看>>
Mysql报错Can‘t create/write to file ‘/tmp/#sql_3a8_0.MYD‘ (Errcode: 28 - No space left on device)
查看>>
MySql报错Deadlock found when trying to get lock; try restarting transaction 的问题解决
查看>>
MySQL报错ERROR 1045 (28000): Access denied for user ‘root‘@‘localhost‘
查看>>
Mysql报错Packet for query is too large问题解决
查看>>
mysql报错级别_更改MySQL日志错误级别记录非法登陆(Access denied)
查看>>
Mysql报错:too many connections
查看>>
MySQL报错:无法启动MySQL服务
查看>>
mysql授权用户,创建用户名密码,授权单个数据库,授权多个数据库
查看>>