矩阵专题(续)

最新推荐文章于 2022-01-29 16:57:38 发布
最新推荐文章于 2022-01-29 16:57:38 发布
阅读量1.1k 收藏
点赞数
分类专栏： 矩阵 文章标签： 矩阵
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/shengtao96/article/details/51849216
版权

哈哈哈，失踪人员回归。最近总是在忙一些乱七八糟的事情，终于有了一段空闲时间可以刷刷题了。学姐给我准备了kuangbin先生的矩阵专题，作为我的矩阵终结专题。（初学者请先看我的矩阵专题前三篇）


第一题 CodeForces 450B

分析：这道题看似是一道矩阵题，但是实际上是一道模拟题。看题目给出的公式 f(i)=f(i-1)+f(i+1),i>=2，那么转换一下得到了f(i)=f(i-1)-f(i-2),i>=3，当然这题可以使用矩阵快速幂来算，但是多此一举。不妨将f(1)=x  f(2)=y代进去计算前几项，那么依次得到 x   y   y-x   -x   -y   x-y   x   y   ......，我们发现这个数列是有循环节的，而且必定为6，那么只需要计算前六个，n取一下模即可得到答案。注意：x  y可能是负数，但是答案一定得是非负整数。

#include <map>
#include <set>
#include <ctime>
#include <stack>
#include <cmath>
#include <queue>
#include <bitset>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <fstream>
#include <cstdlib>
#include <sstream>
#include <cstring>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;
#define   maxn          200000+10
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   clr(x,y)      memset(x,y,sizeof(x))
#define   rep(i,n)      for(int i=0;i<(n);i++)
#define   repf(i,a,b)   for(int i=(a);i<=(b);i++)
#define   pii           pair<int,int>
#define   mp            make_pair
#define   FI            first
#define   SE            second
#define   IT            iterator
#define   PB            push_back
#define   Times         10

typedef   long long     ll;
typedef   unsigned long long ull;
typedef   long double   ld;

const double eps = 1e-10;
const double  pi = acos(-1.0);
const  ll    mod = 1e9+7;
const  int   inf = 0x3f3f3f3f;

inline void RI(int& x)
{
    x=0;
    char c=getchar();
    while(!((c>='0'&&c<='9')||c=='-'))c=getchar();
    bool flag=1;
    if(c=='-')
    {
        flag=0;
        c=getchar();
    }
    while(c<='9'&&c>='0')
    {
        x=x*10+c-'0';
        c=getchar();
    }
    if(!flag)x=-x;
}

//--------------------------------------------------

ll dat[10];
int main(){
	//freopen("d:\\acm\\in.in","r",stdin);
	ll n;
	scanf("%I64d %I64d",&dat[1],&dat[2]);
	dat[1]=(dat[1]+mod)%mod;
	dat[2]=(dat[2]+mod)%mod;
	scanf("%I64d",&n);
	dat[0]=(dat[1]-dat[2]+mod)%mod;
	for(int i=3;i<6;i++)
		dat[i]=(dat[i-1]-dat[i-2]+mod)%mod;
	printf("%I64d\n",dat[n%6]);
    return 0;
}



第二题  HDU 5015

题目大意：给你一个矩阵，矩阵第一行中，第一个数为空，之后的数就是233，2333，23333。。。一直轮下去。然后，矩阵第一列上，第一个数为空（即第一行第一列的那个数为空），之后的所有数会给出。接着，告诉你这个矩阵其他元素都是按照maze[i][j]=maze[i-1][j]+maze[i][j-1]的规律来计算的，让你求出矩阵右下角的数是多少。

分析：一开始没有往矩阵方面想，前段时间数论做多了之后很容易就发现最后的maze[n][m]的数可以通过已知的数乘一个数（满足一个组合数，杨辉三角），然后相加得到。虽然我公式推出来了，但是复杂度太高，而且写起来还麻烦，就被我pass了，有兴趣的朋友可以自己推一推写一写。然后回过去想矩阵做法的话，立刻就想到了（提示n的值很大，但是m的值不大）。将输入的23(第一个23是补充上去的，为了和之后的23333对齐),maze[1][0],maze[2][0],maze[3][0],......作为初始值，根据公式很容易看出转换的矩阵，除了第一个233是前一个23乘10加3，其他的都是maze[i][j]=maze[i-1][j]+maze[i][j-1]。


#include <map>
#include <set>
#include <ctime>
#include <stack>
#include <cmath>
#include <queue>
#include <bitset>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <fstream>
#include <cstdlib>
#include <sstream>
#include <cstring>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;
#define   maxn          10+10
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   clr(x,y)      memset(x,y,sizeof(x))
#define   rep(i,n)      for(int i=0;i<(n);i++)
#define   repf(i,a,b)   for(int i=(a);i<=(b);i++)
#define   pii           pair<int,int>
#define   mp            make_pair
#define   FI            first
#define   SE            second
#define   IT            iterator
#define   PB            push_back
#define   Times         10

typedef   long long     ll;
typedef   unsigned long long ull;
typedef   long double   ld;

const double eps = 1e-10;
const double  pi = acos(-1.0);
const  ll    mod = 1e9+7;
const  int   inf = 0x3f3f3f3f;

inline void RI(int& x)
{
    x=0;
    char c=getchar();
    while(!((c>='0'&&c<='9')||c=='-'))c=getchar();
    bool flag=1;
    if(c=='-')
    {
        flag=0;
        c=getchar();
    }
    while(c<='9'&&c>='0')
    {
        x=x*10+c-'0';
        c=getchar();
    }
    if(!flag)x=-x;
}

//--------------------------------------------------

ll mul(ll a,ll b){
	ll ans=0;
	while(b){
		if(b&1)ans=(ans+a)%mod;
		a=(a+a)%mod;
		b>>=1;
	}
	return ans;
}
struct matrix{
	int n;
	ll maze[maxn][maxn];
	void init(int n){
		this->n=n;
		clr(maze,0);
	}
	matrix operator * (const matrix& rhs){
		matrix ans;
		ans.init(n);
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				for(int k=0;k<n;k++)
					ans.maze[i][j]=(ans.maze[i][j]+maze[i][k]*rhs.maze[k][j])%10000007;
		return ans;
	}
};
matrix qlow(matrix a,int n){
	matrix ans;
	ans.init(a.n);
	for(int i=0;i<a.n;i++)ans.maze[i][i]=1;
	while(n){
		if(n&1)ans=ans*a;
		a=a*a;
		n>>=1;
	}
	return ans;
}
int main(){
	//freopen("d:\\acm\\in.in","r",stdin);
	int n,m;
	while(~scanf("%d %d",&n,&m)){
		matrix ans;
		ans.init(n+2);
		for(int i=1;i<=n;i++)
			scanf("%lld",&ans.maze[0][i]);
		ans.maze[0][0]=23;
		ans.maze[0][n+1]=1;
		matrix ant;
		ant.init(n+2);
		for(int i=0;i<=n;i++)
			ant.maze[0][i]=10,ant.maze[n+1][i]=3;
		ant.maze[n+1][n+1]=1;
		for(int i=1;i<=n;i++)
			for(int j=i;j<=n;j++)
				ant.maze[i][j]=1;
		ant=qlow(ant,m);
		ans=ans*ant;
		printf("%lld\n",ans.maze[0][n]);
	}
    return 0;
}




第三题 HDU 4990
分析：简单来看就是一个数列计算问题。F(0)=0.如果i为奇数，那么F(i)=2*F(i-1)+1；如果i为偶数，那么F(i)=2*F(i-1)。很显然是一道矩阵问题，我们不妨将两步合成一步来做F(i)=4*F(i-2)+2,i为偶数，那么我们按照2个2个来计算，计算到不超过n的最大偶数。最后，如果n为奇数的话，再计算一步即可。 


#include <map>
#include <set>
#include <ctime>
#include <stack>
#include <cmath>
#include <queue>
#include <bitset>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
#include <fstream>
#include <cstdlib>
#include <sstream>
#include <cstring>
#include <iostream>
#include <algorithm>
#pragma comment(linker, "/STACK:1024000000,1024000000")

using namespace std;
#define   maxn         	0+10
#define   lson          l,m,rt<<1
#define   rson          m+1,r,rt<<1|1
#define   clr(x,y)      memset(x,y,sizeof(x))
#define   rep(i,n)      for(int i=0;i<(n);i++)
#define   repf(i,a,b)   for(int i=(a);i<=(b);i++)
#define   pii           pair<int,int>
#define   mp            make_pair
#define   FI            first
#define   SE            second
#define   IT            iterator
#define   PB            push_back
#define   Times         10

typedef   long long     ll;
typedef   unsigned long long ull;
typedef   long double   ld;

const double eps = 1e-10;
const double  pi = acos(-1.0);
const  ll    mod = 1e9+7;
const  int   inf = 0x3f3f3f3f;

inline void RI(int& x)
{
    x=0;
    char c=getchar();
    while(!((c>='0'&&c<='9')||c=='-'))c=getchar();
    bool flag=1;
    if(c=='-')
    {
        flag=0;
        c=getchar();
    }
    while(c<='9'&&c>='0')
    {
        x=x*10+c-'0';
        c=getchar();
    }
    if(!flag)x=-x;
}

//--------------------------------------------------

int m;
struct matrix{
	int n;
	ll maze[maxn][maxn];
	void init(int n){
		this->n=n;
		clr(maze,0);
	}
	matrix operator * (const matrix& rhs){
		matrix ans;
		ans.init(n);
		for(int i=0;i<n;i++)
			for(int j=0;j<n;j++)
				for(int k=0;k<n;k++)
					ans.maze[i][j]=(ans.maze[i][j]+maze[i][k]*rhs.maze[k][j])%m;
		return ans;
	}
};
matrix qlow(matrix a,ll n){
	matrix ans;
	ans.init(a.n);
	for(int i=0;i<a.n;i++)ans.maze[i][i]=1;
	while(n){
		if(n&1)ans=ans*a;
		a=a*a;
		n>>=1;
	}
	return ans;
}
int main(){
	//freopen("d:\\acm\\in.in","r",stdin);
	int n;
	while(~scanf("%d %d",&n,&m)){
		matrix ant;
		ant.init(2);
		ant.maze[0][0]=4;
		ant.maze[1][0]=2;
		ant.maze[1][1]=1;
		ant=qlow(ant,n/2);
		ll ans=ant.maze[1][0];
		if(n&1)ans=(ans*2+1)%m;
		cout<<ans<<endl;
	}
    return 0;
}




第四题 UVA 11651
分析：不会做，看了别人的题解才会的。这里我讲一下思路，作为纯DP题的话，转移方程很好想，但是score太大，无法记忆化搜索；但是作为纯矩阵题的话，一开始完全想不到如何构造矩阵。然后把两者合在一起想，就发现还是矩阵的套路，但是要加上DP的思想（以DP的状态为转移）。将一开始所有可能转移的数据全部算出来作为初始数据，部分数据直接移位就行了，主要就是最后一行的数据（要计算的新数据），需要用到前面所有数据来计算。

我知道读者肯定没有看懂，无妨，请转到这位大神的博客，他写得很好，这里我贴下我的代码

http://www.cnblogs.com/AOQNRMGYXLMV/p/5256508.html

#include <map>
#include <set>
#include <ctime>
#include <stack>
#include <cmath>
#include <queue>
#include <bitset>
#include <string>
#include <vector>
#include <cstdio>
#include <cctype>
最低0.47元/天 解锁文章
shengtao96
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
内存矩阵标签专题图
12-14
本示例聚焦于"内存矩阵标签专题图"，这是一种将数据可视化的重要方法，尤其适用于处理大量地理数据时。下面我们将详细探讨这个主题。 首先，内存矩阵标签专题图是SuperMap iClient for JavaScript中的一个功能，...
const int INF=0x3f3f3f3f;——ACM中的无穷大常量
while_black的博客
04-07 4483
0x3f3f3f3f 0x开头的 是十六进制常数， 等于 十进制 1061109567 等于 二进制： 00111111 00111111 00111111 00111111 const int inf=0x3f3f3f3f -- 声明 inf 是 const int型 , const 表示 inf 一旦有了值以后就不允许（通过赋值 来）改变它的值. 在算法竞赛中，我们常常需要用到设置一个...
INF(infinity) 无穷大 const int INF=0x3f3f3f3f；
这题AC再睡
01-29 665
INF(infinity) 无穷大 const int INF=0x3f3f3f3f;
ACM中关于无穷常数的设置
Kiritow的学园
08-02 1887
无穷大 const int INF = 0x3f3f3f3f; 这样设置的好处是可以通过memset直接对连续的内存进行无穷大初始化，如下 memset(blockptr,0x3f,blocksize); 无穷小 const int NINF = -INF -1; const int NINF = 0xc0c0c0c0; 与无穷大类似，可以使用memset对连续内存进行无穷小初
INF=0x3f3f3f3f背后的知识
Erya Programming learning
09-24 294
我们知道计算机都是按照二进制存储的。在日常练习过程中，我们都会遇到INF无穷大量的问题，那么我们该怎么定义这么一个数字，下面便有一些说明： 0x3f3f3f3f的十进制是1061109567，也就是10^ 9级别的（和0x7fffffff一个数量级），而一般场合下的数据都是小于10^9的，所以它可以作为无穷大使用而不致出现数据大于无穷大的情形。 另一方面，由于一般的数据都不会大于10^9，...
[联合集训6-19] K小数查询 分块+二分答案
DOFYPXY的博客
06-21 331
有一种常数比较小的O(nn−−√log2n)O(nnlog2⁡n)O(n\sqrt n \log^2 n)的做法。 分块，每个块维护一个其中元素排好序之后的数组。修改的时候零散块直接重构，整块打标记。询问的时候先二分答案midmidmid，那么就转化成求小于等于midmidmid的数的个数，对于零散的块重构后暴力数，整块的如果标记≤mid≤mid\le mid答案就是块大小，否则直接在数组上二分...
silverlight实现矩阵标签专题图
12-14
《使用SuperMap iClient实现矩阵标签专题图》 在数据可视化领域，标签专题图是一种有效的方式，用于在地图上清晰地展示数据信息。而矩阵标签专题图则在此基础上增加了更多的表现力，它允许混合使用图片、符号和文字...
JavaScript关联外表做标签矩阵专题图
12-14
本实例主要讲解了如何利用JavaScript结合SuperMap iClient 6R来创建一个标签矩阵专题图，这在地理数据分析和可视化中具有广泛的应用。 首先，我们要理解“关联外表”这个概念。在GIS中，关联外表是指将非空间数据...
天气预报矩阵标签专题图
12-13
本文将深入探讨如何使用SuperMap iClient 6R for JavaScript这一强大的Web GIS开发框架来创建一个天气预报矩阵标签专题图，该图能够高效地展示全国城市的一日天气情况，包括城市名、天气状况以及当日的最高和最低...
python 实现矩阵相加
06-09
# 题目：计算两个矩阵相加。 # 程序分析：创建一个新的矩阵，使用 for 迭代并取出 X 和 Y 矩阵中对应位置的值，相加后放到新矩阵的对应位置中。
生命值区域，全部位置设为空格。然后在原矩阵中取出一个10*30的矩阵，矩阵中的每个位置
07-17
c++应用实例//这是一个虚拟生命游戏。先初始化一个24*80的矩阵区域，全部位置设为空格。然后在原矩阵中取出一个10*30的矩阵，矩阵中的每个位置用X去替换空格， //再显示到屏幕上。然后遵循虚拟生命进化准则：当原来的24*80矩阵中的每个位置前后左右的八个位置，即前、后、左、右、前左、前右、后左、后右, //如果一个位置周围的八个位置中（有的会不足八个），如果最多只有一个X或者多于3个X，那么原来的X会变为空格；如果一个空格位置周围恰好有三个空格， //那么这个位置会变为X。每按一次空格，程序运行一次，并把每次运行后的结果显示在屏幕上。 #include <iostream>：当原来的24*80矩阵中的每个位置前后左右的八个位置，即前、后、左、右、前左、前右、后左、后右, //如果一个位置周围的八个位置中（有的会不足八个），如果最多只有一个X或者多于3个X，那么原来的X会变为空格；如果一个空格位置周围恰好有三个空格， //那么这个位置会变为X。每按一次空格，程序运行一次，并把每次运行后的结果显示在屏幕上。 #include <iostream>
INF=0x3f3f3f3f是 什么意思？
天空的味道
09-06 6335
经常求最短路的时候好多人都把inf 设为0x3f3f3f3f ,今天让我们来看看这个奇妙的无穷大
编程中无穷大常量的设定技巧
Bob Lee's Blog For CS
11-22 4996
本文转自:http://aikilis.tk/2471.html 编程中无穷大常量的设定技巧 如果问题中各数据的范围明确，那么无穷大的设定不是问题，在不明确的情况下，很多程序员都取0x7fffffff作为无穷大，因为这是32-bit int的最大值。如果这个无穷大只用于一般的比较（比如求最小值时min变量的初值），那么0x7fffffff确实是一个完美的选择，但是在更多的情况下，0x7
多机调度问题 贪心策略
The_star_is_at的博客
10-21 1299
const int maxn=100; const int inf=0x3f3f3f3f; int a[maxn],n,m; struct node { int data; bool operator<(const node &a)const { return data>a.data; } }; priority_queueque; int main
INF=0x3f3f3f3f 是什么意思？
sudo-wang的博客
03-16 2004
为什么把无穷大INF设置成0x3f3f3f3f?
xyry
02-16 2662
INF为什么设置成为0x3f3f3f3f
bellman-ford算法模板
yuanba_xs的博客
01-31 318
#include #include #include #define MAXN 1005 const int inf = 0x3f3f3f3f; using namespace std; int n, m; int dis[MAXN],u[MAXN],v[MAXN],w[MAXN]; void bellford(int start) { for (int i = 1; i <= n
Alg5-习题课w1
最新发布
08-03
"本资料涵盖了矩阵理论中的多个重要知识点，包括正交矩阵的证明、向量组的正交化、特征值与特征向量的计算、方阵的特征多项式、矩阵的对角化以及二次型的标准化等。" 在数学的线性代数领域，矩阵是研究的主要对象之...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值