矩阵乘法的应用（hdu1575，hdu1588，poj3070，poj3233，poj3613）

最新推荐文章于 2022-04-19 20:43:25 发布
最新推荐文章于 2022-04-19 20:43:25 发布
阅读量2.5k 收藏
点赞数 2
分类专栏： 矩阵 文章标签： ini matrix struct c io
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/l04205613/article/details/6924306
版权

首先感叹下矩阵的强大，真心强大！

然后推荐下 Matrix67 大神的神作：十个利用矩阵乘法解决的经典题目

大神说的很明白了，我觉得再说啥都多余了，请直接去大神那里看吧，我只是刷了几道简单的入门题，贴下代码，和大家交流下：


hdu1575 Tr A

矩阵乘法 + 快速幂，很简单

代码：

#include<cstdio>

const int mod = 9973;
const int N = 11;

int n;
struct prog{
    int a[N][N],i,j;
    void init(){
        for(i=0;i<N;i++){
            for(j=0;j<N;j++){
                if(i==j) a[i][j]=1;
                else a[i][j]=0;
            }
        }
    }
};

prog matrixmult(prog a,prog b){
    int i,j,k;
    prog c;
    for(i=0;i<n;i++){
        for(j=0;j<n;j++){
            c.a[i][j]=0;
            for(k=0;k<n;k++)
                c.a[i][j]+=((a.a[i][k]*b.a[k][j])%mod);
            c.a[i][j]%=mod;
        }
    }
    return c;
}

prog mult(prog s,int k){
    prog ans;
    ans.init();
    while(k>=1){
        if(k&1)ans=matrixmult(ans,s);
        k=k>>1;
        s=matrixmult(s,s);
    }
    return ans;
}

int main()
{
    int i,j,k,cases;
    prog s;
    scanf("%d",&cases);
    while(cases--)
    {
        scanf("%d%d",&n,&k);
        for(i=0;i<n;i++){
            for(j=0;j<n;j++)scanf("%d",&s.a[i][j]);
        }

        s=mult(s,k);

        int sum=0;
        for(i=0;i<n;i++)
            sum+=s.a[i][i];

        printf("%d\n",sum%mod);
    }
    return 0;
}


hdu1588 Gauss Fibonacci

这个题有点意思，我是没想到怎么处理，这里有篇给力的文章：

http://www.cnblogs.com/Knuth/archive/2009/09/04/1559951.html


他给了两个思路，亚伟用第一个方法过的（代码在这），我用的第二个方法：

/*
A^b*( I + A^k + (A^k)^2 + .... + (A^k)^(N-1) )

设A^k=B

我们来设置这样一个矩阵
B I
O I
其中O是零矩阵，I是单位矩阵

将它乘方，得到
B^2 I+B
O   I

乘三方，得到
B^3 I+B+B^2
O   I

乘四方，得到
B^4 I+B+B^2+B^3
O   I

既然已经转换成矩阵的幂了，继续用我们的二分或者二进制法，直接求出幂就可以了
*/
#include<cstdio>
#include<cstring>

#define INT __int64

struct prog{
	INT a[2][2];
	void init(){
		a[0][0]=a[0][1]=a[1][0]=1;
		a[1][1]=0;
	}
	/*  斐波那契矩阵
	1,1
	1,0
	*/

	void Print(){//打印输出，方便检查
		puts("***********************");
		for(int i=0;i<2;i++){
			for(int j=0;j<2;j++)printf("%d ",a[i][j]);
			puts("");
		}puts("***********************");
	}

};

struct prog_x{//四位构造矩阵
	INT a[4][4];
	void init(){
		memset(a,0,sizeof(a));
		a[0][2]=a[1][3]=a[2][2]=a[3][3]=1;
	}
	/*
	0,0,1,0
	0,0,0,1
	0,0,1,0
	0,0,0,1
	*/
	
	void Print(){//打印输出，方便检查
		puts("***********************");
		for(int i=0;i<4;i++){
			for(int j=0;j<4;j++)printf("%d ",a[i][j]);
			puts("");
		}puts("***********************");
	}
};

int mod;

prog matrixmult(prog a,prog b){
	int i,j,k;
	prog c;
	for(i=0;i<2;i++){
		for(j=0;j<2;j++){
			c.a[i][j]=0;
			for(k=0;k<2;k++)
				c.a[i][j]+=((a.a[i][k]*b.a[k][j])%mod);
			c.a[i][j]%=mod;
		}
	}
	return c;
}

prog mult(prog s,int k){
	prog ans;
	ans.init();
	while(k>=1){
		if(k&1)ans=matrixmult(ans,s);
		k=k>>1;
		s=matrixmult(s,s);
	}
	return ans;
}

prog_x matrixmult_x(prog_x a,prog_x b){
	int i,j,k;
	prog_x c;
	for(i=0;i<4;i++){
		for(j=0;j<4;j++){
			c.a[i][j]=0;
			for(k=0;k<4;k++)
				c.a[i][j]+=((a.a[i][k]*b.a[k][j])%mod);
			c.a[i][j]%=mod;
		}
	}
	return c;
}

prog_x mult_x(prog_x s,int k){
	prog_x ans;
	int i,j;
	for(i=0;i<4;i++){
		for(j=0;j<4;j++)ans.a[i][j]=s.a[i][j];
	}
	while(k>0){
		if(k&1)ans=matrixmult_x(ans,s);
		k=k>>1;
		s=matrixmult_x(s,s);
	}
	return ans;
}

int main()
{
	int k,b,n;
	while(~scanf("%d%d%d%d",&k,&b,&n,&mod))
	{
		prog B,tmp;
		B.init();
		B=mult(B,k-1);					//B=A^k

		/*tmp=A^b,当b==0的时候,tmp=E*/
		if(b==0){
			tmp.a[0][0]=tmp.a[1][1]=1;
			tmp.a[0][1]=tmp.a[1][0]=0;
		}
		else{
			tmp.init();
			tmp=mult(tmp,b-1);			
		}

		prog_x s_x;
		s_x.init();

		s_x.a[0][0]=B.a[0][0];			//将B=A^k填充进构造矩阵的左上角
		s_x.a[0][1]=B.a[0][1];
		s_x.a[1][0]=B.a[1][0];
		s_x.a[1][1]=B.a[1][1];

		//s_x.Print();

		s_x=mult_x(s_x,n-1);

		prog ans;						//取出右上角的2*2的部分
		ans.a[0][0]=s_x.a[0][2];
		ans.a[0][1]=s_x.a[0][3];
		ans.a[1][0]=s_x.a[1][2];
		ans.a[1][1]=s_x.a[1][3];

		ans=matrixmult(ans,tmp);

		printf("%d\n",ans.a[0][1]%mod);
	}
	return 0;
}


poj3070 Fibonacci

最简单的，有人说过是雅礼下至3岁幼童都会的一道题，= =！

#include<cstdio>
struct prog{
	int a[2][2];
	void init(){
		a[1][1]=a[0][1]=a[1][0]=1;
		a[0][0]=0;
	}
};

prog matrixmul(prog a,prog b){
	int i,j,k;
	prog c;
	for(i=0;i<2;i++){
		for(j=0;j<2;j++){
			c.a[i][j]=0;
			for(k=0;k<2;k++)
				c.a[i][j]+=(a.a[i][k]*b.a[k][j]);
			c.a[i][j]%=10000;
		}
	}
	return c;
}

prog mul(prog s,int k){
	prog ans;
	ans.init();
	while(k>=1){
		if(k&1)ans=matrixmul(ans,s);
		k=k>>1;
		s=matrixmul(s,s);
	}
	return ans;
}

int main()
{
	int n;
	while(scanf("%d",&n),n!=-1){
		if(!n){
			printf("0\n");continue;
		}
		prog s;
		s.init();
		s=mul(s,n-1);
		printf("%d\n",s.a[0][1]%10000);
	}
	return 0;
}


poj3233 Matrix Power Series

矩阵的等比数列和，这个用二分实现，表示很经典，很强大，代码不是自己的，贴来Ym，理解：

/*
递归求解
K为偶数的时候不需要计算大括号里面的那一项
S_(K) = A+A^2+A^3+...+A^K 
  = (1+A^(K/2)) * (A+A^2+A^3+...+A^(K/2)) + {A^K}
  = (1+A^(K/2)) * (S_(K/2))
*/

#include<cstdio>
const int N = 30;

int n,mod;
struct prog{
	int a[N][N],i,j;
	void init(){
		for(i=0;i<N;i++){
			for(j=0;j<N;j++){
				if(i==j) a[i][j]=1;
				else a[i][j]=0;
			}
		}
	}
};

prog matrixmult(prog a,prog b){
	int i,j,k;
	prog c;
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			c.a[i][j]=0;
			for(k=0;k<n;k++)
				c.a[i][j]+=((a.a[i][k]*b.a[k][j])%mod);
			c.a[i][j]%=mod;
		}
	}
	return c;
}

prog mult(prog s,int k){
	prog ans;
	ans.init();
	while(k>=1){
		if(k&1)ans=matrixmult(ans,s);
		k=k>>1;
		s=matrixmult(s,s);
	}
	return ans;
}

prog matriplus(prog a,prog b){
	int i,j;
	for(i=0;i<n;i++){
		for(j=0;j<n;j++){
			a.a[i][j]+=b.a[i][j];
			a.a[i][j]%=mod;
		}
	}
	return a;
}

prog dfs(prog s,int k){
	if(k==1)return s;
	prog tmp;
	tmp.init();					//tmp=1
	tmp=matriplus(tmp,mult(s,k>>1));		//tmp = 1+A^(k/2)
	tmp=matrixmult(tmp,dfs(s,k>>1));		//tmp = ( 1+A^(k/2) )*( A+A^2+A^3+...+A^(K/2) )
	if(k&1) tmp=matriplus(tmp,mult(s,k));//奇，偶判断
	return tmp;
}

int main()
{
	int i,j,k;
	prog s;

	scanf("%d%d%d",&n,&k,&mod);

	for(i=0;i<n;i++){
		for(j=0;j<n;j++)
			scanf("%d",&s.a[i][j]);
	}

	s=dfs(s,k);

	for(i=0;i<n;i++){
		printf("%d",s.a[i][0]);
		for(j=1;j<n;j++){
			printf(" %d",s.a[i][j]);
		}puts("");
	}
	return 0;
}

poj3613 Cow Relays

求图中从一个点到另一个点恰好经过K条边的最短路径的权值和

这个题需要改下 矩阵乘法 以及 矩阵快速幂 的函数，把矩阵乘法里面的乘改成类似 floyd 的松弛（因为求的是最短路径权值和），还有要注意的是，这个题的矩阵规模较大，离散处理后也有 200*200 ，所以这个矩阵是不能定义在结构体内部的，用全局变量好了，我就是这么干的。

代码：

#include<cstdio>
#include<cstring>
#include<climits>
#define min(a,b) ((a)<(b))?(a):(b)

const int N = 205;
const int INF = 1000000100;

int m,p[1001];
int ans[N][N];

void matrixmult(int a[N][N],int b[N][N]){
	int i,j,k;
	int c[N][N];
	for(i=0;i<m;i++){
		for(j=0;j<m;j++){
			c[i][j]=INF;
			for(k=0;k<m;k++)
				c[i][j]=min( a[i][k]+b[k][j] , c[i][j] );
		}
	}
	for(i=0;i<m;i++){
		for(j=0;j<m;j++)a[i][j]=c[i][j];
	}
}

void mult(int s[N][N],int k){
	int i,j;
	for(i=0;i<m;i++){
		for(j=0;j<m;j++)
			ans[i][j]=s[i][j];
	}
	k--;
	while(k>0){
		if(k&1) matrixmult(ans,s);
		k=k>>1;
		matrixmult(s,s);
	}
}

int main()
{
	int mat[N][N];
	int i,j,n,t,s,e,from,to,val;
	while(~scanf("%d%d%d%d",&n,&t,&s,&e))
	{
		for(i=0;i<N;i++){
			for(j=0;j<N;j++)mat[i][j]=INF;
		}
		memset(p,-1,sizeof(p));
		m=0;
		for(i=0;i<t;i++){
			scanf("%d%d%d",&val,&from,&to);
			if(p[from]<0) p[from] = m++;
			if(p[to]<0) p[to] = m++;
			if(val<mat[ p[from] ][ p[to] ])
				mat[ p[from] ][ p[to] ] = mat[ p[to] ][ p[from] ] = val;
		}
		s=p[s]; e=p[e];
		mult(mat,n);
		printf("%d\n",ans[s][e]);
	}
	return 0;
}



l04205613
关注
专栏目录
poj3233 hdu 1588 (矩阵)
zylm的专栏
09-27 633
poj3233 做其他的题，要用到矩阵求和的东西嘛，然后
[HDU1575]Tr A（矩阵乘法）
zyf2000
05-22 829
再没有眼泪，再没有恸哭，我只是要靠着你，一如过去的年年月月。
HDU2716 POJ2141 Message Decowding【密码】
海岛Blog
02-19 204
Message Decowding Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 1934 Accepted Submission(s): 834 Problem Description The cows are thrilled because they’ve just learned about encrypting messages.
UVA1566 LA3830 POJ3480 HDU1907 John【博弈】
海岛Blog
04-16 362
John Time Limit: 1000MS Memory Limit: 65536K Total Submissions: 3727 Accepted: 1863 Description Little John is playing very funny game with his younger brother. There is one big box filled with M&am...
最长回文 （hdu3068 && poj3974）Palindrome
寻找&星空の孩子
08-23 1204
最长回文 Time Limit: 4000/2000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 11771    Accepted Submission(s): 4315 Problem Description   给出一个只由小写英文字符a,b,c..
【矩阵加速】POJ3070 & HDU3117 -- 矩阵快速幂数列递推
m0_59701064的博客
04-19 576
矩阵快速幂 给定一个m*m的矩阵A，定义一个结构体以实现矩阵，同时重载矩阵的乘法操作。注意有时需进行取模操作。 const int N=2;//矩阵(方阵)大小 struct Matrix { ll m[N][N]; Matrix() { memset(m,0,sizeof m); } Matrix operator*(Matrix &a) { Matrix ans; memset(ans.m,0,sizeof(ans.m)); for(int i=0;i<N;i
HDU1018 POJ1423 UVALive2697 UVA1185 ZOJ1526 Big Number【阶乘位数】
海岛Blog
11-15 782
Big Number Time Limit: 1000MS   Memory Limit: 65536K Total Submissions: 28524   Accepted: 9075 Description In many applications very large integers numbers are required.
poj 1182,poj 1703,hdu 1829
winhcc的博客
04-26 174
这一类的题目都是n个东西属于不同的集合，集合中的情况要么同时发生，要么同时不发生。 poj 1182 是书上的题目，有详解 poj 1703 和 hdu 1829 都是只有两个集合 poj 1703 #include <cstdio> const int N = 1e5 * 2 + 10; char s[2]; int n,m,x,y,fa[N],t; int find(int x)...
HDU 1402 POJ 2389 BZOJ 2179 大整数乘法 FFT
とある黄鱼のOIブログ
03-20 972
FFT果然还是看算导来得快。。 网上的讲解糊里糊涂的 不过我还是看了2小时QAQ递归形式function recursive-fft(a) { n = a.length // 必为2的幂 if (n == 1) return a wn = Math.pow(Math.e, 2 * Math.PI * Math.i / n) // wn还可以利用复数
博弈问题解析：HDU Online Judge与Poj1067的石子游戏策略
这篇资源主要涉及的是博弈问题的解决方法，特别是关于HDU Online Judge上的几道题目，包括HDU 2516, POJ 1067, HDU 1527, HDU 2177以及HDU 2176。这些题目都是关于石子游戏的，玩家需要在一定的规则下决定如何取石子...
hdu oj poj 题目代码与详解
09-08
【题目代码与详解】是针对在线编程竞赛平台HDU（杭州电子科技大学在线评测系统）和POJ（普林斯顿大学在线判题系统）中的经典题目所编写的解决方案和解析的集合。这些平台是程序员和计算机科学学生提升算法技能、熟悉...
jihe.rar_2289_POJ 3714_poj3714_poj3714 Ra_visual c
09-22
描述中提到的"三道几何题：hdu 1007、hdu 2289、poj 3714"暗示了这是一个关于几何算法的题目集，涉及到三个不同的在线判题系统：HDU（杭州电子科技大学在线评测系统）的1007题和2289题，以及POJ的3714题。...
leetcode下载-algorithm-1:力扣、HDU、ZOJ、POJ
06-29
Interview，ZOJ，POJ 等平台。 欢迎Coders对代码加以指正和提议！ 常见问题总结 两整数求平均值 average = min + (max - min) / 2 防止两整数的和越界 整数乘积对比 1.0 * m * m == num 类似乘积对比， 需转为double...
K8s 很难么？带你从头到尾捋一遍，不信你学不会
民工哥的博客
02-23 9113
点击关注公众号，回复“1024”获取2TB学习资源！为什么要学习 Kubernetes?虽然 Docker 已经很强大了，但是在实际使用上还是有诸多不便，比如集群管理、资源调度、文件管理等...
10天智能锁项目实战第1天（了解单片机STM32F401RET6和C语言基础）
北豼不太皮的博客
02-21 3640
10天智能锁项目实战第1天（了解单片机STM32F401RET6和C语言基础）一、学习目标二、了解单片机STM32F401RET6三、C语言基础 一、学习目标 二、了解单片机STM32F401RET6 4、STM32F401RE特征 三、C语言基础 1.数据类型 常用2的次方： 2^7 = 128 2^8 = 256 2^15 = 32768 2^16= 65536 51单片机常见的数据类型： char: 占内存1字节 取值范围：-128~127 -2^7 ~ 2
Android开发（1） | Fragment 的应用——新闻应用
Jormungand_V的博客
02-21 7132
文章目录 news_item.xml： <TextView xmlns:android="http://schemas.android.com/apk/res/android" android:id="@+id/news_title" android:layout_width="match_parent" android:layout_height="wrap_content" android:lines="1" android:ellipsize="
系统学习 TypeScript（四）——变量声明的初步学习
编程三昧
02-25 1767
认识了 TypeScript 中的基础类型，接下来当然是变量声明的相关学习了，我在这里记录一下我学习过程中的一些总结。
清除浮动的五种方法
weixin_52212950的博客
02-22 2745
为什么要清除浮动？因为往往浮动后的子元素因为脱离了标准流，不能自动撑起父元素的高度，所以会对后续布局产生影响，对此清除浮动不如理解为清除浮动产生的影响更为合理。 例如：我们设置了两个盒子如图所示，粉色为父盒子，只有宽度没有高度，蓝色盒子有宽有高，粉色盒子的高由蓝盒子的高度撑开。 但是给蓝色的子盒子设置了左浮动后，脱离了标准流，无法再撑开粉盒子，可以看到粉盒子高度变成了0； 清除浮动有五种方法： 直接设置父元素的高度 额外标签法 单伪元素法 双伪元素法 ove.
HTML实现学习网站首页
做技术，坚持才是最重要的，一时的热情只是暂时的进步
02-20 4212
项目访问 reset.css /* http://meyerweb.com/eric/tools/css/reset/ v2.0 | 20110126 License: none (public domain) */ html, body, div, span, applet, object, iframe, h1, h2, h3, h4, h5, h6, p, blockquote, pre, a, abbr, acronym, address, big, cite, code, del
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值