Week14 矩阵快速幂

最新推荐文章于 2021-05-20 00:03:01 发布

e青青青

最新推荐文章于 2021-05-20 00:03:01 发布

阅读量201

点赞数

本文链接：https://blog.csdn.net/Hdjjdsbjj/article/details/106373371

版权

Q老师与石头剪刀布

题意：
每一个大人曾经都是一个小孩，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’ 是大小写敏感的。
解题思路：
遍历一遍先把能赢的回合求出来，求出能赢的回合数，再遍历一遍找出剩下的回合，将剩余的出招进行分配，比较能赢得回合和总回合/2即可。
代码实现：

#include <iostream>
#include <string>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;

int main()
{
	int N;
	scanf("%d",&N);
	while(N--)
	{
		int n,a,b,c;
		scanf("%d %d %d %d",&n,&a,&b,&c);
		char s[105],final[105];
		scanf("%s",s+1);
		int cnt=0;
		for (int i=1;i<=n;i++)
			final[i]='0';
		for (int i=1;i<=n;i++)
		{
			if(s[i]=='S'&&a>0)
			{
				final[i]='R';
				a--;
				cnt++;
			}
			
			else if(s[i]=='R'&&b>0)
			{
				final[i]='P';
				b--;
				cnt++;
			}
			
			else if(s[i]=='p'&&c>0)
			{
				final[i]='S';
				c--;
				cnt++;
			}
					
		}
		
		for (int i=1;i<=n;i++)
		{
			if(final[i]=='0')
			{
				if(a>0)
				{
					final[i]='R';
					a--;
				}
				else if(b>0)
				{
					final[i]='P';
					b--;
				}
				else if(c>0)
				{
					final[i]='S';
					c--;
				}
				
			}
		}
		
		if(cnt>=(n/2+n%2))
		{
			printf("YES\n");
			for (int i=1;i<=n;i++)
				printf("%c",final[i]);
			printf("\n");
		}
		else 
			printf("NO\n");
		
	}
	return 0;
}

B - 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， nm 的总和不超过 4 * 10^5
Output
答案输出 q 行, 第 i 行包含一个整数 — 表示第 i 组数据的答案
解题思路：
先横向搜索矩阵，将所有点存储当前行内 ’ . ’ 的数量，然后再纵向搜索，加上之前的数量，时间复杂度为O(n2）。如果当前点是 ’ . ’ ，数量减1。
注意当每个点都重新搜索横向和纵向时复杂度会达到O(n3)，会超时。
还要注意数组的大小，开全局变量会爆空间，可以输入n和m后，再开一个合适大小的数组a[n+10][m+10]。
代码实现：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <string>
using namespace std;

int q,n,m,ans;

int getint()
{
    int x=0,s=1; char ch=' ';
    while(ch<'0' || ch>'9'){ ch=getchar(); if(ch=='-') s=-1;}
    while(ch>='0' && ch<='9'){ x=x*10+ch-'0'; ch=getchar();}
    return x*s;
}

int main()
{
    scanf("%d",&q);
    while(q--)
	{
        scanf("%d %d",&n,&m);
        char ch[n+10][m+10];
        int a[n+10][m+10];
        memset(ch,0,sizeof(ch));
        memset(a,0,sizeof(a));
        ans=INT_MAX/2;
        for (int i=1; i<=n; i++)
            for (int j=1; j<=m; j++)
                scanf(" %c",&ch[i][j]);

        for (int i=1; i<=n; i++){
            int temp=0;
            for (int j=1; j<=m; j++){
                if(ch[i][j]=='.') temp++;
            }
            for (int j=1; j<=m; j++){
                a[i][j]=temp;
            }
        }
        for (int j=1; j<=m; j++){
            int temp=0;
            for (int i=1; i<=n; i++){
                if(ch[i][j]=='.') temp++;
            }
            for (int i=1; i<=n; i++){
                a[i][j]+=temp;
                if(ch[i][j]=='.') a[i][j]--;
            }
        }
        for (int i=1; i<=n; i++){
            for (int j=1; j<=m; j++){
                ans=min(ans,a[i][j]);
            }
        }
        printf("%d\n",ans);
    }

    return 0;
}

C - 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老师的考验吗？
输入输出：
Input
输出文件包含多组测试用例，每组测试用例格式如下：
第一行给定两个正整数 k m。（k < 2e9, m < 1e5）
第二行给定十个整数，分别表示 a0～a9。
Output
对于每一组测试用例输出一行，表示 f(k) % m 的数值大小。
解题思路：
重点在于构造矩阵相乘的递推式，然后利用矩阵快速幂优化递推。
x<10时f(x)=x
x>10时f(x)=a0f(x-1)+a1f(x-2)+…+a9f(x-10)
代码实现：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
using namespace std;
typedef long long ll;
ll k,mod,a[20];
const int N=10;

struct matrix{
	matrix(){memset(a,0,sizeof(a));}
	matrix(const matrix &s){memcpy(a,s.a,sizeof(a));}
	matrix operator*(const matrix &s){
		matrix t;
		for(int i=0;i<N;i++){
			for(int j=0;j<N;j++){
				for(int k=0;k<N;k++){
					t.a[i][j]+=a[i][k]*s.a[k][j]%mod;
					t.a[i][j]%=mod;
				}				
			}	
		}					
		return t;
	}
	ll a[N][N];
};

void solve(matrix x,int n){
	matrix ans;
	for(int i=0;i<10;i++)
		ans.a[i][i]=1;
	while(n){
		if(n&1) ans=ans*x;
		x=x*x;
		n>>=1;
	}
	ll res=0;
	for(int i=0;i<10;i++){
		res+=ans.a[0][i]*(10-i-1)%mod;
	}
	printf("%lld\n",res%mod);
}

int main()
{
	matrix x;
	while(~scanf("%lld%lld",&k,&mod)){
		for(int i=0;i<10;i++)
			scanf("%lld",&a[i]);
		if(k<10){
			printf("%lld\n",k%mod);
			continue;
		}
		memset(x.a,0,sizeof(x.a));
		for(int j=0;j<10;j++)
			x.a[0][j]=a[j];
		for(int i=1;i<10;i++)
			x.a[i][i-1]=1;
		solve(x,k-9);
	}
	return 0;
}

e青青青

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
Week14 矩阵快速幂

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

扫一扫