week14作业矩阵快速幂

week14


A
题目:
每一个大人曾经都是一个小孩,Q老师 也一样。

为了回忆童年,Q老师 和 Monika 玩起了石头剪刀布的游戏,游戏一共 n 轮。无所不知的 Q老师 知道每一轮 Monika 的出招,然而作为限制, Q老师 在这 n 轮游戏中必须恰好出 a 次石头,b 次布和 c 次剪刀。

如果 Q老师 赢了 Monika n/2(上取整) 次,那么 Q老师就赢得了这场游戏,否则 Q老师 就输啦!

Q老师非常想赢,他想知道能否可以赢得这场游戏,如果可以的话,Q老师希望你能告诉他一种可以赢的出招顺序,任意一种都可以。

Input
第一行一个整数 t(1 ≤ t ≤ 100)表示测试数据组数。然后接下来的 t 组数据,每一组都有三个整数:

第一行一个整数 n(1 ≤ n ≤ 100)
第二行包含三个整数 a, b, c(0 ≤ a, b, c ≤ n)。保证 a+b+c=n
第三行包含一个长度为 n 的字符串 s,字符串 s 由且仅由 ‘R’, ‘P’, ‘S’ 这三个字母组成。第 i 个字母 s[i] 表示 Monika 在第 i 轮的出招。字母 ‘R’ 表示石头,字母 ‘P’ 表示布,字母 ‘S’ 表示剪刀
Output
对于每组数据:

如果 Q老师 不能赢,则在第一行输出 “NO”(不含引号)
否则在第一行输出 “YES”(不含引号),在第二行输出 Q老师 的出招序列 t。要求 t 的长度为 n 且仅由 ‘R’, ‘P’, ‘S’ 这三个字母构成。t 中需要正好包含 a 个 ‘R’,b 个 ‘P’ 和 c 个 ‘S’
“YES”/"NO"是大小写不敏感的,但是 ‘R’, ‘P’, ‘S’ 是大小写敏感的。

Example
Input
2
3
1 1 1
RPS
3
3 0 0
RPS
Output
YES
PSR
NO

思路:
此题重点在于将多出来的补位,首先将R用完,其次是P,最后是S,则定义一个函数,用于补位,补R,S,P的顺序,最后按顺序输出即可。

代码:

#include<bits/stdc++.h>
using namespace std;
int t,n,a,b,c,a1,b1,c1,a2,b2,c2,q,w,e;
string s; 
void d()
{
	if(a2>0)
	cout << "R",a2--;
	else if(b2>0)
	cout << "P",b2--;
	else
	cout << "S",c2--;
}
int main()
{
	cin >>t;
	while(t--)
	{
		a1=b1=c1=q=w=e=0;
		cin >> n >> a >>b >> c;
		cin >> s;
		for(int i=0;i<n;i++)
		{
			if(s[i]=='R')
			b1++;
			if(s[i]=='P')
			c1++;
			if(s[i]=='S')
			a1++;
		}
		a2=a-a1;
		a=min(a,a1);
		b2=b-b1;
		b=min(b,b1);
		c2=c-c1;
		c=min(c,c1);
		//cout << a1 <<b1 << c1 <<endl;
		//cout << a2 << b2 << c2 <<endl;		
		if(a+b+c>((n-1)/2))
		{
			cout << "YES" << endl;
			for(int i=0;i<n;i++)
			{
				if(s[i]=='R')
				{
					if(q<b1+b2)
					{
						cout << "P";
						q++;
					}
					else
					{
						d();
					}
				}
				if(s[i]=='P')
				{
					if(w<c1+c2)
					{
						cout << "S";
						w++;
					}
					else
					d();
				}
				if(s[i]=='S')
				{
					if(e<a1+a2)
					{
						cout << "R";
						e++;
					}
					else
					d();
				}
			}
			cout << endl;
		}
		else{
			cout << "NO" << endl;
		}
	}
}

B
题目:
Q老师 得到一张 n 行 m 列的网格图,上面每一个格子要么是白色的要么是黑色的。

Q老师认为失去了 十字叉 的网格图莫得灵魂. 一个十字叉可以用一个数对 x 和 y 来表示, 其中 1 ≤ x ≤ n 并且 1 ≤ y ≤ m, 满足在第 x 行中的所有格子以及在第 y 列的 所有格子都是黑色的

例如下面这5个网格图里都包含十字叉

第四个图有四个十字叉,分别在 (1, 3), (1, 5), (3, 3) 和 (3, 5).

下面的图里没有十字叉

Q老师 得到了一桶黑颜料,他想为这个网格图注入灵魂。 Q老师 每分钟可以选择一个白色的格子并且把它涂黑。现在他想知道要完成这个工作,最少需要几分钟?

Input
第一行包含一个整数 q (1 ≤ q ≤ 5 * 10^4) — 表示测试组数

对于每组数据:
第一行有两个整数 n 和 m (1 ≤ n, m ≤ 5 * 10^4, n * m ≤ 4 * 10^5) — 表示网格图的行数和列数

接下来的 n 行中每一行包含 m 个字符 — ‘.’ 表示这个格子是白色的, ‘*’ 表示这个格子是黑色的

保证 q 组数据中 n 的总和不超过 5 * 10^4, n*m 的总和不超过 4 * 10^5

Output
答案输出 q 行, 第 i 行包含一个整数 — 表示第 i 组数据的答案

Example
Input
9
5 5




3 4


.
.

4 3


*…
*…
*…
5 5


..*


.
…***
1 4


5 5


.
.


5 3

.
.
.*.


..
3 3
.
.
.
.*.
4 4
.*

.*
.*
Output
0
0
0
0
0
4
1
1
2

思路:
此题需要找出一个白格作为答案的一行一列的相交点,这个时候最终答案就要相应减一。先遍历每行每列,找出需要涂黑的最小的数量的行列值,然后判断这些行列是否相较于白格,即可。

代码:

#include<bits/stdc++.h>
using namespace std;
int t,n,m,a,b,c;
bool visn[50002];
bool vism[50002];
string s[50002];
void clear()
{
	//memset(s,0,sizeof(s));
	memset(visn,0,sizeof(visn));
	memset(vism,0,sizeof(vism));
}
int d()
{
		for(int i=0;i<n;i++)
		{
			if(visn[i])
			for(int j=0;j<m;j++)
			{
				if(vism[j]&&s[i][j]=='.')
				return -1; 
			} 
		}
		return 0;
}
int main()
{
	cin >>t;
	while(t--)
	{
		clear();
		b=0;
		c=0;
		scanf("%d%d",&n,&m);
		for(int i=0;i<n;i++)
		{
			a=0;
			cin >> s[i];
			for(int j=0;j<m;j++)
			{
				if(s[i][j]=='*')
				a++;
			}
			b=max(a,b);
		}
		for(int i=0;i<n;i++)
		{
			a=0;
			for(int j=0;j<m;j++)
			{
				if(s[i][j]=='*')
				a++;
			}
			if(a==b)
			{
				visn[i]=1;
			}
		}
		for(int i=0;i<m;i++)
		{
			a=0;
			for(int j=0;j<n;j++)
			{
				if(s[j][i]=='*')
				a++;
			}
			c=max(a,c);
		}
		for(int i=0;i<m;i++)
		{
			a=0;
			for(int j=0;j<n;j++)
			{
				if(s[j][i]=='*')
				a++;
			}
			if(a==c)
			vism[i]=1;
		}
		cout << m+n-b-c+d() << endl;
		
	}
}

C
题目:
Q老师 对数列有一种非同一般的热爱,尤其是优美的斐波那契数列。

这一天,Q老师 为了增强大家对于斐波那契数列的理解,决定在斐波那契的基础上创建一个新的数列 f(x) 来考一考大家。数列 f(x) 定义如下:

当 x < 10 时,f(x) = x;
当 x ≥ 10 时,f(x) = a0 * f(x-1) + a1 * f(x-2) + a2 * f(x-3) + …… + a9 * f(x-10),ai 只能为 0 或 1。

Q老师 将给定 a0~a9,以及两个正整数 k m,询问 f(k) % m 的数值大小。

聪明的你能通过 Q老师 的考验吗?

Input
输出文件包含多组测试用例,每组测试用例格式如下:

第一行给定两个正整数 k m。(k < 2e9, m < 1e5)

第二行给定十个整数,分别表示 a0~a9。

Output
对于每一组测试用例输出一行,表示 f(k) % m 的数值大小。

Sample Input
10 9999
1 1 1 1 1 1 1 1 1 1
20 500
1 0 1 0 1 0 1 0 1 0
Sample Output
45
104

思路:
面对斐波拉契数列类似的表达式,想到用矩阵快速幂,提高效率,表达式需要多少个组成,矩阵的宽度即为多大,快速幂的形式可以参考快速幂的方法,需要在matrix中重载乘法函数。最后判断需要有多少次幂,即可的最终结果。

代码:

#include<bits/stdc++.h>
using namespace std;
int k,p;
int a[10];
long long ans=0;
const int N=10;
struct matrix{
	int x[N][N];
	matrix operator*(const matrix& t) const{
		matrix ret;
		for(int i=0;i<N;i++)
		{
			for(int j=0;j<N;j++)
			{
				ret.x[i][j]=0;
				for(int k=0;k<N;k++)
				{
					ret.x[i][j]+=x[i][k]*t.x[k][j]%p;
					ret.x[i][j]%=p;
				}
			}
		}
		return ret;
	}
	matrix() {memset(x,0,sizeof(x));}
	matrix(const matrix& t) {memcpy(x,t.x,sizeof(x));}
};
void d(matrix &ret)
{
	for(int i=0;i<10;i++)
	{
		for(int j=0;j<10;j++)
		{
			if(i==0)
			ret.x[0][j]=a[j];
			else
			{
				if(j==i-1)
				ret.x[i][j]=1;
				else
				ret.x[i][j]=0;
			}
		}
	}
}
matrix quick_pow(matrix a,int x){
	matrix ret;
	for(int i=0;i<10;i++)
	{
		for(int j=0;j<10;j++)
		{
			if(i==j)
			ret.x[i][j]=1;
			else
			ret.x[i][j]=0;
		}
	}
	while(x){
		if(x&1) ret=ret*a;
		a=a*a;
		x>>=1;
	}
	return ret;
}
int main()
{
	while(scanf("%d%d",&k,&p)!=EOF)
	{
		ans=0;
		for(int i=0;i<10;i++)
		cin >> a[i];
		if(k<10)
		{
			cout << k%p << endl;
			continue;
		}
		matrix a,ret;
		d(a);
		ret=quick_pow(a,k-9);
		for(int i=0;i<10;i++)
		{
			ans+=ret.x[0][i]*(9-i);
		}
		cout << ans%p << endl;
	}
}

D
题目:
衣食无忧的 Q老师 有一天突发奇想,想要去感受一下劳动人民的艰苦生活。

具体工作是这样的,有 N 块砖排成一排染色,每一块砖需要涂上红、蓝、绿、黄这 4 种颜色中的其中 1 种。且当这 N 块砖中红色和绿色的块数均为偶数时,染色效果最佳。

为了使工作效率更高,Q老师 想要知道一共有多少种方案可以使染色效果最佳,你能帮帮他吗?

Input
第一行为 T,代表数据组数。(1 ≤ T ≤ 100)

接下来 T 行每行包括一个数字 N,代表有 N 块砖。(1 ≤ N ≤ 1e9)

Output
输出满足条件的方案数,答案模 10007。

Sample Input
2
1
2
Sample Output
2
6

思路:
直接套用矩阵快速幂即可,其中dp转移方程为:
A[i] = 2 * A[i-1] + C[i-1]
B[i] = 2 * B[i-1] + C[i-1]
C[i] = 2 * A[i-1] + 2 * B[i-1] + 2 * C[i-1]
然后构造矩阵进行快速幂运算。

代码:

#include<iostream>
#include<string.h>
using namespace std;
long long ans=0;
int t,n,p=10007;
const int N=3;
struct matrix{
	int x[N][N];
	matrix operator*(const matrix& t) const{
		matrix ret;
		for(int i=0;i<N;i++)
		{
			for(int j=0;j<N;j++)
			{
				ret.x[i][j]=0;
				for(int k=0;k<N;k++)
				{
					ret.x[i][j]+=x[i][k]*t.x[k][j]%p;
					ret.x[i][j]%=p;
				}
			}
		}
		return ret;
	}
	matrix() {memset(x,0,sizeof(x));}
	matrix(const matrix& t) {memcpy(x,t.x,sizeof(x));}
	//matrix(const matrix& t) {memcpy(x,t.x,sizeof(x));}
};
void init(matrix &ret)
{
	ret.x[0][0]=2;
	ret.x[0][1]=0;
	ret.x[0][2]=1;
	ret.x[1][0]=0;
	ret.x[1][1]=2;
	ret.x[1][2]=1;
	ret.x[2][0]=ret.x[2][1]=ret.x[2][2]=2;
}
matrix quick_pow(matrix a,int x)
{
	matrix ret;
	for(int i=0;i<N;i++)
	{
		for(int j=0;j<N;j++)
		{
			if(i==j)
			ret.x[i][j]=1;
			else
			ret.x[i][j]=0;
		}
	}
	while(x){
		if(x&1) ret=ret*a;
		a=a*a;
		x>>=1;
	}
	return ret;
}
int main()
{
	cin >> t;
	while(t--)
	{
		cin >> n;
		matrix a,ret;
		init(a);
		ret=quick_pow(a,n-1);
		ans=ret.x[0][0]*2+ret.x[0][2]*2;
		cout << ans%p <<endl;
	}
 } 

E
题目:
忙碌了一个学期的 Q老师 决定奖励自己 N 天假期。

假期中不同的穿衣方式会有不同的快乐值。

已知 Q老师 一共有 M 件衬衫,且如果昨天穿的是衬衫 A,今天穿的是衬衫 B,则 Q老师 今天可以获得 f[A][B] 快乐值。

在 N 天假期结束后,Q老师 最多可以获得多少快乐值?

Input
输入文件包含多组测试样例,每组测试样例格式描述如下:

第一行给出两个整数 N M,分别代表假期长度与 Q老师 的衬衫总数。(2 ≤ N ≤ 100000, 1 ≤ M ≤ 100)

接下来 M 行,每行给出 M 个整数,其中第 i 行的第 j 个整数,表示 f[i][j]。(1 ≤ f[i][j] ≤ 1000000)

测试样例组数不会超过 10。

Output
每组测试样例输出一行,表示 Q老师 可以获得的最大快乐值。

Sample Input
3 2
0 1
1 0
4 3
1 2 3
1 2 3
1 2 3
Sample Output
2
9

思路:、
这题用普通的方法会超时,需要用矩阵快速幂进行优化,但由于不满足快速幂的条件,所以需要将矩阵乘法进行修改,改为ret.x[i][j]=max(x[i][k]+t.x[k][j],ret.x[i][j]),而且转换过程
•C[i][j] = ∑(A[i][k] ∗ B[k][j]) ⇒ max (A[i][k]+ B[k][j])
• 将累加换成了max,乘法换成了加法。
最后即可快速求出。

代码:

#include<bits/stdc++.h>
using namespace std;
long long ans=0;
int n,m;
struct matrix{
	long long x[110][110];
	matrix operator*(const matrix& t) const{
		matrix ret;
		for(int i=0;i<n;i++)
		{
			for(int j=0;j<n;j++)
			{
				ret.x[i][j]=0;
				for(int k=0;k<n;k++)
				{
					ret.x[i][j]=max(x[i][k]+t.x[k][j],ret.x[i][j]);
				}
			}
		}
		return ret;
	}
	matrix() {memset(x,0,sizeof(x));}
	matrix(const matrix& t) {memcpy(x,t.x,sizeof(x));}
};
void init(matrix &ret)
{
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<n;j++)
		{
			cin >> ret.x[i][j];
		}
		
	}
}
matrix quick_pow(matrix a,int x)
{
	matrix ret;
	while(x){
		if(x&1) ret=ret*a;
		a=a*a;
		x>>=1;
	}
	return ret;
}
int main()
{
	while(scanf("%d%d",&m,&n)!=EOF)
	{
		ans=0;
		matrix a,ret;
		init(a);
		ret=quick_pow(a,m-1);
		for(int i=0;i<n;i++)
		{
			for(int j=0;j<n;j++)
			{
				ans=max(ans,ret.x[i][j]);
			}
		}
		cout << ans <<endl;
	}
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
完整版:https://download.csdn.net/download/qq_27595745/89522468 【课程大纲】 1-1 什么是java 1-2 认识java语言 1-3 java平台的体系结构 1-4 java SE环境安装和配置 2-1 java程序简介 2-2 计算机中的程序 2-3 java程序 2-4 java类库组织结构和文档 2-5 java虚拟机简介 2-6 java的垃圾回收器 2-7 java上机练习 3-1 java语言基础入门 3-2 数据的分类 3-3 标识符、关键字和常量 3-4 运算符 3-5 表达式 3-6 顺序结构和选择结构 3-7 循环语句 3-8 跳转语句 3-9 MyEclipse工具介绍 3-10 java基础知识章节练习 4-1 一维数组 4-2 数组应用 4-3 多维数组 4-4 排序算法 4-5 增强for循环 4-6 数组和排序算法章节练习 5-0 抽象和封装 5-1 面向过程的设计思想 5-2 面向对象的设计思想 5-3 抽象 5-4 封装 5-5 属性 5-6 方法的定义 5-7 this关键字 5-8 javaBean 5-9 包 package 5-10 抽象和封装章节练习 6-0 继承和多态 6-1 继承 6-2 object类 6-3 多态 6-4 访问修饰符 6-5 static修饰符 6-6 final修饰符 6-7 abstract修饰符 6-8 接口 6-9 继承和多态 章节练习 7-1 面向对象的分析与设计简介 7-2 对象模型建立 7-3 类之间的关系 7-4 软件的可维护与复用设计原则 7-5 面向对象的设计与分析 章节练习 8-1 内部类与包装器 8-2 对象包装器 8-3 装箱和拆箱 8-4 练习题 9-1 常用类介绍 9-2 StringBuffer和String Builder类 9-3 Rintime类的使用 9-4 日期类简介 9-5 java程序国际化的实现 9-6 Random类和Math类 9-7 枚举 9-8 练习题 10-1 java异常处理 10-2 认识异常 10-3 使用try和catch捕获异常 10-4 使用throw和throws引发异常 10-5 finally关键字 10-6 getMessage和printStackTrace方法 10-7 异常分类 10-8 自定义异常类 10-9 练习题 11-1 Java集合框架和泛型机制 11-2 Collection接口 11-3 Set接口实现类 11-4 List接口实现类 11-5 Map接口 11-6 Collections类 11-7 泛型概述 11-8 练习题 12-1 多线程 12-2 线程的生命周期 12-3 线程的调度和优先级 12-4 线程的同步 12-5 集合类的同步问题 12-6 用Timer类调度任务 12-7 练习题 13-1 Java IO 13-2 Java IO原理 13-3 流类的结构 13-4 文件流 13-5 缓冲流 13-6 转换流 13-7 数据流 13-8 打印流 13-9 对象流 13-10 随机存取文件流 13-11 zip文件流 13-12 练习题 14-1 图形用户界面设计 14-2 事件处理机制 14-3 AWT常用组件 14-4 swing简介 14-5 可视化开发swing组件 14-6 声音的播放和处理 14-7 2D图形的绘制 14-8 练习题 15-1 反射 15-2 使用Java反射机制 15-3 反射与动态代理 15-4 练习题 16-1 Java标注 16-2 JDK内置的基本标注类型 16-3 自定义标注类型 16-4 对标注进行标注 16-5 利用反射获取标注信息 16-6 练习题 17-1 顶目实战1-单机版五子棋游戏 17-2 总体设计 17-3 代码实现 17-4 程序的运行与发布 17-5 手动生成可执行JAR文件 17-6 练习题 18-1 Java数据库编程 18-2 JDBC类和接口 18-3 JDBC操作SQL 18-4 JDBC基本示例 18-5 JDBC应用示例 18-6 练习题 19-1 。。。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值