[BZOJ2326] [HNOI2011] 数学作业 (矩阵乘法)

Description

Input

Output

Sample Input

Sample Output

HINT

Source

Solution

  递推式长这样:$f[n]=f[n-1]*10^k+n$

  对于每一段位数个数相同的$n$(如$10\sim99,100\sim999,23333\sim66666,1018701389\sim2147483647$),$k$是个定值

  然后就可以开心地分段矩阵乘法了,剩下的自己推吧

 1 #include <bits/stdc++.h>
 2 using namespace std;
 3 typedef long long ll;
 4 int mod;
 5 struct mat
 6 {
 7     ll a[4][4];
 8     int n, m;
 9  
10     mat()
11     {
12         memset(a, 0, sizeof(a));
13         n = 0, m = 0;
14     }
15  
16     mat(int x, int y)
17     {
18         memset(a, 0, sizeof(a));
19         n = x, m = y;
20     }
21  
22     mat operator* (const mat &rhs) const
23     {
24         mat ans;
25         ans.n = n, ans.m = rhs.m;
26         for(int i = 1; i <= n; ++i)
27             for(int j = 1; j <= rhs.m; ++j)
28                 for(int k = 1; k <= m; ++k)
29                     ans.a[i][j] = (ans.a[i][j] + a[i][k] * rhs.a[k][j]) % mod;
30         return ans;
31     }
32  
33     mat operator^ (ll rhs) const
34     {
35         mat ans(n, n), b = *this;
36         for(int i = 1; i <= n; ++i)
37             ans.a[i][i] = 1;
38         for(; rhs; rhs >>= 1, b = b * b)
39             if(rhs & 1) ans = ans * b;
40         return ans;
41     }
42 };
43  
44 int main()
45 {
46     ll n, c;
47     mat ans(1, 3), b(3, 3);
48     scanf("%lld%d", &n, &mod);
49     ans.a[1][3] = 1;
50     for(int i = 1; i <= 3; ++i)
51         for(int j = 1; j <= i; ++j)
52             b.a[i][j] = 1;
53     for(ll i = 10; ; i *= 10)
54     {
55         b.a[1][1] = i % mod;
56         if(i <= n) c = i / 10 * 9;
57         else c = n - i / 10 + 1;
58         ans = ans * (b ^ c);
59         if(i > n) break;
60     }
61     printf("%lld\n", ans.a[1][1]);
62     return 0;
63 }
View Code

 

转载于:https://www.cnblogs.com/CtrlCV/p/5668799.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/254898.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

HALCON示例程序texture.hdev检测树木

小哥哥小姐姐觉得有用点个赞呗&#xff01; HALCON示例程序texture.hdev检测树木 示例程序源码&#xff08;加注释&#xff09; 关于显示类函数解释 dev_close_window () Interactive : 0 dev_close_window () read_image (MreutHill, ‘mreut_y’) get_image_size (MreutH…

1、python基础速成

基础模块 def prt(age,name):#函数定义 print("%s is %d 年龄 old"%(name,age)) if __name__"__main__":#程序入口 print("Hello World") prt(45,"gaici") 获取输入&#xff1a;使用input()函数 nameinput("you name &#x…

老男孩博客园杨海潮MySQL--MySQL机构逻辑2

转载于:https://blog.51cto.com/yanfeilai528/2103403

法国标致雪铁龙汽车公司采用通快碟片激光器进行焊接

发布日期&#xff1a;2011-10-14 来源&#xff1a;光电新闻网 发布人&#xff1a;星之球科技 摘要&#xff1a;3月11日消息&#xff0c;十一个碟片激光器&#xff08;disk laser&#xff09;将安装在标致雪铁龙集团的工厂&#xff0c;这家法国汽车制造商准备使用4千瓦的激光器…

h.264 rtp打包

(2011-05-27 08:44:13) 转载标签&#xff1a; 杂谈 payload,H.264 RTP payload 格式 on 2011-2-18 in 博文摘选 | 0 Comment 1. 网络抽象层单元类型 (NALU) NALU 头由一个字节组成, 它的语法如下: --------------- |0|1|2|3|4|5|6|7| -------- |F|NRI| Type | --------------…

jquery live hover绑定方法

$(".select_item span").live({mouseenter:function(){$(this).addClass("hover");},mouseleave:function(){$(this).removeClass("hover");} }); 注意&#xff1a;jquery1.9以上版本不支持live&#xff0c;新方法为on 转载于:https://www.cnblo…

HALCON示例程序vessel.hdev血管的分割与测量

小哥哥小姐姐觉得有用点个赞呗&#xff01; HALCON示例程序vessel.hdev血管的分割与测量 示例程序源码&#xff08;加注释&#xff09; 关于显示类函数解释 dev_update_window (‘off’) dev_close_window () dev_open_window (0, 0, 512, 512, ‘black’, WindowID) set_d…

电子凸轮

CAM功能是按照一种人为预先设定的曲线关系(可以在线修改,对SEW的变频/伺服控制器而言)来运动的控制应用。 100%速度前馈的位置控制这个观点偶不敢苟同.典型的一些应用。比如:全自动包装机械上,移动锯,其实大家说的电子齿轮&#xff0c;指的就是一种可以调节主从速度比的同步应用…

浙南联合训练赛20180414

这次题目的代码都不长&#xff0c;CF的一贯风格 A - Game CodeForces - 513A Two players play a simple game. Each player is provided with a box with balls. First players box contains exactly n1 balls and second players box contains exactly n2balls. In one move…

原生JS实现苹果菜单

今天分享下用原生JS实现苹果菜单效果&#xff0c;这个效果的重点有以下几点 图标中心点到鼠标的距离的算法 利用比例计算图标的宽度 代码地址&#xff1a;https://github.com/peng666/blogs/blob/gh-pages/menus/index.html 在线测试地址&#xff1a;http://peng666.github.io/…

Gym 100090D Insomnia

从 n 变到 1&#xff0c;有多少种方案&#xff1f; 打表记忆化。 1 #include <bits/stdc.h>2 3 using namespace std;4 5 int n;6 int dp[1000005];7 int dfs(int n) {8 if(n1)9 return 1; 10 if(dp[n]>0) 11 return dp[n]; 12 int cnt0;…

halcon rectangle1_domain缩减图像域为矩形

目录rectangle1_domain&#xff08;算子&#xff09;描述参数rectangle1_domain&#xff08;算子&#xff09; rectangle1_domain - 将图像的域缩小为矩形。 rectangle1_domain&#xff08;Image&#xff1a;ImageReduced&#xff1a;Row1&#xff0c;Column1&#xff0c;Row…

PC+运动控制卡的控制方案

PC运动控制卡的控制方案&#xff1a; 采用PC&#xff0b;运动控制卡作为上位控制可充分利用计算机资源&#xff0c;用于运动过程、运动轨迹都比较复杂&#xff0c;且柔性比较强的机器和设备。从用户使用的角度来看&#xff0c;基于PC机的运动控制卡主要是硬件接口&#xff08;输…

IP/TCP/UDP/RTP/RTCP 包结构图

IP 包头结构: TCP 包头结构: UDP 包头结构: RTP 包头结构: RTCP 包头结构:

你可能不知道的java、python、JavaScript以及jquary循环语句的区别

一.概述 java循环语句分为四种形式&#xff0c;分别是 while, do/while, for, foreach&#xff1b; python中循环语句有两种&#xff0c;while&#xff0c;for&#xff1b; JavaScript中循环语句有四种&#xff0c;while&#xff0c;do/while&#xff0c;for&#xff0c;for/in…

webservices系列(二)——JAX-WS文件上传下载

新建ImgData类&#xff0c;存放文件javabean DataHandler&#xff1a;使用这个类型存放文件 XmlRootElement(name"ImaData") XmlAccessorType(XmlAccessType.FIELD) public class ImgData {private Integer id;XmlMimeType("application/octet-stream")pri…

halcon sobel边缘检测sobel_amp

目录sobel_amp&#xff08;算子&#xff09;描述参数sobel_amp&#xff08;算子&#xff09; sobel_amp - 使用Sobel算子检测边缘&#xff08;幅度&#xff09;。 sobel_amp&#xff08;图片&#xff1a;边缘图像&#xff1a;滤波器方式&#xff0c;掩膜大小:) 描述 sobel_…

es中的一些知识点记录

1. forcemerge接口 强制段合并&#xff0c;设置为1时&#xff0c;是期望最终只有1个索引段。但实际情况是&#xff0c;合并的结果是段的总数会减少&#xff0c;但仍大于1&#xff0c;可以多次执行强制合并的命令。 设置的的目标值越小。合并消耗的时间会越久。 curl -XPOST htt…

用live555和ffplay搭建流媒体环境

用live555和ffplay搭建流媒体环境 http://bbs.chinavideo.org/viewthread.php?tid12166

如何才能优雅地书写JS代码

第一&#xff1a;关于匿名函数的使用 要避免全局变量泛滥&#xff0c; 可以考虑使用匿名函数&#xff0c; 把不需要在外部访问的变量或者函数限制在一个比较小的范围内。 例如以下代码&#xff1a; <script> function func1(){ var list ["a", "b",…