week14作业矩阵快速幂

活下去

于 2020-06-11 17:03:23 发布

阅读量157

点赞数

本文链接：https://blog.csdn.net/hexiangqin/article/details/106692199

版权

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
评论
week14作业矩阵快速幂

week14ABCA题目：每一个大人曾经都是一个小孩，Q老师也一样。为了回忆童年，Q老师和 Monika 玩起了石头剪刀布的游戏，游戏一共 n 轮。无所不知的 Q老师知道每一轮 Monika 的出招，然而作为限制， Q老师在这 n 轮游戏中必须恰好出 a 次石头，b 次布和 c 次剪刀。如果 Q老师赢了 Monika n/2(上取整) 次，那么 Q老师就赢得了这场游戏，否则 Q老师就输啦！Q老师非常想赢，他想知道能否可以赢得这场游戏，如果可以的话，Q老师希望你能告诉他一种可以赢的出招
复制链接

扫一扫