博客
关于我
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日志
查看>>
mysql日志 事务问题_mysql因为事务日志问题无法启动
查看>>
mysql日志文件
查看>>
mysql日志管理学习笔记
查看>>
mysql日志问题定位实用命令
查看>>
MySQL日期时间函数大全
查看>>
mysql时间相减的问题
查看>>
mysql时间表示和计算
查看>>
MySQL是如何做容器测试的?
查看>>
mysql更改数据库表utf-8_修改mysql数据库为 utf-8
查看>>
mysql更改表引擎INNODB为MyISAM的方法总结
查看>>
mysql更新一个表里的字段等于另一个表某字段的值
查看>>
Mysql更新时间列只改日期为指定日期不更改时间
查看>>
MySQL更新锁(for update)摘要
查看>>
mysql更新频率_MySQL优化之如何了解SQL的执行频率
查看>>
mysql替换表的字段里面内容
查看>>
MySQL最多能有多少连接
查看>>
MySQL最大建议行数 2000w,靠谱吗?
查看>>
MySQL有哪些锁
查看>>
MySQL服务器安装(Linux)
查看>>