week14——作业(矩阵快速幂练习)

Q老师与石头剪刀布 :

问题描述

题目简述

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

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

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

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

输入/输出格式

输入格式:
第一行一个整数 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’ 表示剪刀

输出格式:
对于每组数据:

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

“YES”/"NO"是大小写不敏感的,但是 ‘R’, ‘P’, ‘S’ 是大小写敏感的。

样例

输入样例:
2
3
1 1 1
RPS
3
3 0 0
RPS
输出样例:
YES
PSR
NO

问题分析

解题思路

问Q老师怎么赢,那么能赢的时候就尽可能让Q老师赢就好了。首先根据对手的出拳规律,如果发现现在Q老师还可以出能赢这局的方案,那么就填入该方案,并使可以出此方案的次数减1,并统计上该局为Q老师赢的局。否则,暂时跳过该局。再遍历完毕后如果发现Q老师无法赢得一半及以上的局,那么,Q老师无法取胜,因为这是他最多可能赢得的对局数,如果可以,那么,那些跳过的局就以剩下的可出方案填充即可。

参考代码
#include <iostream>
#include <cstring>
#include <string>
#include <cmath>

using namespace std;

int t;
string cc;
int tt;
int a,b,c;

void solve()
{
	char ans[110];
	for(int i=0;i<110;i++)
	    ans[i]='n';
	int count=0;
	for(int i=0;i<cc.length();i++)
	{
		if(cc[i]=='R')
		{
			if(b!=0) 
			{
				ans[i]='P';
				b--;
				count++;
			} 
		}
		else if(cc[i]=='P')
		{
			if(c!=0) 
			{
				ans[i]='S';
				c--;
				count++;
			} 
		}
		else 
		{
			if(a!=0) 
			{
				ans[i]='R';
				a--;
				count++;
			}
		}
	}
	if(count>=ceil(tt/2.0))
	{
		printf("YES\n");
		for(int i=0;i<tt;i++)
		{
			if(ans[i]!='n') printf("%c",ans[i]);
			else
			{
				if(a!=0)
				{
				    printf("R");
					a--;	
				}
				else if(b!=0)
				{
					printf("P");
					b--;	
				}
				else
				{
					printf("S");
					c--;	
				}
			}
		}
		printf("\n");
	}
	else printf("NO\n");
}

int main()
{
	scanf("%d",&t);
	for(int i=1;i<=t;i++)
	{
		scanf("%d",&tt);
		scanf("%d %d %d",&a,&b,&c); 
	    cin>>cc;
		solve();
	}
	return 0;
}

心得体会

不难,基本是一个贪心的思想。

Q老师与十字叉:

问题描述

题目简述

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老师 每分钟可以选择一个白色的格子并且把它涂黑。现在他想知道要完成这个工作,最少需要几分钟?

输入/输出格式

输入格式:
第一行包含一个整数 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
输出格式:
答案输出 q 行, 第 i 行包含一个整数 — 表示第 i 组数据的答案

样例

输入样例:
1
5 5
00100
00100
11111
00100
00100
输出样例:
0

问题分析

解题思路

十字叉是一行一列均为黑色,那么我们就统计每一行和每一列的白格数。当我们在涂色的时候,涂色的数量为对应的行和列的白格数之和,如果发现行列交叉处的是白格,那么再减去1即可,最终取这些结果的最小值即为所求。

参考代码
#include <iostream>
#include <cstring>
#include <string>
#include <cmath>
#include <algorithm>

using namespace std;

char** array;
int n,m;
int q;

void solve()
{
	int lines[50005];
	int columns[50005];
	memset(lines,0,sizeof(lines));
	memset(columns,0,sizeof(columns));
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			if(array[i][j]=='.')
			{
				lines[i]++;
				columns[j]++;
			}
		}
	}
	int ans=1e8;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			if(array[i][j]=='.')
			{
				ans=min(ans,lines[i]+columns[j]-1);
			}
			else ans=min(ans,lines[i]+columns[j]);
		}
	}
	printf("%d\n",ans);
}

int main()
{
	scanf("%d",&q);
	for(int i=1;i<=q;i++)
	{
		scanf("%d %d",&n,&m);
		array=new char*[n+1];
		for(int j=0;j<=n;j++)
		{
			array[j]=new char[m+1];
		}
		for(int j=1;j<=n;j++)
		{
			getchar();
			//printf("j=%d\n",j);
			for(int k=1;k<=m;k++)
			{
				//printf("k=%d\n",k);
				scanf("%c",&array[j][k]);
			}
		}
		solve();
	}
	return 0;
}

心得体会

也是挺简单的一道题,这里注意就是不要忘记行列交叉处的格子会统计两次,需要减去一次才能得到正确答案,其他的应该很简单。

Q老师的考验:

问题描述

题目简述

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老师 的考验吗?

输入/输出格式

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

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

第二行给定十个整数,分别表示 a0~a9。.
输出格式:
对于每一组测试用例输出一行,表示 f(k) % m 的数值大小。

样例

输入样例:
10 9999
1 1 1 1 1 1 1 1 1 1
20 500
1 0 1 0 1 0 1 0 1 0
输出样例:
45
104

问题分析

解题思路

矩阵快速幂的模板题,题干上已经指明了线性递推式的形式了,这样直接写出转移方程即可。很明显可以发现,这应该是一个10维的向量,转移矩阵是一个10*10的方阵,根据线性递推式写出矩阵中每行和每列的值即可。之后就是矩阵快速幂求出转移矩阵的n-9次方并乘初始向量,即可得到结果向量,取结果向量中的第一维后模m即可。

参考代码
#include <iostream>
#include <cstring>
#include <string>
#include <cmath>

using namespace std;

const int N=10;
int k,m;

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]%m;
					ret.x[i][j]%=m; 
				}
			}
		}
		return ret;
	}
	matrix()
	{
		memset(x,0,sizeof(x));
	}
	matrix(const matrix& t)
	{
		memcpy(x,t.x,sizeof(x));
	}
};

int main()
{
	while(scanf("%d %d",&k,&m)!=EOF)
	{
		if(k<10) 
		{
			printf("%d\n",k);
			continue;
		}
		matrix sp;
		matrix init;
		matrix e;
		for(int i=0;i<10;i++)
		{
			scanf("%d",&sp.x[0][i]);
			if(i+1<10) sp.x[i+1][i]=1;
			init.x[i][0]=10-i-1;
			e.x[i][i]=1;
		} 
		int times=k-9;
		while(times)
		{
			if(times&1) e=e*sp;
			sp=sp*sp;
			times>>=1;
		}
		init=e*init;
		printf("%d\n",init.x[0][0]%m);
	}
	return 0;
}

心得体会

非常直白的一个问题,主要是练习封装矩阵和矩阵快速幂的模板的写法,总体来说比较顺利,也不是很难。

Q老师染砖:

问题描述

题目简述

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

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

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

输入/输出格式

输入格式:
第一行为 T,代表数据组数。(1 ≤ T ≤ 100)

接下来 T 行每行包括一个数字 N,代表有 N 块砖。(1 ≤ N ≤ 1e9)
输出格式:
输出满足条件的方案数,答案模 10007。

样例

输入样例:
2
1
2
输出样例:
2
6

问题分析

解题思路

矩阵快速幂优化dp,这里定义Ai为涂前i块砖后,红色和绿色的砖块数均为偶数的方案数,Bi为涂前i块砖后,红色和绿色的砖块数有且仅有一种为偶数的方案数,Ci为红色和绿色的砖块数均为奇数的方案数。根据以上状态,可以列出三个状态转移方程:
Ai=2Ai-1+Ci-1
Bi=2
Bi-1+Ci-1
Ci=2Ai-1+2Bi-1+2*Ci-1
这三个方程均为线性递推式,因此可以使用矩阵快速幂优化dp求解,即以[Ai,Bi,Ci]为向量,很容易可以写出转移矩阵,之后求解即可。

参考代码
#include <iostream>
#include <cstring>
#include <string>
#include <cmath>

using namespace std;

const int N=3;
int p;
const int m=10007;

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]%m;
					ret.x[i][j]%=m; 
				}
			}
		}
		return ret;
	}
	matrix()
	{
		memset(x,0,sizeof(x));
	}
	matrix(const matrix& t)
	{
		memcpy(x,t.x,sizeof(x));
	}
};

int main()
{
	int tt;
	scanf("%d",&tt);
	for(int m=1;m<=tt;m++)
	{
		scanf("%d",&p);
		if(p==1) 
		{
			printf("%d\n",2);
			continue;
		}
	    matrix sp;
        matrix init;
        matrix e; 
    	sp.x[0][0]=2;
    	sp.x[0][1]=0;
	    sp.x[0][2]=1;
	    sp.x[1][0]=0;
    	sp.x[1][1]=2;
	    sp.x[1][2]=1;
    	sp.x[2][0]=2;
    	sp.x[2][1]=2;
    	sp.x[2][2]=2;
    	for(int i=0;i<3;i++)
	     {
    		e.x[i][i]=1;
    	}
        init.x[0][0]=2;
	    init.x[1][0]=0;
    	init.x[2][0]=2;
		int times=p-1;
		while(times)
		{
			if(times&1) e=e*sp;
			sp=sp*sp;
			times>>=1;
		}
		init=e*init;
		printf("%d\n",init.x[0][0]%10007);
	}
	return 0;
}

心得体会

dp问题还是列出状态非常关键,后面的做法其实非常多。这个题是我第一次接触矩阵快速幂优化dp,做的时候感觉比较生疏,卡的地方有点多,但总体思路感觉还是非常清晰的,也基本理解了算法。

Q老师度假:

问题描述

题目简述

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

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

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

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

输入/输出格式

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

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

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

测试样例组数不会超过 10。
输出格式:
每组测试样例输出一行,表示 Q老师 可以获得的最大快乐值。

样例

输入样例:
3 2
0 1
1 0
4 3
1 2 3
1 2 3
1 2 3
输出样例:
2
9

问题分析

解题思路

矩阵快速幂优化dp的又一种形式,首先根据题目要求,写出状态是非常容易的。设f[i][j]为第i天穿第j件衣服可以获得的最大快乐值,状态转移方程为:f[i][j]=max(f[i-1][k]+h[k][j]),最终结果为max(f[n][k])。但是由于数据范围过大,需要使用矩阵快速幂进行优化。这里使用max和+运算代替原本矩阵乘法中的+和*运算。可以发现,代替后的这两个运算均满足结合律,可以使用快速幂优化。之后的做法和矩阵快速幂相同。

参考代码
#include <iostream>
#include <cstring>
#include <string>
#include <cmath>

using namespace std;

const int N=100;
int m,n;

struct matrix
{
	long long 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]=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));
	}
};

int main()
{
	while(scanf("%d %d",&n,&m)!=EOF)
	{
	    matrix sp;
        matrix init;
        matrix e;
        for(int i=0;i<m;i++)
        {
        	for(int j=0;j<m;j++)
        	{
        		scanf("%lld",&sp.x[i][j]);
			}
		}
		int times=n-1;
		while(times)
		{
			if(times&1) e=e*sp;
			sp=sp*sp;
			times>>=1;
		}
		init=e*init;
		printf("%lld\n",init.x[0][0]);
	}
	return 0;
}

心得体会

自我感觉这个题挺难独立想到做法的,使用矩阵快速幂的优化过程也非常的巧妙,个人感觉是挺难做的一道题。关于矩阵快速幂优化dp还需要进一步的进行落实。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值