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

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

?????????????

?????????????????????????????????????????????????????????????

1. ???????

???????????????????????????????????????? $A$ ? $n$ ??????? $A^n = A \cdot A \cdot \ldots \cdot A$?? $n$ ? $A$??????????? $O(n)$ ??????????????????????????????????????? $O(\log n)$??????????

2. ????

???????????

#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; }};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;}

3. ????

  • ?????????

    ??? mat ????? 10x10 ?????????????????????????????????????

  • ?????

    • init_unit() ????????????? 1???? 0??
    • init() ????????????????????? 1???? 0?
    • set() ????????????????
  • ?????

    fastPow ???????????????????????????????????????????????????

  • ?????

    ?????? k ? m?????? a????????????? a ???????????????

  • 4. ??

    ???????????????????????????????????????????????????????????????????????????

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

    你可能感兴趣的文章
    Openlayers实战:绘制图形,导出geojson文件
    查看>>
    Openlayers实战:绘制图形,导出KML文件
    查看>>
    Openlayers实战:绘制多边形,导出CSV文件
    查看>>
    Openlayers实战:绘制带箭头的线
    查看>>
    Openlayers实战:输入WKT数据,输出GML、Polyline、GeoJSON格式数据
    查看>>
    Openlayers实战:非4326,3857的投影
    查看>>
    Openlayers高级交互(10/20):绘制矩形,截取对应部分的地图并保存
    查看>>
    Openlayers高级交互(11/20):显示带箭头的线段轨迹,箭头居中
    查看>>
    Openlayers高级交互(14/20):汽车移动轨迹动画(开始、暂停、结束)
    查看>>
    Openlayers高级交互(15/20):显示海量多边形,10ms加载完成
    查看>>
    Openlayers高级交互(16/20):两个多边形的交集、差集、并集处理
    查看>>
    Openlayers高级交互(17/20):通过坐标显示多边形,计算出最大幅宽
    查看>>
    Openlayers高级交互(19/20): 地图上点击某处,列表中显示对应位置
    查看>>
    Openlayers高级交互(2/20):清除所有图层的有效方法
    查看>>
    Openlayers高级交互(20/20):超级数据聚合,页面不再混乱
    查看>>
    Openlayers高级交互(3/20):动态添加 layer 到 layerGroup,并动态删除
    查看>>
    Openlayers高级交互(4/20):手绘多边形,导出KML文件,可以自定义name和style
    查看>>
    Openlayers高级交互(5/20):右键点击,获取该点下多个图层的feature信息
    查看>>
    Openlayers高级交互(6/20):绘制某点,判断它是否在一个电子围栏内
    查看>>
    Openlayers高级交互(7/20):点击某点弹出窗口,自动播放视频
    查看>>