Week15-矩阵快速幂&字符串

最新推荐文章于 2022-07-03 23:13:05 发布
最新推荐文章于 2022-07-03 23:13:05 发布
阅读量630 收藏 5
点赞数 4
分类专栏： 程序设计
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_45727515/article/details/118034764
版权

目录

A : 斐波那契数列

题目描述
斐波那契数列的定义如下：
图片1.png

给出 n,pn,p ，求出 f(n)%pf(n)%p 的值。

输入格式
第一行一个数 T(1\le T\le 100)T(1≤T≤100)，表示数据组数。
接下来 TT 行，每行两个整数 n,p(1\le n\le 109,108\le p\le 10^9)n,p(1≤n≤10
9
,10
8
≤p≤10
9
)。

输出格式
输出 TT 行，表示每组数据的答案。

Case 1
Input
6
1 998244353
2 998244353
3 998244353
4 998244353
5 998244353
1000000000 998244353
Output
1
1
2
3
5
990450892

#include<iostream>
using namespace std;
long long mod;

class Mat{
public:
	long long m[2][2];
	Mat(int _a, int _b, int _c, int _d) {
		m[0][0] = _a; m[0][1] = _b;
		m[1][0] = _c; m[1][1] = _d;
	}
};

Mat multi(Mat& a,Mat& b) {
	Mat t(0, 0, 0, 0);
	for (int i = 0; i < 2; i++) {
		for (int j = 0; j < 2; j++) {
			for (int k = 0; k < 2; k++) {
				t.m[i][j] = t.m[i][j] + a.m[i][k] * b.m[k][j];
				t.m[i][j] %= mod;
			}
		}
	}
	return t;
}
long long pow(long long p) {
	Mat t(1, 0, 0, 1);
	Mat c(1, 1, 1, 0);
	while (p) {
		if (p & 1)t = multi(t, c);
		c = multi(c, c);
		p /= 2;
	}
	return t.m[0][1];
}

int main() {
	int T;
	cin >> T;
	while (T--) {
		int n;
		cin >> n >> mod;
		
		long long ans = pow(n);
		cout << ans << endl;
	}
	return 0;
}

在这里插入图片描述

B : 自然数幂和

题目描述
给定 𝑛 和 𝑘，计算 \Sigma{n}_{i=1}ikΣ
i=1
n
​
i
k
对 10^9+710
9
+7 取模的结果。

输入格式
第一行一个数 T(1\le T\le 100)T(1≤T≤100)，表示数据组数。
接下来 TT 行，每行两个整数 n,k(1\le n\le 10^9,1\le k\le 10)n,k(1≤n≤10
9
,1≤k≤10)。

输出格式
输出 TT 行，表示每组数据的答案。

Case 1
Input
6
1 10
2 2
2 3
3 2
3 3
1000000000 9
Output
1
5
9
14
36
12313161

#include<iostream>
#include<string.h>
#include<math.h>
using namespace std;
int n;//矩阵维数
struct Mat {
	long long m[101][101];
	void ini() {
		memset(m, 0, sizeof(m));
	}
	void ini_dw() {
		memset(m, 0, sizeof(m));
		for (int i = 0; i < n; i++)m[i][i] = 1;
	}
};

int mod=1e9+7;//模
Mat a, e;//a是输入矩阵 e是单位矩阵
Mat multi(Mat a, Mat b) {
	Mat c;
	c.ini();
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			for (int k = 0; k < n; k++) {
				c.m[i][j] = ((c.m[i][j] % mod) + (a.m[i][k] * b.m[k][j]) % mod) % mod;
			}
		}
	}
	return c;
}

Mat pow(Mat x, long long y) {
	Mat ans;
	ans.ini_dw();
	//if (y == 0 )return ans;
    //if(y==1)return x;
	while (y) {
		if (y & 1)ans = multi(ans, x);
		x = multi(x, x);
		y /= 2;
	}
	return ans;
}

long long C(int n, int m)//排列函数即C 的个数
{
	if (n < m)return C(m, n);
	long long i, c = 1;
	i = m;
	while (i != 0)
	{
		c *= n; n--; i--;
	}
	while (m != 0)
	{
		c /= m; m--;
	}
	return c;
}

long long fastPower(long long base, long long power) {
	long long result = 1;
	while (power > 0) {
		if (power & 1) {//此处等价于if(power%2==1)
			result = result * base % 1000;
		}
		power >>= 1;//此处等价于power=power/2
		base = (base * base) % 1000;
	}
	return result;
}


int main() {
	int T;
	cin >> T;
	while (T--) {
		int m, k;
		cin >> m >> k;
		n = k + 2;
		Mat t;
        t.ini();
        t.m[0][0]=1;
		for(int i=1;i<n;i++)t.m[0][i]=C(k,i-1);
        for(int i=1;i<n;i++){
            for(int j=0;j<n;j++){
                if(i==j)t.m[i][j]=1;
                else if(i>j)t.m[i][j]=0;
                else if(i<j)t.m[i][j]=C(k+1-i,j-i);
            }
        }
        long long res[11];
		t = pow(t, m - 1);
		long long ans = 0;
		for (int i = 0; i < k + 2; i++) {
			ans = (ans+t.m[0][i])%mod;
            }

		cout <<ans << endl;
	}
	return 0;
}

在这里插入图片描述

C : R??G??B??

题目描述
msy 的“显示器”终于完工啦！但是由于设计出了问题，这个“显示器”存在一个严重的 bug，就是每个像素同时只能显示红色、绿色、蓝色其中的一种颜色。msy 感觉一个学期白忙活了，但是又不能浪费材料，于是他觉得把这个有问题的“显示器”作为装饰。msy 喜欢蓝色和绿色，同时也喜欢偶数，因此他希望“显示器”的每一帧都同时包含偶数个蓝色像素和偶数个绿色像素。此时 msy 想知道，他可以看到多少不同的帧？
由于答案可能很大，你只需输出答案对 998244353998244353 取模的结果即可。

输入格式
第一行一个数 T(1\le T\le 100)T(1≤T≤100)，表示数据组数。
接下来 TT 行，每行一个整数 n(1\le n\le 10^9)n(1≤n≤10
9
)，表示“显示器”上的像素个数。

输出格式
输出 TT 行，表示每组数据的答案。

Case 1
Input
6
1
2
3
4
5
1000000000
Output
1
3
7
21
61
224965630

#include<iostream>
#include<string.h>
#include<math.h>
using namespace std;
int n=3;//矩阵维数
struct Mat {
	long long m[101][101];
	void ini() {
		memset(m, 0, sizeof(m));
	}
	void ini_dw() {
		memset(m, 0, sizeof(m));
		for (int i = 0; i < n; i++)m[i][i] = 1;
	}
};

int mod = 998244353;//模
Mat a, e;//a是输入矩阵 e是单位矩阵
Mat multi(Mat a, Mat b) {
	Mat c;
	c.ini();
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			for (int k = 0; k < n; k++) {
				c.m[i][j] = ((c.m[i][j] % mod) + (a.m[i][k] * b.m[k][j]) % mod) % mod;
			}
		}
	}
	return c;
}

Mat pow(Mat x, long long y) {
	Mat ans;
	ans.ini_dw();
	if (y == 0 )return ans;
	if (y == 1)return x;
	while (y) {
		if (y & 1)ans = multi(ans, x);
		x = multi(x, x);
		y /= 2;
	}
	return ans;
}


int main() {
	int T;
	cin >> T;
	while (T--) {
		int m;
		cin >> m;
		/*if (m == 2) {
			cout << 3 << endl;
			continue;
		}*/
		Mat t;
		t.m[0][0] = 1; t.m[0][1] = 0; t.m[0][2] = 1;
		t.m[1][0] = 0; t.m[1][1] = 1; t.m[1][2] = 1;
		t.m[2][0] = 2; t.m[2][1] = 2; t.m[2][2] = 1;
		t = pow(t, m - 1);
		long long ans = 0;
		ans = (ans%mod + t.m[0][0]%mod + (t.m[0][2] * 2)%mod) % mod;
		cout << ans << endl;
	}
	return 0;
}

在这里插入图片描述

D : 衣柜

题目描述
ZJM 有 mm 件衬衫，现一共有 nn 天。
如果 ZJM 昨天穿衬衫 AA，今天穿衬衫 BB，则他今天可以获得 H[A][B]H[A][B] 快乐值。
询问 nn 天过后，ZJM 最多可以获得多少快乐值？

输入格式
第一行两个数 n,m(2\le n\le 10^8, 1\le m\le 100)n,m(2≤n≤10
8
,1≤m≤100)，表示总天数和衬衫数量。
接下来 mm 行，每行 mm 个数，第 ii 行的第 jj 个数表示 H[i][j](0\le H[i][j]\le 10^6)H[i][j](0≤H[i][j]≤10
6
)，即昨天穿衬衫 ii，今天穿衬衫 jj 可以获得的快乐值。

输出格式
输出一个数表示答案。

Case 1
Input
3 2
0 1
1 0
Output
2

#include<iostream>
#include<string.h>
#include<math.h>
using namespace std;
int n;//矩阵维数
struct Mat {
	long long m[101][101];
	void ini() {
		memset(m, 0, sizeof(m));
	}
	void ini_dw() {
		memset(m, 0, sizeof(m));
		for (int i = 0; i < n; i++)m[i][i] = 1;
	}
};



Mat multi(Mat a, Mat b) {
	Mat c;
	c.ini();
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			for (int k = 0; k < n; k++) {
				c.m[i][j] = max(c.m[i][j],a.m[i][k] + b.m[k][j]);
			}
		}
	}
	return c;
}

Mat pow(Mat x, long long y) {
	Mat ans;
	ans.ini_dw();
	ans = x;
	while (y) {
		if (y & 1)ans = multi(ans, x);
		x = multi(x, x);
		y /= 2;
	}
	return ans;
}

int main() {
	int m;
	cin >> m >> n;
	Mat t;
	t.ini();
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			cin >> t.m[i][j];
		}
	}
	t = pow(t, m-2);
	long long ans = 0;
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			ans = max(ans, t.m[i][j]);
		}
	}
	cout << ans << endl;
}

在这里插入图片描述

思路

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
不会太难的地方，但是有和去年差不多的题目，斐波那契这题网上好多解释，自己改改就行了。
其他题可以看看扩展资料
这里是资料
里面例题能学到一些东西，但是自己做还是有差距的，挺多地方得看着别人的代码想。
主要讲了字典树和KMP还有哈希，都没啥时间去理解了，赶紧准备考试了。希望能做出来两道题。
在这里插入图片描述

正方形的源
关注
  • 4
    点赞
  • 踩
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
山东大学程序设计思维作业15
weixin_43871867的博客
06-04 259
A.ZJM与霍格沃茨 ZJM 为了准备霍格沃兹的期末考试，决心背魔咒词典，一举拿下咒语翻译题 题库格式：[魔咒] 对应功能 背完题库后，ZJM 开始刷题，现共有 N 道题，每道题给出一个字符串，可能是 [魔咒]，也可能是对应功能 ZJM 需要识别这个题目给出的是 [魔咒] 还是对应功能，并写出转换的结果，如果在魔咒词典里找不到，输出 “what?” Input 首先列出魔咒词典中不超过100000条不同的咒语，每条格式为： [魔咒] 对应功能 其中“魔咒”和“对应功能”分别为长度不超过20和80的字符串，字
Week15-Project2-Group4
03-29
加密聊天 目录 应用 连结：XXXXXXXXXXXXXXXXXXX 项目概况 加密货币在过去几年中呈指数增长，并且该应用程序允许用户聚在一起并与他人建立联系。 该应用程序填补了今天所缺少的内容的空白-与其他加密货币迷进行交流...
BZOJ 4002 有意义的字符串(矩阵快速幂，特征方程）
MasterAn的博客
05-02 389
题目链接 分析 这是一个有通项公式的数列。设an=p∗an−1+q∗an−2,an=p∗an−1+q∗an−2,a_n=p*a_{n-1}+q*a_{n-2}, an=(b+d√2)n+(b−d√2)nan=(b+d2)n+(b−d2)na_n={(\frac{b+\sqrt d}{2})}^n+{(\frac{b-\sqrt d}{2})}^n 由特征方程可知，p=b,q=d−b24p=...
矩阵快速幂技巧练习——达标的字符串
Harrison的学习博客
03-28 308
题目 给定一个数N，想象只由0和1两种字符，组成的所有长度为N的字符串。如果某个字符串，任何0字符的左边都有1紧挨着认为这个字符串达标。返回有多少达标的字符串。 package com.harrison.class15; /** * @author Harrison * @create 2022-03-28-13:18 * @motto 众里寻他千百度，蓦然回首，那人却在灯火阑珊处。 */ public class Code03_ZeroLeftOneStringNumber { publ
矩阵快速幂
P19777的博客
09-10 1044
矩阵快速幂其实并不是一个很难的算法，相对来说其实还算是比较简单，但是初看可能会难以理解他到底在干什么，其实有些我们只要动手模拟模拟就好，还有就是公式比较重要，看到有些博客拿斐波那契数列举例子，但是又不给出公式，很是让人头疼。 1.快速幂 在了解矩阵快速幂之前我们需要先了解快速幂。何为快速幂？即加快求幂运算的一种方式。 例如让我们求解2的7次方，我们之前的算法是将2重复的乘7次，也就是循...
HDU - 2604 矩阵快速幂 字符串递推 两种解法
weixin_30438813的博客
02-15 124
记dp[i]为长度i且符合题意的方案数,dp[n]就是解 符合方案的是不含fmf和fff子串的字符串 考虑如何从前面几项递推出后面第i项 (★表示存在生成的非法方案)←其实没啥用处 i=1时 m③ f③ i=2时 mm② mf② fm★② ff★② i=3时 mmm mmf mfm★ mff★ fmm ffm★ i=4时 mmmm① mmmf① mmfm★① mmff★① mfmm① mffm★①...
Week15 字符串和文件操作_1.md
09-10
Week15 字符串和文件操作_1.md
统计作业(week9-15)1
08-03
5. 矩估计 8 . 设 &hellip; 8. 经 量计算 ,
slp-week02-datalab_slp_
10-03
标题 &quot;slp-week02-datalab_slp_&quot; 暗示了这是一个关于系统级编程（System Level Programming）的学习资源，可能是一个课程或工作坊的第二周内容，重点关注数据实验室（datalab）的实践环节。&quot;bits.c&quot; 文件可能是课程...
15-16Week15.pdf
最新发布
10-17
15-16Week15.pdf
【xsy1237】 字符转换 矩阵快速幂
diaoyoutun2652的博客
02-18 140
题目大意：给你两个长度都为n，字符集为{a,b,c}的字符串S和T。 对于字符串S的任意一个字符，我们可以用cost[0]的代价，把字符a变成字符b。用cost[1]的代价，把字符b变成c，用cost[2]的代价，把字符c变成a。 问你在总代价不超过MaxCost的情况下，有多少种方法，使得字符串S与字符串T相同。 两种方法被认为是不同的，当且仅当操作序列的长度是不同的，或者某一步...
山东大学程序设计第九周作业
weixin_43814399的博客
04-26 348
第一题 咕咕东的目录管理器 咕咕东的雪梨电脑的操作系统在上个月受到宇宙射线的影响，时不时发生故障，他受不了了，想要写一个高效易用零bug的操作系统 —— 这工程量太大了，所以他定了一个小目标，从实现一个目录管理器开始。前些日子，东东的电脑终于因为过度收到宇宙射线的影响而宕机，无法写代码。他的好友TT正忙着在B站看猫片，另一位好友瑞神正忙着打守望先锋。现在只有你能帮助东东！ 初始时，咕咕东的硬盘是空...
山东大学计算机科学与技术学院程序设计思维与实践模测
qq_46027075的博客
09-13 2380
20210319 A : R!!G!!B!! 问题描述 msy 的显示器被 yhf 借走了，于是 msy 需要一个新的显示器。他买来了许多 LED 小灯，每个小灯只能发出红、绿、蓝三种颜色光的其中一种。msy 需要三个不同颜色的小灯来拼成一个像素（像素之间并不能共用小灯），但是他并不知道每种颜色的小灯具体有多少个，只知道每个小灯的颜色。msy 想知道他用手头上的小灯可以拼出多少个像素，但是一个个数太麻烦了，他希望你来用程序解决这个问题。 输入格式 第一行输入一个数 nn，表示所有小灯的数量。 第二
山东大学计算机科学与技术学院程序设计思维与实践作业
热门推荐
qq_46027075的博客
09-13 1万+
作业 H1 A : 打酱油 问题描述 小明带着N元钱去买酱油。酱油10块钱一瓶，商家进行促销，每买3瓶送1瓶，或者每买5瓶送2瓶。请问小明最多可以得到多少瓶酱油。 输入格式 输入的第一行包含一个整数N，表示小明可用于买酱油的钱数。N是10的整数倍，N不超过300。 输出格式 输出一个整数，表示小明最多可以得到多少瓶酱油。 样例1 输入 40 输出 5 样例说明 把40元分成30元和10元，分别买3瓶和1瓶，其中3瓶送1瓶，共得到5瓶。 样例2 输入 80 输出
SDU 程序设计思维与实践 Week15 思维题
weixin_43975891的博客
05-29 205
CodeForces - 1165C 新式字符串定义如下： 长度为偶数（可以为空） 每一个奇数位置的字符都不等于它后面的那个字符 例如，字符串 good、string、xyyx 都符合要求，而字符串 bad、aa、aabc 则不符合。 现在你获得了一个字符串 S，你需要从该字符串中删除最少的字符数，使其变成一个新式字符串 代码 #include&lt;iostream&gt; #include&lt;string&gt; #include&lt;algorithm&gt; using namespace s
Vj程序设计作业H15
weixin_53437436的博客
07-03 392
斐波那契数列的定义如下： 给出 n,p ，求出 f(n)%p 的值。第一行一个数 T(1≤T≤100)，表示数据组数。 接下来 T 行，每行两个整数 n,p(1≤n≤109,108≤p≤109)。输出 T 行，表示每组数据的答案。测试样例样例 1输入: 输出: 解答 给定 𝑛 和 𝑘，计算 Σi=1n​ik 对 109+7 取模的结果。第一行一个数 T(1≤T≤100)，表示数据组数。 接下来 T 行，每行两个整数 n,k(1≤n≤109,1≤k≤10)。输出 T 行，表示每组数据的答..
程序设计思维与实践 Week15 实验
qq_44157558的博客
06-05 465
目录A - Q 老师的记录册输入输出输入样例1输出样例1输入样例2输出样例2输入样例3输出样例3思路代码B - ZJM的本领输入输出样例输入1样例输出1样例输入2样例输出2思路代码C - TT 的神秘任务 - XInputOutput样例输入1样例输出1样例输入2样例输出2思路代码D - 瑞瑞爱上字符串InputOutputExample InputExample Output思路代码E - 苟狗之宇宙射线再打击InputOutputExample Input1Example Output1Example
斐波那契数列
Leung2046462的专栏
05-11 314
题目描述： 大家都知道斐波那契数列，现在要求输入一个整数n，请你输出斐波那契数列的第n项。斐波那契数列的定义如下： 输入： 输入可能包含多个测试样例，对于每个测试案例， 输入包括一个整数n(1<=n<=70)。 输出： 对应每个测试案例， 输出第n项斐波那契数列的值。 样例输入： 3 样例输出： 2代码：
程序设计思维与实践 Week4 作业 (3/4/数据班)
明天我能有猫吗的博客
03-20 245
A - DDL 的恐惧 题目 ZJM 有 n 个作业，每个作业都有自己的 DDL，如果 ZJM 没有在 DDL 前做完这个作业，那么老师会扣掉这个作业的全部平时分。 所以 ZJM 想知道如何安排做作业的顺序，才能尽可能少扣一点分。 请你帮帮他吧！ Input 输入包含T个测试用例。输入的第一行是单个整数T，为测试用例的数量。 每个测试用例以一个正整数N开头(1&lt;=N&lt;=1000)，表示...
深度学习week1-RNN心脏病预测
07-27
你好！对于心脏病预测的问题，使用循环神经网络（RNN）是一种常见的方法。RNN适用于处理序列数据，而心电图信号就是一种序列数据。在使用RNN进行心脏病预测时，你可以将心电图信号作为输入序列，然后通过训练RNN模型来预测患者是否患有心脏病。 首先，你需要准备一个合适的数据集，其中包含心电图信号和相应的心脏病标签。可以使用公开的心电图数据集，如PTB数据库或MIT-BIH数据库。然后，你可以对数据进行预处理和特征工程，如数据清洗、滤波、降采样等。 接下来，你可以构建一个RNN模型。RNN模型由一系列循环层组成，每个循环层都会处理一个时间步的输入数据。你可以选择不同类型的RNN单元，如简单循环单元（SimpleRNN）、长短期记忆网络（LSTM）或门控循环单元（GRU）。通过添加适当的全连接层和激活函数，你可以将RNN模型输出映射到二分类问题（有或无心脏病）的结果。 然后，你可以使用训练集对RNN模型进行训练，并使用验证集进行模型调优。在训练过程中，你可以使用适当的损失函数（如交叉熵）和优化算法（如随机梯度下降）来最小化模型的预测误差。 最后，你可以使用测试集对训练好的模型进行评估，并计算模型的性能指标，如准确率、精确率、召回率等。这些指标可以帮助你评估模型的预测能力和泛化能力。 需要注意的是，心脏病预测是一个复杂的医学问题，仅仅使用心电图信号可能不能得到准确的预测结果。通常情况下，还需要结合其他患者的临床信息和医学知识来进行综合评估。因此，在进行心脏病预测时，建议与专业医生合作，并遵循相关的医学准则和规范。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值